The \vec\omega-\vec{D} fluid: general theory with special emphasis on stationary two dimensional flows

In this paper a theory is presented in which the extra stress tensor is allowed to depend not only on the rate of strain tensor but also on the relative vorticity of the fluid, i.e. on the vorticity relative to the rate of rotation of the principal straining directions. This theory has its origin in an expansion of in terms of kinematic tensors in the limit of stationarity in a material sense (constant stretch history flows). For two dimensional flows of an incompressible fluid three tensors suffice to completely specify . The three material functions which appear can depend only on two invariants, namely the second invariant of and on . Using the predictions of an Oldroyd 8 constant fluid in a homogeneous planar flow of constant stretch history, the three material functions are studied in detail. For the special case of a quasi-Newtonian fluid shear thinning and extension thickening can directly be accounted for in the “viscosity” function. Received: September 26, 1996     

A simple integral constitutive equation for polymeric liquids

An integral constitutive equation describing the strain history in terms of the Cauchy-Green and Finger tensors is presented. Its memory functions, which correspond to a flow-dependent continuous spectrum of relaxation times, depend on the invariants of the Cauchy-Green and Finger rate of strain tensors. Theoretical predictions resulting from this equation are compared with experimental data for some types of flow. Using the Laguerre polynomial expansion, the basic equation is generalized.     

Non-swirling axisymmetric flow of the ω-D fluid

In this paper a constitutive equation for the extra stress tensor τ is considered, which can be used for steady axisymmetric flows when u _ϕ=0 (non-swirling). It is explicit with coefficients which, in case of incompressibility, depend only upon three invariants. As such it should prove useful in numerical calculations, especially so if used in its simplest form, namely the quasi-Newtonian fluid. Here, a purely viscous code can be used. The constitutive equation is self-consistent and receives additional justification from the fact that any constitutive equation is bound to result in the form proposed if the flow is a constant stretch history flow. Received: 30 March 1998 Accepted: 6 August 1998     

Generalized Darcy's Law and the Structure of the Phase and Relative Phase Permeabilities for Two-Phase Flows through Anisotropic Porous Media

For two-phase immiscible fluid flows a generalized Darcy's law is written in invariant tensor form for crystallographic point symmetry groups and anisotropic textures. The representation of the phase permeability coefficient tensors and the structure of the expressions for the relative phase permeabilities are analyzed for all symmetry groups. The relation between the phase and absolute permeability coefficient tensors is specified by a fourth-rank tensor with the external symmetry coinciding with external symmetry of the phase permeability tensors. It is shown that the external symmetry of the phase permeability coefficient tensors can differ from the external symmetry of the absolute permeability tensor. For triclinic and monoclinic symmetry groups it is shown that the phase permeability coefficient tensors may not be coaxial with each other and with the absolute permeability tensor; moreover, the directions of the principal axes of the phase permeability coefficient tensors can depend on the saturation.     

Laws of flow with a limiting gradient in anisotropic porous media

The equations of viscoplastic fluid flow through a porous medium are written for all types of anisotropy. It is shown that in anisotropic media the flows with a limiting gradient are characterized by two material tensors: the tensor of permeability (flow resistance) coefficients and the tensor of limiting gradients. A complex of laboratory measurements for determining the tensors of permeability coefficients and limiting gradients is considered for all types of anisotropic media. It is shown that the tensors of permeability coefficients and limiting gradients are coaxial. Conditions of flow onset and fluid flow laws are formulated for media with monoclinic and triclinic symmetries of flow characteristics.     

自旋张量的绝对表示及其在有限变形理论中的应用

基于对一类线性张量方程的一般解法,导出了任一对称张量所对应的自旋张量的绝对表示。该结果可以很自然地用于研究左和右伸长张量的自旋并研讨在连续介质力学中常见到的各种转动率张量间的关系。一个重要的公式,即Hill意义下广义应变的共轭应力和Cauchy应力之间的关系,从功共轭原理建立了起来。尤其是详细讨论了对数应变的时间变率及相应的共轭应力。无疑,上述结果对有限变形条件下本构理论的研究是颇为重要的。     

A new numerical algorithm for viscoelastic fluid flows : The grid-by-grid inversion method

