Positive‐definite q‐families of continuous subcell Darcy‐flux CVD(MPFA) finite‐volume schemes and the mixed finite element method

A new family of locally conservative cell‐centred flux‐continuous schemes is presented for solving the porous media general‐tensor pressure equation. A general geometry‐permeability tensor approximation is introduced that is piecewise constant over the subcells of the control volumes and ensures that the local discrete general tensor is elliptic. A family of control‐volume distributed subcell flux‐continuous schemes are defined in terms of the quadrature parametrization q (Multigrid Methods. Birkhauser: Basel, 1993; Proceedings of the 4th European Conference on the Mathematics of Oil Recovery, Norway, June 1994; Comput. Geosci. 1998; 2 :259–290), where the local position of flux continuity defines the quadrature point and each particular scheme. The subcell tensor approximation ensures that a symmetric positive‐definite (SPD) discretization matrix is obtained for the base member (q=1) of the formulation. The physical‐space schemes are shown to be non‐symmetric for general quadrilateral cells. Conditions for discrete ellipticity of the non‐symmetric schemes are derived with respect to the local symmetric part of the tensor. The relationship with the mixed finite element method is given for both the physical‐space and subcell‐space q‐families of schemes. M‐matrix monotonicity conditions for these schemes are summarized. A numerical convergence study of the schemes shows that while the physical‐space schemes are the most accurate, the subcell tensor approximation reduces solution errors when compared with earlier cell‐wise constant tensor schemes and that subcell tensor approximation using the control‐volume face geometry yields the best SPD scheme results. A particular quadrature point is found to improve numerical convergence of the subcell schemes for the cases tested. Copyright © 2007 John Wiley & Sons, Ltd.     

A least squares finite element method for viscoelastic fluid flow problems

Viscoelastic flows remain a demanding class of problems for approximate analysis, particularly at increasing Weissenberg numbers. Part of the difficulty stems from the convective behavior and in the treatment of the stress field as a primary unknown. This latter aspect has led to the use of higher-order piecewise approximations for the stress approximation spaces in recent finite element research. The computational complexity of the discretized problem is increased significantly by this approach but at present it appears the most viable technique for solving these problems. Motivated by recent success in treating mixed systems and convective problems, we formulate here a least squares finite element method for the viscoelastic flow problem. Numerical experiments are conducted to test the method and examine its strengths and limitations. Some difficulties and open issues are identified through the numerical experiments. We consider the use of high degree elements (p refinement) to improve performance and accuracy.     

An approximate solution for the Couette–Poiseuille flow of the Giesekus model between parallel plates

An approximate analytical solution is derived for the Couette–Poiseuille flow of a nonlinear viscoelastic fluid obeying the Giesekus constitutive equation between parallel plates for the case where the upper plate moves at constant velocity, and the lower one is at rest. Validity of this approximation is examined by comparison to the exact solution during a parametric study. The influence of Deborah number (De) and Giesekus model parameter (α) on the velocity profile, normal stress, and friction factor are investigated. Results show strong effects of viscoelastic parameters on velocity profile and normal stress. In addition, five velocity profile types were obtained for different values of α, De, and the dimensionless pressure gradient (G).     

Resonant waves in elastic structured media: Dynamic homogenisation versus Green’s functions

We address an important issue of dynamic homogenisation in vector elasticity for a doubly periodic mass-spring elastic lattice. The notion of logarithmically growing resonant waves is used in the analysis of star-shaped wave forms induced by an oscillating point force. We note that the dispersion surfaces for Floquet–Bloch waves in the elastic lattice may contain critical points of the saddle type. Based on the local quadratic approximations of a dispersion surface, where the radian frequency is considered as a function of wave vector components, we deduce properties of a transient asymptotic solution associated with the contribution of the point source to the wave form. The notion of local Green’s functions is used to describe localised wave forms corresponding to the resonant frequency. The special feature of the problem is that, at the same resonant frequency, the Taylor quadratic approximations for different groups of the critical points on the dispersion surfaces (and hence different Floquet–Bloch vectors) are different. Thus, it is shown that for the vector case of micro-structured elastic medium there is no uniformly defined dynamic homogenisation procedure for a given resonant frequency. Instead, the continuous approximation of the wave field can be obtained through the asymptotic analysis of the lattice Green’s functions, presented in this paper.     

