Selection of kernel functions used in a new constitutive equation for viscoelastic fluids

A selection of kernel functions is given to be used in a new integral constitutive equation proposed by Piau whereby the deviatoric stress is calculated from the integral of the history of the past intrinsic rate of rotation and rate of deformation tensors through a representation theorem. Piau has demonstrated the objectivity of a frame moving with a given particle whose axis are directed along the eigenvectors of the rate of deformation tensor. The use of such a framework provides a new approach in the attempt to reduce the computational difficulties associated with conventional constitutive equations written in co-deformational or co-rotational reference frames.The shear and primary normal-stress material functions and the extensional (elongational) stress growth function are defined for the proposed integral constitutive equation. These material functions are used to calculate the kernel functions using steady state, stress relaxation and stress growth data of Attané in simple shear flow for monodisperse polystyrene solutions. The shear and extensional stress growth data of Meissner for a polyethylene melt are also used to show the flexibility of the rheological model.The material functions are first written in terms of five monotonically decreasing functions of the time lag between the past and the present time. Then kernel functions are chosen such that when substituted in the new integral constitutive equation they yield the functions used to describe the data. A further condition imposed on the normalized kernel functions is that they be decreasing functions of time lag.     

Nonlinear Oscillations of Viscoelastic Rectangular Plates

Nonlinear oscillations of viscoelastic simply supported rectangular plates are studied by assuming the Voigt–Kelvin constitutive model. Using Hamilton's principle in conjunction with the kinematics associated with Kirchhoff's plate model, the governing equations of motion including the effect of damping are represented in terms of the transversal deflection and a stress function. Utilizing the Bubnov–Galerkin method, the nonlinear partial differential equations are reduced to an ordinary differential equation which is studied geometrically by approximate construction of the Poincaré maps. Explicit expressions are given for periodic solutions.     

Practical comparison of differential viscoelastic constitutive equations in finite element analysis of planar 4:1 contraction flow

Finite element modeling of planar 4:1 contraction flow (isothermal incompressible and creeping) around a sharp entrance corner is performed for favored differential constitutive equations such as the Maxwell, Leonov, Giesekus, FENE-P, Larson, White-Metzner models and the Phan Thien-Tanner model of exponential and linear types. We have implemented the discrete elastic viscous stress splitting and streamline upwinding algorithms in the basic computational scheme in order to augment stability at high flow rate. For each constitutive model, we have obtained the upper limit of the Deborah number under which numerical convergence is guaranteed. All the computational results are analyzed according to consequences of mathematical analyses for constitutive equations from the viewpoint of stability. It is verified that in general the constitutive equations proven globally stable yield convergent numerical solutions for higher Deborah number flows. Therefore one can get solutions for relatively high Deborah number flows when the Leonov, the Phan Thien-Tanner, or the Giesekus constitutive equation is employed as the viscoelastic field equation. The close relationship of numerical convergence with mathematical stability of the model equations is also clearly demonstrated.     

Characteristic solutions of scalar newton equations in one space dimension

Newton equations are dynamical systems on the space of fields. The solutions of a given equation which are curves of characteristic fields for its force are planar and have constant angular momentum. Separable solutions are characteristic with angular momentum equal to zero. A Newton equation is separable if and only if its characteristic equation is homogeneous. Separable equations correspond to invariants of homogeneous ordinary differential equations, and those associated with a given homogenous equation correspond to its generalized dilation symmetries. A Newton equation is compatible with the characteristic condition if and only if its characteristic equation is linear. Such equations correspond to invariants of linear ordinary differential equations. Those associated with a given linear equation correspond to the central force problems on its solution space. Regardless of compatibility, any Newton equation with a plane of characteristic fields has non-separable characteristic solutions.     

A marker-and-cell approach to viscoelastic free surface flows using the PTT model

This work is concerned with the development of a numerical method capable of simulating two-dimensional viscoelastic free surface flows governed by the non-linear constitutive equation PTT (Phan-Thien–Tanner). In particular, we are interested in flows possessing moving free surfaces. The fluid is modelled by a marker-and-cell type method and employs an accurate representation of the fluid surface. Boundary conditions are described in detail and the full free surface stress conditions are considered. The PTT equation is solved by a high order method which requires the calculation of the extra-stress tensor on the mesh contour. The equations describing the numerical technique are solved by the finite difference method on a staggered grid. In order to validate the numerical method fully developed flow in a two-dimensional channel was simulated and the numerical solutions were compared with known analytic solutions. Convergence results were obtained throughout by using mesh refinement. To demonstrate that complex free surface flows using the PTT model can be computed, extrudate swell and a jet flowing onto a rigid plate were simulated.     

