A thermo-viscoelastic–viscoplastic–viscodamage constitutive model for asphaltic materials

A temperature-dependent viscodamage model is proposed and coupled to the temperature-dependent Schapery’s nonlinear viscoelasticity and the temperature-dependent Perzyna’s viscoplasticity constitutive model presented in Abu Al-Rub et al., 2009, Huang et al., in press in order to model the nonlinear constitutive behavior of asphalt mixes. The thermo-viscodamage model is formulated to be a function of temperature, total effective strain, and the damage driving force which is expressed in terms of the stress invariants of the effective stress in the undamaged configuration. This expression for the damage force allows for the distinction between the influence of compression and extension loading conditions on damage nucleation and growth. A systematic procedure for obtaining the thermo-viscodamage model parameters using creep test data at different stress levels and different temperatures is presented. The recursive-iterative and radial return algorithms are used for the numerical implementation of the nonlinear viscoelasticity and viscoplasticity models, respectively, whereas the viscodamage model is implemented using the effective (undamaged) configuration concept. Numerical algorithms are implemented in the well-known finite element code Abaqus via the user material subroutine UMAT. The model is then calibrated and verified by comparing the model predictions with experimental data that include creep-recovery, creep, and uniaxial constant strain rate tests over a range of temperatures, stress levels, and strain rates. It is shown that the presented constitutive model is capable of predicting the nonlinear behavior of asphaltic mixes under different loading conditions.     

A nonlocal phase-field model of Ginzburg–Landau–Korteweg fluids

A thermodynamic model of Korteweg fluids undergoing phase transition and/or phase separation is developed within the framework of weakly nonlocal thermodynamics. Compatibility with second law of thermodynamics is investigated by applying a generalized Liu procedure recently introduced in the literature. Possible forms of the free energy and of the stress tensor, which generalize some earlier ones proposed by several authors in the last decades, are carried out. Owing to the new procedure applied for exploiting the entropy principle, the thermodynamic potentials are allowed to depend on the whole set of variables spanning the state space, including the gradients of the unknown fields, without postulating neither the presence of an energy or entropy extra-flux, nor an additional balance law for microforce.     

A single-camera coupled PTV–LIF technique

 A single-camera coupled particle tracking velocimetry–laser-induced fluorescence (PTV–LIF) technique and validation results from an experiment in a neutrally buoyant turbulent round jet are presented. The single-camera implementation allows the use of a 12-bit 60 frame-per-second 1024 × 1024 pixel digital CCD camera capable of streaming images in real time to hard disk resulting in very accurate PTV and LIF with excellent spatial and temporal resolution. The technique is capable of determining the turbulent scalar flux, as well as the Reynolds stress and mean and fluctuating velocity and concentration fields. Details of dye choice, corrections for attenuation due to dye, particles and water, photobleaching, vignetting, CCD calibration, and illumination power and geometry corrections are presented. Detailed results from the validation experiment confirm the accuracy and resolution of the technique, and in particular, the ability to measure . Bootstrap 95% uncertainty intervals are presented for the calculated statistics. Received: 28 July 2000/Accepted: 8 November 2000     

A Reaction–Diffusion–Advection Equation with Mixed and Free Boundary Conditions

We investigate a reaction–diffusion–advection equation of the form \(u_t-u_{xx}+\beta u_x=f(u)\) \((t>0,\,0<x<h(t))\) with mixed boundary condition at \(x=0\) and Stefan free boundary condition at \(x=h(t)\). Such a model may be applied to describe the dynamical process of a new or invasive species adopting a combination of random movement and advection upward or downward along the resource gradient, with the free boundary representing the expanding front. The goal of this paper is to understand the effect of advection environment and no flux across the left boundary on the dynamics of this species. For the case \(|\beta |<c_0\), we first derive the spreading–vanishing dichotomy and sharp threshold for spreading and vanishing, and then provide a much sharper estimate for the spreading speed of h(t) and the uniform convergence of u(t, x) when spreading happens. For the case \(|\beta |\ge c_0\), some results concerning virtual spreading, vanishing and virtual vanishing are obtained. Here \(c_0\) is the minimal speed of traveling waves of the differential equation.     

A unified model for polystyrene–nanorod and polystyrene–nanoplatelet melt composites

We present a unified constitutive model capable of predicting the steady shear rheology of polystyrene (PS)–nanoparticle melt composites, where particles can be rods, platelets, or any geometry in between, as validated against experimental measurements. The composite model incorporates the rheological properties of the polymer matrix, the aspect ratio and characteristic length scale of the nanoparticles, the orientation of the nanoparticles, hydrodynamic particle–particle interactions, the interaction between the nanoparticles and the polymer, and flow conditions of melt processing. We demonstrate that our constitutive model predicts both the steady rheology of PS–carbon nanofiber composites and the steady rheology of PS–nanoclay composites. Along with presenting the model and validating it against experimental measurements, we evaluate three different closure approximations, an important constitutive assumption in a kinetic theory model, for both polymer–nanoparticle systems. Both composite systems are most accurately modeled with a quadratic closure approximation.     