High-order DG solvers for underresolved turbulent incompressible flows: A comparison of L2 and H(div) methods

The accurate numerical simulation of turbulent incompressible flows is a challenging topic in computational fluid dynamics. For discretisation methods to be robust in the underresolved regime, mass conservation and energy stability are key ingredients to obtain robust and accurate discretisations. Recently, two approaches have been proposed in the context of high-order discontinuous Galerkin (DG) discretisations that address these aspects differently. On the one hand, standard L²-based DG discretisations enforce mass conservation and energy stability weakly by the use of additional stabilisation terms. On the other hand, pointwise divergence-free H(div)-conforming approaches ensure exact mass conservation and energy stability by the use of tailored finite element function spaces. This work raises the question whether and to which extent these two approaches are equivalent when applied to underresolved turbulent flows. This comparative study highlights similarities and differences of these two approaches. The numerical results emphasise that both discretisation strategies are promising for underresolved simulations of turbulent flows due to their inherent dissipation mechanisms.     

A hybrid stress approach for laminated composite plates within the First-order Shear Deformation Theory

This paper presents a hybrid stress approach for the analysis of laminated composite plates. The plate mechanical model is based on the so called First-order Shear Deformation Theory, rationally deduced from the parent three-dimensional theory. Within this framework, a new quadrilateral four-node finite element is developed from a hybrid stress formulation involving, as primary variables, compatible displacements and elementwise equilibrated stress resultants. The element is designed to be simple, stable and locking-free. The displacement interpolation is enhanced by linking the transverse displacement to the nodal rotations and a suitable approximation for stress resultants is selected, ruled by the minimum number of parameters. The transverse stresses through the laminate thickness are reconstructed a posteriori by simply using three-dimensional equilibrium. To improve the results, the stress resultants entering the reconstruction process are first recovered using a superconvergent patch-based procedure called Recovery by Compatibility in Patches, that is properly extended here for laminated plates. This preliminary recovery is very efficient from the computational point of view and generally useful either to accurately evaluate the stress resultants or to estimate the discretization error. Indeed, in the present context, it plays also a key role in effectively predicting the shear stress profiles, since it guarantees the global convergence of the whole reconstruction strategy, that does not need any correction to accommodate equilibrium defects. Actually, this strategy can be adopted together with any plate finite element. Numerical testing demonstrates the excellent performance of both the finite element and the reconstruction strategy.     

Stability of some low‐order approximations for the Stokes problem

Two‐level low‐order finite element approximations are considered for the inhomogeneous Stokes equations. The elements introduced are attractive because of their simplicity and computational efficiency. In this paper, the stability of a Q₁(h)–Q₁(2h) approximation is analysed for general geometries. Using the macroelement technique, we prove the stability condition for both two‐ and three‐dimensional problems. As a result, optimal rates of convergence are found for the velocity and pressure approximations. Numerical results for three test problems are presented. We observe that for the computed examples, the accuracy of the two‐level bilinear approximation is compared favourably with some standard finite elements. Copyright © 2007 John Wiley & Sons, Ltd.     

结构参数大修改时的特征值重分析方法

就结构参数发生大修改的情况提出了两种高精度的特征值重分析方法：Ｐａｄｅ　逼近法和推广的Ｋｉｒｓｃｈ混合法．利用这两种方法,计算了一个具有２０２个结点,３５７个梁单元的平面框架的近似特征值．计算结果表明,所提出的方法是结构参数修改时的特征值重分析的有效方法．     

Generation of the second,fourth, eighth,and subsequent harmonics by a quadratic nonlinear hyperelastic longitudinal plane wave

The perturbation (small-parameter) method is used to obtain the first three approximations for the problem of a harmonic longitudinal plane wave propagating in a quadratic nonlinear hyperelastic material described by the classical Murnaghan model. The subsequent approximations are discussed. The contribution of each approximation to the overall wave pattern is analyzed. It is shown that the third approximation corrects the prediction of the evolution of the initial wave profile. Which of the harmonics dominates depends on the distance traveled by the wave: the second harmonic is generated first, then it transforms into the fourth harmonic, and finally, as the distance increases, the eighth harmonic shows     