Anomalous predictions of pressure drop and heat transfer in ducts of arbitrary cross-section with modified power law fluids

 The apparent viscosities of purely viscous non-Newtonian fluids are shear rate dependent. At low shear rates, many of such fluids exhibit Newtonian behaviour while at higher shear rates non-Newtonian, power law characteristics exist. Between these two ranges, the fluid's viscous properties are neither Newtonian or power law. Utilizing an apparent viscosity constitutive equation called the “Modified Power Law” which accounts for the above behavior, solutions have been obtained for forced convection flows. A shear rate similarity parameter is identified which specifies both the shear rate range for a given fluid and set of operating conditions and the appropriate solution for that range. The results of numerical solutions for the friction factor–Reynolds number product and for the Nusselt number as a function of a dimensionless shear rate parameter have been presented for forced fully developed laminer duct flows of different cross-sections with modified power law fluids. Experimental data is also presented showing the suitability of the “Modified Power Law” constitutive equation to represent the apparent viscosity of various polymer solutions. Received on 21 August 2000     

A phenomenological and extended continuum approach for modelling non-equilibrium flows

This paper presents a new technique that combines Grad’s 13-moment equations (G13) with a phenomenological approach to rarefied gas flows. This combination and the proposed solution technique capture some important non-equilibrium phenomena that appear in the early continuum-transition flow regime. In contrast to the fully coupled 13-moment equation set, a significant advantage of the present solution technique is that it does not require extra boundary conditions explicitly; Grad’s equations for viscous stress and heat flux are used as constitutive relations for the conservation equations instead of being solved as equations of transport. The relative computational cost of this novel technique is low in comparison to other methods, such as fully coupled solutions involving many moments or discrete methods. In this study, the proposed numerical procedure is tested on a planar Couette flow case, and the results are compared to predictions obtained from the direct simulation Monte Carlo method. This test case highlights the presence of normal viscous stresses and tangential heat fluxes that arise from non-equilibrium phenomena, which cannot be captured by the Navier–Stokes–Fourier constitutive equations or phenomenological modifications.      

Decoupled second-order transient schemes for the flow of viscoelastic fluids without a viscous solvent contribution

The simulation of transient flows is relevant in several applications involving viscoelastic fluids. In the last decades, much effort has been spent on deriving time-marching schemes able to efficiently solve the governing equations at low computational cost. In this direction, decoupling schemes, where the global system is split into smaller subsystems, have been particularly successful. However, most of these techniques only work if inertia and/or a large Newtonian solvent contribution is included in the modeling. This is not the case for polymer melts or concentrated polymer solutions.In this work, we propose two second-order time-integration schemes for discretizing the momentum balance as well as the constitutive equation, based on a Gear and a Crank–Nicolson scheme. The solution of the momentum and continuity equations is decoupled from the constitutive one. The stress tensor term in the momentum balance is replaced by its space-continuous but time-discretized form of the constitutive equation through an Euler scheme implicit in the velocity. This adds velocity unknowns in the momentum equation thus an updating of the velocity field is possible even if inertia and solvent viscosity are not included in the model. To further reduce computational costs, the non-linear relaxation term in the constitutive equation is taken explicitly leading to a linear system of equations for each stress component.Four benchmark problems are considered to test the numerical schemes. The results show that a Crank–Nicolson based discretization for the momentum equation produces oscillations when combined with a Crank–Nicolson based scheme for the constitutive equation whereas, if a Gear based scheme is implemented for the constitutive equation, the stability is found to be dependent on the specific problem. However, the Gear based scheme applied to the momentum balance combined with both second-order methods used for the constitutive equation is stable and accurate and performs much better than a first-order Euler scheme. Finally, a numerical proof of the second-order convergence is also carried out.     

Vortex motions of liquid in a narrow channel

Quasi-linear integrodifferential equations that describe vortex flows of an ideal incomparessible liquid in a narrow curved channel in the Eulerian-Lagrangian coordinate system are considered. The necessary and sufficient conditions for hyperbolicity of the system of equations of motion are obtained for flows with a monotonic velocity depth profile. The propagation velocities of the characteristics and the characteristic form of the system are calculated. A particular solution is given in which the system of integrodifferential equations changes type with time. The solution of the Cauchy problem is given for linearized equations. An example of initial data for which the Cauchy problem is ill-posed is constructed. Novosibirsk State University, Novosibirsk 630090. Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 39, No. 4, pp. 38–49, July–August, 1998.     

A generalized variational formulation for convective heat transfer

