Home About us Contact
|
Home About us Contact | ||
|
|
||
Lagrangian
Kinds of Lagrangian Terms modified by Lagrangian Selected AbstractsA thermodynamic approach to the instantaneous non-active powerEUROPEAN TRANSACTIONS ON ELECTRICAL POWER, Issue 6 2001A. P. Morando, Article first published online: 22 MAR 200 Having schematically run through the transition from single-phase to three-phase relationships, the energy balance is formalised using Park vector notation. The imaginary power notation emerges. Leading back, in the sinusoidal case, to the usual reactive power, it generalises its specific contents in the case of variable states, explaining in particular a typical aspect of three-phase networks: the energy bouncing from one phase to another. This aspect can be seen as an index of power quality. At last, these same considerations are obtained by means of Lagrangian and thermodynamic approaches that lend a deeper meaning to the energy related quantities. [source] Scattering of charged tensor bosons in gauge and superstring theoriesFORTSCHRITTE DER PHYSIK/PROGRESS OF PHYSICS, Issue 7-9 2010I. Antoniadis Abstract We calculate the leading-order scattering amplitude of one vector and two tensor gauge bosons in a recently proposed non-Abelian tensor gauge field theory and open superstring theory. The linear in momenta part of the superstring amplitude has identical Lorentz structure with the gauge theory, while its cubic in momenta part can be identified with an effective Lagrangian which is constructed using generalized non-Abelian field strength tensors. [source] Optimised Dirac operators on the lattice: construction, properties and applicationsFORTSCHRITTE DER PHYSIK/PROGRESS OF PHYSICS, Issue 2 20082Article first published online: 29 NOV 200, W. Bietenholz Abstract We review a number of topics related to block variable renormalisation group transformations of quantum fields on the lattice, and to the emerging perfect lattice actions. We first illustrate this procedure by considering scalar fields. Then we proceed to lattice fermions, where we discuss perfect actions for free fields, for the Gross-Neveu model and for a supersymmetric spin model. We also consider the extension to perfect lattice perturbation theory, in particular regarding the axial anomaly and the quark gluon vertex function. Next we deal with properties and applications of truncated perfect fermions, and their chiral correction by means of the overlap formula. This yields a formulation of lattice fermions, which combines exact chiral symmetry with an optimisation of further essential properties. We summarise simulation results for these so-called overlap-hypercube fermions in the two-flavour Schwinger model and in quenched QCD. In the latter framework we establish a link to Chiral Perturbation Theory, both, in the p -regime and in the ,-regime. In particular we present an evaluation of the leading Low Energy Constants of the chiral Lagrangian , the chiral condensate and the pion decay constant , from QCD simulations with extremely light quarks. [source] Arbitrary Lagrangian,Eulerian method for large-strain consolidation problemsINTERNATIONAL JOURNAL FOR NUMERICAL AND ANALYTICAL METHODS IN GEOMECHANICS, Issue 9 2008Majidreza Nazem Abstract In this paper, an arbitrary Lagrangian,Eulerian (ALE) method is generalized to solve consolidation problems involving large deformation. Special issues such as pore-water pressure convection, permeability and void ratio updates due to rotation and convection, mesh refinement and equilibrium checks are discussed. A simple and effective mesh refinement scheme is presented for the ALE method. The ALE method as well as an updated-Lagrangian method is then used to solve some classical consolidation problems involving large deformations with different constitutive laws. The results clearly show the advantage and efficiency of the ALE method for these examples. Copyright © 2007 John Wiley & Sons, Ltd. [source] An operator-split ALE model for large deformation analysis of geomaterialsINTERNATIONAL JOURNAL FOR NUMERICAL AND ANALYTICAL METHODS IN GEOMECHANICS, Issue 12 2007Y. Di Abstract Analysis of large deformation of geomaterials subjected to time-varying load poses a very difficult problem for the geotechnical profession. Conventional finite element schemes using the updated Lagrangian formulation may suffer from serious numerical difficulties when the deformation of geomaterials is significantly large such that the discretized elements are severely distorted. In this paper, an operator-split arbitrary Lagrangian,Eulerian (ALE) finite element model is proposed for large deformation analysis of a soil mass subjected to either static or dynamic loading, where the soil is modelled as a saturated porous material with solid,fluid coupling and strong material non-linearity. Each time step of the operator-split ALE algorithm consists of a Lagrangian step and an Eulerian step. In the Lagrangian step, the equilibrium equation and continuity equation of the saturated soil are solved by the updated Lagrangian method. In the Eulerian step, mesh smoothing is performed for the deformed body and the state variables obtained in the updated Lagrangian step are then transferred to the new mesh system. The accuracy and efficiency of the proposed ALE method are verified by comparison of its results with the results produced by an analytical solution for one-dimensional finite elastic consolidation of a soil column and with the results from the small strain finite element analysis and the updated Lagrangian analysis. Its performance is further illustrated by simulation of a complex problem involving the transient response of an embankment subjected to earthquake loading. Copyright © 2007 John Wiley & Sons, Ltd. [source] Minimum principle and related numerical scheme for simulating initial flow and subsequent propagation of liquefied groundINTERNATIONAL JOURNAL FOR NUMERICAL AND ANALYTICAL METHODS IN GEOMECHANICS, Issue 11 2005Sami Montassar Abstract The problem of predicting the evolution of liquefied ground, modelled as a viscoplastic material, is addressed by combining a minimum principle for the velocity field, which characterizes such an evolution, and a time step integration procedure. Two different numerical schemes are then presented for the finite element implementation of this minimum principle, namely, the regularization technique and the decomposition-co-ordination method by augmented Lagrangian. The second method, which proves more accurate and efficient than the first, is finally applied to simulate the incipient flow failure and subsequent spreading of a liquefied soil embankment subject to gravity. The strong influence of liquefied soil residual shear strength on reducing the maximum amplitude of the ground displacement is particularly emphasized in such an analysis. Copyright © 2005 John Wiley & Sons, Ltd. [source] A new approach to avoid excessive numerical diffusion in Eulerian,Lagrangian methodsINTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN BIOMEDICAL ENGINEERING, Issue 11 2008A. Younes Abstract Lumping is often used to avoid non-physical oscillations for advection,dispersion equations but is known to add numerical diffusion. A new approach is detailed in order to avoid excessive numerical diffusion in Eulerian,Lagrangian methods when several time steps are used. The basic idea of this approach is to keep the same characteristics during all time steps and to interpolate only the concentration variations due to the dispersion process. In this way, numerical diffusion due to the lumping is removed at the end of each time step. The method is combined with the Eulerian,Lagrangian localized adjoint method (ELLAM) which is a mass conservative characteristic method for solving the advection,dispersion equation. Two test problems are modelled to compare the proposed method to the consistent, the full and the selective lumping approaches for linear and non-linear transport equations. Copyright © 2007 John Wiley & Sons, Ltd. [source] Hamiltonian-based error computationsINTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN BIOMEDICAL ENGINEERING, Issue 4 2006Y. L. Kuo Abstract This paper presents two sets of the Hamiltonian for checking errors of approximated solutions. The first set can be applied to those problems having any number of independent and dependent variables. This set of the Hamiltonian can effectively indicate the errors of approximated solutions when requiring a high accuracy. The second set of the Hamiltonian has the invariant property when the Lagrangian is not an explicit function of time, even for non-conservative systems. Both sets can be formulated as error indicators to check errors of approximated solutions. Three illustrative examples demonstrate the error analyses of finite element solutions. Copyright © 2005 John Wiley & Sons, Ltd. [source] Simulation technique for wave generationINTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN BIOMEDICAL ENGINEERING, Issue 5 2003S. Aliabadi Abstract In this paper, we present a new finite element technique for simulation of water waves impacting on floating structures. The emphasis will be on the numerical methods for water wave generation and propagation. In our approach, the governing equations are the Navier,Stokes equations written for two incompressible fluids. An interface function with two distinct values serves as a marker identifying the location of the free-surface. This function is transported throughout the computational domain with a time-dependent advection equation. The stabilized finite element formulations are written and integrated in an arbitrary Lagrangian,Eulerian domain. This allows us to handle the motion of the physical boundaries, such as the wave generator surface by moving the computational nodes. In the mesh-moving scheme, we assume that the computational domain is made of elastic materials. The linear elasticity equations are solved to obtain the displacements for each computational node. The numerical examples include 3D wave generation and wave breaking as they approach the coast, and the waves impacting on near-shore support columns. Copyright © 2003 John Wiley & Sons, Ltd. [source] PC cluster parallel finite element analysis of sloshing problem by earthquake using different network environmentsINTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN BIOMEDICAL ENGINEERING, Issue 10 2002Kazuo Kashiyama Abstract This paper presents a parallel finite element method for the analysis of the sloshing problem caused by earthquakes. The incompressible Navier,Stokes equation based on Arbitrary Lagrangian,Eulerian description is used as the governing equation. The SUPG/PSPG formulation is employed to improve the numerical stability and the accuracy. Parallel implementation of the unstructured grid based formulation was carried out on a PC cluster. The present method was applied to analyse the sloshing problem of a rectangular tank and an actual reservoir. The effect of parallelization on the efficiency of the computations was examined using a number of different network environments. Copyright © 2002 John Wiley & Sons, Ltd. [source] ODDLS: A new unstructured mesh finite element method for the analysis of free surface flow problemsINTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING, Issue 9 2008Julio Garcia-Espinosa Abstract This paper introduces a new stabilized finite element method based on the finite calculus (Comput. Methods Appl. Mech. Eng. 1998; 151:233,267) and arbitrary Lagrangian,Eulerian techniques (Comput. Methods Appl. Mech. Eng. 1998; 155:235,249) for the solution to free surface problems. The main innovation of this method is the application of an overlapping domain decomposition concept in the statement of the problem. The aim is to increase the accuracy in the capture of the free surface as well as in the resolution of the governing equations in the interface between the two fluids. Free surface capturing is based on the solution to a level set equation. The Navier,Stokes equations are solved using an iterative monolithic predictor,corrector algorithm (Encyclopedia of Computational Mechanics. Wiley: New York, 2004), where the correction step is based on imposing the divergence-free condition in the velocity field by means of the solution to a scalar equation for the pressure. Examples of application of the ODDLS formulation (for overlapping domain decomposition level set) to the analysis of different free surface flow problems are presented. Copyright © 2008 John Wiley & Sons, Ltd. [source] Stabilized updated Lagrangian corrected SPH for explicit dynamic problemsINTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING, Issue 13 2007Y. Vidal Abstract Smooth particle hydrodynamics with a total Lagrangian formulation are, in general, more robust than finite elements for large distortion problems. Nevertheless, updating the reference configuration may still be necessary in some problems involving extremely large distortions. However, as discussed here, a standard updated formulation suffers the presence of zero-energy modes that are activated and may completely spoil the solution. It is important to note that, unlike an Eulerian formulation, the updated Lagrangian does not present tension instability but only zero-energy modes. Here a stabilization technique is incorporated to the updated formulation to obtain an improved method without any mechanisms and which is capable to solve problems with extremely large distortions. Copyright © 2006 John Wiley & Sons, Ltd. [source] Prediction of the non-ideal detonation performance of commercial explosives using the DeNE and JWL++ codesINTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING, Issue 14 2005S. Esen Abstract The non-ideal detonation performance of two commercial explosives is determined using the DeNE and JWL++ codes. These two codes differ in that DeNE is based on a pseudo-one-dimensional theory which is valid on the central stream-tube and capable of predicting the non-ideal detonation characteristics of commercial explosives as a function of the explosive type, rock properties and blasthole diameter. On the other hand, JWL++ is a hydrocode running in a 2-D arbitrary Lagrangian,Eulerian code with CALE-like properties and can determine the flow properties in all stream lines within the reaction zone. The key flow properties (detonation velocity, pressure, specific volume, extent of reaction and reaction zone length) at the sonic locus on the charge axis have been compared. In general, it is shown that the flow parameters determined using both codes agree well. The pressure contours determined using the JWL++ are analysed in detail for two explosives at 165 mm blastholes confined in limestone and kimberlite with a view to further investigate the explosive/rock interface. The DeNE and JWL++ codes have been validated using the measured in-hole detonation velocity data. Copyright © 2005 John Wiley & Sons, Ltd. [source] Reduced modified quadratures for quadratic membrane finite elementsINTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING, Issue 6 2004Craig S. Long Abstract Reduced integration is frequently used in evaluating the element stiffness matrix of quadratically interpolated finite elements. Typical examples are the serendipity (Q8) and Lagrangian (Q9) membrane finite elements, for which a reduced 2 × 2 Gauss,Legendre integration rule is frequently used, as opposed to full 3 × 3 Gauss,Legendre integration. This ,softens' these element, thereby increasing accuracy, albeit at the introduction of spurious zero energy modes on the element level. This is in general not considered problematic for the ,hourglass' mode common to Q8 and Q9 elements, since this spurious mode is non-communicable. The remaining two zero energy modes occurring in the Q9 element are indeed communicable. However, in topology optimization for instance, conditions may arise where the non-communicable spurious mode associated with the elements becomes activated. To effectively suppress these modes altogether in elements employing quadratic interpolation fields, two modified quadratures are employed herein. For the Q8 and Q9 membrane elements, the respective rules are a five and an eight point rule. As compared to fully integrated elements, the new rules enhance element accuracy due to the introduction of soft, higher-order deformation modes. A number of standard test problems reveal that element accuracy remains comparable to that of the under-integrated counterparts. Copyright © 2004 John Wiley & Sons, Ltd. [source] Numerical analysis of Augmented Lagrangian algorithms in complementary elastoplasticityINTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING, Issue 14 2004L. Contrafatto The main subject of the paper is the investigation of Augmented Lagrangian algorithms and update formulas in the solution of elastoplastic problems. A stress rate formulation for elastoplastic models with internal variables and its finite increment form is employed to state the mechanical problem. In this formulation the Augmented Lagrangian is used to enforce the constraint of plastic admissibility directly on the stresses and thermodynamic forces. This is not a limitation of the Augmented Lagrangian approach, and the same framework can be built on more classical displacement formulations as well. The meaning and the derivation of various first and second order Lagrangian multipliers update formulas and iterative schemes is shown. A new diagonal iteration algorithm and the introduction of a scale factor for the Augmented Lagrangian term are proposed. Numerical examples compare the efficiency of several forms of Augmented Lagrangian algorithms and illustrate the influence of the scale factor and of the penalty parameter. Copyright © 2004 John Wiley & Sons, Ltd. [source] An arbitrary Lagrangian,Eulerian finite element method for finite strain plasticityINTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING, Issue 4 2003Francisco Armero Abstract This paper presents a new arbitrary Lagrangian,Eulerian (ALE) finite element formulation for finite strain plasticity in non-linear solid mechanics. We consider the models of finite strain plasticity defined by the multiplicative decomposition of the deformation gradient in an elastic and a plastic part (F = FeFp), with the stresses given by a hyperelastic relation. In contrast with more classical ALE approaches based on plastic models of the hypoelastic type, the ALE formulation presented herein considers the direct interpolation of the motion of the material with respect to the reference mesh together with the motion of the spatial mesh with respect to this same reference mesh. This aspect is shown to be crucial for a simple treatment of the advection of the plastic internal variables and dynamic variables. In fact, this advection is carried out exactly through a particle tracking in the reference mesh, a calculation that can be accomplished very efficiently with the use of the connectivity graph of the fixed reference mesh. A staggered scheme defined by three steps (the smoothing, the advection and the Lagrangian steps) leads to an efficient method for the solution of the resulting equations. We present several representative numerical simulations that illustrate the performance of the newly proposed methods. Both quasi-static and dynamic conditions are considered in these model examples. Copyright © 2003 John Wiley & Sons, Ltd. [source] Flow-induced vibrations of non-linear cables.INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING, Issue 5 2002Part 1: Models, algorithms Abstract In this paper, we develop governing equations for non-linear cables as well as a formulation for the coupled flow-structure problem. The structure is discretized with second-order accuracy while the flow is discretized using spectral/hp elements in the context of the arbitrary Lagrangian,Eulerian formulation (ALE). Several benchmark problems are considered and the computational implementation is detailed. In the second part of this work large-scale simulation examples are presented. Copyright © 2002 John Wiley & Sons, Ltd. [source] Free vibration of sandwich plates with laminated facesINTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING, Issue 2 2002W. X. Yuan Description is given of the development of a spline finite strip method for predicting the natural frequencies and modes of conventional rectangular sandwich plates. The faceplates are treated as being classically thin and may be of composite laminated construction. The core is modelled as a three-dimensional body. Finite strip stiffness and mass properties are based on a displacement field which represents eight fundamental through-thickness displacements as a series of products of longitudinal B-spline functions and crosswise Lagrangian or Hermitian polynominal shape functions. The solution procedure utilizes the efficient superstrip concept in conjunction with the extended Sturm sequence-bisection approach. A variety of applications of the developed analysis capability is described which demonstrates the nature of the convergence of the finite strip predictions of natural frequencies and the close comparison of these predictions with available results in the literature, and also the use of the capability in parametric studies. Copyright © 2002 John Wiley & Sons, Ltd. [source] Use of the tangent derivative boundary integral equations for the efficient computation of stresses and error indicatorsINTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING, Issue 4 2002K. H. Muci-Küchler Abstract In this work, a new global reanalysis technique for the efficient computation of stresses and error indicators in two-dimensional elastostatic problems is presented. In the context of the boundary element method, the global reanalysis technique can be viewed as a post-processing activity that is carried out once an analysis using Lagrangian elements has been performed. To do the reanalysis, the functional representation for the displacements is changed from Lagrangian to Hermite, introducing the nodal values of the tangential derivatives of those quantities as additional degrees of freedom. Next, assuming that the nodal values of the displacements and the tractions remain practically unchanged from the ones obtained in the analysis using Lagrangian elements, the tangent derivative boundary integral equations are collocated at each functional node in order to determine the additional degrees of freedom that were introduced. Under this scheme, a second system of equations is generated and, once it is solved, the nodal values of the tangential derivatives of the displacements are obtained. This approach gives more accurate results for the stresses at the nodes since it avoids the need to differentiate the shape functions in order to obtain the normal strain in the tangential direction. When compared with the use of Hermite elements, the global reanalysis technique has the attraction that the user does not have to give as input data the additional information required by this type of elements. Another important feature of the proposed approach is that an efficient error indicator for the values of the stresses can also be obtained comparing the values for the stresses obtained through the use of Lagrangian elements and the global reanalysis technique. Copyright © 2001 John Wiley & Sons, Ltd. [source] Remarks on tension instability of Eulerian and Lagrangian corrected smooth particle hydrodynamics (CSPH) methodsINTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING, Issue 11 2001Javier Bonet Abstract The paper discusses the problem of tension instability of particle-based methods such as smooth particle hydrodynamics (SPH) or corrected SPH (CSPH). It is shown that tension instability is a property of a continuum where the stress tensor is isotropic and the value of the pressure is a function of the density or volume ratio. The paper will show that, for this material model, the non-linear continuum equations fail to satisfy the stability condition in the presence of tension. Consequently, any discretization of this continuum will result in negative eigenvalues in the tangent stiffness matrix that will lead to instabilities in the time integration process. An important exception is the 1-D case where the continuum becomes stable but SPH or CSPH can still exhibit negative eigenvalues. The paper will show that these negative eigenvalues can be eliminated if a Lagrangian formulation is used whereby all derivatives are referred to a fixed reference configuration. The resulting formulation maintains the momentum preservation properties of its Eulerian equivalent. Finally a simple 1-D wave propagation example will be used to demonstrate that a stable solution can be obtained using Lagrangian CSPH without the need for any artificial viscosity. Copyright © 2001 John Wiley & Sons, Ltd. [source] An arbitrary Lagrangian,Eulerian method based on the adaptive Riemann solvers for general equations of stateINTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS, Issue 11 2009Baolin Tian Abstract Approximate or exact Riemann solvers play a key role in Godunov-type methods. In this paper, three approximate Riemann solvers, the MFCAV, DKWZ and weak wave approximation method schemes, are investigated through numerical experiments, and their numerical features, such as the resolution for shock and contact waves, are analyzed and compared. Based on the analysis, two new adaptive Riemann solvers for general equations of state are proposed, which can resolve both shock and contact waves well. As a result, an ALE method based on the adaptive Riemann solvers is formulated. A number of numerical experiments show good performance of the adaptive solvers in resolving both shock waves and contact discontinuities. Copyright © 2008 John Wiley & Sons, Ltd. [source] Sigma transformation and ALE formulation for three-dimensional free surface flowsINTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS, Issue 4 2009A. Decoene Abstract In this paper we establish a link between the sigma transformation approach and the arbitrary Lagrangian,Eulerian (ALE) approach. For that purpose we introduce the ALE-sigma (ALES) approach, which consists in an ALE interpretation of the sigma transformation. Taking advantage of this new approach, we propose a general ALES transformation, allowing for a great adaptability of the vertical discretization and therefore overcoming some drawbacks of the classical sigma transformation. Numerical results are presented, showing the advantages of this general coordinate system, as, for example, a better representation of horizontal stratifications. Copyright © 2008 John Wiley & Sons, Ltd. [source] A new modification of the immersed-boundary method for simulating flows with complex moving boundariesINTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS, Issue 11 2006Jian Deng Abstract In this paper, a new immersed-boundary method for simulating flows over complex immersed, moving boundaries is presented. The flow is computed on a fixed Cartesian mesh and the solid boundaries are allowed to move freely through the mesh. The present method is based on a finite-difference approach on a staggered mesh together with a fractional-step method. It must be noted that the immersed boundary is generally not coincident with the position of the solution variables on the grid, therefore, an appropriate strategy is needed to construct a relationship between the curved boundary and the grid points nearby. Furthermore, a momentum forcing is added on the body boundaries and also inside the body to satisfy the no-slip boundary condition. The immersed boundary is represented by a series of interfacial markers, and the markers are also used as Lagrangian forcing points. A linear interpolation is then used to scale the Lagrangian forcing from the interfacial markers to the corresponding grid points nearby. This treatment of the immersed-boundary is used to simulate several problems, which have been validated with previous experimental results in the open literature, verifying the accuracy of the present method. Copyright © 2006 John Wiley & Sons, Ltd. [source] A particle finite element method applied to long wave run-upINTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS, Issue 3 2006J. Birknes Abstract This paper presents a Lagrangian,Eulerian finite element formulation for solving fluid dynamics problems with moving boundaries and employs the method to long wave run-up. The method is based on a set of Lagrangian particles which serve as moving nodes for the finite element mesh. Nodes at the moving shoreline are identified by the alpha shape concept which utilizes the distance from neighbouring nodes in different directions. An efficient triangulation technique is then used for the mesh generation at each time step. In order to validate the numerical method the code has been compared with analytical solutions and a preexisting finite difference model. The main focus of our investigation is to assess the numerical method through simulations of three-dimensional dam break and long wave run-up on curved beaches. Particularly the method is put to test for cases where different shoreline segments connect and produce a computational domain surrounding dry regions. Copyright © 2006 John Wiley & Sons, Ltd. [source] Flow simulation on moving boundary-fitted grids and application to fluid,structure interaction problemsINTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS, Issue 4 2006Martin Engel Abstract We present a method for the parallel numerical simulation of transient three-dimensional fluid,structure interaction problems. Here, we consider the interaction of incompressible flow in the fluid domain and linear elastic deformation in the solid domain. The coupled problem is tackled by an approach based on the classical alternating Schwarz method with non-overlapping subdomains, the subproblems are solved alternatingly and the coupling conditions are realized via the exchange of boundary conditions. The elasticity problem is solved by a standard linear finite element method. A main issue is that the flow solver has to be able to handle time-dependent domains. To this end, we present a technique to solve the incompressible Navier,Stokes equation in three-dimensional domains with moving boundaries. This numerical method is a generalization of a finite volume discretization using curvilinear coordinates to time-dependent coordinate transformations. It corresponds to a discretization of the arbitrary Lagrangian,Eulerian formulation of the Navier,Stokes equations. Here the grid velocity is treated in such a way that the so-called Geometric Conservation Law is implicitly satisfied. Altogether, our approach results in a scheme which is an extension of the well-known MAC-method to a staggered mesh in moving boundary-fitted coordinates which uses grid-dependent velocity components as the primary variables. To validate our method, we present some numerical results which show that second-order convergence in space is obtained on moving grids. Finally, we give the results of a fully coupled fluid,structure interaction problem. It turns out that already a simple explicit coupling with one iteration of the Schwarz method, i.e. one solution of the fluid problem and one solution of the elasticity problem per time step, yields a convergent, simple, yet efficient overall method for fluid,structure interaction problems. Copyright © 2005 John Wiley & Sons, Ltd. [source] Generation of Arbitrary Lagrangian,Eulerian (ALE) velocities, based on monitor functions, for the solution of compressible fluid equationsINTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS, Issue 10-11 2005B. V. Wells Abstract A moving mesh method is outlined based on the use of monitor functions. The method is developed from a weak conservation principle. From this principle a conservation law for the mesh position is derived. Using the Helmholtz decomposition theorem, this conservation law can be converted into an elliptic equation for a mesh velocity potential. The moving mesh method is discretized using standard finite elements. Once the mesh velocities are obtained an arbitrary Lagrangian,Eulerian (ALE) (Journal of Computational Physics 1974; 14:227) fluid solver is used to update the solution on the adaptive mesh. Results are shown for the compressible Euler equations of gas dynamics in one and two spatial dimensions. Two monitor functions are used, the fluid density (which corresponds to a Lagrangian description), and a function which includes the density gradient. A variety of test problems are considered. Copyright © 2005 John Wiley & Sons, Ltd. [source] An implicit edge-based ALE method for the incompressible Navier,Stokes equations,INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS, Issue 3 2003Richard W. Smith Abstract A new finite volume method for the incompressible Navier,Stokes equations, expressed in arbitrary Lagrangian,Eulerian (ALE) form, is presented. The method uses a staggered storage arrangement for the pressure and velocity variables and adopts an edge-based data structure and assembly procedure which is valid for arbitrary n-sided polygonal meshes. Edge formulas are presented for assembling the ALE form of the momentum and pressure equations. An implicit multi-stage time integrator is constructed that is geometrically conservative to the precision of the arithmetic used in the computation. The method is shown to be second-order-accurate in time and space for general time-dependent polygonal meshes. The method is first evaluated using several well-known unsteady incompressible Navier,Stokes problems before being applied to a periodically forced aeroelastic problem and a transient free surface problem. Published in 2003 by John Wiley & Sons, Ltd. [source] Numerical simulation of disperse multiphase flows with an application in power engineeringINTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS, Issue 12 2003K. Bernert Abstract This paper deals with the numerical simulation of two-phase flows based on the solution of the Navier,Stokes equations with a k,,turbulence model for the gas phase and a particle tracking model of the disperse phase fulfilling the framework of the Eulerian,Lagrangian (PSI-Cell) approach. The numerical procedures for the two phases are based on the domain decomposition method applied to a block-structured grid. The complete code is parallelized for computers of MIMD architecture. The paper gives a description of the numerical methods with special attention to the parallelization. Some test calculations demonstrate the performance of the code. The numerical simulation of a flow splitter from the field of power engineering is presented as an example for a real world application of the method. Copyright © 2003 John Wiley & Sons, Ltd. [source] Three-dimensional transient free-surface flow of viscous fluids inside cavities of arbitrary shapeINTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS, Issue 10 2003Kyu-Tae Kim Abstract The three-dimensional transient free-surface flow inside cavities of arbitrary shape is examined in this study. An adaptive (Lagrangian) boundary-element approach is proposed for the general three-dimensional simulation of confined free-surface flow of viscous incompressible fluids. The method is stable as it includes remeshing capabilities of the deforming free-surface, and thus can handle large deformations. A simple algorithm is developed for mesh refinement of the deforming free-surface mesh. Smooth transition between large and small elements is achieved without significant degradation of the aspect ratio of the elements in the mesh. The method is used to determine the flow field and free-surface evolution inside cubic, rectangular and cylindrical containers. These problems illustrate the transient nature of the flow during the mixing process. Surface tension effects are also explored. Copyright © 2003 John Wiley & Sons, Ltd. [source] Interface reconstruction with least-square fit and split Eulerian,Lagrangian advectionINTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS, Issue 3 2003Ruben Scardovelli Abstract Two new volume-of-fluid (VOF) reconstruction algorithms, which are based on a least-square fit technique, are presented. Their performance is tested for several standard shapes and is compared to a few other VOF/PLIC reconstruction techniques, showing in general a better convergence rate. The geometric nature of Lagrangian and Eulerian split advection algorithms is investigated in detail and a new mixed split Eulerian implicit,Lagrangian explicit (EI,LE) scheme is presented. This method conserves the mass to machine error, performs better than split Eulerian and Lagrangian algorithms, and it is only slightly worse than unsplit schemes. However, the combination of the interface reconstruction with the least-square fit and its advection with the EI,LE scheme appears superior to other existing approaches. Copyright © 2003 John Wiley & Sons, Ltd. [source] | |||||||||||||||||||