On the numerical stability of mixed finite-element methods for viscoelastic flows governed by differential constitutive equations

Mixed finite-element methods for computation of viscoelastic flows governed by differential constitutive equations vary by the polynomial approximations used for the velocity, pressure, and stress fields, and by the weighted residual methods used to discretize the momentum, continuity, and constitutive equations. This paper focuses on computation of the linear stability of the planar Couette flow as a test of the numerical stability for solution of the upper-convected Maxwell model. Previous theoretical results prove this inertialess flow to be always stable, but that accurate calculation is difficult at high De because eigenvalues with fine spatial structure and high temporal frequency approach neutral stability. Computations with the much used biquadratic finite-element approximations for velocity and deviatoric stress and bilinear interpolation for pressure demonstrate numerical instability beyond a critical value of De for either the explicitly elliptic momentum equation (EEME) or elastic-viscous split-stress (EVSS) formulations, applying Galerkin's method for solution of the momentum and continuity equations, and using streamline upwind Petrov-Galerkin (SUPG) method for solution of the hyperbolic constitutive equation. The disturbance that causes the instability is concentrated near the stationary streamline of the base flow. The removal of this instability in a slightly modified form of the EEME formulation suggests that the instability results from coupling the approximations to the variables. A new mixed finite-element method, EVSS-G, is presented that includes smooth interpolation of the velocity gradients in the constitutive equation that is compatible with bilinear interpolation of the stress field. This formulation is tested with SUPG, streamline upwinding (SU), and Galerkin least squares (GLS) discretization of the constitutive equation. The EVSS-G/SUPG and EVSS-G/SU do not have the numerical instability described above; linear stability calculations for planar Couette flow are stable to values of De in excess of 50 and converge with mesh and time step. Calculations for the steady-state flow and its linear stability for a sphere falling in a tube demonstrate the appearance of linear instability to a time-periodic instability simultaneously with the apparent loss of existence of the steady-state solution. The instability appears as finely structured secondary cells that move from the front to the back of the sphere.Financial support for this research was given by the National Science Foundation, the Office of Naval Research, and the Defense Research Projects Agency. Computational resources were supplied by a grant from the Pittsburgh National Supercomputer Center and by the MIT Supercomputer Facility.     

Viscoelasticity: Fractal parameters studied on mammalian erythrocytes under shear stress

In this work, experimental determinations are carried out using a home-made device called an erythrodeformeter, which has been developed and constructed for rheological measurements on red blood cells subjected to definite fluid shear stress. A numerical method formulated on the basis of the fractal approximation for ordinary and fractionary Brownian motion¹ is proposed to evaluate the viscoelastic behavior of mammalian erythrocyte membranes. The diffraction pattern, which is circular when the mammalian erythrocyte membranes are at rest, becomes elliptical when the cells undergo shear stress. Photometric readings of light intensity variation along the major axis of the elliptical diffraction pattern are recorded during the creep and recovery process. These data series are used to calculate, fractal rheological parameters of self-affine Brownian motion on the erythrocytes, averaged over several millions of cells. Three different parameters over the time dependent process could be obtained, which are: correlation coefficient <C(t)>, correlation integral, andK ₂-entropy, and very different results were obtained.     

Woven fabric composite material model with material nonlinearity for nonlinear finite element simulation

The objective of the current investigation is to develop a simple, yet generalized, model which considers the two-dimensional extent of woven fabric, and to have an interface with nonlinear finite element codes. A micromechanical composite material model for woven fabric with nonlinear stress-strain relations is developed and implemented in ABAQUS for nonlinear finite element structural analysis. Within the model a representative volume cell is assumed. Using the iso-stress and iso-strain assumptions the constitutive equations are averaged along the thickness direction. The cell is then divided into many subcells and an averaging is performed again by assuming uniform stress distribution in each subcell to obtain the effective stress–strain relations of the subcell. The stresses and strains within the subcells are combined to yield the effective stresses and strains in the representative cell. Then this information is passed to the finite element code at each material point of the shell element. In this manner structural analysis of woven composites can be performed. Also, at each load increment global stresses and strains are communicated to the representative cell and subsequently distributed to each subcell. Once stresses and strains are associated to a subcell they can be distributed to each constituent of the subcell i.e. fill, warp, and resin. Consequently micro-failure criteria (MFC) can be defined for each constituent of a subcell and the proper stiffness degradation can be modeled if desired. This material model is suitable for implicit and could be modified for explicit finite element codes to deal with problems such as crashworthiness, impact, and failure analysis under static loads.     

