On frequency–amplitude dependences for surface and internal standing waves

The analytic approach proposed by Sekerzh-Zenkovich [On the theory of standing waves of finite amplitude, Dokl. Akad. Nauk USSR 58 (1947) 551–554] is developed in the present study of standing waves. Generalizing the solution method, a set of standing wave problems are solved, namely, the infinite- and finite-depth surface standing waves and the infinite- and finite-depth internal standing waves. Two-dimensional wave motion of an irrotational incompressible fluid in a rectangular domain is considered to study weakly nonlinear surface and internal standing waves. The Lagrangian formulation of the problems is used and the fifth-order perturbation solutions are determined. Since most of the approximate analytic solutions to these problems were obtained using the Eulerian formulation, the comparison of the results, as an example the analytic frequency–amplitude dependences, obtained in Lagrangian variables with the corresponding ones known in Eulerian variables has been carried out in the paper. The analytic frequency–amplitude dependences are in complete agreement with previous results known in the literature. Computer algebra procedures were written for the construction of asymptotic solutions. The application of the model constructed in Lagrangian formulation to a set of different problems shows the ability to correctly reproduce and predict a wide range of situations with different characteristics and some advantages of Lagrangian particle models (for example, the bigger radius of convergence of an expansion parameter than in Eulerian variables, simplification of the boundary conditions, parametrization of a free boundary).     

Numerical Investigation of Dense Turbulent Sprays with an Eulerian Approach and DQMOM