The Scope of this paper is to develop the basic equations for a variational formulation which can be used to solve problems related to convection and/or diffusion dominated flows. The formulation is based on the introduction of a generalized quantity defined as the hear displacement. The governing equation is expressed in terms of this quantity and a variational formulation is developed which leads to a system of equations similar in form to Lagrange's equations of mechanics. These equations can be used for obtaining approximate solutions, though they are of particular interest for application of the finite element method. As an example of the formulation two finite element models are derived for solving convectiondiffusion boundary value problems. The performance of the two models is investigated and numerical results are given for different cases of convection and diffusion with two types of boundary conditions. The applications of the developed formulations are not limited to convection-diffusion problems but can also be applied to other types of problems such as mass transfer, hydrodynamics and wave propagation.     

A method for implementing time-integral constitutive equations in commercial CFD packages

A method is described for obtaining viscoelastic flow solutions, based upon time-integral constitutive equations, using a general purpose CFD package. The method is general enough to be applied to any available software that has rudimentary input and output facilities and can solve a Stoke's flow problem. From this basis, flexibility of choice of constitutive equation and computational techniques is available. The method is presented in a form appropriate for solving both planar and axisymmetric flows. Delaunay triangulation is used to reconstruct a mesh for external code, and stress computation procedures are performed on this mesh. The method has only two particular requirements for the CFD package used – it must be able to output nodal values (of position and velocity) to file and it must be able to read body-forces from a file. Two methods of velocity adjustment were compared: an incremental method and a method whereby the viscoelastic stresses were incorporated directly as body-forces. Results and convergence from the two methods were found to be essentially identical, hence the direct body-force method (which is considerably easier to implement) is described in detail. The method is applied to a well-known flow problem of LDPE melt through a 4 : 1 abrupt contraction axisymmetric die. Convergence was obtained up to nearly the highest value of apparent shear rate for which published simulation results are available. Quantitative results for vortex strength, vortex opening angle and Couette correction are presented which are compared with earlier work on the problem using other methods. Agreement is generally good, giving confidence in the method. Simulations of planar flows of the same melt are performed: a decreasing corner vortex was observed. This phenomenon has been observed experimentally for flows of other substances, but is not expected for flows of LDPE melt. A parametric study of a critical strain hardening parameter is conducted to help explain the cause of the results.     

Plane Strain Bending of Cylindrical Sectors of Admissible Compressible Hyperelastic Materials

In the theory of nonlinear elasticity of rubber-like materials, if a homogeneous isotropic compressible material is described by a strain–energy function that is a homogeneous function of the principal stretches, then the equations of equilibrium for axisymmetric deformations reduce to a separable first-order ordinary differential equation. For a particular class of such strain–energy functions, this property is used to obtain a general parametric solution to the equilibrium equation for plane strain bending of cylindrical sectors. Specification of the arbitrary function that appears in such strain–energy functions yields some parametric solutions. In some cases, the parameter can be eliminated to yield closed-form solutions in implicit or explicit form. Other possible forms for the arbitrary constitutive function that are likely to yield such solutions are also indicated.     

Comparison of experimental data with the numerical simulation of planar entry flow: Role of the constitutive equation

The goal of this research was to determine whether there is any interaction between the type of constitutive equation used and the degree of mesh refinement, as well as how the type of constitutive equation might affect the convergence and quality of the solution, for a planar 4:1 contraction in the finite eiement method. Five constitutive equations were used in this work: the Phan-Thien–Tanner (PTT), Johnson–Segalman (JS), White–Metzner (WM), Leonov-like and upper convected Maxwell (UCM) models. A penalty Galerkin finite element technique was used to solve the system of non-linear differential equations. The constitutive equations were fitted to the steady shear viscosity and normal stress data for a polystyrene melt. In general it was found that the convergence limit based on the Deborah number De and the Weissenberg number We varied from model to model and from mesh to mesh. From a practical point of view it was observed that the wall shear stress in the downstream region should also be indicated at the point where convergence is lost, since this parameter reflects the throughput conditions. Because of the dependence of convergence on the combination of mesh size and constitutive equation, predictions of the computations were compared with birefringence data obtained for the same polystyrene melt flowing through a 4:1 planar contraction. Refinement in the mesh led to better agreement between the predictions using the PTT model and flow birefringence, but the oscillations became worse in the corner region as the mesh was further refined, eventually leading to the loss of convergence of the numerical algorithm. In comparing results using different models at the same wall shear stress conditions and on the same mesh, it was found that the PTT model gave less overshoot of the stresses at the re-entrant corner. Away from the corner there were very small differences between the quality of the solutions obtained using different models. All the models predicted solutions with oscillations. However, the values of the solutions oscillated around the experimental birefringence data, even when the numerical algorithm would not converge. Whereas the stresses are predicted to oscillate, the streamlines and velocity field remained smooth. Predictions for the existence of vortices as well as for the entrance pressure loss (ΔP_ent) varied from model to model. The UCM and WM models predicted negative values for ΔP_ent.     

