An arbitrary Lagrangian–Eulerian finite element approach to non‐steady state turbulent fluid flow with application to mould filling in casting

This paper presents a two‐dimensional Lagrangian–Eulerian finite element approach of non‐steady state turbulent fluid flows with free surfaces. The proposed model is based on a velocity–pressure finite element Navier–Stokes solver, including an augmented Lagrangian technique and an iterative resolution of Uzawa type. Turbulent effects are taken into account with the k–ε two‐equation statistical model. Mesh updating is carried out through an arbitrary Lagrangian–Eulerian (ALE) method in order to describe properly the free surface evolution. Three comparisons between experimental and numerical results illustrate the efficiency of the method. The first one is turbulent flow in an academic geometry, the second one is a mould filling in effective casting conditions and the third one is a precise confrontation to a water model. Copyright © 2000 John Wiley & Sons, Ltd.     

Unsteady circular Couette flow of a Bingham plastic with the Augmented Lagrangian Method

Numerical simulations are undertaken for unsteady flows of an ideal Bingham fluid in a circular Couette viscometer. The main difficulties in such simulations are caused by the non-differentiability of the constitutive equation and the need to determine the position and shape of the yield surface separating the yielded zones from the unyielded ones. In this work, these difficulties are overcome by using a numerical method based on variational inequalities, i.e. the augmented Lagrangian/Uzawa method. The start-up and cessation of circular Couette flows of a Bingham fluid are solved numerically assuming that only one of the cylinders is rotating. An improved theoretical upper bound for the stopping time in the case of cessation is derived. The numerical estimates for the stopping time compare well with the theoretical bounds. Moreover, with the adopted method the evolution of the velocity profiles and the locations of yielded/unyielded surfaces are accurately calculated. In flow cessation, we observe an interesting effect, namely the appearance of a small unyielded region adjoined to the outer cylinder shortly before cessation.     

Numerical simulations of cessation flows of a Bingham plastic with the augmented Lagrangian method

The augmented Lagrangian/Uzawa method has been used to study benchmark one-dimensional cessation flow problems of a Bingham fluid, such as the plane Couette flow, and the plane, round, and annular Poiseuille flows. The calculated stopping times agree well with available theoretical upper bounds for the whole range of Bingham numbers and with previous numerical results. The applied method allows for easy determination of the yielded and unyielded regions. The evolution of the rigid zones in these unsteady flows is presented. It is demonstrated that the appearance of an unyielded zone near the wall occurs for any non-zero Bingham number not only in the case of a round tube but also in the case of an annular tube of small radii ratio. The advantages of using the present method instead of regularizing the constitutive equation are also discussed.     

Kinematic shakedown by the Norton–Hoff–Friaa regularising method and augmented LagrangianAdaptation cinématique par la methode de régularisation de Norton–Hoff–Friaa et du Lagrangien augmenté

The shakedown analysis of elastic perfectly plastic structures is formulated as a discrete nonlinear mathematical programming problem by means of the finite element technique. The kinematical problem is regularized through the introduction of the Norton–Hoff viscoplastic material to overcome the non-differentiability of the objective function, and can be solved numerically by the augmented Lagrangian technique. To cite this article: M.A. Hamadouche, C. R. Mecanique 330 (2002) 305–311.     

EXACT AUGMENTED LAGRANGIAN FUNCTION FOR NONLINEAR PROGRAMMING PROBLEMS WITH INEQUALITY CONSTRAINTS

An exact augmented Lagrangian function for the nonlinear nonconvex programming problems with inequality constraints was discussed. Under suitable hypotheses, the relationship was established between the local unconstrained minimizers of the augmented Lagrangian function on the space of problem variables and the local minimizers of the original constrained problem. Furthermore, under some assumptions, the relationship was also established between the global solutions of the augmented Lagrangian function on some compact subset of the space of problem variables and the global solutions of the constrained problem. Therefore, from the theoretical point of view, a solution of the inequality constrained problem and the corresponding values of the Lagrange multipliers can be found by the well-known method of multipliers which resort to the unconstrained minimization of the augmented Lagrangian function presented.     