We present a new algorithm for solving viscoelastic flows with a general constitutive equation. In our approach the hyperbolic constitutive equation is split such that the term for the convective transport of stress tensor is treated as a source. This allows the stress tensor at each grid point to be expressed mainly in terms of the velocity gradient tensor at the same point. Then, the set of six stress tensor components is found after inverting a six by six matrix at each grid point. Thus we call this algorithm the grid-by-grid inversion method. The convective transport of stress tensor in the constitutive equation, which has been treated as a source, is updated iteratively. The present algorithm can be combined with finite volume method, finite element method or the spectral methods. To corroborate the accuracy and robustness of the present algorithm we consider viscoelastic flow past a cylinder placed at the center between two plates, which has served as a benchmark problem. Also considered is the investigation of the pattern and strength of the secondary flows in the viscoelastic flows through a rectangular pipe. It is found that the present method yields accurate results even for large relaxation times.     

The Giesekus fluid in ω–D form for steady two-dimensional flows

Using the fact that for simple fluids the most general constitutive equation in constant stretch history flows for the extra stress tensor τ is known in an explicit form, the Giesekus fluid model is cast into this (ω–D) form for two-dimensional flows. The three material functions needed to characterize τ are listed. The explicit results for simple shear and planar elongation reveal that the parameter α should be restricted to values less than 0.5. It is demonstrated that in this explicit form the constitutive equation is free from thermodynamic objections and can thus be used as a starting point for numerical calculations of general, but steady, two-dimensional flows. Received: 9 November 1998 Accepted: 20 May 1999     

A priori evaluations and least-squares optimizations of turbulence models for fully developed rotating turbulent channel flow

The present study involves a priori tests of pressure-strain and dissipation rate tensor models using data from direct numerical simulations (DNS) of fully developed turbulent channel flow with and without spanwise system rotation. Three different pressure-strain rate models are tested ranging from a simple quasi-linear model to a realizable fourth order model. The evaluations demonstrate the difficulties of developing RANS-models that accurately describe the flow for a wide range of rotation numbers. Furthermore, least-squares based tensor representations of the exact pressure-strain and dissipation rate tensors are derived pointwise in space. The relation obtained for the rapid pressure-strain rate is exact for general 2D mean flows. Hence, the corresponding distribution of the optimized coefficients show the ideal behaviour. The corresponding representations for the slow pressure-strain and dissipation rate tensors are incomplete but still optimal in a least-squares sense. On basis of the least-squares analysis it is argued that the part of the representation that is tensorially linear in the Reynolds stress anisotropy is the most important for these parts.     

The invariant representation of spins with applications in the theory of finite deformations

Based on the general solution given to a kind of linear tensor equations, the spin of a symmetric tensor is derived in an invariant form. The result is applied to find the spins of the left and the right stretch tensors and the relation among different rotation rate tensors has been discussed. According to work conjugacy, the relations between Cauchy stress and the stresses conjugate to Hill's generalized strains are obtained. Particularly, the logarithmic strain, its time rate and the conjugate stress have been discussed in detail. These results are important in modeling the constitutive relations for finite deformations in continuum mechanics. The project is supported by the National Natural Science Foundation of China and the Chinese Academy of Sciences (No. 87-52).     

Invariant Properties for Finding Distance in Space of?Elasticity Tensors

Using orthogonal projections, we investigate distance of a given elasticity tensor to classes of elasticity tensors exhibiting particular material symmetries. These projections depend on the orientation of the elasticity tensor; hence the distance is obtained as the minimization of corresponding expressions with respect to the action of the orthogonal group. These expressions are stated in terms of the eigenvalues of both the given tensor and the projected one. The process of minimization is facilitated by the fact that, as we prove, the traces of the corresponding Voigt and dilatation tensors are invariant under these orthogonal projections. For isotropy, cubic symmetry and transverse isotropy, we formulate algorithms to find both the orientation and the eigenvalues of the elasticity tensor endowed with a particular symmetry and closest to the given elasticity tensor.      

On the consequences of material frame-indifference in algebraic stress models