A discontinuous Galerkin‐based sharp‐interface method to simulate three‐dimensional compressible two‐phase flow

A numerical method for the simulation of compressible two‐phase flows is presented in this paper. The sharp‐interface approach consists of several components: a discontinuous Galerkin solver for compressible fluid flow, a level‐set tracking algorithm to follow the movement of the interface and a coupling of both by a ghost‐fluid approach with use of a local Riemann solver at the interface. There are several novel techniques used: the discontinuous Galerkin scheme allows locally a subcell resolution to enhance the interface resolution and an interior finite volume Total Variation Diminishing (TVD) approximation at the interface. The level‐set equation is solved by the same discontinuous Galerkin scheme. To obtain a very good approximation of the interface curvature, the accuracy of the level‐set field is improved and smoothed by an additional P_NP_M‐reconstruction. The capabilities of the method for the simulation of compressible two‐phase flow are demonstrated for a droplet at equilibrium, an oscillating ellipsoidal droplet, and a shock‐droplet interaction problem at Mach 3. Copyright © 2015 John Wiley & Sons, Ltd.     

An approximation method to relate the linear viscoelastic functions

A method based on rational approximations is presented to interpolate the data from sinusoidal experiments in linear viscoelasticity. Bounds to the corresponding dynamical function and a discrete approximation to the spectrum are established. From this approximation the related viscoelastic functions can be computed. The method is checked by considering two theoretical models of physical interest and a satisfactory accuracy is achieved.     

A combined approach of the Laplace transform and Padé approximation solving viscoelasticity problems

In this paper the Laplace transform method is combined with Padé approximations to solve linear viscoelastic problems. This approach allows to avoid the usual difficulties of original function determination. An algorithm is given to find solution with arbitrary precision. As an example the solution for problem of viscoelastic orthotropic half-plane stress state under concentrated normal force is given.     

Accurate and straightforward symplectic approach for fracture analysis of fractional viscoelastic media

An accurate and straightforward symplectic method is presented for the fracture analysis of fractional two-dimensional(2D) viscoelastic media. The fractional Kelvin-Zener constitutive model is used to describe the time-dependent behavior of viscoelastic materials. Within the framework of symplectic elasticity, the governing equations in the Hamiltonian form for the frequency domain(s-domain) can be directly and rigorously calculated. In the s-domain, the analytical solutions of the displacement ...     

The distinction of viscoelastic phases from entangled wormlike micelles and of densely packed multilamellar vesicles on the basis of rheological measurements

It is shown in this work how two viscoelastic surfactant systems that are both shear thinning but differ in their morphology can be distinguished on the basis of rheological measurements. The measurements were carried out on the novel surfactant system cetyltrimethylammonium 2-hydroxy-1-naphthoate. The phases in this system are produced by mixing cetyltrimethylammonium hydroxide and 2-hydroxy-1-naphthoic acid. With increasing counterion surfactant ratio X, the system has two viscoelastic regions that are separated by a two-phase region. It is shown by cryo-transmission electron microscopy and by small angle neutron scattering that the first viscoelastic region which exists between X=0.5 and X=0.75 contains wormlike micelles, while the second viscoelastic region that exists between X=0.9 and X=1.4 contains multilamellar vesicles. Both phases look alike, are highly viscoelastic, have similar storage modulus values, and are shear thinning. The phases and the properties of the phases for the studied system are very similar as the phases for the system CTA-3-hydroxy-2-naphthoate that has been studied before (see Hassan et al. Langmuir, 12:4350–4357, 1996; Horbaschek et al. J Colloid Interface Sci, 206:439–456, 1998). The two viscoelastic phases show the same shear-thinning behavior, but differ in other rheological results. The phases can most easily be distinguished with the help of normal stress measurements. The wormlike viscoelastic solutions show large normal stresses that give rise to a large Weissenberg effect while the vesicle phases show no Weissenberg effect.     