A substructure‐based iterative inner solver coupled with Uzawa's algorithm for the Stokes problem

A domain decomposition method with Lagrange multipliers for the Stokes problem is developed and analysed. A common approach to solve the Stokes problem, termed the Uzawa algorithm, is to decouple the velocity and the pressure. This approach yields the Schur complement system for the pressure Lagrange multiplier which is solved with an iterative solver. Each outer iteration of the Uzawa procedure involves the inversion of a Laplacian in each spatial direction. The objective of this paper is to effectively solve this inner system (the vector Laplacian system) by applying the finite‐element tearing and interconnecting (FETI) method. Previously calculated search directions for the FETI solver are reused in subsequent outer Uzawa iterations. The advantage of the approach proposed in this paper is that pressure is continuous across the entire computational domain. Numerical tests are performed by solving the driven cavity problem. An analysis of the number of outer Uzawa iterations and inner FETI iterations is reported. Results show that the total number of inner iterations is almost numerically scalable since it grows asymptotically with the mesh size and the number of subdomains. Copyright © 2003 John Wiley & Sons, Ltd.     

An efficient operator splitting scheme for three-dimensional hydrodynamic computations

A new efficient numerical method for three-dimensional hydrodynamic computations is presented and discussed in this paper. The method is based on the operator splitting method and combined with Eulerian–Lagrangian method, finite element method and finite difference method. To increase the efficiency and stability of the numerical solutions, the operator splitting method is employed to partition the momentum equations into three parts, according to physical phenomena. A time step is divided into three time substeps. In the first substep, advection and Coriolis force are solved using the explicit Eulerian–Lagrangian method. In the second substep, horizontal diffusion is approximated by implicit FEM in each horizontal layer. In the last substep, the continuity equation is solved by implicit FEM, and vertical diffusion and pressure gradient are discretized by implicit FDM in each nodal column. The stability analysis shows that this method is unconditionally stable. A number of numerical experiments have been performed. The results simulated by the present scheme agree well with analytical solutions and the other documented model results. The method is efficient for 3D shallow water flow computations and fully fits complicated configurations. © 1998 John Wiley & Sons, Ltd.     

A staggered grid‐based explicit jump immersed interface method for two‐dimensional Stokes flows

In this paper the explicit jump immersed interface method (EJIIM) is applied to stationary Stokes flows. The boundary value problem in a general, non‐grid aligned domain is reduced by the EJIIM to a sequence of problems in a rectangular domain, where staggered grid‐based finite differences for velocity and pressure variables are used. Each of these subproblems is solved by the fast Stokes solver, consisting of the pressure equation (known also as conjugate gradient Uzawa) method and a fast Fourier transform‐based Poisson solver. This results in an effective algorithm with second‐order convergence for the velocity and first order for the pressure. In contrast to the earlier versions of the EJIIM, the Dirichlét boundary value problem is solved very efficiently also in the case when the computational domain is not simply connected. Copyright © 2007 John Wiley & Sons, Ltd.     

A meshless finite point method for three‐dimensional analysis of compressible flow problems involving moving boundaries and adaptivity

A finite point method for solving compressible flow problems involving moving boundaries and adaptivity is presented. The numerical methodology is based on an upwind‐biased discretization of the Euler equations, written in arbitrary Lagrangian–Eulerian form and integrated in time by means of a dual‐time steeping technique. In order to exploit the meshless potential of the method, a domain deformation approach based on the spring network analogy is implemented, and h‐adaptivity is also employed in the computations. Typical movable boundary problems in transonic flow regime are solved to assess the performance of the proposed technique. In addition, an application to a fluid–structure interaction problem involving static aeroelasticity illustrates the capability of the method to deal with practical engineering analyses. The computational cost and multi‐core performance of the proposed technique is also discussed through the examples provided. Copyright © 2013 John Wiley & Sons, Ltd.     

Efficient meshless SPH method for the numerical modeling of thick shell structures undergoing large deformations