The spray characteristics of the injected liquid fuel predominantly influence the combustion and emissions in IC-engines and gas turbines. This is predetermined by a dense spray region in which droplet-droplet interaction processes play a significant role. In order to accurately describe and control the dense spray behavior in modern engines, an appropriate numerical modeling tool is needed. This contribution aims at including droplet-droplet interactions into an Eulerian approach coupled to the Direct Quadrature Method of Moments (DQMOM) in order to describe evaporating droplet polydispersity and analyze dense turbulent sprays phenomena. Among the advantages of the Eulerian approach are a lower computational cost through optimal parallel computing and a straightforward liquid-gas phase coupling. To assess the designed tool, numerical results are compared to Phase Doppler Anemometry (PDA) measurements of a hollow-cone water spray. The experiment provides comprehensive validation data that include gas velocities, droplet size distribution, droplet mass fluxes and droplet velocities. Turbulence is captured by two different k-ε based models. Preliminary results show that the designed tool is able to capture the process under study in a satisfactory way. (© 2011 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

A numerical study on the effect of cavitation erosion in a diesel injector

The consequences of geometry alterations in a diesel injector caused by cavitation erosion are investigated with numerical simulations. The differences in the results between the nominal design geometry and the eroded one are analyzed for the internal injector flow and spray formation. The flow in the injector is modeled with a three-phase Eulerian approach using a compressible pressure-based multiphase flow solver. Cavitation is simulated with a nonequilibrium mass transfer rate model based on the simplified form of the Rayleigh–Plesset equation. Slip velocity between the liquid-vapor mixture and air is included in the model by solving two separate momentum conservation equations. The eroded injector is found to result in a loss in the rate of injection but also lower cavitation volume fraction inside the nozzle. The injected sprays are then simulated with a Lagrangian method considering as initial conditions the predicted flow characteristics at the exit of the nozzle. The results obtained show wider spray dispersion for the eroded injector and shorter spray tip penetration.     

A Lagrangian formulation for steady three-dimensional Newton-Busemann flow

A new Lagrangian formulation for steady three dimensional inviscid flow over rigid bodies is developed. First, the continuity equation is eliminated by the use of two stream functions. This is followed by a transformation to new independent variables, two of which are these stream functions and the third one is a Lagrangian time distinct from the Eulerian time. This Lagrangian formulation with the use of Lagrangian time requires only three independent variables and allows the free boundary problem of flow with shock wave to be rendered a fixed boundary one thereby making it easier to solve. In the Newtonian limit the governing equations reduce to the geodesic equations of the body surface. For flow past two-dimensional bodies, bodies of revolution, and conical bodies and wings, the problems are solved to within quadrature. All known Newtonian flow solutions are found to be special cases of the present theory.     

Formulation and survey of ALE method in nonlinear solid mechanics

This paper investigates the applicability and accuracy of existing formulation methods in general purpose finite element programs to the finite strain deformation problems. The basic shortcomings in using such programs in these applications are then pointed out and the need for a different type of formulation is discussed. An arbitrary Lagrangian-Eulerian (ALE) method is proposed and a concise survey of ALE formulation is given. A consistent and complete ALE formulation is derived from the virtual work equation transformed to arbitrary computational reference configurations. Differences between the proposed formulations and similar ones in the literature are discussed. The proposed formulation presents a general approach to ALE method. It includes load correction terms and is suitable for rate-dependent and rate-independent material constitutive law. The proposed formulation reduces to both updated Lagrangian and Eulerian formulations as special cases.     

Perturbation analysis of short-crested waves in Lagrangian coordinates

A new second-order asymptotic solution that describes short-crested waves is derived in Lagrangian coordinates. The analytical Lagrangian solution that is uniformly valid satisfies the irrotational condition and there being zero pressure at the free surface, in contrast with the Eulerian solution, in which there is residual pressure at the free surface. The explicit parametric solution highlights the trajectory of a water particle and the wave kinematics above the mean water level. The mass transport velocity and Lagrangian mean level associated with particle displacement can also be obtained directly. In particular, the mean level of the particle motion in a Lagrangian form differs that of the Eulerian form. The new formulation reduces to second-order standing or progressive wave solutions in Lagrangian coordinates at the limiting angles of approach. Expressions for kinematic quantities are also presented.     

A family of ELLAM schemes for advection-diffusion-reaction equations and their convergence analyses

A family of ELLAM (Eulerian–Lagrangian localized adjoint method) schemes is developed and analyzed for linear advection-diffusion-reaction transport partial differential equations with any combination of inflow and outflow Dirichlet, Neumann, or flux boundary conditions. The formulation uses space-time finite elements, with edges oriented along Lagrangian flow paths, in a time–stepping procedure, where space-time test functions are chosen to satisfy a local adjoint condition. This allows Eulerian–Lagrangian concepts to be applied in a systematic mass-conservative manner, yielding numerical schemes defined at each discrete time level. Optimal-order error estimates and superconvergence results are derived. Numerical experiments are performed to verify the theoretical estimates. © 1998 John Wiley & Sons, Inc. Numer Methods Partial Differential Eq 14: 739–780, 1998     

A numerical study on statistical temporal scales in inertia particle dispersion

The statistical temporal scales involved in inertia particle dispersion are analyzed numerically. The numerical method of large eddy simulation, solving a filtered Navier-Stokes equation, is utilized to calculate fully developed turbulent channel flows with Reynolds numbers of 180 and 640, and the particle Lagrangian trajectory method is employed to track inertia particles released into the flow fields. The Lagrangian and Eulerian temporal scales are obtained statistically for fluid tracer particles and three different inertia particles with Stokes numbers of 1, 10 and 100. The Eulerian temporal scales, decreasing with the velocity of advection from the wall to the channel central plane, are smaller than the Lagrangian ones. The Lagrangian temporal scales of inertia particles increase with the particle Stokes number. The Lagrangian temporal scales of the fluid phase ‘seen’ by inertia particles are separate from those of the fluid phase, where inertia particles travel in turbulent vortices, due to the particle inertia and particle trajectory crossing effects. The effects of the Reynolds number on the integral temporal scales are also discussed. The results are worthy of use in examining and developing engineering prediction models of particle dispersion.     

On the Eulerian–Lagrangian formulation of balance equations in porous media

The article presents a general approach to modeling the transport of extensive quantities in the case of flow of multiple multicomponent fluid phases in a deformable porous medium domain under nonisothermal conditions. The models are written in a modified Eulerian–Lagrangian formulation. In this modified formulation, the material derivatives are written in terms of modified velocities. These are the velocities at which the various phase and component variables propagate in the domain, along their respective characteristic curves. It is shown that these velocities depend on the heterogeneity of various solid matrix and fluid properties. The advantage of this formulation, with respect to the usually employed Eulerian one, is that numerical dispersion, associated with the advective fluxes of extensive quantities, are eliminated. The methodology presented in the article shows how the Eulerian–Lagrangian formulation is written in terms of the relatively small number of primary variables of a transport problem. © 1997 John Wiley & Sons, Inc. Numer Methods Partial Differential Eq 13: 505–530, 1997     

A comparison of Eulerian and updated Lagrangian finite element algorithms for simulating film casting

This paper describes and compares Eulerian (E) and Updated Lagrangian (UL) finite element algorithms for simulating film casting. It is found that the E algorithm is more accurate, faster and simpler to implement for one-dimensional film casting of a viscous fluid. However, the UL algorithm provides a natural and intuitive framework for accommodating viscoelasticity and for investigating the stability of the film.     

Comparison between Eulerian and Lagrangian semi-analytical models to simulate the pollutant dispersion in the PBL

In this work we display a numerical comparison, under statistical and computational point of view, between semi-analytical Eulerian and Lagrangian dispersion models to simulate the ground-level concentration values of a passive pollutant released from a low height source. The Eulerian approach is based on the solution of the advection–diffusion equation by the Laplace transform technique. The Lagrangian approach is based on solution of the Langevin equation through the Picard’s Iterative Method. Turbulence inputs are calculated according to a parameterization capable of generating continuous values in all stability conditions and in all heights of the Planetary Boundary Layer (PBL). Numerical simulations and comparisons show a good agreement between predicted and observed concentrations values. The comparison reveals the main advantages and disadvantages between the models.     

Predictions of the gas–liquid flow in wet electrostatic precipitators

Based on Computational Fluid Dynamics (CFD), the present paper aims to simulate several important phenomena in a wet type ESP from the liquid spray generation to gas-droplet flow in electric field. A single passage between the adjacent plates is considered for the simulation domain. Firstly, the electric field intensity and ion charge density are solved locally around a corona emitter of a barbed wire electrode, which are applied to the entire ESP using periodic conditions. Next, the Euler–Lagrange method is used to simulate the gas-droplet flow. Water droplets are tracked statistically along their trajectories, together with evaporation and particle charging. Finally, the deposition density on the plate is taken as the input for the liquid film model. The liquid film is simulated separately using the homogenous Eulerian approach in ANSYS-CFX. In the current case, since the free surface of the thin water film is difficult to resolve, a special method is devised to determine the film thickness.As parametric study, the variables considered include the nozzle pressure, initial spray spreading patterns (solid versus hollow spray) and plate wettability. The droplet emission rate and film thickness distribution are the results of interest. Main findings: electric field has strong effect on the droplet trajectories. Hollow spray is preferred to solid spray for its lower droplet emission. The liquid film uniformity is sensitive to the plate wettability.     

Adjoint-based optimal control using meshfree discretizations

The paper at hand presents a combination of optimal control approaches for PDEs with meshless discretizations. Applying a classical Lagrangian type particle method to optimization problems with hyperbolic constraints, several adjoint-based strategies differing in the sequential order of optimization and discretization of the Lagrangian or Eulerian problem formulation are proposed and compared. The numerical results confirm the theoretically predicted independence principle of the optimization approaches and show the expected convergence behavior. Moreover, they exemplify the superiority of meshless methods over the conventional mesh-based approaches for the numerical handling and optimization of problems with time-dependent geometries and freely moving boundaries.     

An Eulerian‐Lagrangian substructuring domain decomposition method for unsteady‐state advection‐diffusion equations

We develop an Eulerian‐Lagrangian substructuring domain decomposition method for the solution of unsteady‐state advection‐diffusion transport equations. This method reduces to an Eulerian‐Lagrangian scheme within each subdomain and to a type of Dirichlet‐Neumann algorithm at subdomain interfaces. The method generates accurate and stable solutions that are free of artifacts even if large time‐steps are used in the simulation. Numerical experiments are presented to show the strong potential of the method. © 2001 John Wiley & Sons, Inc. Numer Methods Partial Differential Eq 17:565–583, 2001     

Thermo-mechanical computation of inelastic pavement structures under rolling tire loads

A realistic and numerically efficient computation of the tire-pavement interaction is essential for the investigation of the structural behavior of pavements under rolling tire load as base of the development of more durable pavement structures. The paper presents a thermo-mechanical coupled simulation model based on an Arbitrary Lagrangian Eulerian (ALE) formulation that considers the inelastic and temperature dependent material properties of asphalt. (© 2015 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Multilevel numerical solutions of convection‐dominated diffusion problems by spline wavelets

In this article, we utilize spline wavelets to establish an adaptive multilevel numerical scheme for time‐dependent convection‐dominated diffusion problems within the frameworks of Galerkin formulation and Eulerian‐Lagrangian localized adjoint methods (ELLAM). In particular, we shall use linear Chui‐Quak semi‐orthogonal wavelets, which have explicit expressions and compact supports. Therefore, both the diffusion term and boundary conditions in the convection‐diffusion problems can be readily handled. Strategies for efficiently implementing the scheme are discussed and numerical results are interpreted from the viewpoint of nonlinear approximation. © 2005 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2006     

Convergence of a finite element/ALE method for the Stokes equations in a domain depending on time

We consider the approximation of the unsteady Stokes equations in a time dependent domain when the motion of the domain is given. More precisely, we apply the finite element method to an Arbitrary Lagrangian Eulerian (ALE) formulation of the system. Our main results state the convergence of the solutions of the semi-discretized (with respect to the space variable) and of the fully-discrete problems towards the solutions of the Stokes system.     

Positive and Entropy Stable Godunov-type Schemes for Gas Dynamics and MHD Equations in Lagrangian or Eulerian Coordinates

Summary. Using the concept of simple Riemann solvers, we present entropic and positive Godunov-type schemes preserving contact discontinuities for both Lagrangian and Eulerian systems of gas dynamics and magnetohydrodynamics (MHD). On the one hand, for the Lagrangian form, we develop positive and entropic Riemann solvers which can be considered as a natural extension of Roe's solvers in which the sound speed is relaxed. On the other hand, for the Eulerian form, we are able to construct by two ways Godunov-type schemes based on Lagrangian simple Riemann solvers. The first method establishes a relation between the jump of the intermediate states and the second one between the intermediate states themselves.     

An Eulerian‐Lagrangian Concept for the Numerical Simulation of Flat Steady‐State Hot Rolling Processes

For the finite element analysis of stationary flat hot rolling processes, a new and efficient mathematical model was developed. The method is based on an intermediary Eulerian‐Lagrangian concept, where an Eulerian coordinate is employed in the rolling direction, while Lagrangian coordinates are used in the direction of the thickness and width of the strip. This approach yields an efficient algorithm, where the time is eliminated as an independent variable in the steady‐state case. Further, the vector of independent field variables consists of a velocity component in Eulerian and of displacement components in Lagrangian directions. Due to this concept, the free surface deformations can be accounted for directly and the problems encountered with pure Eulerian or Lagrangian models now appear with reduced complexity and can thus be tackled more easily. The general formalism was applied to different practical hot rolling situations, ranging from thick slabs to ultra‐thin hot strips. (© 2004 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

An Eulerian‐Lagrangian control volume scheme for two‐dimensional unsteady advection‐diffusion problems

We develop a mass conservative Eulerian‐Lagrangian control volume scheme (ELCVS) for the solution of the transient advection‐diffusion equations in two space dimensions. This method uses finite volume test functions over the space‐time domain defined by the characteristics within the framework of the class of Eulerian‐Lagrangian localized adjoint characteristic methods (ELLAM). It, therefore, maintains the advantages of characteristic methods in general, and of this class in particular, which include global mass conservation as well as a natural treatment of all types of boundary conditions. However, it differs from other methods in that class in the treatment of the mass storage integrals at the previous time step defined on deformed Lagrangian regions. This treatment is especially attractive for orthogonal rectangular Eulerian grids composed of block elements. In the algorithm, each deformed region is approximated by an eight‐node region with sides drawn parallel to the Eulerian grid, which significantly simplifies the integration compared with the approach used in finite volume ELLAM methods, based on backward tracking, while retaining local mass conservation at no additional expenses in terms of accuracy or CPU consumption. This is verified by numerical tests which show that ELCVS performs as well as standard finite volume ELLAM methods, which have previously been shown to outperform many other well‐received classes of numerical methods for the equations considered. © 2011 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2012