Computer-extended series for a source/sink driven gas centrifuge

We have reformulated the general problem of internal flow in a modern, high speed gas centrifuge with sources and sinks in such a way as to obtain new, simple, rigorous closed form analytical solutions. Both symmetric and antisymmetric drives lead us to an ordinary differential equation in place of the usual inhomogeneous Onsager partial differential equation. Owing to the difficulties of exactly solving this sixth order, inhomogeneous, variable coefficient ordinary differential equation we appeal to the power of perturbation theory and techniques. Two extreme parameter regimes are identified, the so-called semi-long bowl approximation and a new short bowl approximation. Only the former class of problems is treated here. The long bowl solution for axial drive is the correct leading order term, just as for pure thermal drive. New O(1) results are derived for radial, drag and heat drives in two dimensions. Then regular asymptotic, even ordered power series expansions for the flow field are carried out on the computer to O(ε⁴) using MACSYMA. These approximations are valid for values of ε near unity. In the spirit of Van Dyke, one can carry out this expansion process, in theory, to apparently arbitrary order for arbitrary but finite decay length ratio. Curiously, the flows induced by axial and radial forces are proportional for asymptotically large source scale heights, x^*. Corresponding isotope separation integral parameters will be given in a companion paper.     

A physically based flux limiter for QUICK calculations of advective scalar transport

Transient, advective transport of a contaminant into a clean domain will exhibit a moving sharp front that separates contaminated and clean regions. Due to ‘numerical diffusion’—the combined effects of ‘cross‐wind diffusion’ and ‘artificial dispersion’—a numerical solution based on a first‐order (upwind) treatment will smear out the sharp front. The use of higher‐order schemes, e.g. QUICK (quadratic upwinding) reduces the smearing but can introduce non‐physical oscillations in the solution. A common approach to reduce numerical diffusion without oscillations is to use a scheme that blends low‐order and high‐order approximations of the advective transport. Typically, the blending is based on a parameter that measures the local monotonicity in the predicted scalar field. In this paper, an alternative approach is proposed for use in scalar transport problems where physical bounds C_Low?C?C_High on the scalar are known a priori. For this class of problems, the proposed scheme switches from a QUICK approximation to an upwind approximation whenever the predicted upwind nodal value falls outside of the physical range [C_Low, C_High]. On two‐dimensional steady‐state and one‐dimensional transient test problems predictions obtained with the proposed scheme are essentially indistinguishable from those obtained with monotonic flux‐limiter schemes. An analysis of the modified equation explains the observed performance of first‐ and second‐order time‐stepping schemes in predicting the advective transport of a step. In application to the transient two‐dimensional problem of contaminate transport into a streambed, predictions obtained with the proposed flux‐limiter scheme agree with those obtained with a scheme from the literature. Copyright © 2007 John Wiley & Sons, Ltd.     

Drag of a strongly heated sphere at small Reynolds numbers

Using an asymptotic small-perturbation method, the flow around a strongly heated sphere at small Reynolds numbers Re ≪ 1 is considered with account for thermal stresses in the gas in the higher-order approximations, beyond the Stokes one. It is assumed that the value of the Prandtl number Pr is arbitrary and the temperature dependence of the viscosity is described by a power law with an arbitrary exponent. In the O(Re²) and O(Re³ ln(Re)) approximations, the drag force of a heated sphere is found over a wide range of the ratios of sphere’s temperature to the gas free-stream temperature T _W /T _∞. The limits of applicability of the first (in Re) approximation are investigated, including the negative-drag effect, attributable to the action of the thermal stresses. The results are compared with numerical calculations of the flow around a hot sphere. The limits of applicability of the approximations found are examined. Similar results are obtained for the standard Navier-Stokes equations in which the thermal stresses are neglected.