The objective of this paper is to present an extension of the Lagrangian Smoothed Particle Hydrodynamics (SPH) method to solve three-dimensional shell-like structures undergoing large deformations. The present method is an enhancement of the classical stabilized SPH commonly used for 3D continua, by introducing a Reissner–Mindlin shell formulation, allowing the modeling of moderately thin structure using only one layer of particles in the shell mid-surface. The proposed Shell-based SPH method is efficient and very fast compared to the classical continuum SPH method. The Total Lagrangian Formulation valid for large deformations is adopted using a strong formulation of the differential equilibrium equations based on the principle of collocation. The resulting non-linear dynamic problem is solved incrementally using the explicit time integration scheme, suited to highly dynamic applications. To validate the reliability and accuracy of the proposed Shell-based SPH method in solving shell-like structure problems, several numerical applications including geometrically non-linear behavior are performed and the results are compared with analytical solutions when available and also with numerical reference solutions available in the literature or obtained using the Finite Element method by means of ABAQUS^© commercial software.     

Partitioned solution to fluid–structure interaction problem in application to free-surface flows

In this work we discuss a way to compute the impact of free-surface flow on nonlinear structures. The approach chosen relies on a partitioned strategy that allows us to solve the strongly coupled fluid–structure interaction problem. It is then possible to re-use the existing and validated strategy for each sub-problem. The structure is formulated in a Lagrangian way and solved by the finite element method. The free-surface flow approach considers a Volume-Of-Fluid (VOF) strategy formulated in an Arbitrary Lagrangian–Eulerian (ALE) framework, and the finite volume is used to discrete and solve this problem. The software coupling is ensured in an efficient way using the Communication Template Library (CTL). Numerical examples presented herein concern the 2D validation case but also 3D problems with a large number of equations to be solved.     

Viscoelastic flow between eccentric rotating cylinders: unstructured control volume method

This paper reports a convergent numerical algorithm for the Upper-Convected Maxwell (UCM) fluid between two eccentric cylinders at various eccentricity ratios (?); the outer cylinder is stationary, and the inner one rotating. The problem is solved by an unstructured control volume method (UCV), which is designed for a general viscoelastic flow problem with an arbitrary computational domain. A self-consistent false diffusion technique and an iteration scheme are used in combination to solve the problem. The computations of the UCM fluid using the numerical algorithm are carried out to a higher value of the Deborah number (De) at each eccentricity tested than hitherto possible with previous numerical simulations. The solutions are compared with previous numerical results, confirming the effectiveness of the UCV method as a general technique for solving viscoelastic flow problems.     

An implicit method for frictional contact, impact and rolling

In this paper an implicit method for frictional contact, impact and rolling is suggested. A nonclassical formulation of a two-dimensional hyperelastic body unilaterally constrained to rigid supports is proposed by following the ideas of Moreau and Jean. A total Lagrangian formulation of the system is given. The elastic properties are defined by coupling the second Piola–Kirchhoff stress to the Green–Lagrange strain via the Kirchhoff–St. Venant law. The equation of motion is written in the spirit of Moreau by using the mean value impulses introduced by Jean. The mean value impulses appear explicitly in the equation of motion. In such manner the treatment of nonconstant kinematic transformation matrices becomes straightforward. The rigid supports are described by smooth functions. By utilizing these functions and the mean value impulses, new contact/impact laws of Signorini and Coulomb type are formulated. The governing equations are solved by a nonsmooth Newton method. This is performed by following the augmented Lagrangian approach and deriving the consistent stiffness matrix as well as the contact stiffness matrices. Three two-dimensional examples are solved by the method: a contact problem, an impact problem and a rolling contact problem.     

Global Continuous Riemann Solver for Nonlinear Elasticity

We consider in this paper an isothermal model of nonlinear elasticity. This model is described by two conservation laws that define a problem of mixed type, both elliptic and hyperbolic. We restrict ourselves to the linearly degenerate case, and consider Riemann data that lies in the hyperbolic regions. The lack of uniqueness of the Riemann problem is solved by the introduction of a so-called kinetic relation, used to narrow the set of admissible subsonic phase transitions. In this situation, we consider the Riemann problem for any data lying in the hyperbolic region, using either explicit computations or geometric arguments. This construction allows us to give sufficient conditions on the kinetic relation in order that the generated Riemann solver possesses properties of uniqueness, globality, and continuous dependence on the initial data in the L ¹ distance. Accepted October 1, 2000?Published online January 22, 2001     