A Beale–Kato–Majda Criterion for Three Dimensional Compressible Viscous Heat-Conductive Flows

We prove a blow-up criterion in terms of the upper bound of (ρ, ρ ⁻¹, θ) for a strong solution to three dimensional compressible viscous heat-conductive flows. The main ingredient of the proof is an a priori estimate for a quantity independently introduced in Haspot (Regularity of weak solutions of the compressible isentropic Navier–Stokes equation, arXiv:1001.1581, 2010) and Sun et al. (J Math Pure Appl 95:36–47, 2011), whose divergence can be viewed as the effective viscous flux.     

A space-time lower–upper symmetric Gauss–Seidel scheme for the time-spectral method

The time-spectral method (TSM) offers the advantage of increased order of accuracy compared to methods using finite-difference in time for periodic unsteady flow problems. Explicit Runge–Kutta pseudo-time marching and implicit schemes have been developed to solve iteratively the space-time coupled nonlinear equations resulting from TSM. Convergence of the explicit schemes is slow because of the stringent time-step limit. Many implicit methods have been developed for TSM. Their computational efficiency is, however, still limited in practice because of delayed implicit temporal coupling, multiple iterative loops, costly matrix operations, or lack of strong diagonal dominance of the implicit operator matrix. To overcome these shortcomings, an efficient space-time lower–upper symmetric Gauss–Seidel (ST-LU-SGS) implicit scheme with multigrid acceleration is presented. In this scheme, the implicit temporal coupling term is split as one additional dimension of space in the LU-SGS sweeps. To improve numerical stability for periodic flows with high frequency, a modification to the ST-LU-SGS scheme is proposed. Numerical results show that fast convergence is achieved using large or even infinite Courant–Friedrichs–Lewy (CFL) numbers for unsteady flow problems with moderately high frequency and with the use of moderately high numbers of time intervals. The ST-LU-SGS implicit scheme is also found to work well in calculating periodic flow problems where the frequency is not known a priori and needed to be determined by using a combined Fourier analysis and gradient-based search algorithm.     

A new class of micro–macro models for elastic–viscoplastic heterogeneous materials

The determination of the effective behavior of heterogeneous materials from the properties of the components and the microstructure constitutes a major task in the design of new materials and the modeling of their mechanical behavior. In real heterogeneous materials, the simultaneous presence of instantaneous mechanisms (elasticity) and time dependent ones (non-linear viscoplasticity) leads to a complex space–time coupling between the mechanical fields, difficult to represent in a simple and efficient way. In this work, a new self-consistent model is proposed, starting from the integral equation for a translated strain rate field. The chosen translated field is the (compatible) viscoplastic strain rate of the (fictitious) viscoplastic heterogeneous medium submitted to a uniform (unknown) boundary condition. The self-consistency condition allows to define these boundary conditions so that a relative simple and compact strain rate concentration equation is obtained. This equation is explained in terms of interactions between an inclusion and a matrix, which lead to interesting conclusions. The model is first applied to the case of two-phase composites with isotropic, linear and incompressible viscoelastic properties. In that case, an exact self-consistent solution using the Laplace–Carson transform is available. The agreement between both approaches appears quite good. Results for elastic–viscoplastic BCC polycrystals are also presented and compared with results obtained from Kröner–Weng's and Paquin et al. (Arch. Appl. Mech. 69 (1999) 14)'s model.     

A phase-field model for quasi-incompressible solid–liquid transitions

According to a unified thermodynamic scheme, we derive the general kinetic equation ruling the phase-field evolution in a binary quasi-incompressible mixture for both transition and separation phenomena. When diffusion effects are negligible in comparison with source and production terms, a solid–liquid phase transition induced by temperature and pressure variations is obtained. In particular, we recover the explicit expression of the liquid–pressure curve separating the solid from the liquid stability regions in the pressure–temperature plane. Consistently with physical evidence, its slope is positive (negative) for substances which compress (expand) during the freezing process.     

A note on the Mooney–Rivlin material model

In finite elasticity, the Mooney–Rivlin material model for the Cauchy stress tensor T in terms of the left Cauchy–Green strain tensor B is given by $$T = -pI + s_1 B + s_2 B^{-1},$$ where p is the pressure and s ₁, s ₂ are two material constants. It is usually assumed that s ₁ >?0 and s ₂?≤ 0, known as E-inequalities, based on the assumption that the free energy function be positive definite for any deformation. In this note, we shall relax this assumption and with a thermodynamic stability analysis, prove that s ₂ need not be negative so that some typical behavior of materials under contraction can also be modeled.     

A void–crack nucleation model for ductile metals