The principle of material frame-indifference (MFI) is a fundamental and controversial principle of continuum mechanics that has been invoked recently to derive nonlinear algebraic models for stresses of viscoelastic liquids. The purpose of the present study is to identify regions of a flow field where MFI should be considered. Such regions are identified by computing the angular velocity of the principal directions of the rate-of-deformation tensor in order to obtain an Euclidean objective vorticity tensor. An analysis is carried out for uniform shear and extensional flows, and for a Couette flow. The method is then applied to the planar flow through an abrupt 4:1 contraction and to the two-dimensional stream past a circular cylinder. The main results are: (1) MFI should be taken into account in regions characterized by the transition between two different kinematics and a significant velocity magnitude, and (2) MFI can be safely ignored in regions of pure viscometric behaviour as well as in recirculation regions. The consequences of MFI being taken into account are then examined upon using the Euclidean objective vorticity tensor in a simple algebraic constitutive law for viscoelastic fluids.     

Stress and strain-rate formulations for fabric evolution in polar ice

Re-orientation of individual crystal glide planes, as isotropic surface ice is deformed during its passage to depth in an ice sheet, creates a fabric and associated anisotropy. We re-examine an orthotropic viscous law which was developed to reflect the induced anisotropy arising from the mean rotation of crystal axes during deformation. This expresses the deviatoric stress, the stress formulation, in terms of the strain-rate, strain, and three structure tensors based on the principal stretch axes, and involves two fabric response coefficient functions which determine the strength of the anisotropy. A validity condition implicitly relates the two response functions, so the model law has only one independent fabric response function. A modified formulation is now presented in which the two fabric response coefficients are expressed as functions of different invariant arguments, and the validity condition becomes an explicit algebraic relation between the two functions. The response can therefore be described explicitly in terms of a single fabric response function. An analogous orthotropic viscous law for the strain-rate, the strain-rate formulation, akin to the conventional “flow law” for isotropic ice, expressed in terms of the deviatoric stresss, strain and the three structure tensors, is also constructed. Correlations with complete (idealised) uni-axial compression and shearing responses are made for the stress formulation, to determine the fabric response function which would yield these responses. Received January 30, 2002 / Published online October 15, 2002 RID="*" ID="*" On leave from the Institute of Hydroengineering, Polish Academy of Sciences, ul. Waryńskiego 17, 71-310 Szczecin, Poland Communicated by Kolumban Hutter, Darmstadt     

Constitutive modeling of temperature and strain rate dependent elastoplastic hardening materials using a corotational rate associated with the plastic deformation

In this paper, a constitutive model with a temperature and strain rate dependent flow stress (Bergstrom hardening rule) and modified Armstrong-Frederick kinematic evolution equation for elastoplastic hardening materials is introduced. Based on the multiplicative decomposition of the deformation gradient,new kinematic relations for the elastic and plastic left stretch tensors as well as the plastic deformation-dependent spin tensor are proposed. Also, a closed-form solution has been obtained for the elastic and plastic left stretch tensors for the simple shear problem.To evaluate model validity, results are compared with known experimental data for SUS 304 stainless steel, which shows a good agreement with the results of the proposed theoretical model.Finally, the stress-deformation curve, as predicted by the model, is plotted for the simple shear problem at room and elevated temperatures using the same material properties for AA5754-O aluminium alloy.     

Equivalent Diffusion in a Thin Film Flow with a Non-One-Dimensional Velocity Field and Anisotropic Diffusion Tensor

For describing the mass transfer processes in channels, Taylor's dispersion theory is widely used. This theory makes it possible, with asymptotic rigor, to replace the complete diffusion (heat conduction) equation with a convective term that depends on the coordinate transverse to the flow by an effective diffusion (dispersion) equation with constant coefficients, averaged over the channel cross-section. In numerous subsequent studies, Taylor's theory was generalized to include more complex situations, and novel algorithms for constructing the dispersion equations were proposed. For thin film flows a theory similar to Taylor's leads to a matrix of dispersion coefficients.In this study, Taylor's theory is extended to film flows with a non-one-dimensional velocity field and anisotropic diffusion tensor. These characteristics also depend to a considerable extent on the spatial coordinates and time. The dispersion equations obtained can be simplified in regions in which the effective diffusion coefficient tensor changes sharply.     

An ‘extended’ volumetric/deviatoric formulation of anisotropic damage based on a pseudo-log rate