A preconditioned conjugate gradient Uzawa-type method for the solution of the stokes problem by mixed Q1–P0 stabilized finite elements

We study the behaviour of a conjugate gradient Uzawa-type method for a stabilized finite element approximation of the Stokes problem. Many variants of the Uzawa algorithm have been described for different finite elements satisfying the well-known Inf-Sup condition of Babu?ka and Brezzi, but it is surprising that developments for unstable ‘low-order’ discretizations with stabilization procedures are still missing. Our paper is presented in this context for the popular (so-called) Q1–P0 element. First we show that a simple stabilization technique for this element permits us to retain the property of a convergence factor bounded independently of the discretization mesh size. The second contribution of this work deals with the construction of a less costly preconditioner taking full advantages of the block diagonal structure of the stabilization matrix. Its efficiency is supported by 2D and. 3D numerical results.     

Optimum Suction Distribution for Transition Control

The optimum suction distribution which gives the longest laminar region for a given total suction is computed. The goal here is to provide the designer with a method of finding the best suction distribution subject to some overall constraints applied to the suction. We formulate the problem using the Lagrangian multiplier method with constraints. The resulting nonlinear system of equations is solved using the Newton–Raphson technique. The computations are performed for a Blasius boundary layer on flat-plate and crossflow cases. For the Blasius boundary layer, the optimum suction distribution peaks upstream of the maximum growth rate region and remains flat in the middle before it decreases to zero at the end of the transition point. For the stationary and travelling crossflow instability, the optimum suction peaks upstream of the maximum growth rate region and decreases gradually to zero. Received 8 May 1997 and accepted 5 November 1998     

多层横观各向同性圆柱壳的三维弹性理论解

本文从三维弹性理论出发,用特征函数法研究多层横观各向同性圆柱壳的轴对称问题.把位移和应力分量的齐次解表达成特征函数展开式,并把特解部分用Fourier级数表示.以多层圆柱壳的内、外柱面作为齐次边界,同时考虑层间的连续条件,推导出问题的特征方程并用Muller法求解.文中运用传递矩阵技术处理多层问题,并用边界型最小二乘配点法处理端部边界条件.作为实例,对双层圆柱壳作了数值计算.     

A marker‐and‐cell approach to free surface 2‐D multiphase flows

This work describes a methodology to simulate free surface incompressible multiphase flows. This novel methodology allows the simulation of multiphase flows with an arbitrary number of phases, each of them having different densities and viscosities. Surface and interfacial tension effects are also included. The numerical technique is based on the GENSMAC front‐tracking method. The velocity field is computed using a finite‐difference discretization of a modification of the Navier–Stokes equations. These equations together with the continuity equation are solved for the two‐dimensional multiphase flows, with different densities and viscosities in the different phases. The governing equations are solved on a regular Eulerian grid, and a Lagrangian mesh is employed to track free surfaces and interfaces. The method is validated by comparing numerical with analytic results for a number of simple problems; it was also employed to simulate complex problems for which no analytic solutions are available. The method presented in this paper has been shown to be robust and computationally efficient. Copyright © 2012 John Wiley & Sons, Ltd.     

The Stress–Strain State of Geometrically and Physically Nonlinear Plates with One-Sided Corrosive Wear

A technique is proposed to investigate one-sided corrosive wear. The problem is solved with regard for geometric and physical nonlinearity. Two, Dolinskii's and Gutman's corrosion models are considered. The quasistatic problem is solved by the method of variational iterations, which reduce ordinary differential equations to a system of nonlinear equations with approximation o(h ²) to be solved by Newton's method. At each step, to allow for physical nonlinearity, the method of variable elastic parameters is used. Also a technique is developed to consider various boundary conditions and ⁱ(e ⁱ) diagrams. Specific numerical results are presented.