Computation of Viscoelastic flow in a cylindrical tank with a rotating lid

An orthogonal collocation method is used to compute steady flows of viscoelastic fluids in a cylindrical tank covered by a rotating disk. The stability of these flows with respect to small disturbances is also analyzed. The Criminale-Ericksen-Filbey constitutive equation is used, since the primary flow is viscometric and the secondary flow is small. The solutions satisfy the equation of continuity and the boundary conditions exactly, including the velocity discontinuity at the edge of the disk. The computed flows for aqueous solutions of Separan 30 exhibit single or double vortices, according to the concentration and the rotation speed. Reasonable agreement is found with the data of Hill [7,8].     

Development of an artificial compressibility methodology using flux vector splitting

An implicit, upwind arithmetic scheme that is efficient for the solution of laminar, steady, incompressible, two-dimensional flow fields in a generalised co-ordinate system is presented in this paper. The developed algorithm is based on the extended flux-vector-splitting (FVS) method for solving incompressible flow fields. As in the case of compressible flows, the FVS method consists of the decomposition of the convective fluxes into positive and negative parts that transmit information from the upstream and downstream flow field respectively. The extension of this method to the solution of incompressible flows is achieved by the method of artificial compressibility, whereby an artificial time derivative of the pressure is added to the continuity equation. In this way the incompressible equations take on a hyperbolic character with pseudopressure waves propagating with finite speed. In such problems the ‘information’ inside the field is transmitted along its characteristic curves. In this sense, we can use upwind schemes to represent the finite volume scheme of the problem's governing equations. For the representation of the problem variables at the cell faces, upwind schemes up to third order of accuracy are used, while for the development of a time-iterative procedure a first-order-accurate Euler backward-time difference scheme is used and a second-order central differencing for the shear stresses is presented. The discretized Navier–Stokes equations are solved by an implicit unfactored method using Newton iterations and Gauss–Siedel relaxation. To validate the derived arithmetical results against experimental data and other numerical solutions, various laminar flows with known behaviour from the literature are examined. © 1997 John Wiley & Sons, Ltd.     

Construction of approximate solutions to two-dimensional potential problems of electrohydrodynamics

After transformation to new variables the system of equations describing planar potential electrohydrodynamic flows with a small interaction parameter is converted to a single equation. The particular solution of this equation, which is the electrohydrodynamic analog of Hamel's solution in the dynamics of a viscous liquid, is found. Two types of flows, described by simplified equations, can be distinguished when certain constraints are imposed on the manner in which the electrical parameters vary along the coordinate lines and the terms of the equation correspondingly estimated. These flows are the jet and quasione-dimensional flows of the charged component in a curvilinear electrostatic field and a supper-posed two-dimensional potential flow of the carrier medium. Solutions to the approximate equations are obtained for certain particular cases.Kiev. Transklated from Izvestiya Akademii Nauk, Mekhanika Zhidkosti i Gaza, No. 4, pp. 140–147, July–August, 1972.     

粘弹性固体的精细积分有限元算法

粘弹性固体本构方程的数学表达式分为微分型和积分型两种，其数值求解主要是时域上离散计算。文中从微分型表达式出发导出其状态空间方程的数学表达式，通过严格推导论证了它与微、积分型表达式的等价性；引入状态空间方程，从而利用精细积分格式来求解粘弹性固体本构方程；给出了粘弹性固体本构方程的精细积分有限元算法，为求解粘弹性固体本构方程的数值解提供了一个新的途径，具有计算简便，求解精度高等优点。     

Some properties of three-dimensional Biltrami flows

The investigation of Beltrami flows is important for the research on the mechanism of turbulent structure. In this paper the general solutions of the Beltrami flows are given, which depend explicitly on the solutions of three independent Helmholtz equations with scalar unknowns. Velocity fields of Beltrami flows can then be obtained explicitly after the application of some curl operations on the solutions of Helmholtz equations. On the basis of the exact solutions of Euler equations given above, we obtain one kind of exact solutions of non-steady Navier-Stokes equations which are also the Beltrami flows. Some interesting examples of Beltrami flows other than “ABC flows”, “Kolmogolov flows”, “Rayleigh-Bernard flows”, “Q-flows” are given. The detailed analytic results of these examples will be published in the near future.