Following a framework of elastic degradation and damage previously proposed by the authors, an ‘extended’ formulation of orthotropic damage in initially isotropic materials, based on volumetric/deviatoric decomposition, is presented. The formulation is founded on the concept of energy equivalence and makes use of second-order symmetric tensor damage variables. It is characterized by fourth-order damage-effect tensors (relating nominal to effective stresses and strains) built from the underlying second-order damage tensors and decomposed in product-form in isotropic and anisotropic parts. The formulation is developed in two steps. First, secant relations are established. In the isotropic case, the model embeds a path parameter allowing to range between pure volumetric to pure deviatoric damage. With the two undamaged material constants this makes a total of three constant parameters plus an evolving scalar damage variable, giving rise to a four-parameter model with two varying isotropic material coefficients. In the anisotropic case, the model is still characterized by the same three material constants plus three evolving variables which are the principal values of a second-order damage tensor. This leads to a six-parameter restricted form of orthotropic damage. In the second step, damage evolution rules are formulated in terms of a pseudo-logarithmic rate of damage. This allows to define meaningful conjugate forces that constitute a feasible space in which loading functions and damage evolution rules can be defined. The present ‘extended’ formulation is closed by the derivation of the tangent stiffness.     

n the fourth order tensor valued function of the stress in return map algorithm

针对各向同性材料,基于一组相互正交的基张量,建立了一套有 效的相关运算方法. 基张量中的两个分别是归一化的二阶单位张量和偏应力张量,另一个则 使用应力的各向同性二阶张量值函数经过归一化构造所得,三者共主轴. 根据张量函数表示 定理,本构方程和返回映射算法中所涉及到的应力的二阶、四阶张量值函数及其逆都由这组 基所表示. 推演结果表明：这些张量之间的运算,表现为对应系数矩阵之间的简单 关系. 其中,四阶张量求逆归结为对应的3\times3系数矩阵求逆,它对二阶张量的变换 则表现为该矩阵对3times 1列阵的变换. 最后,对这些变换关系应用于返回映 射算法的迭代格式进行了相关讨论.     

General invariant representations of the constitutive equations for isotropic nonlinearly elastic materials

This paper develops general invariant representations of the constitutive equations for isotropic nonlinearly elastic materials. Different sets of mutually orthogonal unit tensor bases are constructed from the strain argument tensor by using the representation theorem and corresponding irreducible invariants are defined. Their relations and geometrical interpretations are established in three dimensional principal space. It is shown that the constitutive law linking the stress and strain tensors is revealed to be a simple relationship between two vectors in the principal space. Relative to two different sets of the basis tensors, the constitutive equations are transformed according to the transformation rule of vectors. When a potential function is assumed to exist, the vector associated with the stress tensor is expressed in terms of its gradient with respect to the vector associated with the strain tensor. The Hill’s stability condition is shown to be that the scalar product of the increment of those two vectors must be positive. When potential function exists, it becomes to be that the 3 × 3 constitutive matrix derived from its second order derivative with respect to the vector associated with the strain must be positive definite. By decomposing the second order symmetric tensor space into the direct sum of a coaxial tensor subspace and another one orthogonal to it, the closed form representations for the fourth order tangent operator and its inversion are derived in an extremely simple way.     

Flow of a quasi-Newtonian fluid through a planar contraction

Flows through abrupt contractions are dominated by the rapid extension experienced in passing through the contraction. Thus, it is useful to employ a fluid model which considers the extensional viscosity explicitly in its constitutive equation. In this paper, the quasi-Newtonian fluid model, which admits shear thinning and extension thickening of the viscosity depending on the local type of flow as proposed by Schunk and Scriven [P. Schunk, L. Scriven, J, Rheol 34 (1990) 1085], is applied to the numerical simulation of the flow of a dilute polyacrylamide solution through a planar 4 : 1 contraction. In this theory the extra stress tensor does not only depend on the rate of strain tensor but also on the relative rate of rotation of the fluid. The material function – the viscosity function – is allowed to depend on the invariants of these two kinematic tensors yielding a local distinction between extensional, shear or rotation dominated flow. The governing equations are discretized using a finite volume method. Different model parameters are varied and the simulation results are compared with the generalized Newtonian fluid and experimental data.