A phenomenological void–crack nucleation model for ductile metals with secondphases is described which is motivated from fracture mechanics and microscale physicalobservations. The void–crack nucleation model is a function of the fracture toughness of theaggregate material, length scale parameter (taken to be the average size of the second phaseparticles in the examples shown in this writing) , the volume fraction of the second phase, strainlevel, and stress state. These parameters are varied to explore their effects upon the nucleationand damage rates. Examples of correlating the void–crack nucleation model to tension data in theliterature illustrate the utility of the model for several ductile metals. Furthermore, compression,tension, and torsion experiments on a cast Al–Si–Mg alloy were conducted to determinevoid–crack nucleation rates under different loading conditions. The nucleation model was thencorrelated to the cast Al–Si–Mg data as well.     

A Variational Approach for the Navier–Stokes System

We introduce a variational approach to treat the regularity of the Navier–Stokes equations both in dimensions 2 and 3. Though the method allows the full treatment in dimension 2, we seek to precisely stress where it breaks down for dimension 3. The basic feature of the procedure is to look directly for strong solutions, by minimizing a suitable error functional that measures the departure of feasible fields from being a solution of the problem. By considering the divergence-free property as part of feasibility, we are able to avoid the explicit analysis of the pressure. Two main points in our analysis are:

A Whitham–Boussinesq long-wave model for variable topography

We study the propagation of water waves in a channel of variable depth using the long-wave asymptotic regime. We use the Hamiltonian formulation of the problem in which the non-local Dirichlet–Neumann operator appears explicitly in the Hamiltonian, and propose a Hamiltonian model for bidirectional wave propagation in shallow water that involves pseudo-differential operators that simplify the variable-depth Dirichlet–Neumann operator. The model generalizes the Boussinesq system, as it includes the exact dispersion relation in the case of constant depth. Analogous models were proposed by Whitham for unidirectional wave propagation. We first present results for the normal modes and eigenfrequencies of the linearized problem. We see that variable depth introduces effects such as a steepening of the normal modes with the increase in depth variation, and a modulation of the normal mode amplitude. Numerical integration also suggests that the constant depth nonlocal Boussinesq model can capture qualitative features of the evolution obtained with higher order approximations of the Dirichlet–Neumann operator. In the case of variable depth we observe that wave-crests have variable speeds that depend on the depth. We also study the evolutions of Stokes waves initial conditions and observe certain oscillations in width of the crest and also some interesting textures and details in the evolution of wave-crests during the passage over obstacles.     

A spectral solver for the Navier–Stokes equations in the velocity–vorticity formulation

A numerical algorithm intended for the study of flows in a cylindrical container under laminar flow conditions is proposed. High resolution of the flow field, governed by the Navier–Stokes equations in velocity–vorticity formulation relative to a cylindrical frame of reference, is achieved through spatial discretisation by means of the spectral method. This method is based on a Fourier expansion in the azimuthal direction and an expansion in Chebyshev polynomials in the (nonperiodic) radial and axial directions. Several regularity constraints are used to take care of the coordinate singularity. These constraints are implemented, together with the boundary conditions at the top, bottom and mantle of the cylinder, via the tau method. The a priori unknown boundary values of the vorticity are evaluated by means of the influence-matrix technique. The compatibility between the mathematical and numerical formulation of the Navier–Stokes equations is established through a tau-correction procedure. The resolved flow field exhibits high-precision satisfaction of the incompressibility constraints for velocity and vorticity and the definition of the vorticity. The performance of the solver is illustrated by resolution of several configurations representative of generic three-dimensional laminar flows.     

A LBM–DEM solver for fast discrete particle simulation of particle–fluid flows

The lattice Boltzmann method (LBM) for simulating fluid phases was coupled with the discrete element method (DEM) for studying solid phases to formulate a novel solver for fast discrete particle simulation (DPS) of particle–fluid flows. The fluid hydrodynamics was obtained by solving LBM equations instead of solving the Navier–Stokes equation by the finite volume method (FVM). Interparticle and particle–wall collisions were determined by DEM. The new DPS solver was validated by simulating a three-dimensional gas–solid bubbling fluidized bed. The new solver was found to yield results faster than its FVM–DEM counterpart, with the increase in the domain-averaged gas volume fraction. Additionally, the scalability of the LBM–DEM DPS solver was superior to that of the FVM–DEM DPS solver in parallel computing. Thus, the LBM–DEM DPS solver is highly suitable for use in simulating dilute and large-scale particle–fluid flows.     

A Weak-L p Prodi–Serrin Type Regularity Criterion for the Navier–Stokes Equations

We give simple proofs that a weak solution u of the Navier–Stokes equations with H ¹ initial data remains strong on the time interval [0, T] if it satisfies the Prodi–Serrin type condition u ∈ L ^s (0, T;L ^r,∞(Ω)) or if its L ^s,∞(0, T;L ^r,∞(Ω)) norm is sufficiently small, where 3 < r ≤ ∞ and (3/r) + (2/s) = 1.