Forced pitch motion of wind turbines

NASA Astrophysics Data System (ADS)

Leble, V.; Barakos, G.

2016-09-01

The possibility of a wind turbine entering vortex ring state during pitching oscillations is explored in this paper. The aerodynamic performance of the rotor was computed using the Helicopter Multi-Block flow solver. This code solves the Navier-Stokes equations in integral form using the arbitrary Lagrangian-Eulerian formulation for time-dependent domains with moving boundaries. A 10-MW wind turbine was put to perform yawing and pitching oscillations suggesting the partial vortex ring state during pitching motion. The results also show the strong effect of the frequency and amplitude of oscillations on the wind turbine performance.

An Operator-Integration-Factor Splitting (OIFS) method for Incompressible Flows in Moving Domains

DOE Office of Scientific and Technical Information (OSTI.GOV)

Patel, Saumil S.; Fischer, Paul F.; Min, Misun

In this paper, we present a characteristic-based numerical procedure for simulating incompressible flows in domains with moving boundaries. Our approach utilizes an operator-integration-factor splitting technique to help produce an effcient and stable numerical scheme. Using the spectral element method and an arbitrary Lagrangian-Eulerian formulation, we investigate flows where the convective acceleration effects are non-negligible. Several examples, ranging from laminar to turbulent flows, are considered. Comparisons with a standard, semi-implicit time-stepping procedure illustrate the improved performance of the scheme.

CELFE: Coupled Eulerian-Lagrangian Finite Element program for high velocity impact. Part 1: Theory and formulation. [hydroelasto-viscoplastic model

NASA Technical Reports Server (NTRS)

Lee, C. H.

1978-01-01

A 3-D finite element program capable of simulating the dynamic behavior in the vicinity of the impact point, together with predicting the dynamic response in the remaining part of the structural component subjected to high velocity impact is discussed. The finite algorithm is formulated in a general moving coordinate system. In the vicinity of the impact point contained by a moving failure front, the relative velocity of the coordinate system will approach the material particle velocity. The dynamic behavior inside the region is described by Eulerian formulation based on a hydroelasto-viscoplastic model. The failure front which can be regarded as the boundary of the impact zone is described by a transition layer. The layer changes the representation from the Eulerian mode to the Lagrangian mode outside the failure front by varying the relative velocity of the coordinate system to zero. The dynamic response in the remaining part of the structure described by the Lagrangian formulation is treated using advanced structural analysis. An interfacing algorithm for coupling CELFE with NASTRAN is constructed to provide computational capabilities for large structures.

Constrained optimization framework for interface-aware sub-scale dynamics models for voids closure in Lagrangian hydrodynamics

DOE PAGES

Barlow, Andrew; Klima, Matej; Shashkov, Mikhail

2018-04-02

In hydrocodes, voids are used to represent vacuum and model free boundaries between vacuum and real materials. We give a systematic description of a new treatment of void closure in the framework of the multimaterial arbitrary Lagrangian–Eulerian (ALE) methods. This includes a new formulation of the interface-aware sub-scale-dynamics (IA-SSD) closure model for multimaterial cells with voids, which is used in the Lagrangian stage of our indirect ALE scheme. The results of the comprehensive testing of the new model are presented for one- and two-dimensional multimaterial calculations in the presence of voids. Finally, we also present a sneak peek of amore » realistic shaped charge calculation in the presence of voids and solids.« less

Constrained optimization framework for interface-aware sub-scale dynamics models for voids closure in Lagrangian hydrodynamics

DOE Office of Scientific and Technical Information (OSTI.GOV)

Barlow, Andrew; Klima, Matej; Shashkov, Mikhail

In hydrocodes, voids are used to represent vacuum and model free boundaries between vacuum and real materials. We give a systematic description of a new treatment of void closure in the framework of the multimaterial arbitrary Lagrangian–Eulerian (ALE) methods. This includes a new formulation of the interface-aware sub-scale-dynamics (IA-SSD) closure model for multimaterial cells with voids, which is used in the Lagrangian stage of our indirect ALE scheme. The results of the comprehensive testing of the new model are presented for one- and two-dimensional multimaterial calculations in the presence of voids. Finally, we also present a sneak peek of amore » realistic shaped charge calculation in the presence of voids and solids.« less

A hybridizable discontinuous Galerkin method for modeling fluid-structure interaction

NASA Astrophysics Data System (ADS)

Sheldon, Jason P.; Miller, Scott T.; Pitt, Jonathan S.

2016-12-01

This work presents a novel application of the hybridizable discontinuous Galerkin (HDG) finite element method to the multi-physics simulation of coupled fluid-structure interaction (FSI) problems. Recent applications of the HDG method have primarily been for single-physics problems including both solids and fluids, which are necessary building blocks for FSI modeling. Utilizing these established models, HDG formulations for linear elastostatics, a nonlinear elastodynamic model, and arbitrary Lagrangian-Eulerian Navier-Stokes are derived. The elasticity formulations are written in a Lagrangian reference frame, with the nonlinear formulation restricted to hyperelastic materials. With these individual solid and fluid formulations, the remaining challenge in FSI modeling is coupling together their disparate mathematics on the fluid-solid interface. This coupling is presented, along with the resultant HDG FSI formulation. Verification of the component models, through the method of manufactured solutions, is performed and each model is shown to converge at the expected rate. The individual components, along with the complete FSI model, are then compared to the benchmark problems proposed by Turek and Hron [1]. The solutions from the HDG formulation presented in this work trend towards the benchmark as the spatial polynomial order and the temporal order of integration are increased.

Modeling shockwaves and impact phenomena with Eulerian peridynamics

DOE PAGES

Silling, Stewart A.; Parks, Michael L.; Kamm, James R.; ...

2017-05-09

Most previous development of the peridynamic theory has assumed a Lagrangian formulation, in which the material model refers to an undeformed reference configuration. Here, an Eulerian form of material modeling is developed, in which bond forces depend only on the positions of material points in the deformed configuration. The formulation is consistent with the thermodynamic form of the peridynamic model and is derivable from a suitable expression for the free energy of a material. We show that the resulting formulation of peridynamic material models can be used to simulate strong shock waves and fluid response in which very large deformationsmore » make the Lagrangian form unsuitable. The Eulerian capability is demonstrated in numerical simulations of ejecta from a wavy free surface on a metal subjected to strong shock wave loading. The Eulerian and Lagrangian contributions to bond force can be combined in a single material model, allowing strength and fracture under tensile or shear loading to be modeled consistently with high compressive stresses. Furthermore, we demonstrate this capability in numerical simulation of bird strike against an aircraft, in which both tensile fracture and high pressure response are important.« less

Modeling shockwaves and impact phenomena with Eulerian peridynamics

DOE Office of Scientific and Technical Information (OSTI.GOV)

Silling, Stewart A.; Parks, Michael L.; Kamm, James R.

Most previous development of the peridynamic theory has assumed a Lagrangian formulation, in which the material model refers to an undeformed reference configuration. Here, an Eulerian form of material modeling is developed, in which bond forces depend only on the positions of material points in the deformed configuration. The formulation is consistent with the thermodynamic form of the peridynamic model and is derivable from a suitable expression for the free energy of a material. We show that the resulting formulation of peridynamic material models can be used to simulate strong shock waves and fluid response in which very large deformationsmore » make the Lagrangian form unsuitable. The Eulerian capability is demonstrated in numerical simulations of ejecta from a wavy free surface on a metal subjected to strong shock wave loading. The Eulerian and Lagrangian contributions to bond force can be combined in a single material model, allowing strength and fracture under tensile or shear loading to be modeled consistently with high compressive stresses. Furthermore, we demonstrate this capability in numerical simulation of bird strike against an aircraft, in which both tensile fracture and high pressure response are important.« less

Verification of transport equations in a general purpose commercial CFD code.

NASA Astrophysics Data System (ADS)

Melot, Matthieu; Nennemann, Bernd; Deschênes, Claire

2016-11-01

In this paper, the Verification and Validation methodology is presented. This method aims to increase the reliability and the trust that can be placed into complex CFD simulations. The first step of this methodology, the code verification is presented in greater details. The CFD transport equations in steady state, transient and Arbitrary Eulerian Lagrangian (ALE, used for transient moving mesh) formulations in Ansys CFX are verified. It is shown that the expected spatial and temporal order of convergence are achieved for the steady state and the transient formulations. Unfortunately this is not completely the case for the ALE formulation. As for a lot of other commercial and in-house CFD codes, the temporal convergence of the velocity is limited to a first order where a second order would have been expected.

Adjoint of the global Eulerian-Lagrangian coupled atmospheric transport model (A-GELCA v1.0): development and validation

NASA Astrophysics Data System (ADS)

Belikov, Dmitry A.; Maksyutov, Shamil; Yaremchuk, Alexey; Ganshin, Alexander; Kaminski, Thomas; Blessing, Simon; Sasakawa, Motoki; Gomez-Pelaez, Angel J.; Starchenko, Alexander

2016-02-01

We present the development of the Adjoint of the Global Eulerian-Lagrangian Coupled Atmospheric (A-GELCA) model that consists of the National Institute for Environmental Studies (NIES) model as an Eulerian three-dimensional transport model (TM), and FLEXPART (FLEXible PARTicle dispersion model) as the Lagrangian Particle Dispersion Model (LPDM). The forward tangent linear and adjoint components of the Eulerian model were constructed directly from the original NIES TM code using an automatic differentiation tool known as TAF (Transformation of Algorithms in Fortran; http://www.FastOpt.com, with additional manual pre- and post-processing aimed at improving transparency and clarity of the code and optimizing the performance of the computing, including MPI (Message Passing Interface). The Lagrangian component did not require any code modification, as LPDMs are self-adjoint and track a significant number of particles backward in time in order to calculate the sensitivity of the observations to the neighboring emission areas. The constructed Eulerian adjoint was coupled with the Lagrangian component at a time boundary in the global domain. The simulations presented in this work were performed using the A-GELCA model in forward and adjoint modes. The forward simulation shows that the coupled model improves reproduction of the seasonal cycle and short-term variability of CO2. Mean bias and standard deviation for five of the six Siberian sites considered decrease roughly by 1 ppm when using the coupled model. The adjoint of the Eulerian model was shown, through several numerical tests, to be very accurate (within machine epsilon with mismatch around to ±6 e-14) compared to direct forward sensitivity calculations. The developed adjoint of the coupled model combines the flux conservation and stability of an Eulerian discrete adjoint formulation with the flexibility, accuracy, and high resolution of a Lagrangian backward trajectory formulation. A-GELCA will be incorporated into a variational inversion system designed to optimize surface fluxes of greenhouse gases.

Implementing a Loosely Coupled Fluid Structure Interaction Finite Element Model in PHASTA

NASA Astrophysics Data System (ADS)

Pope, David

Fluid Structure Interaction problems are an important multi-physics phenomenon in the design of aerospace vehicles and other engineering applications. A variety of computational fluid dynamics solvers capable of resolving the fluid dynamics exist. PHASTA is one such computational fluid dynamics solver. Enhancing the capability of PHASTA to resolve Fluid-Structure Interaction first requires implementing a structural dynamics solver. The implementation also requires a correction of the mesh used to solve the fluid equations to account for the deformation of the structure. This results in mesh motion and causes the need for an Arbitrary Lagrangian-Eulerian modification to the fluid dynamics equations currently implemented in PHASTA. With the implementation of both structural dynamics physics, mesh correction, and the Arbitrary Lagrangian-Eulerian modification of the fluid dynamics equations, PHASTA is made capable of solving Fluid-Structure Interaction problems.

Parallel computing using a Lagrangian formulation

NASA Technical Reports Server (NTRS)

Liou, May-Fun; Loh, Ching Yuen

1991-01-01

A new Lagrangian formulation of the Euler equation is adopted for the calculation of 2-D supersonic steady flow. The Lagrangian formulation represents the inherent parallelism of the flow field better than the common Eulerian formulation and offers a competitive alternative on parallel computers. The implementation of the Lagrangian formulation on the Thinking Machines Corporation CM-2 Computer is described. The program uses a finite volume, first-order Godunov scheme and exhibits high accuracy in dealing with multidimensional discontinuities (slip-line and shock). By using this formulation, a better than six times speed-up was achieved on a 8192-processor CM-2 over a single processor of a CRAY-2.

Parallel computing using a Lagrangian formulation

NASA Technical Reports Server (NTRS)

Liou, May-Fun; Loh, Ching-Yuen

1992-01-01

This paper adopts a new Lagrangian formulation of the Euler equation for the calculation of two dimensional supersonic steady flow. The Lagrangian formulation represents the inherent parallelism of the flow field better than the common Eulerian formulation and offers a competitive alternative on parallel computers. The implementation of the Lagrangian formulation on the Thinking Machines Corporation CM-2 Computer is described. The program uses a finite volume, first-order Godunov scheme and exhibits high accuracy in dealing with multidimensional discontinuities (slip-line and shock). By using this formulation, we have achieved better than six times speed-up on a 8192-processor CM-2 over a single processor of a CRAY-2.

A hybridized discontinuous Galerkin framework for high-order particle-mesh operator splitting of the incompressible Navier-Stokes equations

NASA Astrophysics Data System (ADS)

Maljaars, Jakob M.; Labeur, Robert Jan; Möller, Matthias

2018-04-01

A generic particle-mesh method using a hybridized discontinuous Galerkin (HDG) framework is presented and validated for the solution of the incompressible Navier-Stokes equations. Building upon particle-in-cell concepts, the method is formulated in terms of an operator splitting technique in which Lagrangian particles are used to discretize an advection operator, and an Eulerian mesh-based HDG method is employed for the constitutive modeling to account for the inter-particle interactions. Key to the method is the variational framework provided by the HDG method. This allows to formulate the projections between the Lagrangian particle space and the Eulerian finite element space in terms of local (i.e. cellwise) ℓ2-projections efficiently. Furthermore, exploiting the HDG framework for solving the constitutive equations results in velocity fields which excellently approach the incompressibility constraint in a local sense. By advecting the particles through these velocity fields, the particle distribution remains uniform over time, obviating the need for additional quality control. The presented methodology allows for a straightforward extension to arbitrary-order spatial accuracy on general meshes. A range of numerical examples shows that optimal convergence rates are obtained in space and, given the particular time stepping strategy, second-order accuracy is obtained in time. The model capabilities are further demonstrated by presenting results for the flow over a backward facing step and for the flow around a cylinder.

Lagrangian Perturbation Approach to the Formation of Large-scale Structure

NASA Astrophysics Data System (ADS)

Buchert, Thomas

The present lecture notes address three columns on which the Lagrangian perturbation approach to cosmological dynamics is based: 1. the formulation of a Lagrangian theory of self-gravitating flows in which the dynamics is described in terms of a single field variable; 2. the procedure, how to obtain the dynamics of Eulerian fields from the Lagrangian picture, and 3. a precise definition of a Newtonian cosmology framework in which Lagrangian perturbation solutions can be studied. While the first is a discussion of the basic equations obtained by transforming the Eulerian evolution and field equations to the Lagrangian picture, the second exemplifies how the Lagrangian theory determines the evolution of Eulerian fields including kinematical variables like expansion, vorticity, as well as the shear and tidal tensors. The third column is based on a specification of initial and boundary conditions, and in particular on the identification of the average flow of an inhomogeneous cosmology with a `Hubble-flow'. Here, we also look at the limits of the Lagrangian perturbation approach as inferred from comparisons with N-body simulations and illustrate some striking properties of the solutions.

Full-Scale Direct Numerical Simulation of Two- and Three-Dimensional Instabilities and Rivulet Formulation in Heated Falling Films

NASA Technical Reports Server (NTRS)

Krishnamoorthy, S.; Ramaswamy, B.; Joo, S. W.

1995-01-01

A thin film draining on an inclined plate has been studied numerically using finite element method. Three-dimensional governing equations of continuity, momentum and energy with a moving boundary are integrated in an arbitrary Lagrangian Eulerian frame of reference. Kinematic equation is solved to precisely update interface location. Rivulet formation based on instability mechanism has been simulated using full-scale computation. Comparisons with long-wave theory are made to validate the numerical scheme. Detailed analysis of two- and three-dimensional nonlinear wave formation and spontaneous rupture forming rivulets under the influence of combined thermocapillary and surface-wave instabilities is performed.

Numerical simulation of fluid flow through simplified blade cascade with prescribed harmonic motion using discontinuous Galerkin method

NASA Astrophysics Data System (ADS)

Vimmr, Jan; Bublík, Ondřej; Prausová, Helena; Hála, Jindřich; Pešek, Luděk

2018-06-01

This paper deals with a numerical simulation of compressible viscous fluid flow around three flat plates with prescribed harmonic motion. This arrangement presents a simplified blade cascade with forward wave motion. The aim of this simulation is to determine the aerodynamic forces acting on the flat plates. The mathematical model describing this problem is formed by Favre-averaged system of Navier-Stokes equations in arbitrary Lagrangian-Eulerian (ALE) formulation completed by one-equation Spalart-Allmaras turbulence model. The simulation was performed using the developed in-house CFD software based on discontinuous Galerkin method, which offers high order of accuracy.

Eulerian-Lagrangian Simulations of Transonic Flutter Instabilities

NASA Technical Reports Server (NTRS)

Bendiksen, Oddvar O.

1994-01-01

This paper presents an overview of recent applications of Eulerian-Lagrangian computational schemes in simulating transonic flutter instabilities. This approach, the fluid-structure system is treated as a single continuum dynamics problem, by switching from an Eulerian to a Lagrangian formulation at the fluid-structure boundary. This computational approach effectively eliminates the phase integration errors associated with previous methods, where the fluid and structure are integrated sequentially using different schemes. The formulation is based on Hamilton's Principle in mixed coordinates, and both finite volume and finite element discretization schemes are considered. Results from numerical simulations of transonic flutter instabilities are presented for isolated wings, thin panels, and turbomachinery blades. The results suggest that the method is capable of reproducing the energy exchange between the fluid and the structure with significantly less error than existing methods. Localized flutter modes and panel flutter modes involving traveling waves can also be simulated effectively with no a priori knowledge of the type of instability involved.

Hyperviscosity for unstructured ALE meshes

NASA Astrophysics Data System (ADS)

Cook, Andrew W.; Ulitsky, Mark S.; Miller, Douglas S.

2013-01-01

An artificial viscosity, originally designed for Eulerian schemes, is adapted for use in arbitrary Lagrangian-Eulerian simulations. Changes to the Eulerian model (dubbed 'hyperviscosity') are discussed, which enable it to work within a Lagrangian framework. New features include a velocity-weighted grid scale and a generalised filtering procedure, applicable to either structured or unstructured grids. The model employs an artificial shear viscosity for treating small-scale vorticity and an artificial bulk viscosity for shock capturing. The model is based on the Navier-Stokes form of the viscous stress tensor, including the diagonal rate-of-expansion tensor. A second-order version of the model is presented, in which Laplacian operators act on the velocity divergence and the grid-weighted strain-rate magnitude to ensure that the velocity field remains smooth at the grid scale. Unlike sound-speed-based artificial viscosities, the hyperviscosity model is compatible with the low Mach number limit. The new model outperforms a commonly used Lagrangian artificial viscosity on a variety of test problems.

An adaptive reconstruction for Lagrangian, direct-forcing, immersed-boundary methods

NASA Astrophysics Data System (ADS)

Posa, Antonio; Vanella, Marcos; Balaras, Elias

2017-12-01

Lagrangian, direct-forcing, immersed boundary (IB) methods have been receiving increased attention due to their robustness in complex fluid-structure interaction problems. They are very sensitive, however, on the selection of the Lagrangian grid, which is typically used to define a solid or flexible body immersed in a fluid flow. In the present work we propose a cost-efficient solution to this problem without compromising accuracy. Central to our approach is the use of isoparametric mapping to bridge the relative resolution requirements of Lagrangian IB, and Eulerian grids. With this approach, the density of surface Lagrangian markers, which is essential to properly enforce boundary conditions, is adapted dynamically based on the characteristics of the underlying Eulerian grid. The markers are not stored and the Lagrangian data-structure is not modified. The proposed scheme is implemented in the framework of a moving least squares reconstruction formulation, but it can be adapted to any Lagrangian, direct-forcing formulation. The accuracy and robustness of the approach is demonstrated in a variety of test cases of increasing complexity.

Drekar v.2.0

DOE Office of Scientific and Technical Information (OSTI.GOV)

Seefeldt, Ben; Sondak, David; Hensinger, David M.

Drekar is an application code that solves partial differential equations for fluids that can be optionally coupled to electromagnetics. Drekar solves low-mach compressible and incompressible computational fluid dynamics (CFD), compressible and incompressible resistive magnetohydrodynamics (MHD), and multiple species plasmas interacting with electromagnetic fields. Drekar discretization technology includes continuous and discontinuous finite element formulations, stabilized finite element formulations, mixed integration finite element bases (nodal, edge, face, volume) and an initial arbitrary Lagrangian Eulerian (ALE) capability. Drekar contains the implementation of the discretized physics and leverages the open source Trilinos project for both parallel solver capabilities and general finite element discretization tools.more » The code will be released open source under a BSD license. The code is used for fundamental research for simulation of fluids and plasmas on high performance computing environments.« less

An arbitrary Lagrangian–Eulerian finite element formulation for a poroelasticity problem stemming from mixture theory

DOE PAGES

Costanzo, Francesco; Miller, Scott T.

2017-05-22

In this paper, a finite element formulation is developed for a poroelastic medium consisting of an incompressible hyperelastic skeleton saturated by an incompressible fluid. The governing equations stem from mixture theory and the application is motivated by the study of interstitial fluid flow in brain tissue. The formulation is based on the adoption of an arbitrary Lagrangian–Eulerian (ALE) perspective. We focus on a flow regime in which inertia forces are negligible. Finally, the stability and convergence of the formulation is discussed, and numerical results demonstrate agreement with the theory.

An arbitrary Lagrangian–Eulerian finite element formulation for a poroelasticity problem stemming from mixture theory

DOE Office of Scientific and Technical Information (OSTI.GOV)

Costanzo, Francesco; Miller, Scott T.

In this paper, a finite element formulation is developed for a poroelastic medium consisting of an incompressible hyperelastic skeleton saturated by an incompressible fluid. The governing equations stem from mixture theory and the application is motivated by the study of interstitial fluid flow in brain tissue. The formulation is based on the adoption of an arbitrary Lagrangian–Eulerian (ALE) perspective. We focus on a flow regime in which inertia forces are negligible. Finally, the stability and convergence of the formulation is discussed, and numerical results demonstrate agreement with the theory.

Lagrangian continuum dynamics in ALEGRA.

DOE Office of Scientific and Technical Information (OSTI.GOV)

Wong, Michael K. W.; Love, Edward

Alegra is an ALE (Arbitrary Lagrangian-Eulerian) multi-material finite element code that emphasizes large deformations and strong shock physics. The Lagrangian continuum dynamics package in Alegra uses a Galerkin finite element spatial discretization and an explicit central-difference stepping method in time. The goal of this report is to describe in detail the characteristics of this algorithm, including the conservation and stability properties. The details provided should help both researchers and analysts understand the underlying theory and numerical implementation of the Alegra continuum hydrodynamics algorithm.

Eulerian-Lagrangian solution of the convection-dispersion equation in natural coordinates

USGS Publications Warehouse

Cheng, Ralph T.; Casulli, Vincenzo; Milford, S. Nevil

1984-01-01

The vast majority of numerical investigations of transport phenomena use an Eulerian formulation for the convenience that the computational grids are fixed in space. An Eulerian-Lagrangian method (ELM) of solution for the convection-dispersion equation is discussed and analyzed. The ELM uses the Lagrangian concept in an Eulerian computational grid system. The values of the dependent variable off the grid are calculated by interpolation. When a linear interpolation is used, the method is a slight improvement over the upwind difference method. At this level of approximation both the ELM and the upwind difference method suffer from large numerical dispersion. However, if second-order Lagrangian polynomials are used in the interpolation, the ELM is proven to be free of artificial numerical dispersion for the convection-dispersion equation. The concept of the ELM is extended for treatment of anisotropic dispersion in natural coordinates. In this approach the anisotropic properties of dispersion can be conveniently related to the properties of the flow field. Several numerical examples are given to further substantiate the results of the present analysis.

Discrete-element simulation of sea-ice mechanics: Contact mechanics and granular jamming

NASA Astrophysics Data System (ADS)

Damsgaard, A.; Adcroft, A.; Sergienko, O. V.; Stern, A. A.

2017-12-01

Lagrangian models of sea-ice dynamics offer several advantages to Eulerian continuum methods. Spatial discretization on the ice-floe scale is natural for Lagrangian models, which additionally offer the convenience of being able to handle arbitrary sea-ice concentrations. This is likely to improve model performance in ice-marginal zones with strong advection. Furthermore, phase transitions in granular rheology around the jamming limit, such as observed when sea ice moves through geometric confinements, includes sharp thresholds in effective viscosity which are typically ignored in Eulerian models. Granular jamming is a stochastic process dependent on having the right grains in the right place at the right time, and the jamming likelihood over time can be described by a probabilistic model. Difficult to parameterize in continuum formulations, jamming occurs naturally in dense granular systems simulated in a Lagrangian framework, and is a very relevant process controlling sea-ice transport through narrow straits. We construct a flexible discrete-element framework for simulating Lagrangian sea-ice dynamics at the ice-floe scale, forced by ocean and atmosphere velocity fields. Using this framework, we demonstrate that frictionless contact models based on compressive stiffness alone are unlikely to jam, and describe two different approaches based on friction and tensile strength which both result in increased bulk shear strength of the granular assemblage. The frictionless but cohesive contact model, with certain tensile strength values, can display jamming behavior which on the large scale is very similar to a more complex and realistic model with contact friction and ice-floe rotation.

Modeling of combustion processes of stick propellants via combined Eulerian-Lagrangian approach

NASA Technical Reports Server (NTRS)

Kuo, K. K.; Hsieh, K. C.; Athavale, M. M.

1988-01-01

This research is motivated by the improved ballistic performance of large-caliber guns using stick propellant charges. A comprehensive theoretical model for predicting the flame spreading, combustion, and grain deformation phenomena of long, unslotted stick propellants is presented. The formulation is based upon a combined Eulerian-Lagrangian approach to simulate special characteristics of the two phase combustion process in a cartridge loaded with a bundle of sticks. The model considers five separate regions consisting of the internal perforation, the solid phase, the external interstitial gas phase, and two lumped parameter regions at either end of the stick bundle. For the external gas phase region, a set of transient one-dimensional fluid-dynamic equations using the Eulerian approach is obtained; governing equations for the stick propellants are formulated using the Lagrangian approach. The motion of a representative stick is derived by considering the forces acting on the entire propellant stick. The instantaneous temperature and stress fields in the stick propellant are modeled by considering the transient axisymmetric heat conduction equation and dynamic structural analysis.

Well-balanced Arbitrary-Lagrangian-Eulerian finite volume schemes on moving nonconforming meshes for the Euler equations of gas dynamics with gravity

NASA Astrophysics Data System (ADS)

Gaburro, Elena; Castro, Manuel J.; Dumbser, Michael

2018-06-01

In this work, we present a novel second-order accurate well-balanced arbitrary Lagrangian-Eulerian (ALE) finite volume scheme on moving nonconforming meshes for the Euler equations of compressible gas dynamics with gravity in cylindrical coordinates. The main feature of the proposed algorithm is the capability of preserving many of the physical properties of the system exactly also on the discrete level: besides being conservative for mass, momentum and total energy, also any known steady equilibrium between pressure gradient, centrifugal force, and gravity force can be exactly maintained up to machine precision. Perturbations around such equilibrium solutions are resolved with high accuracy and with minimal dissipation on moving contact discontinuities even for very long computational times. This is achieved by the novel combination of well-balanced path-conservative finite volume schemes, which are expressly designed to deal with source terms written via non-conservative products, with ALE schemes on moving grids, which exhibit only very little numerical dissipation on moving contact waves. In particular, we have formulated a new HLL-type and a novel Osher-type flux that are both able to guarantee the well balancing in a gas cloud rotating around a central object. Moreover, to maintain a high level of quality of the moving mesh, we have adopted a nonconforming treatment of the sliding interfaces that appear due to the differential rotation. A large set of numerical tests has been carried out in order to check the accuracy of the method close and far away from the equilibrium, both, in one- and two-space dimensions.

Finite Element Simulation of a Space Shuttle Solid Rocket Booster Aft Skirt Splashdown Using an Arbitrary Lagrangian-Eulerian Approach

NASA Astrophysics Data System (ADS)

Melis, Matthew E.

2003-01-01

Explicit finite element techniques employing an Arbitrary Lagrangian-Eulerian (ALE) methodology, within the transient dynamic code LS-DYNA, are used to predict splashdown loads on a proposed replacement/upgrade of the hydrazine tanks on the thrust vector control system housed within the aft skirt of a Space Shuttle Solid Rocket Booster. Two preliminary studies are performed prior to the full aft skirt analysis: An analysis of the proposed tank impacting water without supporting aft skirt structure, and an analysis of space capsule water drop tests conducted at NASA's Langley Research Center. Results from the preliminary studies provide confidence that useful predictions can be made by applying the ALE methodology to a detailed analysis of a 26-degree section of the skirt with proposed tank attached. Results for all three studies are presented and compared to limited experimental data. The challenges of using the LS-DYNA ALE capability for this type of analysis are discussed.

Finite Element Simulation of a Space Shuttle Solid Rocket Booster Aft Skirt Splashdown Using an Arbitrary Lagrangian-eulerian Approach

NASA Technical Reports Server (NTRS)

Melis, Matthew E.

2003-01-01

Explicit finite element techniques employing an Arbitrary Lagrangian-Eulerian (ALE) methodology, within the transient dynamic code LS-DYNA, are used to predict splashdown loads on a proposed replacement/upgrade of the hydrazine tanks on the thrust vector control system housed within the aft skirt of a Space Shuttle Solid Rocket Booster. Two preliminary studies are performed prior to the full aft skirt analysis: An analysis of the proposed tank impacting water without supporting aft skirt structure, and an analysis of space capsule water drop tests conducted at NASA's Langley Research Center. Results from the preliminary studies provide confidence that useful predictions can be made by applying the ALE methodology to a detailed analysis of a 26-degree section of the skirt with proposed tank attached. Results for all three studies are presented and compared to limited experimental data. The challenges of using the LS-DYNA ALE capability for this type of analysis are discussed.

Vaporization of irradiated droplets

NASA Astrophysics Data System (ADS)

Armstrong, R. L.; O'Rourke, P. J.; Zardecki, A.

1986-11-01

The vaporization of a spherically symmetric liquid droplet subject to a high-intensity laser flux is investigated on the basis of a hydrodynamic description of the system composed of the vapor and ambient gas. In the limit of the convective vaporization, the boundary conditions at the fluid-gas interface are formulated by using the notion of a Knudsen layer in which translational equilibrium is established. This leads to approximate jump conditions at the interface. For homogeneous energy deposition, the hydrodynamic equations are solved numerically with the aid of the CON1D computer code (``CON1D: A computer program for calculating spherically symmetric droplet combustion,'' Los Alamos National Laboratory Report No. LA-10269-MS, December, 1984), based on the implict continuous-fluid Eulerian (ICE) [J. Comput. Phys. 8, 197 (1971)] and arbitrary Lagrangian-Eulerian (ALE) [J. Comput. Phys. 14, 1227 (1974)] numerical mehtods. The solutions exhibit the existence of two shock waves propagating in opposite directions with respect to the contact discontinuity surface that separates the ambient gas and vapor.

Coupling fluid-structure interaction with phase-field fracture

NASA Astrophysics Data System (ADS)

Wick, Thomas

2016-12-01

In this work, a concept for coupling fluid-structure interaction with brittle fracture in elasticity is proposed. The fluid-structure interaction problem is modeled in terms of the arbitrary Lagrangian-Eulerian technique and couples the isothermal, incompressible Navier-Stokes equations with nonlinear elastodynamics using the Saint-Venant Kirchhoff solid model. The brittle fracture model is based on a phase-field approach for cracks in elasticity and pressurized elastic solids. In order to derive a common framework, the phase-field approach is re-formulated in Lagrangian coordinates to combine it with fluid-structure interaction. A crack irreversibility condition, that is mathematically characterized as an inequality constraint in time, is enforced with the help of an augmented Lagrangian iteration. The resulting problem is highly nonlinear and solved with a modified Newton method (e.g., error-oriented) that specifically allows for a temporary increase of the residuals. The proposed framework is substantiated with several numerical tests. In these examples, computational stability in space and time is shown for several goal functionals, which demonstrates reliability of numerical modeling and algorithmic techniques. But also current limitations such as the necessity of using solid damping are addressed.

Asymptotic-preserving Lagrangian approach for modeling anisotropic transport in magnetized plasmas for arbitrary magnetic fields

NASA Astrophysics Data System (ADS)

Chacon, Luis; Del-Castillo-Negrete, Diego; Hauck, Cory

2012-10-01

Modeling electron transport in magnetized plasmas is extremely challenging due to the extreme anisotropy between parallel (to the magnetic field) and perpendicular directions (χ/χ˜10^10 in fusion plasmas). Recently, a Lagrangian Green's function approach, developed for the purely parallel transport case,footnotetextD. del-Castillo-Negrete, L. Chac'on, PRL, 106, 195004 (2011)^,footnotetextD. del-Castillo-Negrete, L. Chac'on, Phys. Plasmas, 19, 056112 (2012) has been extended to the anisotropic transport case in the tokamak-ordering limit with constant density.footnotetextL. Chac'on, D. del-Castillo-Negrete, C. Hauck, JCP, submitted (2012) An operator-split algorithm is proposed that allows one to treat Eulerian and Lagrangian components separately. The approach is shown to feature bounded numerical errors for arbitrary χ/χ ratios, which renders it asymptotic-preserving. In this poster, we will present the generalization of the Lagrangian approach to arbitrary magnetic fields. We will demonstrate the potential of the approach with various challenging configurations, including the case of transport across a magnetic island in cylindrical geometry.

Edge remap for solids

DOE Office of Scientific and Technical Information (OSTI.GOV)

Kamm, James R.; Love, Edward; Robinson, Allen C.

We review the edge element formulation for describing the kinematics of hyperelastic solids. This approach is used to frame the problem of remapping the inverse deformation gradient for Arbitrary Lagrangian-Eulerian (ALE) simulations of solid dynamics. For hyperelastic materials, the stress state is completely determined by the deformation gradient, so remapping this quantity effectively updates the stress state of the material. A method, inspired by the constrained transport remap in electromagnetics, is reviewed, according to which the zero-curl constraint on the inverse deformation gradient is implicitly satisfied. Open issues related to the accuracy of this approach are identified. An optimization-based approachmore » is implemented to enforce positivity of the determinant of the deformation gradient. The efficacy of this approach is illustrated with numerical examples.« less

A Lagrangian effective field theory

DOE Office of Scientific and Technical Information (OSTI.GOV)

Vlah, Zvonimir; White, Martin; Aviles, Alejandro

We have continued the development of Lagrangian, cosmological perturbation theory for the low-order correlators of the matter density field. We provide a new route to understanding how the effective field theory (EFT) of large-scale structure can be formulated in the Lagrandian framework and a new resummation scheme, comparing our results to earlier work and to a series of high-resolution N-body simulations in both Fourier and configuration space. The `new' terms arising from EFT serve to tame the dependence of perturbation theory on small-scale physics and improve agreement with simulations (though with an additional free parameter). We find that all ofmore » our models fare well on scales larger than about two to three times the non-linear scale, but fail as the non-linear scale is approached. This is slightly less reach than has been seen previously. At low redshift the Lagrangian model fares as well as EFT in its Eulerian formulation, but at higher z the Eulerian EFT fits the data to smaller scales than resummed, Lagrangian EFT. Furthermore, all the perturbative models fare better than linear theory.« less

A Lagrangian effective field theory

DOE Office of Scientific and Technical Information (OSTI.GOV)

Vlah, Zvonimir; White, Martin; Aviles, Alejandro, E-mail: zvlah@stanford.edu, E-mail: mwhite@berkeley.edu, E-mail: aviles@berkeley.edu

We have continued the development of Lagrangian, cosmological perturbation theory for the low-order correlators of the matter density field. We provide a new route to understanding how the effective field theory (EFT) of large-scale structure can be formulated in the Lagrandian framework and a new resummation scheme, comparing our results to earlier work and to a series of high-resolution N-body simulations in both Fourier and configuration space. The 'new' terms arising from EFT serve to tame the dependence of perturbation theory on small-scale physics and improve agreement with simulations (though with an additional free parameter). We find that all ofmore » our models fare well on scales larger than about two to three times the non-linear scale, but fail as the non-linear scale is approached. This is slightly less reach than has been seen previously. At low redshift the Lagrangian model fares as well as EFT in its Eulerian formulation, but at higher z the Eulerian EFT fits the data to smaller scales than resummed, Lagrangian EFT. All the perturbative models fare better than linear theory.« less

A Lagrangian effective field theory

DOE PAGES

Vlah, Zvonimir; White, Martin; Aviles, Alejandro

2015-09-02

We have continued the development of Lagrangian, cosmological perturbation theory for the low-order correlators of the matter density field. We provide a new route to understanding how the effective field theory (EFT) of large-scale structure can be formulated in the Lagrandian framework and a new resummation scheme, comparing our results to earlier work and to a series of high-resolution N-body simulations in both Fourier and configuration space. The `new' terms arising from EFT serve to tame the dependence of perturbation theory on small-scale physics and improve agreement with simulations (though with an additional free parameter). We find that all ofmore » our models fare well on scales larger than about two to three times the non-linear scale, but fail as the non-linear scale is approached. This is slightly less reach than has been seen previously. At low redshift the Lagrangian model fares as well as EFT in its Eulerian formulation, but at higher z the Eulerian EFT fits the data to smaller scales than resummed, Lagrangian EFT. Furthermore, all the perturbative models fare better than linear theory.« less

Developing a Learning Algorithm-Generated Empirical Relaxer

DOE Office of Scientific and Technical Information (OSTI.GOV)

Mitchell, Wayne; Kallman, Josh; Toreja, Allen

2016-03-30

One of the main difficulties when running Arbitrary Lagrangian-Eulerian (ALE) simulations is determining how much to relax the mesh during the Eulerian step. This determination is currently made by the user on a simulation-by-simulation basis. We present a Learning Algorithm-Generated Empirical Relaxer (LAGER) which uses a regressive random forest algorithm to automate this decision process. We also demonstrate that LAGER successfully relaxes a variety of test problems, maintains simulation accuracy, and has the potential to significantly decrease both the person-hours and computational hours needed to run a successful ALE simulation.

An Eulerian/Lagrangian coupling procedure for three-dimensional vortical flows

NASA Technical Reports Server (NTRS)

Felici, Helene M.; Drela, Mark

1993-01-01

A coupled Eulerian/Lagrangian method is presented for the reduction of numerical diffusion observed in solutions of 3D vortical flows using standard Eulerian finite-volume time-marching procedures. A Lagrangian particle tracking method, added to the Eulerian time-marching procedure, provides a correction of the Eulerian solution. In turn, the Eulerian solution is used to integrate the Lagrangian state-vector along the particles trajectories. While the Eulerian solution ensures the conservation of mass and sets the pressure field, the particle markers describe accurately the convection properties and enhance the vorticity and entropy capturing capabilities of the Eulerian solver. The Eulerian/Lagrangian coupling strategies are discussed and the combined scheme is tested on a constant stagnation pressure flow in a 90 deg bend and on a swirling pipe flow. As the numerical diffusion is reduced when using the Lagrangian correction, a vorticity gradient augmentation is identified as a basic problem of this inviscid calculation.

Compatible, energy conserving, bounds preserving remap of hydrodynamic fields for an extended ALE scheme

NASA Astrophysics Data System (ADS)

Burton, D. E.; Morgan, N. R.; Charest, M. R. J.; Kenamond, M. A.; Fung, J.

2018-02-01

From the very origins of numerical hydrodynamics in the Lagrangian work of von Neumann and Richtmyer [83], the issue of total energy conservation as well as entropy production has been problematic. Because of well known problems with mesh deformation, Lagrangian schemes have evolved into Arbitrary Lagrangian-Eulerian (ALE) methods [39] that combine the best properties of Lagrangian and Eulerian methods. Energy issues have persisted for this class of methods. We believe that fundamental issues of energy conservation and entropy production in ALE require further examination. The context of the paper is an ALE scheme that is extended in the sense that it permits cyclic or periodic remap of data between grids of the same or differing connectivity. The principal design goals for a remap method then consist of total energy conservation, bounded internal energy, and compatibility of kinetic energy and momentum. We also have secondary objectives of limiting velocity and stress in a non-directional manner, keeping primitive variables monotone, and providing a higher than second order reconstruction of remapped variables. In particular, the new contributions fall into three categories associated with: energy conservation and entropy production, reconstruction and bounds preservation of scalar and tensor fields, and conservative remap of nonlinear fields. The paper presents a derivation of the methods, details of implementation, and numerical results for a number of test problems. The methods requires volume integration of polynomial functions in polytopal cells with planar facets, and the requisite expressions are derived for arbitrary order.

Stability analysis of Eulerian-Lagrangian methods for the one-dimensional shallow-water equations

USGS Publications Warehouse

Casulli, V.; Cheng, R.T.

1990-01-01

In this paper stability and error analyses are discussed for some finite difference methods when applied to the one-dimensional shallow-water equations. Two finite difference formulations, which are based on a combined Eulerian-Lagrangian approach, are discussed. In the first part of this paper the results of numerical analyses for an explicit Eulerian-Lagrangian method (ELM) have shown that the method is unconditionally stable. This method, which is a generalized fixed grid method of characteristics, covers the Courant-Isaacson-Rees method as a special case. Some artificial viscosity is introduced by this scheme. However, because the method is unconditionally stable, the artificial viscosity can be brought under control either by reducing the spatial increment or by increasing the size of time step. The second part of the paper discusses a class of semi-implicit finite difference methods for the one-dimensional shallow-water equations. This method, when the Eulerian-Lagrangian approach is used for the convective terms, is also unconditionally stable and highly accurate for small space increments or large time steps. The semi-implicit methods seem to be more computationally efficient than the explicit ELM; at each time step a single tridiagonal system of linear equations is solved. The combined explicit and implicit ELM is best used in formulating a solution strategy for solving a network of interconnected channels. The explicit ELM is used at channel junctions for each time step. The semi-implicit method is then applied to the interior points in each channel segment. Following this solution strategy, the channel network problem can be reduced to a set of independent one-dimensional open-channel flow problems. Numerical results support properties given by the stability and error analyses. ?? 1990.

Application of an Arbitrary Lagrangian Eulerian Method to Describe High Velocity Gas-Particle Flow Behavior

DTIC Science & Technology

2011-09-01

applied in this work was based on some of the standard definitions of soil constitutive properties as found in, e.g., Chen and Baladi [20], and...Livermore, CA. [20] Chen, W. F., and Baladi , G. Y., 1985. Soil Plasticity: Theory and Implementation. Elsevier Science, New York. [21] Zimmerman, H.D

Blast Fragmentation Modeling and Analysis

DTIC Science & Technology

2010-10-31

weapons device containing a multiphase blast explosive (MBX). 1. INTRODUCTION The ARL Survivability Lethality and Analysis Directorate (SLAD) is...velocity. In order to simulate the highly complex phenomenon, the exploding cylinder is modeled with the hydrodynamics code ALE3D , an arbitrary...Lagrangian-Eulerian multiphysics code, developed at Lawrence Livermore National Laboratory. ALE3D includes physical properties, constitutive models for

A hybridizable discontinuous Galerkin method for modeling fluid–structure interaction

DOE PAGES

Sheldon, Jason P.; Miller, Scott T.; Pitt, Jonathan S.

2016-08-31

This study presents a novel application of the hybridizable discontinuous Galerkin (HDG) finite element method to the multi-physics simulation of coupled fluid–structure interaction (FSI) problems. Recent applications of the HDG method have primarily been for single-physics problems including both solids and fluids, which are necessary building blocks for FSI modeling. Utilizing these established models, HDG formulations for linear elastostatics, a nonlinear elastodynamic model, and arbitrary Lagrangian–Eulerian Navier–Stokes are derived. The elasticity formulations are written in a Lagrangian reference frame, with the nonlinear formulation restricted to hyperelastic materials. With these individual solid and fluid formulations, the remaining challenge in FSI modelingmore » is coupling together their disparate mathematics on the fluid–solid interface. This coupling is presented, along with the resultant HDG FSI formulation. Verification of the component models, through the method of manufactured solutions, is performed and each model is shown to converge at the expected rate. The individual components, along with the complete FSI model, are then compared to the benchmark problems proposed by Turek and Hron [1]. The solutions from the HDG formulation presented in this work trend towards the benchmark as the spatial polynomial order and the temporal order of integration are increased.« less

Reduction of numerical diffusion in three-dimensional vortical flows using a coupled Eulerian/Lagrangian solution procedure

NASA Technical Reports Server (NTRS)

Felici, Helene M.; Drela, Mark

1993-01-01

A new approach based on the coupling of an Eulerian and a Lagrangian solver, aimed at reducing the numerical diffusion errors of standard Eulerian time-marching finite-volume solvers, is presented. The approach is applied to the computation of the secondary flow in two bent pipes and the flow around a 3D wing. Using convective point markers the Lagrangian approach provides a correction of the basic Eulerian solution. The Eulerian flow in turn integrates in time the Lagrangian state-vector. A comparison of coarse and fine grid Eulerian solutions makes it possible to identify numerical diffusion. It is shown that the Eulerian/Lagrangian approach is an effective method for reducing numerical diffusion errors.

Imposing a Lagrangian Particle Framework on an Eulerian Hydrodynamics Infrastructure in Flash

NASA Technical Reports Server (NTRS)

Dubey, A.; Daley, C.; ZuHone, J.; Ricker, P. M.; Weide, K.; Graziani, C.

2012-01-01

In many astrophysical simulations, both Eulerian and Lagrangian quantities are of interest. For example, in a galaxy cluster merger simulation, the intracluster gas can have Eulerian discretization, while dark matter can be modeled using particles. FLASH, a component-based scientific simulation code, superimposes a Lagrangian framework atop an adaptive mesh refinement Eulerian framework to enable such simulations. The discretization of the field variables is Eulerian, while the Lagrangian entities occur in many different forms including tracer particles, massive particles, charged particles in particle-in-cell mode, and Lagrangian markers to model fluid structure interactions. These widely varying roles for Lagrangian entities are possible because of the highly modular, flexible, and extensible architecture of the Lagrangian framework. In this paper, we describe the Lagrangian framework in FLASH in the context of two very different applications, Type Ia supernovae and galaxy cluster mergers, which use the Lagrangian entities in fundamentally different ways.

Imposing a Lagrangian Particle Framework on an Eulerian Hydrodynamics Infrastructure in FLASH

NASA Astrophysics Data System (ADS)

Dubey, A.; Daley, C.; ZuHone, J.; Ricker, P. M.; Weide, K.; Graziani, C.

2012-08-01

In many astrophysical simulations, both Eulerian and Lagrangian quantities are of interest. For example, in a galaxy cluster merger simulation, the intracluster gas can have Eulerian discretization, while dark matter can be modeled using particles. FLASH, a component-based scientific simulation code, superimposes a Lagrangian framework atop an adaptive mesh refinement Eulerian framework to enable such simulations. The discretization of the field variables is Eulerian, while the Lagrangian entities occur in many different forms including tracer particles, massive particles, charged particles in particle-in-cell mode, and Lagrangian markers to model fluid-structure interactions. These widely varying roles for Lagrangian entities are possible because of the highly modular, flexible, and extensible architecture of the Lagrangian framework. In this paper, we describe the Lagrangian framework in FLASH in the context of two very different applications, Type Ia supernovae and galaxy cluster mergers, which use the Lagrangian entities in fundamentally different ways.

Dynamic analysis of a needle insertion for soft materials: Arbitrary Lagrangian-Eulerian-based three-dimensional finite element analysis.

PubMed

Yamaguchi, Satoshi; Tsutsui, Kihei; Satake, Koji; Morikawa, Shigehiro; Shirai, Yoshiaki; Tanaka, Hiromi T

2014-10-01

Our goal was to develop a three-dimensional finite element model that enables dynamic analysis of needle insertion for soft materials. To demonstrate large deformation and fracture, we used the arbitrary Lagrangian-Eulerian (ALE) method for fluid analysis. We performed ALE-based finite element analysis for 3% agar gel and three types of copper needle with bevel tips. To evaluate simulation results, we compared the needle deflection and insertion force with corresponding experimental results acquired with a uniaxial manipulator. We studied the shear stress distribution of agar gel on various time scales. For 30°, 45°, and 60°, differences in deflections of each needle between both sets of results were 2.424, 2.981, and 3.737mm, respectively. For the insertion force, there was no significant difference for mismatching area error (p<0.05) between simulation and experimental results. Our results have the potential to be a stepping stone to develop pre-operative surgical planning to estimate an optimal needle insertion path for MR image-guided microwave coagulation therapy and for analyzing large deformation and fracture in biological tissues. Copyright © 2014 Elsevier Ltd. All rights reserved.

A point-centered arbitrary Lagrangian Eulerian hydrodynamic approach for tetrahedral meshes

DOE PAGES

Morgan, Nathaniel R.; Waltz, Jacob I.; Burton, Donald E.; ...

2015-02-24

We present a three dimensional (3D) arbitrary Lagrangian Eulerian (ALE) hydrodynamic scheme suitable for modeling complex compressible flows on tetrahedral meshes. The new approach stores the conserved variables (mass, momentum, and total energy) at the nodes of the mesh and solves the conservation equations on a control volume surrounding the point. This type of an approach is termed a point-centered hydrodynamic (PCH) method. The conservation equations are discretized using an edge-based finite element (FE) approach with linear basis functions. All fluxes in the new approach are calculated at the center of each tetrahedron. A multidirectional Riemann-like problem is solved atmore » the center of the tetrahedron. The advective fluxes are calculated by solving a 1D Riemann problem on each face of the nodal control volume. A 2-stage Runge–Kutta method is used to evolve the solution forward in time, where the advective fluxes are part of the temporal integration. The mesh velocity is smoothed by solving a Laplacian equation. The details of the new ALE hydrodynamic scheme are discussed. Results from a range of numerical test problems are presented.« less

TESS: A RELATIVISTIC HYDRODYNAMICS CODE ON A MOVING VORONOI MESH

DOE Office of Scientific and Technical Information (OSTI.GOV)

Duffell, Paul C.; MacFadyen, Andrew I., E-mail: pcd233@nyu.edu, E-mail: macfadyen@nyu.edu

2011-12-01

We have generalized a method for the numerical solution of hyperbolic systems of equations using a dynamic Voronoi tessellation of the computational domain. The Voronoi tessellation is used to generate moving computational meshes for the solution of multidimensional systems of conservation laws in finite-volume form. The mesh-generating points are free to move with arbitrary velocity, with the choice of zero velocity resulting in an Eulerian formulation. Moving the points at the local fluid velocity makes the formulation effectively Lagrangian. We have written the TESS code to solve the equations of compressible hydrodynamics and magnetohydrodynamics for both relativistic and non-relativistic fluidsmore » on a dynamic Voronoi mesh. When run in Lagrangian mode, TESS is significantly less diffusive than fixed mesh codes and thus preserves contact discontinuities to high precision while also accurately capturing strong shock waves. TESS is written for Cartesian, spherical, and cylindrical coordinates and is modular so that auxiliary physics solvers are readily integrated into the TESS framework and so that this can be readily adapted to solve general systems of equations. We present results from a series of test problems to demonstrate the performance of TESS and to highlight some of the advantages of the dynamic tessellation method for solving challenging problems in astrophysical fluid dynamics.« less

A direct Arbitrary-Lagrangian-Eulerian ADER-WENO finite volume scheme on unstructured tetrahedral meshes for conservative and non-conservative hyperbolic systems in 3D

NASA Astrophysics Data System (ADS)

Boscheri, Walter; Dumbser, Michael

2014-10-01

In this paper we present a new family of high order accurate Arbitrary-Lagrangian-Eulerian (ALE) one-step ADER-WENO finite volume schemes for the solution of nonlinear systems of conservative and non-conservative hyperbolic partial differential equations with stiff source terms on moving tetrahedral meshes in three space dimensions. A WENO reconstruction technique is used to achieve high order of accuracy in space, while an element-local space-time Discontinuous Galerkin finite element predictor on moving curved meshes is used to obtain a high order accurate one-step time discretization. Within the space-time predictor the physical element is mapped onto a reference element using a high order isoparametric approach, where the space-time basis and test functions are given by the Lagrange interpolation polynomials passing through a predefined set of space-time nodes. Since our algorithm is cell-centered, the final mesh motion is computed by using a suitable node solver algorithm. A rezoning step as well as a flattener strategy are used in some of the test problems to avoid mesh tangling or excessive element deformations that may occur when the computation involves strong shocks or shear waves. The ALE algorithm presented in this article belongs to the so-called direct ALE methods because the final Lagrangian finite volume scheme is based directly on a space-time conservation formulation of the governing PDE system, with the rezoned geometry taken already into account during the computation of the fluxes. We apply our new high order unstructured ALE schemes to the 3D Euler equations of compressible gas dynamics, for which a set of classical numerical test problems has been solved and for which convergence rates up to sixth order of accuracy in space and time have been obtained. We furthermore consider the equations of classical ideal magnetohydrodynamics (MHD) as well as the non-conservative seven-equation Baer-Nunziato model of compressible multi-phase flows with stiff relaxation source terms.

An extended Lagrangian method for subsonic flows

NASA Technical Reports Server (NTRS)

Liou, Meng-Sing; Loh, Ching Y.

1992-01-01

It is well known that fluid motion can be specified by either the Eulerian of Lagrangian description. Most of Computational Fluid Dynamics (CFD) developments over the last three decades have been based on the Eulerian description and considerable progress has been made. In particular, the upwind methods, inspired and guided by the work of Gudonov, have met with many successes in dealing with complex flows, especially where discontinuities exist. However, this shock capturing property has proven to be accurate only when the discontinuity is aligned with one of the grid lines since most upwind methods are strictly formulated in 1-D framework and only formally extended to multi-dimensions. Consequently, the attractive property of crisp resolution of these discontinuities is lost and research on genuine multi-dimensional approach has just been undertaken by several leading researchers. Nevertheless they are still based on the Eulerian description.

A coupled Eulerian/Lagrangian method for the solution of three-dimensional vortical flows

NASA Technical Reports Server (NTRS)

Felici, Helene Marie

1992-01-01

A coupled Eulerian/Lagrangian method is presented for the reduction of numerical diffusion observed in solutions of three-dimensional rotational flows using standard Eulerian finite-volume time-marching procedures. A Lagrangian particle tracking method using particle markers is added to the Eulerian time-marching procedure and provides a correction of the Eulerian solution. In turn, the Eulerian solutions is used to integrate the Lagrangian state-vector along the particles trajectories. The Lagrangian correction technique does not require any a-priori information on the structure or position of the vortical regions. While the Eulerian solution ensures the conservation of mass and sets the pressure field, the particle markers, used as 'accuracy boosters,' take advantage of the accurate convection description of the Lagrangian solution and enhance the vorticity and entropy capturing capabilities of standard Eulerian finite-volume methods. The combined solution procedures is tested in several applications. The convection of a Lamb vortex in a straight channel is used as an unsteady compressible flow preservation test case. The other test cases concern steady incompressible flow calculations and include the preservation of turbulent inlet velocity profile, the swirling flow in a pipe, and the constant stagnation pressure flow and secondary flow calculations in bends. The last application deals with the external flow past a wing with emphasis on the trailing vortex solution. The improvement due to the addition of the Lagrangian correction technique is measured by comparison with analytical solutions when available or with Eulerian solutions on finer grids. The use of the combined Eulerian/Lagrangian scheme results in substantially lower grid resolution requirements than the standard Eulerian scheme for a given solution accuracy.

Mine Blast Loading: Experiments and Simulations

DTIC Science & Technology

2010-04-01

plates by approximately 50%. We investigated the root cause for this discrepancy. The simulations calculate a turbulent-like flow field characterized...Toussaint [19] evaluated two numerical methods, Smooth Particle Hydrodynamics ( SPH ) and Arbitrary Lagrangian Eulerian (ALE), to simulate a mine blast on...That is, the mine blast products were not flowing along the solid plate boundary in the simulations as freely as they should. 6 In particular, the V

Lagrangian and Eulerian statistics obtained from direct numerical simulations of homogeneous turbulence

NASA Technical Reports Server (NTRS)

Squires, Kyle D.; Eaton, John K.

1991-01-01

Direct numerical simulation is used to study dispersion in decaying isotropic turbulence and homogeneous shear flow. Both Lagrangian and Eulerian data are presented allowing direct comparison, but at fairly low Reynolds number. The quantities presented include properties of the dispersion tensor, isoprobability contours of particle displacement, Lagrangian and Eulerian velocity autocorrelations and time scale ratios, and the eddy diffusivity tensor. The Lagrangian time microscale is found to be consistently larger than the Eulerian microscale, presumably due to the advection of the small scales by the large scales in the Eulerian reference frame.

Compatible, energy conserving, bounds preserving remap of hydrodynamic fields for an extended ALE scheme

DOE PAGES

Burton, Donald E.; Morgan, Nathaniel Ray; Charest, Marc Robert Joseph; ...

2017-11-22

From the very origins of numerical hydrodynamics in the Lagrangian work of von Neumann and Richtmyer [83], the issue of total energy conservation as well as entropy production has been problematic. Because of well known problems with mesh deformation, Lagrangian schemes have evolved into Arbitrary Lagrangian–Eulerian (ALE) methods [39] that combine the best properties of Lagrangian and Eulerian methods. Energy issues have persisted for this class of methods. We believe that fundamental issues of energy conservation and entropy production in ALE require further examination. The context of the paper is an ALE scheme that is extended in the sense thatmore » it permits cyclic or periodic remap of data between grids of the same or differing connectivity. The principal design goals for a remap method then consist of total energy conservation, bounded internal energy, and compatibility of kinetic energy and momentum. We also have secondary objectives of limiting velocity and stress in a non-directional manner, keeping primitive variables monotone, and providing a higher than second order reconstruction of remapped variables. Particularly, the new contributions fall into three categories associated with: energy conservation and entropy production, reconstruction and bounds preservation of scalar and tensor fields, and conservative remap of nonlinear fields. Our paper presents a derivation of the methods, details of implementation, and numerical results for a number of test problems. The methods requires volume integration of polynomial functions in polytopal cells with planar facets, and the requisite expressions are derived for arbitrary order.« less

Compatible, energy conserving, bounds preserving remap of hydrodynamic fields for an extended ALE scheme

DOE Office of Scientific and Technical Information (OSTI.GOV)

Burton, Donald E.; Morgan, Nathaniel Ray; Charest, Marc Robert Joseph

From the very origins of numerical hydrodynamics in the Lagrangian work of von Neumann and Richtmyer [83], the issue of total energy conservation as well as entropy production has been problematic. Because of well known problems with mesh deformation, Lagrangian schemes have evolved into Arbitrary Lagrangian–Eulerian (ALE) methods [39] that combine the best properties of Lagrangian and Eulerian methods. Energy issues have persisted for this class of methods. We believe that fundamental issues of energy conservation and entropy production in ALE require further examination. The context of the paper is an ALE scheme that is extended in the sense thatmore » it permits cyclic or periodic remap of data between grids of the same or differing connectivity. The principal design goals for a remap method then consist of total energy conservation, bounded internal energy, and compatibility of kinetic energy and momentum. We also have secondary objectives of limiting velocity and stress in a non-directional manner, keeping primitive variables monotone, and providing a higher than second order reconstruction of remapped variables. Particularly, the new contributions fall into three categories associated with: energy conservation and entropy production, reconstruction and bounds preservation of scalar and tensor fields, and conservative remap of nonlinear fields. Our paper presents a derivation of the methods, details of implementation, and numerical results for a number of test problems. The methods requires volume integration of polynomial functions in polytopal cells with planar facets, and the requisite expressions are derived for arbitrary order.« less

A robust and efficient polyhedron subdivision and intersection algorithm for three-dimensional MMALE remapping

NASA Astrophysics Data System (ADS)

Chen, Xiang; Zhang, Xiong; Jia, Zupeng

2017-06-01

The Multi-Material Arbitrary Lagrangian Eulerian (MMALE) method is an effective way to simulate the multi-material flow with severe surface deformation. Comparing with the traditional Arbitrary Lagrangian Eulerian (ALE) method, the MMALE method allows for multiple materials in a single cell which overcomes the difficulties in grid refinement process. In recent decades, many researches have been conducted for the Lagrangian, rezoning and surface reconstruction phases, but less attention has been paid to the multi-material remapping phase especially for the three-dimensional problems due to two complex geometric problems: the polyhedron subdivision and the polyhedron intersection. In this paper, we propose a ;Clipping and Projecting; algorithm for polyhedron intersection whose basic idea comes from the commonly used method by Grandy (1999) [29] and Jia et al. (2013) [34]. Our new algorithm solves the geometric problem by an incremental modification of the topology based on segment-plane intersections. A comparison with Jia et al. (2013) [34] shows our new method improves the efficiency by 55% to 65% when calculating polyhedron intersections. Moreover, the instability caused by the geometric degeneracy can be thoroughly avoided because the geometry integrity is preserved in the new algorithm. We also focus on the polyhedron subdivision process and describe an algorithm which could automatically and precisely tackle the various situations including convex, non-convex and multiple subdivisions. Numerical studies indicate that by using our polyhedron subdivision and intersection algorithm, the volume conversation of the remapping phase can be exactly preserved in the MMALE simulation.

Variational Lagrangian data assimilation in open channel networks

NASA Astrophysics Data System (ADS)

Wu, Qingfang; Tinka, Andrew; Weekly, Kevin; Beard, Jonathan; Bayen, Alexandre M.

2015-04-01

This article presents a data assimilation method in a tidal system, where data from both Lagrangian drifters and Eulerian flow sensors were fused to estimate water velocity. The system is modeled by first-order, hyperbolic partial differential equations subject to periodic forcing. The estimation problem can then be formulated as the minimization of the difference between the observed variables and model outputs, and eventually provide the velocity and water stage of the hydrodynamic system. The governing equations are linearized and discretized using an implicit discretization scheme, resulting in linear equality constraints in the optimization program. Thus, the flow estimation can be formed as an optimization problem and efficiently solved. The effectiveness of the proposed method was substantiated by a large-scale field experiment in the Sacramento-San Joaquin River Delta in California. A fleet of 100 sensors developed at the University of California, Berkeley, were deployed in Walnut Grove, CA, to collect a set of Lagrangian data, a time series of positions as the sensors moved through the water. Measurements were also taken from Eulerian sensors in the region, provided by the United States Geological Survey. It is shown that the proposed method can effectively integrate Lagrangian and Eulerian measurement data, resulting in a suited estimation of the flow variables within the hydraulic system.

A Numerical Investigation of Two-Different Drosophila Forward Flight Modes

NASA Astrophysics Data System (ADS)

Sahin, Mehmet; Dilek, Ezgi; Erzincanli, Belkis

2016-11-01

The parallel large-scale unstructured finite volume method based on an Arbitrary Lagrangian-Eulerian (ALE) formulation has been applied in order to investigate the near wake structure of Drosophila in forward flight. DISTENE MeshGems-Hexa algorithm based on the octree method is used to generate the all hexahedral mesh for the wing-body combination. The mesh deformation algorithm is based on the indirect radial basis function (RBF) method at each time level while avoiding remeshing in order to enhance numerical robustness. The large-scale numerical simulations are carried out for a flapping Drosophila in forward flight. In the first case, the wing tip-path plane is tilted forward to generate forward force. In the second case, paddling wing motion is used to generate the forward fore. The λ2-criterion proposed by Jeong and Hussain (1995) is used for investigating the time variation of the Eulerian coherent structures in the near wake. The present simulations reveal highly detailed near wake topology for a hovering Drosophila. This is very useful in terms of understanding physics in biological flights which can provide a very useful tool for designing bio-inspired MAVs.

Structure of sheared and rotating turbulence: Multiscale statistics of Lagrangian and Eulerian accelerations and passive scalar dynamics.

PubMed

Jacobitz, Frank G; Schneider, Kai; Bos, Wouter J T; Farge, Marie

2016-01-01

The acceleration statistics of sheared and rotating homogeneous turbulence are studied using direct numerical simulation results. The statistical properties of Lagrangian and Eulerian accelerations are considered together with the influence of the rotation to shear ratio, as well as the scale dependence of their statistics. The probability density functions (pdfs) of both Lagrangian and Eulerian accelerations show a strong and similar dependence on the rotation to shear ratio. The variance and flatness of both accelerations are analyzed and the extreme values of the Eulerian acceleration are observed to be above those of the Lagrangian acceleration. For strong rotation it is observed that flatness yields values close to three, corresponding to Gaussian-like behavior, and for moderate and vanishing rotation the flatness increases. Furthermore, the Lagrangian and Eulerian accelerations are shown to be strongly correlated for strong rotation due to a reduced nonlinear term in this case. A wavelet-based scale-dependent analysis shows that the flatness of both Eulerian and Lagrangian accelerations increases as scale decreases, which provides evidence for intermittent behavior. For strong rotation the Eulerian acceleration is even more intermittent than the Lagrangian acceleration, while the opposite result is obtained for moderate rotation. Moreover, the dynamics of a passive scalar with gradient production in the direction of the mean velocity gradient is analyzed and the influence of the rotation to shear ratio is studied. Concerning the concentration of a passive scalar spread by the flow, the pdf of its Eulerian time rate of change presents higher extreme values than those of its Lagrangian time rate of change. This suggests that the Eulerian time rate of change of scalar concentration is mainly due to advection, while its Lagrangian counterpart is only due to gradient production and viscous dissipation.

An unstructured mesh arbitrary Lagrangian-Eulerian unsteady incompressible flow solver and its application to insect flight aerodynamics

NASA Astrophysics Data System (ADS)

Su, Xiaohui; Cao, Yuanwei; Zhao, Yong

2016-06-01

In this paper, an unstructured mesh Arbitrary Lagrangian-Eulerian (ALE) incompressible flow solver is developed to investigate the aerodynamics of insect hovering flight. The proposed finite-volume ALE Navier-Stokes solver is based on the artificial compressibility method (ACM) with a high-resolution method of characteristics-based scheme on unstructured grids. The present ALE model is validated and assessed through flow passing over an oscillating cylinder. Good agreements with experimental results and other numerical solutions are obtained, which demonstrates the accuracy and the capability of the present model. The lift generation mechanisms of 2D wing in hovering motion, including wake capture, delayed stall, rapid pitch, as well as clap and fling are then studied and illustrated using the current ALE model. Moreover, the optimized angular amplitude in symmetry model, 45°, is firstly reported in details using averaged lift and the energy power method. Besides, the lift generation of complete cyclic clap and fling motion, which is simulated by few researchers using the ALE method due to large deformation, is studied and clarified for the first time. The present ALE model is found to be a useful tool to investigate lift force generation mechanism for insect wing flight.

Bayesian Nonlinear Assimilation of Eulerian and Lagrangian Coastal Flow Data

DTIC Science & Technology

2015-09-30

Lagrangian Coastal Flow Data Dr. Pierre F.J. Lermusiaux Department of Mechanical Engineering Center for Ocean Science and Engineering Massachusetts...Develop and apply theory, schemes and computational systems for rigorous Bayesian nonlinear assimilation of Eulerian and Lagrangian coastal flow data...coastal ocean fields, both in Eulerian and Lagrangian forms. - Further develop and implement our GMM-DO schemes for robust Bayesian nonlinear estimation

Structured population dynamics: continuous size and discontinuous stage structures.

PubMed

Buffoni, Giuseppe; Pasquali, Sara

2007-04-01

A nonlinear stochastic model for the dynamics of a population with either a continuous size structure or a discontinuous stage structure is formulated in the Eulerian formalism. It takes into account dispersion effects due to stochastic variability of the development process of the individuals. The discrete equations of the numerical approximation are derived, and an analysis of the existence and stability of the equilibrium states is performed. An application to a copepod population is illustrated; numerical results of Eulerian and Lagrangian models are compared.

Chemical Continuous Time Random Walks

NASA Astrophysics Data System (ADS)

Aquino, T.; Dentz, M.

2017-12-01

Traditional methods for modeling solute transport through heterogeneous media employ Eulerian schemes to solve for solute concentration. More recently, Lagrangian methods have removed the need for spatial discretization through the use of Monte Carlo implementations of Langevin equations for solute particle motions. While there have been recent advances in modeling chemically reactive transport with recourse to Lagrangian methods, these remain less developed than their Eulerian counterparts, and many open problems such as efficient convergence and reconstruction of the concentration field remain. We explore a different avenue and consider the question: In heterogeneous chemically reactive systems, is it possible to describe the evolution of macroscopic reactant concentrations without explicitly resolving the spatial transport? Traditional Kinetic Monte Carlo methods, such as the Gillespie algorithm, model chemical reactions as random walks in particle number space, without the introduction of spatial coordinates. The inter-reaction times are exponentially distributed under the assumption that the system is well mixed. In real systems, transport limitations lead to incomplete mixing and decreased reaction efficiency. We introduce an arbitrary inter-reaction time distribution, which may account for the impact of incomplete mixing. This process defines an inhomogeneous continuous time random walk in particle number space, from which we derive a generalized chemical Master equation and formulate a generalized Gillespie algorithm. We then determine the modified chemical rate laws for different inter-reaction time distributions. We trace Michaelis-Menten-type kinetics back to finite-mean delay times, and predict time-nonlocal macroscopic reaction kinetics as a consequence of broadly distributed delays. Non-Markovian kinetics exhibit weak ergodicity breaking and show key features of reactions under local non-equilibrium.

A cavitation model based on Eulerian stochastic fields

NASA Astrophysics Data System (ADS)

Magagnato, F.; Dumond, J.

2013-12-01

Non-linear phenomena can often be described using probability density functions (pdf) and pdf transport models. Traditionally the simulation of pdf transport requires Monte-Carlo codes based on Lagrangian "particles" or prescribed pdf assumptions including binning techniques. Recently, in the field of combustion, a novel formulation called the stochastic-field method solving pdf transport based on Eulerian fields has been proposed which eliminates the necessity to mix Eulerian and Lagrangian techniques or prescribed pdf assumptions. In the present work, for the first time the stochastic-field method is applied to multi-phase flow and in particular to cavitating flow. To validate the proposed stochastic-field cavitation model, two applications are considered. Firstly, sheet cavitation is simulated in a Venturi-type nozzle. The second application is an innovative fluidic diode which exhibits coolant flashing. Agreement with experimental results is obtained for both applications with a fixed set of model constants. The stochastic-field cavitation model captures the wide range of pdf shapes present at different locations.

Eulerian and Lagrangian Plasma Jet Modeling for the Plasma Liner Experiment

NASA Astrophysics Data System (ADS)

Hatcher, Richard; Cassibry, Jason; Stanic, Milos; Loverich, John; Hakim, Ammar

2011-10-01

The Plasma Liner Experiment (PLX) aims to demonstrate the feasibility of using spherically-convergent plasma jets to from an imploding plasma liner. Our group has modified two hydrodynamic simulation codes to include radiative loss, tabular equations of state (EOS), and thermal transport. Nautilus, created by TechX Corporation, is a finite-difference Eulerian code which solves the MHD equations formulated as systems of hyperbolic conservation laws. The other is SPHC, a smoothed particle hydrodynamics code produced by Stellingwerf Consulting. Use of the Lagrangian fluid particle approach of SPH is motivated by the ability to accurately track jet interfaces, the plasma vacuum boundary, and mixing of various layers, but Eulerian codes have been in development for much longer and have better shock capturing. We validate these codes against experimental measurements of jet propagation, expansion, and merging of two jets. Precursor jets are observed to form at the jet interface. Conditions that govern evolution of two and more merging jets are explored.

An updated Lagrangian particle hydrodynamics (ULPH) for Newtonian fluids

NASA Astrophysics Data System (ADS)

Tu, Qingsong; Li, Shaofan

2017-11-01

In this work, we have developed an updated Lagrangian particle hydrodynamics (ULPH) for Newtonian fluid. Unlike the smoothed particle hydrodynamics, the non-local particle hydrodynamics formulation proposed here is consistent and convergence. Unlike the state-based peridynamics, the discrete particle dynamics proposed here has no internal material bond between particles, and it is not formulated with respect to initial or a fixed referential configuration. In specific, we have shown that (1) the non-local update Lagrangian particle hydrodynamics formulation converges to the conventional local fluid mechanics formulation; (2) the non-local updated Lagrangian particle hydrodynamics can capture arbitrary flow discontinuities without any changes in the formulation, and (3) the proposed non-local particle hydrodynamics is computationally efficient and robust.

A Fluid Structure Algorithm with Lagrange Multipliers to Model Free Swimming

NASA Astrophysics Data System (ADS)

Sahin, Mehmet; Dilek, Ezgi

2017-11-01

A new monolithic approach is prosed to solve the fluid-structure interaction (FSI) problem with Lagrange multipliers in order to model free swimming/flying. In the present approach, the fluid domain is modeled by the incompressible Navier-Stokes equations and discretized using an Arbitrary Lagrangian-Eulerian (ALE) formulation based on the stable side-centered unstructured finite volume method. The solid domain is modeled by the constitutive laws for the nonlinear Saint Venant-Kirchhoff material and the classical Galerkin finite element method is used to discretize the governing equations in a Lagrangian frame. In order to impose the body motion/deformation, the distance between the constraint pair nodes is imposed using the Lagrange multipliers, which is independent from the frame of reference. The resulting algebraic linear equations are solved in a fully coupled manner using a dual approach (null space method). The present numerical algorithm is initially validated for the classical FSI benchmark problems and then applied to the free swimming of three linked ellipses. The authors are grateful for the use of the computing resources provided by the National Center for High Performance Computing (UYBHM) under Grant Number 10752009 and the computing facilities at TUBITAK-ULAKBIM, High Performance and Grid Computing Center.

Effects of Helicity on Lagrangian and Eulerian Time Correlations in Turbulence

NASA Technical Reports Server (NTRS)

Rubinstein, Robert; Zhou, Ye

1998-01-01

Taylor series expansions of turbulent time correlation functions are applied to show that helicity influences Eulerian time correlations more strongly than Lagrangian time correlations: to second order in time, the helicity effect on Lagrangian time correlations vanishes, but the helicity effect on Eulerian time correlations is nonzero. Fourier analysis shows that the helicity effect on Eulerian time correlations is confined to the largest inertial range scales. Some implications for sound radiation by swirling flows are discussed.

On hydrostatic flows in isentropic coordinates

NASA Astrophysics Data System (ADS)

Bokhove, Onno

2000-01-01

The hydrostatic primitive equations of motion which have been used in large-scale weather prediction and climate modelling over the last few decades are analysed with variational methods in an isentropic Eulerian framework. The use of material isentropic coordinates for the Eulerian hydrostatic equations is known to have distinct conceptual advantages since fluid motion is, under inviscid and statically stable circumstances, confined to take place on quasi-horizontal isentropic surfaces. First, an Eulerian isentropic Hamilton's principle, expressed in terms of fluid parcel variables, is therefore derived by transformation of a Lagrangian Hamilton's principle to an Eulerian one. This Eulerian principle explicitly describes the boundary dynamics of the time-dependent domain in terms of advection of boundary isentropes sB; these are the values the isentropes have at their intersection with the (lower) boundary. A partial Legendre transform for only the interior variables yields an Eulerian ‘action’ principle. Secondly, Noether's theorem is used to derive energy and potential vorticity conservation from the Eulerian Hamilton's principle. Thirdly, these conservation laws are used to derive a wave-activity invariant which is second-order in terms of small-amplitude disturbances relative to a resting or moving basic state. Linear stability criteria are derived but only for resting basic states. In mid-latitudes a time- scale separation between gravity and vortical modes occurs. Finally, this time-scale separation suggests that conservative geostrophic and ageostrophic approximations can be made to the Eulerian action principle for hydrostatic flows. Approximations to Eulerian variational principles may be more advantageous than approximations to Lagrangian ones because non-dimensionalization and scaling tend to be based on Eulerian estimates of the characteristic scales involved. These approximations to the stratified hydrostatic formulation extend previous approximations to the shallow- water equations. An explicit variational derivation is given of an isentropic version of Hoskins & Bretherton's model for atmospheric fronts.

Using Underwater Explosion and Cylinder Expansion Tests to Calibrate Afterburn Models for Aluminized Explosives

NASA Astrophysics Data System (ADS)

Wedberg, Rasmus

2017-06-01

The study explores the combined use of underwater performance tests and cylinder expansion tests in order to parameterize detonation models for aluminized explosives which exhibit afterburning. The approach is suggested to be used in conjunction with thermochemical computation. A formulation containing RDX and aluminum powder is considered and several charges with varying masses are submerged and detonated. Pressure gauges are employed at horizontal distances scaling with the charge diameter, and the specific shock wave energy is shown to increase with charge mass. This is attributed to the combustion of aluminum particles after the Chapman-Jouguet plane. Cylinder expansion tests are carried out using Photon Doppler Velocimetry to register the wall expansion velocity. The tests are modeled using a multi-material arbitrary Lagrangian-Eulerian approach with the Guirguis-Miller model describing detonation with afterburning. The equation of state and afterburn rate law parameters are adjusted such that the model reproduces the results from the cylinder expansion and underwater tests. The approach seems promising, and might be valuable for aluminized explosive formulations intended to be used in a variety of confinement conditions. Swedish Armed Forces.

An extended Lagrangian method

NASA Technical Reports Server (NTRS)

Liou, Meng-Sing

1992-01-01

A unique formulation of describing fluid motion is presented. The method, referred to as 'extended Lagrangian method', is interesting from both theoretical and numerical points of view. The formulation offers accuracy in numerical solution by avoiding numerical diffusion resulting from mixing of fluxes in the Eulerian description. Meanwhile, it also avoids the inaccuracy incurred due to geometry and variable interpolations used by the previous Lagrangian methods. Unlike the Lagrangian method previously imposed which is valid only for supersonic flows, the present method is general and capable of treating subsonic flows as well as supersonic flows. The method proposed in this paper is robust and stable. It automatically adapts to flow features without resorting to clustering, thereby maintaining rather uniform grid spacing throughout and large time step. Moreover, the method is shown to resolve multi-dimensional discontinuities with a high level of accuracy, similar to that found in one-dimensional problems.

A Shock-Adaptive Godunov Scheme Based on the Generalised Lagrangian Formulation

NASA Astrophysics Data System (ADS)

Lepage, C. Y.; Hui, W. H.

1995-12-01

Application of the Godunov scheme to the Euler equations of gas dynamics based on the Eulerian formulation of flow smears discontinuities, sliplines especially, over several computational cells, while the accuracy in the smooth flow region is of the order O( h), where h is the cell width. Based on the generalised Lagrangian formulation (GLF) of Hui et al., the Godunov scheme yields superior accuracy. By the use of coordinate streamlines in the GLF, the slipline—itself a streamline—is resolved crisply. Infinite shock resolution is achieved through the splitting of shock-cells. An improved entropy-conservation formulation of the governing equations is also proposed for computations in smooth flow regions. Finally, the use of the GLF substantially simplifies the programming logic resulting in a very robust, accurate, and efficient scheme.

An extended Lagrangian method

NASA Technical Reports Server (NTRS)

Liou, Meng-Sing

1993-01-01

A unique formulation of describing fluid motion is presented. The method, referred to as 'extended Lagrangian method', is interesting from both theoretical and numerical points of view. The formulation offers accuracy in numerical solution by avoiding numerical diffusion resulting from mixing of fluxes in the Eulerian description. Meanwhile, it also avoids the inaccuracy incurred due to geometry and variable interpolations used by the previous Lagrangian methods. The present method is general and capable of treating subsonic flows as well as supersonic flows. The method proposed in this paper is robust and stable. It automatically adapts to flow features without resorting to clustering, thereby maintaining rather uniform grid spacing throughout and large time step. Moreover, the method is shown to resolve multidimensional discontinuities with a high level of accuracy, similar to that found in 1D problems.

General Lagrangian formulation for higher spin fields with arbitrary index symmetry. 2. Fermionic fields

NASA Astrophysics Data System (ADS)

Reshetnyak, A.

2013-04-01

We continue the construction of a Lagrangian description of irreducible half-integer higher-spin representations of the Poincare group with an arbitrary Young tableaux having k rows, on a basis of the BRST-BFV approach suggested for bosonic fields in our first article [I.L. Buchbinder, A. Reshetnyak, Nucl. Phys. B 862 (2012) 270, arXiv:1110.5044 [hep-th

A new Lagrangian random choice method for steady two-dimensional supersonic/hypersonic flow

NASA Technical Reports Server (NTRS)

Loh, C. Y.; Hui, W. H.

1991-01-01

Glimm's (1965) random choice method has been successfully applied to compute steady two-dimensional supersonic/hypersonic flow using a new Lagrangian formulation. The method is easy to program, fast to execute, yet it is very accurate and robust. It requires no grid generation, resolves slipline and shock discontinuities crisply, can handle boundary conditions most easily, and is applicable to hypersonic as well as supersonic flow. It represents an accurate and fast alternative to the existing Eulerian methods. Many computed examples are given.

A Eulerian-Lagrangian Model to Simulate Two-Phase/Particulate Flows

NASA Technical Reports Server (NTRS)

Apte, S. V.; Mahesh, K.; Lundgren, T.

2003-01-01

Figure 1 shows a snapshot of liquid fuel spray coming out of an injector nozzle in a realistic gas-turbine combustor. Here the spray atomization was simulated using a stochastic secondary breakup model (Apte et al. 2003a) with point-particle approximation for the droplets. Very close to the injector, it is observed that the spray density is large and the droplets cannot be treated as point-particles. The volume displaced by the liquid in this region is significant and can alter the gas-phase ow and spray evolution. In order to address this issue, one can compute the dense spray regime by an Eulerian-Lagrangian technique using advanced interface tracking/level-set methods (Sussman et al. 1994; Tryggvason et al. 2001; Herrmann 2003). This, however, is computationally intensive and may not be viable in realistic complex configurations. We therefore plan to develop a methodology based on Eulerian-Lagrangian technique which will allow us to capture the essential features of primary atomization using models to capture interactions between the fluid and droplets and which can be directly applied to the standard atomization models used in practice. The numerical scheme for unstructured grids developed by Mahesh et al. (2003) for incompressible flows is modified to take into account the droplet volume fraction. The numerical framework is directly applicable to realistic combustor geometries. Our main objectives in this work are: Develop a numerical formulation based on Eulerian-Lagrangian techniques with models for interaction terms between the fluid and particles to capture the Kelvin- Helmholtz type instabilities observed during primary atomization. Validate this technique for various two-phase and particulate flows. Assess its applicability to capture primary atomization of liquid jets in conjunction with secondary atomization models.

Sensitivity of Particle Size in Discrete Element Method to Particle Gas Method (DEM_PGM) Coupling in Underbody Blast Simulations

DTIC Science & Technology

2016-06-12

Particle Size in Discrete Element Method to Particle Gas Method (DEM_PGM) Coupling in Underbody Blast Simulations Venkatesh Babu, Kumar Kulkarni, Sanjay...buried in soil viz., (1) coupled discrete element & particle gas methods (DEM-PGM) and (2) Arbitrary Lagrangian-Eulerian (ALE), are investigated. The...DEM_PGM and identify the limitations/strengths compared to the ALE method. Discrete Element Method (DEM) can model individual particle directly, and

Numerical modeling of pulsed laser-material interaction and of laser plume dynamics

DOE Office of Scientific and Technical Information (OSTI.GOV)

Zhao, Qiang; Shi, Yina

2015-03-10

We have developed two-dimensional Arbitrary Lagrangian Eulerian (ALE) code which is used to study the physical processes, the plasma absorption, the crater profile, and the temperature distribution on metallic target and below the surface. The ALE method overcomes problems with Lagrangian moving mesh distortion by mesh smoothing and conservative quantities remapping from Lagrangian mesh to smoothed one. A new second order accurate diffusion solver has been implemented for the thermal conduction and radiation transport on distorted mesh. The results of numerical simulation of pulsed laser ablation are presented. The influences of different processes, such as time evolution of the surfacemore » temperature, interspecies interactions (elastic collisions, recombination-dissociation reaction), interaction with an ambient gas are examined. The study presents particular interest for the analysis of experimental results obtained during pulsed laser ablation.« less

Differential geometry based solvation model II: Lagrangian formulation.

PubMed

Chen, Zhan; Baker, Nathan A; Wei, G W

2011-12-01

Solvation is an elementary process in nature and is of paramount importance to more sophisticated chemical, biological and biomolecular processes. The understanding of solvation is an essential prerequisite for the quantitative description and analysis of biomolecular systems. This work presents a Lagrangian formulation of our differential geometry based solvation models. The Lagrangian representation of biomolecular surfaces has a few utilities/advantages. First, it provides an essential basis for biomolecular visualization, surface electrostatic potential map and visual perception of biomolecules. Additionally, it is consistent with the conventional setting of implicit solvent theories and thus, many existing theoretical algorithms and computational software packages can be directly employed. Finally, the Lagrangian representation does not need to resort to artificially enlarged van der Waals radii as often required by the Eulerian representation in solvation analysis. The main goal of the present work is to analyze the connection, similarity and difference between the Eulerian and Lagrangian formalisms of the solvation model. Such analysis is important to the understanding of the differential geometry based solvation model. The present model extends the scaled particle theory of nonpolar solvation model with a solvent-solute interaction potential. The nonpolar solvation model is completed with a Poisson-Boltzmann (PB) theory based polar solvation model. The differential geometry theory of surfaces is employed to provide a natural description of solvent-solute interfaces. The optimization of the total free energy functional, which encompasses the polar and nonpolar contributions, leads to coupled potential driven geometric flow and PB equations. Due to the development of singularities and nonsmooth manifolds in the Lagrangian representation, the resulting potential-driven geometric flow equation is embedded into the Eulerian representation for the purpose of computation, thanks to the equivalence of the Laplace-Beltrami operator in the two representations. The coupled partial differential equations (PDEs) are solved with an iterative procedure to reach a steady state, which delivers desired solvent-solute interface and electrostatic potential for problems of interest. These quantities are utilized to evaluate the solvation free energies and protein-protein binding affinities. A number of computational methods and algorithms are described for the interconversion of Lagrangian and Eulerian representations, and for the solution of the coupled PDE system. The proposed approaches have been extensively validated. We also verify that the mean curvature flow indeed gives rise to the minimal molecular surface and the proposed variational procedure indeed offers minimal total free energy. Solvation analysis and applications are considered for a set of 17 small compounds and a set of 23 proteins. The salt effect on protein-protein binding affinity is investigated with two protein complexes by using the present model. Numerical results are compared to the experimental measurements and to those obtained by using other theoretical methods in the literature. © Springer-Verlag 2011

Differential geometry based solvation model II: Lagrangian formulation

PubMed Central

Chen, Zhan; Baker, Nathan A.; Wei, G. W.

2010-01-01

Solvation is an elementary process in nature and is of paramount importance to more sophisticated chemical, biological and biomolecular processes. The understanding of solvation is an essential prerequisite for the quantitative description and analysis of biomolecular systems. This work presents a Lagrangian formulation of our differential geometry based solvation model. The Lagrangian representation of biomolecular surfaces has a few utilities/advantages. First, it provides an essential basis for biomolecular visualization, surface electrostatic potential map and visual perception of biomolecules. Additionally, it is consistent with the conventional setting of implicit solvent theories and thus, many existing theoretical algorithms and computational software packages can be directly employed. Finally, the Lagrangian representation does not need to resort to artificially enlarged van der Waals radii as often required by the Eulerian representation in solvation analysis. The main goal of the present work is to analyze the connection, similarity and difference between the Eulerian and Lagrangian formalisms of the solvation model. Such analysis is important to the understanding of the differential geometry based solvation model. The present model extends the scaled particle theory (SPT) of nonpolar solvation model with a solvent-solute interaction potential. The nonpolar solvation model is completed with a Poisson-Boltzmann (PB) theory based polar solvation model. The differential geometry theory of surfaces is employed to provide a natural description of solvent-solute interfaces. The minimization of the total free energy functional, which encompasses the polar and nonpolar contributions, leads to coupled potential driven geometric flow and Poisson-Boltzmann equations. Due to the development of singularities and nonsmooth manifolds in the Lagrangian representation, the resulting potential-driven geometric flow equation is embedded into the Eulerian representation for the purpose of computation, thanks to the equivalence of the Laplace-Beltrami operator in the two representations. The coupled partial differential equations (PDEs) are solved with an iterative procedure to reach a steady state, which delivers desired solvent-solute interface and electrostatic potential for problems of interest. These quantities are utilized to evaluate the solvation free energies and protein-protein binding affinities. A number of computational methods and algorithms are described for the interconversion of Lagrangian and Eulerian representations, and for the solution of the coupled PDE system. The proposed approaches have been extensively validated. We also verify that the mean curvature flow indeed gives rise to the minimal molecular surface (MMS) and the proposed variational procedure indeed offers minimal total free energy. Solvation analysis and applications are considered for a set of 17 small compounds and a set of 23 proteins. The salt effect on protein-protein binding affinity is investigated with two protein complexes by using the present model. Numerical results are compared to the experimental measurements and to those obtained by using other theoretical methods in the literature. PMID:21279359

Solution of mixed convection heat transfer from isothermal in-line fins

NASA Technical Reports Server (NTRS)

Khalilollahi, Amir

1993-01-01

Transient and steady state combined natural and forced convective flows over two in-line finite thickness fins (louvers) in a vertical channel are numerically solved using two methods. The first method of solution is based on the 'Simple Arbitrary Lagrangian Eulerian' (SALE) technique which incorporates mainly two computational phases: (1) a Lagrangian phase in which the velocity field is updated by the effects of all forces, and (2) an Eulerian phase that executes all advective fluxes of mass, momentum and energy. The second method of solution uses the finite element code entitled FIDAP. In the first part, comparison of the results by FIDAP, SALE, and available experimental work were done and discussed for steady state forced convection over louvered fins. Good agreements were deduced between the three sets of results especially for the flow over a single fin. In the second part and in the absence of experimental literature, the numerical predictions were extended to the transient transports and to the opposing flow where pressure drop is reversed. Results are presented and discussed for heat transfer and pressure drop in assisting and opposing mixed convection flows.

Action principle for Coulomb collisions in plasmas

DOE Office of Scientific and Technical Information (OSTI.GOV)

Hirvijoki, Eero

In this study, an action principle for Coulomb collisions in plasmas is proposed. Although no natural Lagrangian exists for the Landau-Fokker-Planck equation, an Eulerian variational formulation is found considering the system of partial differential equations that couple the distribution function and the Rosenbluth-MacDonald-Judd potentials. Conservation laws are derived after generalizing the energy-momentum stress tensor for second order Lagrangians and, in the case of a test-particle population in a given plasma background, the action principle is shown to correspond to the Langevin equation for individual particles.

Action principle for Coulomb collisions in plasmas

DOE PAGES

Hirvijoki, Eero

2016-09-14

In this study, an action principle for Coulomb collisions in plasmas is proposed. Although no natural Lagrangian exists for the Landau-Fokker-Planck equation, an Eulerian variational formulation is found considering the system of partial differential equations that couple the distribution function and the Rosenbluth-MacDonald-Judd potentials. Conservation laws are derived after generalizing the energy-momentum stress tensor for second order Lagrangians and, in the case of a test-particle population in a given plasma background, the action principle is shown to correspond to the Langevin equation for individual particles.

Numerical simulation of the fluid-structure interaction between air blast waves and soil structure

NASA Astrophysics Data System (ADS)

Umar, S.; Risby, M. S.; Albert, A. Luthfi; Norazman, M.; Ariffin, I.; Alias, Y. Muhamad

2014-03-01

Normally, an explosion threat on free field especially from high explosives is very dangerous due to the ground shocks generated that have high impulsive load. Nowadays, explosion threats do not only occur in the battlefield, but also in industries and urban areas. In industries such as oil and gas, explosion threats may occur on logistic transportation, maintenance, production, and distribution pipeline that are located underground to supply crude oil. Therefore, the appropriate blast resistances are a priority requirement that can be obtained through an assessment on the structural response, material strength and impact pattern of material due to ground shock. A highly impulsive load from ground shocks is a dynamic load due to its loading time which is faster than ground response time. Of late, almost all blast studies consider and analyze the ground shock in the fluid-structure interaction (FSI) because of its influence on the propagation and interaction of ground shock. Furthermore, analysis in the FSI integrates action of ground shock and reaction of ground on calculations of velocity, pressure and force. Therefore, this integration of the FSI has the capability to deliver the ground shock analysis on simulation to be closer to experimental investigation results. In this study, the FSI was implemented on AUTODYN computer code by using Euler-Godunov and the arbitrary Lagrangian-Eulerian (ALE). Euler-Godunov has the capability to deliver a structural computation on a 3D analysis, while ALE delivers an arbitrary calculation that is appropriate for a FSI analysis. In addition, ALE scheme delivers fine approach on little deformation analysis with an arbitrary motion, while the Euler-Godunov scheme delivers fine approach on a large deformation analysis. An integrated scheme based on Euler-Godunov and the arbitrary Lagrangian-Eulerian allows us to analyze the blast propagation waves and structural interaction simultaneously.

Duality for massive spin two theories in arbitrary dimensions

NASA Astrophysics Data System (ADS)

González, B.; Khoudeir, A.; Montemayor, R.; Urrutia, L. F.

2008-09-01

Using the parent Lagrangian approach we construct a dual formulation, in the sense originally proposed by Curtright and Freund, of a massive spin two Fierz-Pauli theory in arbitrary dimensions D. This is achieved in terms of a mixed symmetry tensor TA[B1B2...BD-2], without the need of auxiliary fields. The relation of this method with an alternative formulation based on a gauge symmetry principle proposed by Zinoviev is elucidated. We show that the latter formulation in four dimensions, with a given gauge fixing together with a definite sequence of auxiliary fields elimination via their equations of motion, leads to the parent Lagrangian already considered by West completed by a Fierz-Pauli mass term, which in turns yields the Curtright-Freund action. This motivates our generalization to arbitrary dimensions leading to the corresponding extension of the four dimensional result. We identify the transverse true degrees of freedom of the dual theory and verify that their number is in accordance with those of the massive Fierz-Pauli field.

Multi-Phase CFD Modeling of Solid Sorbent Carbon Capture System

DOE Office of Scientific and Technical Information (OSTI.GOV)

Ryan, Emily M.; DeCroix, David; Breault, Ronald W.

2013-07-30

Computational fluid dynamics (CFD) simulations are used to investigate a low temperature post-combustion carbon capture reactor. The CFD models are based on a small scale solid sorbent carbon capture reactor design from ADA-ES and Southern Company. The reactor is a fluidized bed design based on a silica-supported amine sorbent. CFD models using both Eulerian-Eulerian and Eulerian-Lagrangian multi-phase modeling methods are developed to investigate the hydrodynamics and adsorption of carbon dioxide in the reactor. Models developed in both FLUENT® and BARRACUDA are presented to explore the strengths and weaknesses of state of the art CFD codes for modeling multi-phase carbon capturemore » reactors. The results of the simulations show that the FLUENT® Eulerian-Lagrangian simulations (DDPM) are unstable for the given reactor design; while the BARRACUDA Eulerian-Lagrangian model is able to simulate the system given appropriate simplifying assumptions. FLUENT® Eulerian-Eulerian simulations also provide a stable solution for the carbon capture reactor given the appropriate simplifying assumptions.« less

Edgeworth streaming model for redshift space distortions

NASA Astrophysics Data System (ADS)

Uhlemann, Cora; Kopp, Michael; Haugg, Thomas

2015-09-01

We derive the Edgeworth streaming model (ESM) for the redshift space correlation function starting from an arbitrary distribution function for biased tracers of dark matter by considering its two-point statistics and show that it reduces to the Gaussian streaming model (GSM) when neglecting non-Gaussianities. We test the accuracy of the GSM and ESM independent of perturbation theory using the Horizon Run 2 N -body halo catalog. While the monopole of the redshift space halo correlation function is well described by the GSM, higher multipoles improve upon including the leading order non-Gaussian correction in the ESM: the GSM quadrupole breaks down on scales below 30 Mpc /h whereas the ESM stays accurate to 2% within statistical errors down to 10 Mpc /h . To predict the scale-dependent functions entering the streaming model we employ convolution Lagrangian perturbation theory (CLPT) based on the dust model and local Lagrangian bias. Since dark matter halos carry an intrinsic length scale given by their Lagrangian radius, we extend CLPT to the coarse-grained dust model and consider two different smoothing approaches operating in Eulerian and Lagrangian space, respectively. The coarse graining in Eulerian space features modified fluid dynamics different from dust while the coarse graining in Lagrangian space is performed in the initial conditions with subsequent single-streaming dust dynamics, implemented by smoothing the initial power spectrum in the spirit of the truncated Zel'dovich approximation. Finally, we compare the predictions of the different coarse-grained models for the streaming model ingredients to N -body measurements and comment on the proper choice of both the tracer distribution function and the smoothing scale. Since the perturbative methods we considered are not yet accurate enough on small scales, the GSM is sufficient when applied to perturbation theory.

AN EULERIAN-LAGRANGIAN LOCALIZED ADJOINT METHOD FOR THE ADVECTION-DIFFUSION EQUATION

EPA Science Inventory

Many numerical methods use characteristic analysis to accommodate the advective component of transport. Such characteristic methods include Eulerian-Lagrangian methods (ELM), modified method of characteristics (MMOC), and operator splitting methods. A generalization of characteri...

A LES-based Eulerian-Lagrangian approach to predict the dynamics of bubble plumes

NASA Astrophysics Data System (ADS)

Fraga, Bruño; Stoesser, Thorsten; Lai, Chris C. K.; Socolofsky, Scott A.

2016-01-01

An approach for Eulerian-Lagrangian large-eddy simulation of bubble plume dynamics is presented and its performance evaluated. The main numerical novelties consist in defining the gas-liquid coupling based on the bubble size to mesh resolution ratio (Dp/Δx) and the interpolation between Eulerian and Lagrangian frameworks through the use of delta functions. The model's performance is thoroughly validated for a bubble plume in a cubic tank in initially quiescent water using experimental data obtained from high-resolution ADV and PIV measurements. The predicted time-averaged velocities and second-order statistics show good agreement with the measurements, including the reproduction of the anisotropic nature of the plume's turbulence. Further, the predicted Eulerian and Lagrangian velocity fields, second-order turbulence statistics and interfacial gas-liquid forces are quantified and discussed as well as the visualization of the time-averaged primary and secondary flow structure in the tank.

Symmetry- and essentially-bound-preserving flux-corrected remapping of momentum in staggered ALE hydrodynamics

NASA Astrophysics Data System (ADS)

Velechovský, J.; Kuchařík, M.; Liska, R.; Shashkov, M.; Váchal, P.

2013-12-01

We present a new flux-corrected approach for remapping of velocity in the framework of staggered arbitrary Lagrangian-Eulerian methods. The main focus of the paper is the definition and preservation of coordinate invariant local bounds for velocity vector and development of momentum remapping method such that the radial symmetry of the radially symmetric flows is preserved when remapping from one equiangular polar mesh to another. The properties of this new method are demonstrated on a set of selected numerical cyclic remapping tests and a full hydrodynamic example.

Simulation of Particle Size Effect on Dynamic Properties and Fracture of PTFE-W-Al Composites

NASA Astrophysics Data System (ADS)

Herbold, E. B.; Cai, J.; Benson, D. J.; Nesterenko, V. F.

2007-12-01

Recent investigations of the dynamic compressive strength of cold isostatically pressed composites of polytetrafluoroethylene (PTFE), tungsten (W) and aluminum (Al) powders show significant differences depending on the size of metallic particles. The addition of W increases the density and changes the overall strength of the sample depending on the size of W particles. To investigate relatively large deformations, multi-material Eulerian and arbitrary Lagrangian-Eulerian methods, which have the ability to efficiently handle the formation of free surfaces, were used. The calculations indicate that the increased sample strength with fine metallic particles is due to the dynamic formation of force chains. This phenomenon occurs for samples with a higher porosity of the PTFE matrix compared to samples with larger particle size of W and a higher density PTFE matrix.

Unsteady Cascade Aerodynamic Response Using a Multiphysics Simulation Code

NASA Technical Reports Server (NTRS)

Lawrence, C.; Reddy, T. S. R.; Spyropoulos, E.

2000-01-01

The multiphysics code Spectrum(TM) is applied to calculate the unsteady aerodynamic pressures of oscillating cascade of airfoils representing a blade row of a turbomachinery component. Multiphysics simulation is based on a single computational framework for the modeling of multiple interacting physical phenomena, in the present case being between fluids and structures. Interaction constraints are enforced in a fully coupled manner using the augmented-Lagrangian method. The arbitrary Lagrangian-Eulerian method is utilized to account for deformable fluid domains resulting from blade motions. Unsteady pressures are calculated for a cascade designated as the tenth standard, and undergoing plunging and pitching oscillations. The predicted unsteady pressures are compared with those obtained from an unsteady Euler co-de refer-red in the literature. The Spectrum(TM) code predictions showed good correlation for the cases considered.

Topology of two-dimensional turbulent flows of dust and gas

NASA Astrophysics Data System (ADS)

Mitra, Dhrubaditya; Perlekar, Prasad

2018-04-01

We perform direct numerical simulations (DNS) of passive heavy inertial particles (dust) in homogeneous and isotropic two-dimensional turbulent flows (gas) for a range of Stokes number, St<1 . We solve for the particles using both a Lagrangian and an Eulerian approach (with a shock-capturing scheme). In the latter, the particles are described by a dust-density field and a dust-velocity field. We find the following: the dust-density field in our Eulerian simulations has the same correlation dimension d2 as obtained from the clustering of particles in the Lagrangian simulations for St<1 ; the cumulative probability distribution function of the dust density coarse grained over a scale r , in the inertial range, has a left tail with a power-law falloff indicating the presence of voids; the energy spectrum of the dust velocity has a power-law range with an exponent that is the same as the gas-velocity spectrum except at very high Fourier modes; the compressibility of the dust-velocity field is proportional to St2. We quantify the topological properties of the dust velocity and the gas velocity through their gradient matrices, called A and B , respectively. Our DNS confirms that the statistics of topological properties of B are the same in Eulerian and Lagrangian frames only if the Eulerian data are weighed by the dust density. We use this correspondence to study the statistics of topological properties of A in the Lagrangian frame from our Eulerian simulations by calculating density-weighted probability distribution functions. We further find that in the Lagrangian frame, the mean value of the trace of A is negative and its magnitude increases with St approximately as exp(-C /St) with a constant C ≈0.1 . The statistical distribution of different topological structures that appear in the dust flow is different in Eulerian and Lagrangian (density-weighted Eulerian) cases, particularly for St close to unity. In both of these cases, for small St the topological structures have close to zero divergence and are either vortical (elliptic) or strain dominated (hyperbolic, saddle). As St increases, the contribution to negative divergence comes mostly from saddles and the contribution to positive divergence comes from both vortices and saddles. Compared to the Eulerian case, the Lagrangian (density-weighted Eulerian) case has less outward spirals and more converging saddles. Inward spirals are the least probable topological structures in both cases.

Modeling and Numerical Challenges in Eulerian-Lagrangian Computations of Shock-driven Multiphase Flows

NASA Astrophysics Data System (ADS)

Diggs, Angela; Balachandar, Sivaramakrishnan

2015-06-01

The present work addresses the numerical methods required for particle-gas and particle-particle interactions in Eulerian-Lagrangian simulations of multiphase flow. Local volume fraction as seen by each particle is the quantity of foremost importance in modeling and evaluating such interactions. We consider a general multiphase flow with a distribution of particles inside a fluid flow discretized on an Eulerian grid. Particle volume fraction is needed both as a Lagrangian quantity associated with each particle and also as an Eulerian quantity associated with the flow. In Eulerian Projection (EP) methods, the volume fraction is first obtained within each cell as an Eulerian quantity and then interpolated to each particle. In Lagrangian Projection (LP) methods, the particle volume fraction is obtained at each particle and then projected onto the Eulerian grid. Traditionally, EP methods are used in multiphase flow, but sub-grid resolution can be obtained through use of LP methods. By evaluating the total error and its components we compare the performance of EP and LP methods. The standard von Neumann error analysis technique has been adapted for rigorous evaluation of rate of convergence. The methods presented can be extended to obtain accurate field representations of other Lagrangian quantities. Most importantly, we will show that such careful attention to numerical methodologies is needed in order to capture complex shock interaction with a bed of particles. Supported by U.S. Department of Defense SMART Program and the U.S. Department of Energy PSAAP-II program under Contract No. DE-NA0002378.

Hybrid Eulerian and Lagrangian Simulation of Steep and Breaking Waves and Surface Fluxes in High Winds

DTIC Science & Technology

2010-09-30

simulating violent free - surface flows , and show the importance of wave breaking in energy transport...using Eulerian simulation . 3 IMPACT/APPLICATION This project aims at developing an advanced simulation tool for multi-fluids free - surface flows that...several Eulerian and Lagrangian methods for free - surface turbulence and wave simulation . The WIND–SNOW is used to simulate 1 Report

Computing the stresses and deformations of the human eye components due to a high explosive detonation using fluid-structure interaction model.

PubMed

Karimi, Alireza; Razaghi, Reza; Navidbakhsh, Mahdi; Sera, Toshihiro; Kudo, Susumu

2016-05-01

In spite the fact that a very small human body surface area is comprised by the eye, its wounds due to detonation have recently been dramatically amplified. Although many efforts have been devoted to measure injury of the globe, there is still a lack of knowledge on the injury mechanism due to Primary Blast Wave (PBW). The goal of this study was to determine the stresses and deformations of the human eye components, including the cornea, aqueous, iris, ciliary body, lens, vitreous, retina, sclera, optic nerve, and muscles, attributed to PBW induced by trinitrotoluene (TNT) explosion via a Lagrangian-Eulerian computational coupling model. Magnetic Resonance Imaging (MRI) was employed to establish a Finite Element (FE) model of the human eye according to a normal human eye. The solid components of the eye were modelled as Lagrangian mesh, while an explosive TNT, air domain, and aqueous were modelled using Arbitrary Lagrangian-Eulerian (ALE) mesh. Nonlinear dynamic FE simulations were accomplished using the explicit FE code, namely LS-DYNA. In order to simulate the blast wave generation, propagation, and interaction with the eye, the ALE formulation with Jones-Wilkins-Lee (JWL) equation defining the explosive material were employed. The results revealed a peak stress of 135.70kPa brought about by detonation upsurge on the cornea at the distance of 25cm. The highest von Mises stresses were observed on the sclera (267.3kPa), whereas the lowest one was seen on the vitreous body (0.002kPa). The results also showed a relatively high resultant displacement for the macula as well as a high variation for the radius of curvature for the cornea and lens, which can result in both macular holes, optic nerve damage and, consequently, vision loss. These results may have implications not only for understanding the value of stresses and strains in the human eye components but also giving an outlook about the process of PBW triggers damage to the eye. Copyright © 2016 Elsevier Ltd. All rights reserved.

Forms of null Lagrangians in field theories of continuum mechanics

NASA Astrophysics Data System (ADS)

Kovalev, V. A.; Radaev, Yu. N.

2012-02-01

The divergence representation of a null Lagrangian that is regular in a star-shaped domain is used to obtain its general expression containing field gradients of order ≤ 1 in the case of spacetime of arbitrary dimension. It is shown that for a static three-component field in the three-dimensional space, a null Lagrangian can contain up to 15 independent elements in total. The general form of a null Lagrangian in the four-dimensional Minkowski spacetime is obtained (the number of physical field variables is assumed arbitrary). A complete theory of the null Lagrangian for the n-dimensional spacetime manifold (including the four-dimensional Minkowski spacetime as a special case) is given. Null Lagrangians are then used as a basis for solving an important variational problem of an integrating factor. This problem involves searching for factors that depend on the spacetime variables, field variables, and their gradients and, for a given system of partial differential equations, ensure the equality between the scalar product of a vector multiplier by the system vector and some divergence expression for arbitrary field variables and, hence, allow one to formulate a divergence conservation law on solutions to the system.

Stochastic partial differential fluid equations as a diffusive limit of deterministic Lagrangian multi-time dynamics.

PubMed

Cotter, C J; Gottwald, G A; Holm, D D

2017-09-01

In Holm (Holm 2015 Proc. R. Soc. A 471 , 20140963. (doi:10.1098/rspa.2014.0963)), stochastic fluid equations were derived by employing a variational principle with an assumed stochastic Lagrangian particle dynamics. Here we show that the same stochastic Lagrangian dynamics naturally arises in a multi-scale decomposition of the deterministic Lagrangian flow map into a slow large-scale mean and a rapidly fluctuating small-scale map. We employ homogenization theory to derive effective slow stochastic particle dynamics for the resolved mean part, thereby obtaining stochastic fluid partial equations in the Eulerian formulation. To justify the application of rigorous homogenization theory, we assume mildly chaotic fast small-scale dynamics, as well as a centring condition. The latter requires that the mean of the fluctuating deviations is small, when pulled back to the mean flow.

A Lagrangian Transport Eulerian Reaction Spatial (LATERS) Markov Model for Prediction of Effective Bimolecular Reactive Transport

NASA Astrophysics Data System (ADS)

Sund, Nicole; Porta, Giovanni; Bolster, Diogo; Parashar, Rishi

2017-11-01

Prediction of effective transport for mixing-driven reactive systems at larger scales, requires accurate representation of mixing at small scales, which poses a significant upscaling challenge. Depending on the problem at hand, there can be benefits to using a Lagrangian framework, while in others an Eulerian might have advantages. Here we propose and test a novel hybrid model which attempts to leverage benefits of each. Specifically, our framework provides a Lagrangian closure required for a volume-averaging procedure of the advection diffusion reaction equation. This hybrid model is a LAgrangian Transport Eulerian Reaction Spatial Markov model (LATERS Markov model), which extends previous implementations of the Lagrangian Spatial Markov model and maps concentrations to an Eulerian grid to quantify closure terms required to calculate the volume-averaged reaction terms. The advantage of this approach is that the Spatial Markov model is known to provide accurate predictions of transport, particularly at preasymptotic early times, when assumptions required by traditional volume-averaging closures are least likely to hold; likewise, the Eulerian reaction method is efficient, because it does not require calculation of distances between particles. This manuscript introduces the LATERS Markov model and demonstrates by example its ability to accurately predict bimolecular reactive transport in a simple benchmark 2-D porous medium.

Dual domain material point method for multiphase flows

NASA Astrophysics Data System (ADS)

Zhang, Duan

2017-11-01

Although the particle-in-cell method was first invented in the 60's for fluid computations, one of its later versions, the material point method, is mostly used for solid calculations. Recent development of the multi-velocity formulations for multiphase flows and fluid-structure interactions requires the Lagrangian capability of the method be combined with Eulerian calculations for fluids. Because of different numerical representations of the materials, additional numerical schemes are needed to ensure continuity of the materials. New applications of the method to compute fluid motions have revealed numerical difficulties in various versions of the method. To resolve these difficulties, the dual domain material point method is introduced and improved. Unlike other particle based methods, the material point method uses both Lagrangian particles and Eulerian mesh, therefore it avoids direct communication between particles. With this unique property and the Lagrangian capability of the method, it is shown that a multiscale numerical scheme can be efficiently built based on the dual domain material point method. In this talk, the theoretical foundation of the method will be introduced. Numerical examples will be shown. Work sponsored by the next generation code project of LANL.

Stochastic-field cavitation model

NASA Astrophysics Data System (ADS)

Dumond, J.; Magagnato, F.; Class, A.

2013-07-01

Nonlinear phenomena can often be well described using probability density functions (pdf) and pdf transport models. Traditionally, the simulation of pdf transport requires Monte-Carlo codes based on Lagrangian "particles" or prescribed pdf assumptions including binning techniques. Recently, in the field of combustion, a novel formulation called the stochastic-field method solving pdf transport based on Eulerian fields has been proposed which eliminates the necessity to mix Eulerian and Lagrangian techniques or prescribed pdf assumptions. In the present work, for the first time the stochastic-field method is applied to multi-phase flow and, in particular, to cavitating flow. To validate the proposed stochastic-field cavitation model, two applications are considered. First, sheet cavitation is simulated in a Venturi-type nozzle. The second application is an innovative fluidic diode which exhibits coolant flashing. Agreement with experimental results is obtained for both applications with a fixed set of model constants. The stochastic-field cavitation model captures the wide range of pdf shapes present at different locations.

Stochastic-field cavitation model

DOE Office of Scientific and Technical Information (OSTI.GOV)

Dumond, J., E-mail: julien.dumond@areva.com; AREVA GmbH, Erlangen, Paul-Gossen-Strasse 100, D-91052 Erlangen; Magagnato, F.

2013-07-15

Nonlinear phenomena can often be well described using probability density functions (pdf) and pdf transport models. Traditionally, the simulation of pdf transport requires Monte-Carlo codes based on Lagrangian “particles” or prescribed pdf assumptions including binning techniques. Recently, in the field of combustion, a novel formulation called the stochastic-field method solving pdf transport based on Eulerian fields has been proposed which eliminates the necessity to mix Eulerian and Lagrangian techniques or prescribed pdf assumptions. In the present work, for the first time the stochastic-field method is applied to multi-phase flow and, in particular, to cavitating flow. To validate the proposed stochastic-fieldmore » cavitation model, two applications are considered. First, sheet cavitation is simulated in a Venturi-type nozzle. The second application is an innovative fluidic diode which exhibits coolant flashing. Agreement with experimental results is obtained for both applications with a fixed set of model constants. The stochastic-field cavitation model captures the wide range of pdf shapes present at different locations.« less

A second-order shock-adaptive Godunov scheme based on the generalized Lagrangian formulation

NASA Astrophysics Data System (ADS)

Lepage, Claude

Application of the Godunov scheme to the Euler equations of gas dynamics, based on the Eulerian formulation of flow, smears discontinuities (especially sliplines) over several computational cells, while the accuracy in the smooth flow regions is of the order of a function of the cell width. Based on the generalized Lagrangian formulation (GLF), the Godunov scheme yields far superior results. By the use of coordinate streamlines in the GLF, the slipline (itself a streamline) is resolved crisply. Infinite shock resolution is achieved through the splitting of shock cells, while the accuracy in the smooth flow regions is improved using a nonconservative formulation of the governing equations coupled to a second order extension of the Godunov scheme. Furthermore, GLF requires no grid generation for boundary value problems and the simple structure of the solution to the Riemann problem in the GLF is exploited in the numerical implementation of the shock adaptive scheme. Numerical experiments reveal high efficiency and unprecedented resolution of shock and slipline discontinuities.

Examination of Eulerian and Lagrangian Coordinate Systems.

ERIC Educational Resources Information Center

Remillard, Wilfred J.

1978-01-01

Studies the relationship between Eulerian and Lagrangian coordinate systems with the help of computer plots of variables such as density and particle displacement. Gives examples which illustrate the differences in the shape of a traveling wave as seen by observers in the two systems. (Author/GA)

The Trapping Index: How to integrate the Eulerian and the Lagrangian approach for the computation of the transport time scales of semi-enclosed basins.

PubMed

Cucco, Andrea; Umgiesser, Georg

2015-09-15

In this work, we investigated if the Eulerian and the Lagrangian approaches for the computation of the Transport Time Scales (TTS) of semi-enclosed water bodies can be used univocally to define the spatial variability of basin flushing features. The Eulerian and Lagrangian TTS were computed for both simplified test cases and a realistic domain: the Venice Lagoon. The results confirmed the two approaches cannot be adopted univocally and that the spatial variability of the water renewal capacity can be investigated only through the computation of both the TTS. A specific analysis, based on the computation of a so-called Trapping Index, was then suggested to integrate the information provided by the two different approaches. The obtained results proved the Trapping Index to be useful to avoid any misleading interpretation due to the evaluation of the basin renewal features just from an Eulerian only or from a Lagrangian only perspective. Copyright © 2015 Elsevier Ltd. All rights reserved.

Parallel implementation of a Lagrangian-based model on an adaptive mesh in C++: Application to sea-ice

NASA Astrophysics Data System (ADS)

Samaké, Abdoulaye; Rampal, Pierre; Bouillon, Sylvain; Ólason, Einar

2017-12-01

We present a parallel implementation framework for a new dynamic/thermodynamic sea-ice model, called neXtSIM, based on the Elasto-Brittle rheology and using an adaptive mesh. The spatial discretisation of the model is done using the finite-element method. The temporal discretisation is semi-implicit and the advection is achieved using either a pure Lagrangian scheme or an Arbitrary Lagrangian Eulerian scheme (ALE). The parallel implementation presented here focuses on the distributed-memory approach using the message-passing library MPI. The efficiency and the scalability of the parallel algorithms are illustrated by the numerical experiments performed using up to 500 processor cores of a cluster computing system. The performance obtained by the proposed parallel implementation of the neXtSIM code is shown being sufficient to perform simulations for state-of-the-art sea ice forecasting and geophysical process studies over geographical domain of several millions squared kilometers like the Arctic region.

Coupled Eulerian-Lagrangian transport of large debris by tsunamis

NASA Astrophysics Data System (ADS)

Conde, Daniel A. S.; Ferreira, Rui M. L.; Sousa Oliveira, Carlos

2016-04-01

Tsunamis are notorious for the large disruption they can cause on coastal environments, not only due to the imparted momentum of the incoming wave but also due to its capacity to transport large quantities of solid debris, either from natural or human-made sources, over great distances. A 2DH numerical model under development at CERIS-IST (Ferreira et al., 2009; Conde, 2013) - STAV2D - capable of simulating solid transport in both Eulerian and Lagrangian paradigms will be used to assess the relevance of Lagrangian-Eulerian coupling when modelling the transport of solid debris by tsunamis. The model has been previously validated and applied to tsunami scenarios (Conde, 2013), being well-suited for overland tsunami propagation and capable of handling morphodynamic changes in estuaries and seashores. The discretization scheme is an explicit Finite Volume technique employing flux-vector splitting and a reviewed Roe-Riemann solver. Source term formulations are employed in a semi-implicit way, including the two-way coupling of the Lagrangian and Eulerian solvers by means of conservative mass and momentum transfers between fluid and solid phases. The model was applied to Sines Port, a major commercial port in Portugal, where two tsunamigenic scenarios are considered: an 8.5 Mw scenario, consistent with the Great Lisbon Earthquake and Tsunami of the 1st November 1755 (Baptista, 2009), and an hypothetical 9.5 Mw worst-case scenario based on the same historical event. Open-ocean propagation of these scenarios were simulated with GeoClaw model from ClawPack (Leveque, 2011). Following previous efforts on the modelling of debris transport by tsunamis in seaports (Conde, 2015), this work discusses the sensitivity of the obtained results with respect to the phenomenological detail of the employed Eulerian-Lagrangian formulation and the resolution of the mesh used in the Eulerian solver. The results have shown that the fluid to debris mass ratio is the key parameter regarding the conservativeness of the model. This way, in highly resolved meshes and high quantities of debris, the model approaches full conservativeness only if the two-way coupling feature is present, an effect that is attenuated in coarse meshes or with small debris quantities. Aknownledgements: This work was partially funded by FEDER, program COMPETE, and by national funds through the Portuguese Foundation for Science and Technology (FCT) with project RECI/ECM-HID/0371/2012. References: Baptista M.A. & Miranda, J.M. (2009) Revision of the Portuguese catalog of tsunamis. Nat. Hazards Earth Syst. Sci., 9, 25-42. Conde, D. A. S.; Baptista, M. A. V.; Sousa Oliveira, C. & Ferreira, R. M. L. (2013) A shallow-flow model for the propagation of tsunamis over complex geometries and mobile beds, Nat. Hazards Earth Syst. Sci., 13, 2533-2542. Conde, D. A. S.; Baptista, M. A. V.; Sousa Oliveira, C. & Ferreira, R. M. L. (2015) Mathematical modelling of tsunami impacts on critical infrastructures: exposure and severity associated with debris transport at Sines port. EGU General Assembly 2015, Vienna, Austria. Ferreira, R. M. L.; Franca, M. J.; Leal, J. G. & Cardoso, A. H. (2009) Mathematical modelling of shallow flows: Closure models drawn from grain-scale mechanics of sediment transport and flow hydrodynamics, Can. J. Civil. Eng., 36, 1604-1621. LeVeque, R. J., George, D. L., & Berger, M. J. (2011) Tsunami modelling with adaptively refined finite volume methods, Acta Numerica, pp. 211-289.

Multiphase Fluid Dynamics for Spacecraft Applications

NASA Astrophysics Data System (ADS)

Shyy, W.; Sim, J.

2011-09-01

Multiphase flows involving moving interfaces between different fluids/phases are observed in nature as well as in a wide range of engineering applications. With the recent development of high fidelity computational techniques, a number of challenging multiphase flow problems can now be computed. We introduce the basic notion of the main categories of multiphase flow computation; Lagrangian, Eulerian, and Eulerian-Lagrangian techniques to represent and follow interface, and sharp and continuous interface methods to model interfacial dynamics. The marker-based adaptive Eulerian-Lagrangian method, which is one of the most popular methods, is highlighted with microgravity and space applications including droplet collision and spacecraft liquid fuel tank surface stability.

A coupled ALE-AMR method for shock hydrodynamics

DOE PAGES

Waltz, J.; Bakosi, J.

2018-03-05

We present a numerical method combining adaptive mesh refinement (AMR) with arbitrary Lagrangian-Eulerian (ALE) mesh motion for the simulation of shock hydrodynamics on unstructured grids. The primary goal of the coupled method is to use AMR to reduce numerical error in ALE simulations at reduced computational expense relative to uniform fine mesh calculations, in the same manner that AMR has been used in Eulerian simulations. We also identify deficiencies with ALE methods that AMR is able to mitigate, and discuss the unique coupling challenges. The coupled method is demonstrated using three-dimensional unstructured meshes of up to O(10 7) tetrahedral cells.more » Convergence of ALE-AMR solutions towards both uniform fine mesh ALE results and analytic solutions is demonstrated. Speed-ups of 5-10× for a given level of error are observed relative to uniform fine mesh calculations.« less

A coupled ALE-AMR method for shock hydrodynamics

DOE Office of Scientific and Technical Information (OSTI.GOV)

Waltz, J.; Bakosi, J.

We present a numerical method combining adaptive mesh refinement (AMR) with arbitrary Lagrangian-Eulerian (ALE) mesh motion for the simulation of shock hydrodynamics on unstructured grids. The primary goal of the coupled method is to use AMR to reduce numerical error in ALE simulations at reduced computational expense relative to uniform fine mesh calculations, in the same manner that AMR has been used in Eulerian simulations. We also identify deficiencies with ALE methods that AMR is able to mitigate, and discuss the unique coupling challenges. The coupled method is demonstrated using three-dimensional unstructured meshes of up to O(10 7) tetrahedral cells.more » Convergence of ALE-AMR solutions towards both uniform fine mesh ALE results and analytic solutions is demonstrated. Speed-ups of 5-10× for a given level of error are observed relative to uniform fine mesh calculations.« less

Simulation of Particle Size Effect on Dynamic Properties and Fracture of PTFE-W-Al Composites

NASA Astrophysics Data System (ADS)

Herbold, Eric; Cai, Jing; Benson, David; Nesterenko, Vitali

2007-06-01

Recent investigations of the dynamic compressive strength of cold isostatically pressed (CIP) composites of polytetrafluoroethylene (PTFE), tungsten and aluminum powders show significant differences depending on the size of metallic particles. PTFE and aluminum mixtures are known to be energetic under dynamic and thermal loading. The addition of tungsten increases density and overall strength of the sample. Multi-material Eulerian and arbitrary Lagrangian-Eulerian methods were used for the investigation due to the complexity of the microstructure, relatively large deformations and the ability to handle the formation of free surfaces in a natural manner. The calculations indicate that the observed dependence of sample strength on particle size is due to the formation of force chains under dynamic loading in samples with small particle sizes even at larger porosity in comparison with samples with large grain size and larger density.

Stochastic partial differential fluid equations as a diffusive limit of deterministic Lagrangian multi-time dynamics

PubMed Central

Cotter, C. J.

2017-01-01

In Holm (Holm 2015 Proc. R. Soc. A 471, 20140963. (doi:10.1098/rspa.2014.0963)), stochastic fluid equations were derived by employing a variational principle with an assumed stochastic Lagrangian particle dynamics. Here we show that the same stochastic Lagrangian dynamics naturally arises in a multi-scale decomposition of the deterministic Lagrangian flow map into a slow large-scale mean and a rapidly fluctuating small-scale map. We employ homogenization theory to derive effective slow stochastic particle dynamics for the resolved mean part, thereby obtaining stochastic fluid partial equations in the Eulerian formulation. To justify the application of rigorous homogenization theory, we assume mildly chaotic fast small-scale dynamics, as well as a centring condition. The latter requires that the mean of the fluctuating deviations is small, when pulled back to the mean flow. PMID:28989316

Partitioned fluid-solid coupling for cardiovascular blood flow: left-ventricular fluid mechanics.

PubMed

Krittian, Sebastian; Janoske, Uwe; Oertel, Herbert; Böhlke, Thomas

2010-04-01

We present a 3D code-coupling approach which has been specialized towards cardiovascular blood flow. For the first time, the prescribed geometry movement of the cardiovascular flow model KaHMo (Karlsruhe Heart Model) has been replaced by a myocardial composite model. Deformation is driven by fluid forces and myocardial response, i.e., both its contractile and constitutive behavior. Whereas the arbitrary Lagrangian-Eulerian formulation (ALE) of the Navier-Stokes equations is discretized by finite volumes (FVM), the solid mechanical finite elasticity equations are discretized by a finite element (FEM) approach. Taking advantage of specialized numerical solution strategies for non-matching fluid and solid domain meshes, an iterative data-exchange guarantees the interface equilibrium of the underlying governing equations. The focus of this work is on left-ventricular fluid-structure interaction based on patient-specific magnetic resonance imaging datasets. Multi-physical phenomena are described by temporal visualization and characteristic FSI numbers. The results gained show flow patterns that are in good agreement with previous observations. A deeper understanding of cavity deformation, blood flow, and their vital interaction can help to improve surgical treatment and clinical therapy planning.

A Finite Element Method for Simulation of Compressible Cavitating Flows

NASA Astrophysics Data System (ADS)

Shams, Ehsan; Yang, Fan; Zhang, Yu; Sahni, Onkar; Shephard, Mark; Oberai, Assad

2016-11-01

This work focuses on a novel approach for finite element simulations of multi-phase flows which involve evolving interface with phase change. Modeling problems, such as cavitation, requires addressing multiple challenges, including compressibility of the vapor phase, interface physics caused by mass, momentum and energy fluxes. We have developed a mathematically consistent and robust computational approach to address these problems. We use stabilized finite element methods on unstructured meshes to solve for the compressible Navier-Stokes equations. Arbitrary Lagrangian-Eulerian formulation is used to handle the interface motions. Our method uses a mesh adaptation strategy to preserve the quality of the volumetric mesh, while the interface mesh moves along with the interface. The interface jump conditions are accurately represented using a discontinuous Galerkin method on the conservation laws. Condensation and evaporation rates at the interface are thermodynamically modeled to determine the interface velocity. We will present initial results on bubble cavitation the behavior of an attached cavitation zone in a separated boundary layer. We acknowledge the support from Army Research Office (ARO) under ARO Grant W911NF-14-1-0301.

Computational comparison of aortic root stresses in presence of stentless and stented aortic valve bio-prostheses.

PubMed

Nestola, M G C; Faggiano, E; Vergara, C; Lancellotti, R M; Ippolito, S; Antona, C; Filippi, S; Quarteroni, A; Scrofani, R

2017-02-01

We provide a computational comparison of the performance of stentless and stented aortic prostheses, in terms of aortic root displacements and internal stresses. To this aim, we consider three real patients; for each of them, we draw the two prostheses configurations, which are characterized by different mechanical properties and we also consider the native configuration. For each of these scenarios, we solve the fluid-structure interaction problem arising between blood and aortic root, through Finite Elements. In particular, the Arbitrary Lagrangian-Eulerian formulation is used for the numerical solution of the fluid-dynamic equations and a hyperelastic material model is adopted to predict the mechanical response of the aortic wall and the two prostheses. The computational results are analyzed in terms of aortic flow, internal wall stresses and aortic wall/prosthesis displacements; a quantitative comparison of the mechanical behavior of the three scenarios is reported. The numerical results highlight a good agreement between stentless and native displacements and internal wall stresses, whereas higher/non-physiological stresses are found for the stented case.

Interpretation of the lime column penetration test

NASA Astrophysics Data System (ADS)

Liyanapathirana, D. S.; Kelly, R. B.

2010-06-01

Dry soil mix (DSM) columns are used to reduce the settlement and to improve the stability of embankments constructed on soft clays. During construction the shear strength of the columns needs to be confirmed for compliance with technical assumptions. A specialized blade shaped penetrometer known as the lime column probe, has been developed for testing DSM columns. This test can be carried out as a pull out resistance test (PORT) or a push in resistance test (PIRT). The test is considered to be more representative of average column shear strength than methods that test only a limited area of the column. Both PORT and PIRT tests require empirical correlations of measured resistance to an absolute measure of shear strength, in a similar manner to the cone penetration test. In this paper, finite element method is used to assess the probe factor, N, for the PORT test. Due to the large soil deformations around the probe, an Arbitrary Lagrangian Eulerian (ALE) based finite element formulation has been used. Variation of N with rigidity index and the friction at the probe-soil interface are investigated to establish a range for the probe factor.

Fluid-structure interaction analysis on the effect of vessel wall hypertrophy and stiffness on the blood flow in carotid artery bifurcation

NASA Astrophysics Data System (ADS)

Lee, Sang Hoon; Choi, Hyoung Gwon; Yoo, Jung Yul

2012-11-01

The effect of artery wall hypertrophy and stiffness on the flow field is investigated using three-dimensional finite element method for simulating the blood flow. To avoid the complexity due to the necessity of additional mechanical constraints, we use the combined formulation which includes both the fluid and structural equations of motion into single coupled variational equation. A P2P1 Galerkin finite element method is used to solve the Navier-Stokes equations for fluid flow and arbitrary Lagrangian-Eulerian formulation is used to achieve mesh movement. The Newmark method is employed for solving the dynamic equilibrium equations for linear elastic solid mechanics. The pulsatile, incompressible flows of Newtonian fluids constrained in the flexible wall are analyzed with Womersley velocity profile at the inlet and constant pressure at the outlet. The study shows that the stiffness of carotid artery wall affects significantly the flow phenomena during the pulse cycle. Similarly, it is found that the flow field is also strongly influenced by wall hypertrophy. This work was supported by Mid-career Researcher Program and Priority Research Centers Program through the National Research Foundation of Korea (NRF) funded by the Ministry of Education, Science and Technology (2009-0079936 & 2011-0029613).

Mass and tracer transport within oceanic Lagrangian coherent vortices as diagnosed in a global mesoscale eddying climate model

NASA Astrophysics Data System (ADS)

Tarshish, Nathaniel; Abernathey, Ryan; Dufour, Carolina; Frenger, Ivy; Griffies, Stephen

2017-04-01

Transient ocean mesoscale fluctuations play a central role in the global climate system, transporting climate relevant tracers such as heat and carbon. In satellite observations and numerical simulations, mesoscale vortices feature prominently as collectively rotating regions that remain visibly coherent. Prior studies on transport from ocean vortices typically rely on Eulerian identification methods, in which vortices are identified by selecting closed contours of Eulerian fields (e.g. sea surface height, or the Okubo-Weiss parameter) that satisfy geometric criteria and anomaly thresholds. In contrast, recent studies employ Lagrangian analysis of virtual particle trajectories initialized within the selected Eulerian contours, revealing significant discrepancies between the advection of the contour's material interior and the evolution of the Eulerian field contour. This work investigates the global mass and tracer transport associated with materially coherent surface ocean vortices. Further, it addresses differences between Eulerian and Lagrangian analyses for the detection of vortices. To do so, we use GFDL's CM2.6 coupled climate model with 5-10km horizontal grid spacing. We identify coherent vortices in CM2.6 by implementing the Rotationally Coherent Lagrangian Vortex (RCLV) framework, which recently emerged from dynamical systems theory. This approach involves the numerical advection of millions of Lagrangian particles and guarantees material coherence by construction. We compute the statistics, spatial distribution, and lifetimes of coherent vortices in addition to calculating the associated mass and tracer transports. We offer compelling evidence that Eulerian vortex methods are poorly suited to answer questions of mass and tracer transport.

Hybrid finite difference/finite element immersed boundary method.

PubMed

E Griffith, Boyce; Luo, Xiaoyu

2017-12-01

The immersed boundary method is an approach to fluid-structure interaction that uses a Lagrangian description of the structural deformations, stresses, and forces along with an Eulerian description of the momentum, viscosity, and incompressibility of the fluid-structure system. The original immersed boundary methods described immersed elastic structures using systems of flexible fibers, and even now, most immersed boundary methods still require Lagrangian meshes that are finer than the Eulerian grid. This work introduces a coupling scheme for the immersed boundary method to link the Lagrangian and Eulerian variables that facilitates independent spatial discretizations for the structure and background grid. This approach uses a finite element discretization of the structure while retaining a finite difference scheme for the Eulerian variables. We apply this method to benchmark problems involving elastic, rigid, and actively contracting structures, including an idealized model of the left ventricle of the heart. Our tests include cases in which, for a fixed Eulerian grid spacing, coarser Lagrangian structural meshes yield discretization errors that are as much as several orders of magnitude smaller than errors obtained using finer structural meshes. The Lagrangian-Eulerian coupling approach developed in this work enables the effective use of these coarse structural meshes with the immersed boundary method. This work also contrasts two different weak forms of the equations, one of which is demonstrated to be more effective for the coarse structural discretizations facilitated by our coupling approach. © 2017 The Authors International  Journal  for  Numerical  Methods  in  Biomedical  Engineering Published by John Wiley & Sons Ltd.

Laser Ray Tracing in a Parallel Arbitrary Lagrangian-Eulerian Adaptive Mesh Refinement Hydrocode

DOE Office of Scientific and Technical Information (OSTI.GOV)

Masters, N D; Kaiser, T B; Anderson, R W

2009-09-28

ALE-AMR is a new hydrocode that we are developing as a predictive modeling tool for debris and shrapnel formation in high-energy laser experiments. In this paper we present our approach to implementing laser ray-tracing in ALE-AMR. We present the equations of laser ray tracing, our approach to efficient traversal of the adaptive mesh hierarchy in which we propagate computational rays through a virtual composite mesh consisting of the finest resolution representation of the modeled space, and anticipate simulations that will be compared to experiments for code validation.

On tide-induced Lagrangian residual current and residual transport: 1. Lagrangian residual current

USGS Publications Warehouse

Feng, Shizuo; Cheng, Ralph T.; Pangen, Xi

1986-01-01

Residual currents in tidal estuaries and coastal embayments have been recognized as fundamental factors which affect the long-term transport processes. It has been pointed out by previous studies that it is more relevant to use a Lagrangian mean velocity than an Eulerian mean velocity to determine the movements of water masses. Under weakly nonlinear approximation, the parameter k, which is the ratio of the net displacement of a labeled water mass in one tidal cycle to the tidal excursion, is assumed to be small. Solutions for tides, tidal current, and residual current have been considered for two-dimensional, barotropic estuaries and coastal seas. Particular attention has been paid to the distinction between the Lagrangian and Eulerian residual currents. When k is small, the first-order Lagrangian residual is shown to be the sum of the Eulerian residual current and the Stokes drift. The Lagrangian residual drift velocity or the second-order Lagrangian residual current has been shown to be dependent on the phase of tidal current. The Lagrangian drift velocity is induced by nonlinear interactions between tides, tidal currents, and the first-order residual currents, and it takes the form of an ellipse on a hodograph plane. Several examples are given to further demonstrate the unique properties of the Lagrangian residual current.

Two-dimensional simulation by regularization of free surface viscoplastic flows with Drucker-Prager yield stress and application to granular collapse

NASA Astrophysics Data System (ADS)

Lusso, Christelle; Ern, Alexandre; Bouchut, François; Mangeney, Anne; Farin, Maxime; Roche, Olivier

2017-03-01

This work is devoted to numerical modeling and simulation of granular flows relevant to geophysical flows such as avalanches and debris flows. We consider an incompressible viscoplastic fluid, described by a rheology with pressure-dependent yield stress, in a 2D setting with a free surface. We implement a regularization method to deal with the singularity of the rheological law, using a mixed finite element approximation of the momentum and incompressibility equations, and an arbitrary Lagrangian Eulerian (ALE) formulation for the displacement of the domain. The free surface is evolved by taking care of its deposition onto the bottom and of preventing it from folding over itself. Several tests are performed to assess the efficiency of our method. The first test is dedicated to verify its accuracy and cost on a one-dimensional simple shear plug flow. On this configuration we setup rules for the choice of the numerical parameters. The second test aims to compare the results of our numerical method to those predicted by an augmented Lagrangian formulation in the case of the collapse and spreading of a granular column over a horizontal rigid bed. Finally we show the reliability of our method by comparing numerical predictions to data from experiments of granular collapse of both trapezoidal and rectangular columns over horizontal rigid or erodible granular bed made of the same material. We compare the evolution of the free surface, the velocity profiles, and the static-flowing interface. The results show the ability of our method to deal numerically with the front behavior of granular collapses over an erodible bed.

Quasiconservation laws for compressible three-dimensional Navier-Stokes flow.

PubMed

Gibbon, J D; Holm, D D

2012-10-01

We formulate the quasi-Lagrangian fluid transport dynamics of mass density ρ and the projection q=ω·∇ρ of the vorticity ω onto the density gradient, as determined by the three-dimensional compressible Navier-Stokes equations for an ideal gas, although the results apply for an arbitrary equation of state. It turns out that the quasi-Lagrangian transport of q cannot cross a level set of ρ. That is, in this formulation, level sets of ρ (isopycnals) are impermeable to the transport of the projection q.

ALEGRA -- A massively parallel h-adaptive code for solid dynamics

DOE Office of Scientific and Technical Information (OSTI.GOV)

Summers, R.M.; Wong, M.K.; Boucheron, E.A.

1997-12-31

ALEGRA is a multi-material, arbitrary-Lagrangian-Eulerian (ALE) code for solid dynamics designed to run on massively parallel (MP) computers. It combines the features of modern Eulerian shock codes, such as CTH, with modern Lagrangian structural analysis codes using an unstructured grid. ALEGRA is being developed for use on the teraflop supercomputers to conduct advanced three-dimensional (3D) simulations of shock phenomena important to a variety of systems. ALEGRA was designed with the Single Program Multiple Data (SPMD) paradigm, in which the mesh is decomposed into sub-meshes so that each processor gets a single sub-mesh with approximately the same number of elements. Usingmore » this approach the authors have been able to produce a single code that can scale from one processor to thousands of processors. A current major effort is to develop efficient, high precision simulation capabilities for ALEGRA, without the computational cost of using a global highly resolved mesh, through flexible, robust h-adaptivity of finite elements. H-adaptivity is the dynamic refinement of the mesh by subdividing elements, thus changing the characteristic element size and reducing numerical error. The authors are working on several major technical challenges that must be met to make effective use of HAMMER on MP computers.« less

Simulating Space Capsule Water Landing with Explicit Finite Element Method

NASA Technical Reports Server (NTRS)

Wang, John T.; Lyle, Karen H.

2007-01-01

A study of using an explicit nonlinear dynamic finite element code for simulating the water landing of a space capsule was performed. The finite element model contains Lagrangian shell elements for the space capsule and Eulerian solid elements for the water and air. An Arbitrary Lagrangian Eulerian (ALE) solver and a penalty coupling method were used for predicting the fluid and structure interaction forces. The space capsule was first assumed to be rigid, so the numerical results could be correlated with closed form solutions. The water and air meshes were continuously refined until the solution was converged. The converged maximum deceleration predicted is bounded by the classical von Karman and Wagner solutions and is considered to be an adequate solution. The refined water and air meshes were then used in the models for simulating the water landing of a capsule model that has a flexible bottom. For small pitch angle cases, the maximum deceleration from the flexible capsule model was found to be significantly greater than the maximum deceleration obtained from the corresponding rigid model. For large pitch angle cases, the difference between the maximum deceleration of the flexible model and that of its corresponding rigid model is smaller. Test data of Apollo space capsules with a flexible heat shield qualitatively support the findings presented in this paper.

A Combined Eulerian-Lagrangian Data Representation for Large-Scale Applications.

PubMed

Sauer, Franz; Xie, Jinrong; Ma, Kwan-Liu

2017-10-01

The Eulerian and Lagrangian reference frames each provide a unique perspective when studying and visualizing results from scientific systems. As a result, many large-scale simulations produce data in both formats, and analysis tasks that simultaneously utilize information from both representations are becoming increasingly popular. However, due to their fundamentally different nature, drawing correlations between these data formats is a computationally difficult task, especially in a large-scale setting. In this work, we present a new data representation which combines both reference frames into a joint Eulerian-Lagrangian format. By reorganizing Lagrangian information according to the Eulerian simulation grid into a "unit cell" based approach, we can provide an efficient out-of-core means of sampling, querying, and operating with both representations simultaneously. We also extend this design to generate multi-resolution subsets of the full data to suit the viewer's needs and provide a fast flow-aware trajectory construction scheme. We demonstrate the effectiveness of our method using three large-scale real world scientific datasets and provide insight into the types of performance gains that can be achieved.

Eulerian-Lagrangian numerical scheme for simulating advection, dispersion, and transient storage in streams and a comparison of numerical methods

USGS Publications Warehouse

Cox, T.J.; Runkel, R.L.

2008-01-01

Past applications of one-dimensional advection, dispersion, and transient storage zone models have almost exclusively relied on a central differencing, Eulerian numerical approximation to the nonconservative form of the fundamental equation. However, there are scenarios where this approach generates unacceptable error. A new numerical scheme for this type of modeling is presented here that is based on tracking Lagrangian control volumes across a fixed (Eulerian) grid. Numerical tests are used to provide a direct comparison of the new scheme versus nonconservative Eulerian numerical methods, in terms of both accuracy and mass conservation. Key characteristics of systems for which the Lagrangian scheme performs better than the Eulerian scheme include: nonuniform flow fields, steep gradient plume fronts, and pulse and steady point source loadings in advection-dominated systems. A new analytical derivation is presented that provides insight into the loss of mass conservation in the nonconservative Eulerian scheme. This derivation shows that loss of mass conservation in the vicinity of spatial flow changes is directly proportional to the lateral inflow rate and the change in stream concentration due to the inflow. While the nonconservative Eulerian scheme has clearly worked well for past published applications, it is important for users to be aware of the scheme's limitations. ?? 2008 ASCE.

Eulerian and Lagrangian Parameterization of the Oceanic Mixed Layer using Large Eddy Simulation and MPAS-Ocean

DOE Office of Scientific and Technical Information (OSTI.GOV)

Van Roekel, Luke

We have conducted a suite of Large Eddy Simulation (LES) to form the basis of a multi-model comparison (left). The results have led to proposed model improvements. We have verified that Eulerian-Lagrangian effective diffusivity estimates of mesoscale mixing are consistent with traditional particle statistics metrics (right). LES and Lagrangian particles will be utilized to better represent the movement of water into and out of the mixed layer.

Bayesian Lagrangian Data Assimilation and Drifter Deployment Strategies

NASA Astrophysics Data System (ADS)

Dutt, A.; Lermusiaux, P. F. J.

2017-12-01

Ocean currents transport a variety of natural (e.g. water masses, phytoplankton, zooplankton, sediments, etc.) and man-made materials and other objects (e.g. pollutants, floating debris, search and rescue, etc.). Lagrangian Coherent Structures (LCSs) or the most influential/persistent material lines in a flow, provide a robust approach to characterize such Lagrangian transports and organize classic trajectories. Using the flow-map stochastic advection and a dynamically-orthogonal decomposition, we develop uncertainty prediction schemes for both Eulerian and Lagrangian variables. We then extend our Bayesian Gaussian Mixture Model (GMM)-DO filter to a joint Eulerian-Lagrangian Bayesian data assimilation scheme. The resulting nonlinear filter allows the simultaneous non-Gaussian estimation of Eulerian variables (e.g. velocity, temperature, salinity, etc.) and Lagrangian variables (e.g. drifter/float positions, trajectories, LCSs, etc.). Its results are showcased using a double-gyre flow with a random frequency, a stochastic flow past a cylinder, and realistic ocean examples. We further show how our Bayesian mutual information and adaptive sampling equations provide a rigorous efficient methodology to plan optimal drifter deployment strategies and predict the optimal times, locations, and types of measurements to be collected.

Assimilating Eulerian and Lagrangian data in traffic-flow models

NASA Astrophysics Data System (ADS)

Xia, Chao; Cochrane, Courtney; DeGuire, Joseph; Fan, Gaoyang; Holmes, Emma; McGuirl, Melissa; Murphy, Patrick; Palmer, Jenna; Carter, Paul; Slivinski, Laura; Sandstede, Björn

2017-05-01

Data assimilation of traffic flow remains a challenging problem. One difficulty is that data come from different sources ranging from stationary sensors and camera data to GPS and cell phone data from moving cars. Sensors and cameras give information about traffic density, while GPS data provide information about the positions and velocities of individual cars. Previous methods for assimilating Lagrangian data collected from individual cars relied on specific properties of the underlying computational model or its reformulation in Lagrangian coordinates. These approaches make it hard to assimilate both Eulerian density and Lagrangian positional data simultaneously. In this paper, we propose an alternative approach that allows us to assimilate both Eulerian and Lagrangian data. We show that the proposed algorithm is accurate and works well in different traffic scenarios and regardless of whether ensemble Kalman or particle filters are used. We also show that the algorithm is capable of estimating parameters and assimilating real traffic observations and synthetic observations obtained from microscopic models.

The Case for Including Eulerian Kinematics in Undergraduate Dynamics.

ERIC Educational Resources Information Center

Uram, Earl M.

A Eulerian framework is proposed as an alternative to the Lagrangian framework usually used in undergraduate dynamics courses. An attempt to introduce Eulerian kinematics into a dynamics course is discussed. (LMH)

A finite-volume Eulerian-Lagrangian Localized Adjoint Method for solution of the advection-dispersion equation

USGS Publications Warehouse

Healy, R.W.; Russell, T.F.

1993-01-01

A new mass-conservative method for solution of the one-dimensional advection-dispersion equation is derived and discussed. Test results demonstrate that the finite-volume Eulerian-Lagrangian localized adjoint method (FVELLAM) outperforms standard finite-difference methods, in terms of accuracy and efficiency, for solute transport problems that are dominated by advection. For dispersion-dominated problems, the performance of the method is similar to that of standard methods. Like previous ELLAM formulations, FVELLAM systematically conserves mass globally with all types of boundary conditions. FVELLAM differs from other ELLAM approaches in that integrated finite differences, instead of finite elements, are used to approximate the governing equation. This approach, in conjunction with a forward tracking scheme, greatly facilitates mass conservation. The mass storage integral is numerically evaluated at the current time level, and quadrature points are then tracked forward in time to the next level. Forward tracking permits straightforward treatment of inflow boundaries, thus avoiding the inherent problem in backtracking, as used by most characteristic methods, of characteristic lines intersecting inflow boundaries. FVELLAM extends previous ELLAM results by obtaining mass conservation locally on Lagrangian space-time elements. Details of the integration, tracking, and boundary algorithms are presented. Test results are given for problems in Cartesian and radial coordinates.

SANTA BARBARA CLUSTER COMPARISON TEST WITH DISPH

DOE Office of Scientific and Technical Information (OSTI.GOV)

Saitoh, Takayuki R.; Makino, Junichiro, E-mail: saitoh@elsi.jp

2016-06-01

The Santa Barbara cluster comparison project revealed that there is a systematic difference between entropy profiles of clusters of galaxies obtained by Eulerian mesh and Lagrangian smoothed particle hydrodynamics (SPH) codes: mesh codes gave a core with a constant entropy, whereas SPH codes did not. One possible reason for this difference is that mesh codes are not Galilean invariant. Another possible reason is the problem of the SPH method, which might give too much “protection” to cold clumps because of the unphysical surface tension induced at contact discontinuities. In this paper, we apply the density-independent formulation of SPH (DISPH), whichmore » can handle contact discontinuities accurately, to simulations of a cluster of galaxies and compare the results with those with the standard SPH. We obtained the entropy core when we adopt DISPH. The size of the core is, however, significantly smaller than those obtained with mesh simulations and is comparable to those obtained with quasi-Lagrangian schemes such as “moving mesh” and “mesh free” schemes. We conclude that both the standard SPH without artificial conductivity and Eulerian mesh codes have serious problems even with such an idealized simulation, while DISPH, SPH with artificial conductivity, and quasi-Lagrangian schemes have sufficient capability to deal with it.« less

Hydrodynamic interactions for complex-shaped nanocarriers in targeted drug delivery

NASA Astrophysics Data System (ADS)

Wang, Yaohong; Eckmann, David; Radhakrishnan, Ravi; Ayyaswamy, Portonovo

2014-11-01

Nanocarrier motion in a blood vessel involves hydrodynamic and Brownian interactions, which collectively dictate the efficacy in targeted drug delivery. The shape of nanocarriers plays a crucial role in drug delivery. In order to quantify the flow and association properties of elliptical nanoparticles, we have developed an arbitrary Lagrangian-Eulerian framework with capabilities to simulate the hydrodynamic motion of nanoparticles of arbitrary shapes. We introduce the quaternions for rotational motion, and two collision models, namely, (a) an impulse-based model for wall-particle collision, and (b) the short-range repulsive Gay-Berne potential for particle-particle collision. We also study the red blood cell and nanocarrier (such as ellipsoid) interactions. We compare our results with those obtained for a hard sphere model for both RBCs and nanocarriers. Supported by NIH through grant U01-EB016027.

Computing eddy-driven effective diffusivity using Lagrangian particles

DOE PAGES

Wolfram, Phillip J.; Ringler, Todd D.

2017-08-14

A novel method to derive effective diffusivity from Lagrangian particle trajectory data sets is developed and then analyzed relative to particle-derived meridional diffusivity for eddy-driven mixing in an idealized circumpolar current. Quantitative standard dispersion- and transport-based mixing diagnostics are defined, compared and contrasted to motivate the computation and use of effective diffusivity derived from Lagrangian particles. We compute the effective diffusivity by first performing scalar transport on Lagrangian control areas using stored trajectories computed from online Lagrangian In-situ Global High-performance particle Tracking (LIGHT) using the Model for Prediction Across Scales Ocean (MPAS-O). Furthermore, the Lagrangian scalar transport scheme is comparedmore » against an Eulerian scalar transport scheme. Spatially-variable effective diffusivities are computed from resulting time-varying cumulative concentrations that vary as a function of cumulative area. The transport-based Eulerian and Lagrangian effective diffusivity diagnostics are found to be qualitatively consistent with the dispersion-based diffusivity. All diffusivity estimates show a region of increased subsurface diffusivity within the core of an idealized circumpolar current and results are within a factor of two of each other. The Eulerian and Lagrangian effective diffusivities are most similar; smaller and more spatially diffused values are obtained with the dispersion-based diffusivity computed with particle clusters.« less

Computing eddy-driven effective diffusivity using Lagrangian particles

DOE Office of Scientific and Technical Information (OSTI.GOV)

Wolfram, Phillip J.; Ringler, Todd D.

A novel method to derive effective diffusivity from Lagrangian particle trajectory data sets is developed and then analyzed relative to particle-derived meridional diffusivity for eddy-driven mixing in an idealized circumpolar current. Quantitative standard dispersion- and transport-based mixing diagnostics are defined, compared and contrasted to motivate the computation and use of effective diffusivity derived from Lagrangian particles. We compute the effective diffusivity by first performing scalar transport on Lagrangian control areas using stored trajectories computed from online Lagrangian In-situ Global High-performance particle Tracking (LIGHT) using the Model for Prediction Across Scales Ocean (MPAS-O). Furthermore, the Lagrangian scalar transport scheme is comparedmore » against an Eulerian scalar transport scheme. Spatially-variable effective diffusivities are computed from resulting time-varying cumulative concentrations that vary as a function of cumulative area. The transport-based Eulerian and Lagrangian effective diffusivity diagnostics are found to be qualitatively consistent with the dispersion-based diffusivity. All diffusivity estimates show a region of increased subsurface diffusivity within the core of an idealized circumpolar current and results are within a factor of two of each other. The Eulerian and Lagrangian effective diffusivities are most similar; smaller and more spatially diffused values are obtained with the dispersion-based diffusivity computed with particle clusters.« less

Uncertainty quantification in Eulerian-Lagrangian models for particle-laden flows

NASA Astrophysics Data System (ADS)

Fountoulakis, Vasileios; Jacobs, Gustaaf; Udaykumar, Hs

2017-11-01

A common approach to ameliorate the computational burden in simulations of particle-laden flows is to use a point-particle based Eulerian-Lagrangian model, which traces individual particles in their Lagrangian frame and models particles as mathematical points. The particle motion is determined by Stokes drag law, which is empirically corrected for Reynolds number, Mach number and other parameters. The empirical corrections are subject to uncertainty. Treating them as random variables renders the coupled system of PDEs and ODEs stochastic. An approach to quantify the propagation of this parametric uncertainty to the particle solution variables is proposed. The approach is based on averaging of the governing equations and allows for estimation of the first moments of the quantities of interest. We demonstrate the feasibility of our proposed methodology of uncertainty quantification of particle-laden flows on one-dimensional linear and nonlinear Eulerian-Lagrangian systems. This research is supported by AFOSR under Grant FA9550-16-1-0008.

The Role of Multiphysics Simulation in Multidisciplinary Analysis

NASA Technical Reports Server (NTRS)

Rifai, Steven M.; Ferencz, Robert M.; Wang, Wen-Ping; Spyropoulos, Evangelos T.; Lawrence, Charles; Melis, Matthew E.

1998-01-01

This article describes the applications of the Spectrum(Tm) Solver in Multidisciplinary Analysis (MDA). Spectrum, a multiphysics simulation software based on the finite element method, addresses compressible and incompressible fluid flow, structural, and thermal modeling as well as the interaction between these disciplines. Multiphysics simulation is based on a single computational framework for the modeling of multiple interacting physical phenomena. Interaction constraints are enforced in a fully-coupled manner using the augmented-Lagrangian method. Within the multiphysics framework, the finite element treatment of fluids is based on Galerkin-Least-Squares (GLS) method with discontinuity capturing operators. The arbitrary-Lagrangian-Eulerian method is utilized to account for deformable fluid domains. The finite element treatment of solids and structures is based on the Hu-Washizu variational principle. The multiphysics architecture lends itself naturally to high-performance parallel computing. Aeroelastic, propulsion, thermal management and manufacturing applications are presented.

A volume-of-fluid method for simulation of compressible axisymmetric multi-material flow

NASA Astrophysics Data System (ADS)

de Niem, D.; Kührt, E.; Motschmann, U.

2007-02-01

A two-dimensional Eulerian hydrodynamic method for the numerical simulation of inviscid compressible axisymmetric multi-material flow in external force fields for the situation of pure fluids separated by macroscopic interfaces is presented. The method combines an implicit Lagrangian step with an explicit Eulerian advection step. Individual materials obey separate energy equations, fulfill general equations of state, and may possess different temperatures. Material volume is tracked using a piecewise linear volume-of-fluid method. An overshoot-free logically simple and economic material advection algorithm for cylinder coordinates is derived, in an algebraic formulation. New aspects arising in the case of more than two materials such as the material ordering strategy during transport are presented. One- and two-dimensional numerical examples are given.

Computations of spray, fuel-air mixing, and combustion in a lean-premixed-prevaporized combustor

NASA Technical Reports Server (NTRS)

Dasgupta, A.; Li, Z.; Shih, T. I.-P.; Kundu, K.; Deur, J. M.

1993-01-01

A code was developed for computing the multidimensional flow, spray, combustion, and pollutant formation inside gas turbine combustors. The code developed is based on a Lagrangian-Eulerian formulation and utilizes an implicit finite-volume method. The focus of this paper is on the spray part of the code (both formulation and algorithm), and a number of issues related to the computation of sprays and fuel-air mixing in a lean-premixed-prevaporized combustor. The issues addressed include: (1) how grid spacings affect the diffusion of evaporated fuel, and (2) how spurious modes can arise through modelling of the spray in the Lagrangian computations. An upwind interpolation scheme is proposed to account for some effects of grid spacing on the artificial diffusion of the evaporated fuel. Also, some guidelines are presented to minimize errors associated with the spurious modes.

Biomechanical responses of a pig head under blast loading: a computational simulation.

PubMed

Zhu, Feng; Skelton, Paul; Chou, Cliff C; Mao, Haojie; Yang, King H; King, Albert I

2013-03-01

A series of computational studies were performed to investigate the biomechanical responses of the pig head under a specific shock tube environment. A finite element model of the head of a 50-kg Yorkshire pig was developed with sufficient details, based on the Lagrangian formulation, and a shock tube model was developed using the multimaterial arbitrary Lagrangian-Eulerian (MMALE) approach. These two models were integrated and a fluid/solid coupling algorithm was used to simulate the interaction of the shock wave with the pig's head. The finite element model-predicted incident and intracranial pressure traces were in reasonable agreement with those obtained experimentally. Using the verified numerical model of the shock tube and pig head, further investigations were carried out to study the spatial and temporal distributions of pressure, shear stress, and principal strain within the head. Pressure enhancement was found in the skull, which is believed to be caused by shock wave reflection at the interface of the materials with distinct wave impedances. Brain tissue has a shock attenuation effect and larger pressures were observed in the frontal and occipital regions, suggesting a greater possibility of coup and contrecoup contusion. Shear stresses in the brain and deflection in the skull remained at a low level. Higher principal strains were observed in the brain near the foramen magnum, suggesting that there is a greater chance of cellular or vascular injuries in the brainstem region. Copyright © 2012 John Wiley & Sons, Ltd.

Ice Accretion Modeling using an Eulerian Approach for Droplet Impingement

NASA Technical Reports Server (NTRS)

Kim, Joe Woong; Garza, Dennis P.; Sankar, Lakshmi N.; Kreeger, Richard E.

2012-01-01

A three-dimensional Eulerian analysis has been developed for modeling droplet impingement on lifting bodes. The Eulerian model solves the conservation equations of mass and momentum to obtain the droplet flow field properties on the same mesh used in CFD simulations. For complex configurations such as a full rotorcraft, the Eulerian approach is more efficient because the Lagrangian approach would require a significant amount of seeding for accurate estimates of collection efficiency. Simulations are done for various benchmark cases such as NACA0012 airfoil, MS317 airfoil and oscillating SC2110 airfoil to illustrate its use. The present results are compared with results from the Lagrangian approach used in an industry standard analysis called LEWICE.

A combined Eulerian-volume of fraction-Lagrangian method for atomization simulation

NASA Technical Reports Server (NTRS)

Seung, S. P.; Chen, C. P.; Ziebarth, John P.

1994-01-01

The tracking of free surfaces between liquid and gas phases and analysis of the interfacial phenomena between the two during the atomization and breakup process of a liquid fuel jet is modeled. Numerical modeling of liquid-jet atomization requires the resolution of different conservation equations. Detailed formulation and validation are presented for the confined dam broken problem, the water surface problem, the single droplet problem, a jet breakup problem, and the liquid column instability problem.

Dynamic Deployment Simulations of Inflatable Space Structures

NASA Technical Reports Server (NTRS)

Wang, John T.

2005-01-01

The feasibility of using Control Volume (CV) method and the Arbitrary Lagrangian Eulerian (ALE) method in LSDYNA to simulate the dynamic deployment of inflatable space structures is investigated. The CV and ALE methods were used to predict the inflation deployments of three folded tube configurations. The CV method was found to be a simple and computationally efficient method that may be adequate for modeling slow inflation deployment sine the inertia of the inflation gas can be neglected. The ALE method was found to be very computationally intensive since it involves the solving of three conservative equations of fluid as well as dealing with complex fluid structure interactions.

Lagrangian averages, averaged Lagrangians, and the mean effects of fluctuations in fluid dynamics.

PubMed

Holm, Darryl D.

2002-06-01

We begin by placing the generalized Lagrangian mean (GLM) equations for a compressible adiabatic fluid into the Euler-Poincare (EP) variational framework of fluid dynamics, for an averaged Lagrangian. This is the Lagrangian averaged Euler-Poincare (LAEP) theorem. Next, we derive a set of approximate small amplitude GLM equations (glm equations) at second order in the fluctuating displacement of a Lagrangian trajectory from its mean position. These equations express the linear and nonlinear back-reaction effects on the Eulerian mean fluid quantities by the fluctuating displacements of the Lagrangian trajectories in terms of their Eulerian second moments. The derivation of the glm equations uses the linearized relations between Eulerian and Lagrangian fluctuations, in the tradition of Lagrangian stability analysis for fluids. The glm derivation also uses the method of averaged Lagrangians, in the tradition of wave, mean flow interaction. Next, the new glm EP motion equations for incompressible ideal fluids are compared with the Euler-alpha turbulence closure equations. An alpha model is a GLM (or glm) fluid theory with a Taylor hypothesis closure. Such closures are based on the linearized fluctuation relations that determine the dynamics of the Lagrangian statistical quantities in the Euler-alpha equations. Thus, by using the LAEP theorem, we bridge between the GLM equations and the Euler-alpha closure equations, through the small-amplitude glm approximation in the EP variational framework. We conclude by highlighting a new application of the GLM, glm, and alpha-model results for Lagrangian averaged ideal magnetohydrodynamics. (c) 2002 American Institute of Physics.

Modeling NIF experimental designs with adaptive mesh refinement and Lagrangian hydrodynamics

NASA Astrophysics Data System (ADS)

Koniges, A. E.; Anderson, R. W.; Wang, P.; Gunney, B. T. N.; Becker, R.; Eder, D. C.; MacGowan, B. J.; Schneider, M. B.

2006-06-01

Incorporation of adaptive mesh refinement (AMR) into Lagrangian hydrodynamics algorithms allows for the creation of a highly powerful simulation tool effective for complex target designs with three-dimensional structure. We are developing an advanced modeling tool that includes AMR and traditional arbitrary Lagrangian-Eulerian (ALE) techniques. Our goal is the accurate prediction of vaporization, disintegration and fragmentation in National Ignition Facility (NIF) experimental target elements. Although our focus is on minimizing the generation of shrapnel in target designs and protecting the optics, the general techniques are applicable to modern advanced targets that include three-dimensional effects such as those associated with capsule fill tubes. Several essential computations in ordinary radiation hydrodynamics need to be redesigned in order to allow for AMR to work well with ALE, including algorithms associated with radiation transport. Additionally, for our goal of predicting fragmentation, we include elastic/plastic flow into our computations. We discuss the integration of these effects into a new ALE-AMR simulation code. Applications of this newly developed modeling tool as well as traditional ALE simulations in two and three dimensions are applied to NIF early-light target designs.

On inter-tidal transport equation

USGS Publications Warehouse

Cheng, Ralph T.; Feng, Shizuo; Pangen, Xi

1989-01-01

The transports of solutes, sediments, nutrients, and other tracers are fundamental to the interactive physical, chemical, and biological processes in estuaries. The characteristic time scales for most estuarine biological and chemical processes are on the order of several tidal cycles or longer. To address the long-term transport mechanism meaningfully, the formulation of an inter-tidal conservation equation is the main subject of this paper. The commonly used inter-tidal conservation equation takes the form of a convection-dispersion equation in which the convection is represented by the Eulerian residual current, and the dispersion terms are due to the introduction of a Fickian hypothesis, unfortunately, the physical significance of this equation is not clear, and the introduction of a Fickian hypothesis is at best an ad hoc approximation. Some recent research results on the Lagrangian residual current suggest that the long-term transport problem is more closely related to the Lagrangian residual current than to the Eulerian residual current. With the aid of additional insight of residual current, the inter-tidal transport equation has been reformulated in this paper using a small perturbation method for a weakly nonlinear tidal system. When tidal flows can be represented by an M2 system, the new intertidal transport equation also takes the form of a convective-dispersion equation without the introduction of a Fickian hypothesis. The convective velocity turns out to be the first order Lagrangian residual current (the sum of the Eulerian residual current and the Stokes’ drift), and the correlation terms take the form of convection with the Stokes’ drift as the convective velocity. The remaining dispersion terms are perturbations of lower order solution to higher order solutions due to shear effect and turbulent mixing.

An incompressible two-dimensional multiphase particle-in-cell model for dense particle flows

DOE Office of Scientific and Technical Information (OSTI.GOV)

Snider, D.M.; O`Rourke, P.J.; Andrews, M.J.

1997-06-01

A two-dimensional, incompressible, multiphase particle-in-cell (MP-PIC) method is presented for dense particle flows. The numerical technique solves the governing equations of the fluid phase using a continuum model and those of the particle phase using a Lagrangian model. Difficulties associated with calculating interparticle interactions for dense particle flows with volume fractions above 5% have been eliminated by mapping particle properties to a Eulerian grid and then mapping back computed stress tensors to particle positions. This approach utilizes the best of Eulerian/Eulerian continuum models and Eulerian/Lagrangian discrete models. The solution scheme allows for distributions of types, sizes, and density of particles,more » with no numerical diffusion from the Lagrangian particle calculations. The computational method is implicit with respect to pressure, velocity, and volume fraction in the continuum solution thus avoiding courant limits on computational time advancement. MP-PIC simulations are compared with one-dimensional problems that have analytical solutions and with two-dimensional problems for which there are experimental data.« less

Acoustic Radiation Pressure

NASA Technical Reports Server (NTRS)

Cantrell, John H.

2018-01-01

The theoretical foundation of acoustic radiation pressure in plane wave beams is reexamined. It is shown from finite deformation theory and the Boltzmann-Ehrenfest Adiabatic Principle that the Brillouin stress tensor (BST) is the radiation stress in Lagrangian coordinates (not Eulerian coordinates) and that the terms in the BST are not the momentum flux density and mean excess Eulerian stress but are simply contributions to the variation in the wave oscillation period resulting from changes in path length and true wave velocity, respectively, from virtual variations in the strain. It is shown that the radiation stress in Eulerian coordinates is the mean Cauchy stress (not the momentum flux density, as commonly assumed) and that Langevin's second relation does not yield an assessment of the mean Eulerian pressure, since the enthalpy used in the traditional derivations is a function of the thermodynamic tensions - not the Eulerian pressure. It is shown that the transformation between Lagrangian and Eulerian quantities cannot be obtained from the commonly-used expansion of one of the quantities in terms of the particle displacement, since the expansion provides only the difference between the value of the quantity at two different points in Cartesian space separated by the displacement. The proper transformation is obtained only by employing the transformation coefficients of finite deformation theory, which are defined in terms of the displacement gradients. Finite deformation theory leads to the result that for laterally unconfined, plane waves the Lagrangian and Eulerian radiation pressures are equal with the value (1/4)(2K) along the direction of wave propagation, where (K) is the mean kinetic energy density, and zero in directions normal to the propagation direction. This is contrary to the Langevin result that the Lagrangian radiation pressure in the propagation direction is equal to (2K) and the BST result that the Eulerian radiation pressure in that direction is the momentum flux density.

Advances in Quantum Trajectory Approaches to Dynamics

NASA Astrophysics Data System (ADS)

Askar, Attila

2001-03-01

The quantum fluid dynamics (QFD) formulation is based on the separation of the amplitude and phase of the complex wave function in Schrodinger's equation. The approach leads to conservation laws for an equivalent "gas continuum". The Lagrangian [1] representation corresponds to following the particles of the fluid continuum, i. e. calculating "quantum trajectories". The Eulerian [2] representation on the other hand, amounts to observing the dynamics of the gas continuum at the points of a fixed coordinate frame. The combination of several factors leads to a most encouraging computational efficiency. QFD enables the numerical analysis to deal with near monotonic amplitude and phase functions. The Lagrangian description concentrates the computation effort to regions of highest probability as an optimal adaptive grid. The Eulerian representation allows the study of multi-coordinate problems as a set of one-dimensional problems within an alternating direction methodology. An explicit time integrator limits the increase in computational effort with the number of discrete points to linear. Discretization of the space via local finite elements [1,2] and global radial functions [3] will be discussed. Applications include wave packets in four-dimensional quadratic potentials and two coordinate photo-dissociation problems for NOCl and NO2. [1] "Quantum fluid dynamics (QFD) in the Lagrangian representation with applications to photo-dissociation problems", F. Sales, A. Askar and H. A. Rabitz, J. Chem. Phys. 11, 2423 (1999) [2] "Multidimensional wave-packet dynamics within the fluid dynamical formulation of the Schrodinger equation", B. Dey, A. Askar and H. A. Rabitz, J. Chem. Phys. 109, 8770 (1998) [3] "Solution of the quantum fluid dynamics equations with radial basis function interpolation", Xu-Guang Hu, Tak-San Ho, H. A. Rabitz and A. Askar, Phys. Rev. E. 61, 5967 (2000)

Generalized Lee-Wick formulation from higher derivative field theories

DOE Office of Scientific and Technical Information (OSTI.GOV)

Cho, Inyong; Kwon, O-Kab; Department of Physics, BK21 Physics Research Division, Institute of Basic Science, Sungkyunkwan University, Suwon 440-746

2010-07-15

We study a higher derivative (HD) field theory with an arbitrary order of derivative for a real scalar field. The degree of freedom for the HD field can be converted to multiple fields with canonical kinetic terms up to the overall sign. The Lagrangian describing the dynamics of the multiple fields is known as the Lee-Wick (LW) form. The first step to obtain the LW form for a given HD Lagrangian is to find an auxiliary field (AF) Lagrangian which is equivalent to the original HD Lagrangian up to the quantum level. Until now, the AF Lagrangian has been studiedmore » only for N=2 and 3 cases, where N is the number of poles of the two-point function of the HD scalar field. We construct the AF Lagrangian for arbitrary N. By the linear combinations of AF fields, we also obtain the corresponding LW form. We find the explicit mapping matrices among the HD fields, the AF fields, and the LW fields. As an exercise of our construction, we calculate the relations among parameters and mapping matrices for N=2, 3, and 4 cases.« less

Extension of rezoned Eulerian-Lagrangian method to astrophysical plasma applications

NASA Technical Reports Server (NTRS)

Song, M. T.; Wu, S. T.; Dryer, Murray

1993-01-01

The rezoned Eulerian-Lagrangian procedure developed by Brackbill and Pracht (1973), which is limited to simple configurations of the magnetic fields, is modified in order to make it applicable to astrophysical plasma. For this purpose, two specific methods are introduced, which make it possible to determine the initial field topology for which no analytical expressions are available. Numerical examples illustrating these methods are presented.

Large-Deformation Displacement Transfer Functions for Shape Predictions of Highly Flexible Slender Aerospace Structures

NASA Technical Reports Server (NTRS)

Ko, William L.; Fleischer, Van Tran

2013-01-01

Large deformation displacement transfer functions were formulated for deformed shape predictions of highly flexible slender structures like aircraft wings. In the formulation, the embedded beam (depth wise cross section of structure along the surface strain sensing line) was first evenly discretized into multiple small domains, with surface strain sensing stations located at the domain junctures. Thus, the surface strain (bending strains) variation within each domain could be expressed with linear of nonlinear function. Such piecewise approach enabled piecewise integrations of the embedded beam curvature equations [classical (Eulerian), physical (Lagrangian), and shifted curvature equations] to yield closed form slope and deflection equations in recursive forms.

A general higher-order remap algorithm for ALE calculations

DOE Office of Scientific and Technical Information (OSTI.GOV)

Chiravalle, Vincent P

2011-01-05

A numerical technique for solving the equations of fluid dynamics with arbitrary mesh motion is presented. The three phases of the Arbitrary Lagrangian Eulerian (ALE) methodology are outlined: the Lagrangian phase, grid relaxation phase and remap phase. The Lagrangian phase follows a well known approach from the HEMP code; in addition the strain rate andflow divergence are calculated in a consistent manner according to Margolin. A donor cell method from the SALE code forms the basis of the remap step, but unlike SALE a higher order correction based on monotone gradients is also added to the remap. Four test problemsmore » were explored to evaluate the fidelity of these numerical techniques, as implemented in a simple test code, written in the C programming language, called Cercion. Novel cell-centered data structures are used in Cercion to reduce the complexity of the programming and maximize the efficiency of memory usage. The locations of the shock and contact discontinuity in the Riemann shock tube problem are well captured. Cercion demonstrates a high degree of symmetry when calculating the Sedov blast wave solution, with a peak density at the shock front that is similar to the value determined by the RAGE code. For a flyer plate test problem both Cercion and FLAG give virtually the same velocity temporal profile at the target-vacuum interface. When calculating a cylindrical implosion of a steel shell, Cercion and FLAG agree well and the Cercion results are insensitive to the use of ALE.« less

A volume-filtered formulation to capture particle-shock interactions in multiphase compressible flows

NASA Astrophysics Data System (ADS)

Shallcross, Gregory; Capecelatro, Jesse

2017-11-01

Compressible particle-laden flows are common in engineering systems. Applications include but are not limited to water injection in high-speed jet flows for noise suppression, rocket-plume surface interactions during planetary landing, and explosions during coal mining operations. Numerically, it is challenging to capture these interactions due to the wide range of length and time scales. Additionally, there are many forms of the multiphase compressible flow equations with volume fraction effects, some of which are conflicting in nature. The purpose of this presentation is to develop the capability to accurately capture particle-shock interactions in systems with a large number of particles from dense to dilute regimes. A thorough derivation of the volume filtered equations is presented. The volume filtered equations are then implemented in a high-order, energy-stable Eulerian-Lagrangian framework. We show this framework is capable of decoupling the fluid mesh from the particle size, enabling arbitrary particle size distributions in the presence of shocks. The proposed method is then assessed against particle-laden shock tube data. Quantities of interest include fluid-phase pressure profiles and particle spreading rates. The effect of collisions in 2D and 3D are also evaluated.

Unsteady Flow in Different Atmospheric Boundary Layer Regimes and Its Impact on Wind-Turbine Performance

NASA Astrophysics Data System (ADS)

Gohari, Iman; Korobenko, Artem; Yan, Jinhui; Bazilevs, Yuri; Sarkar, Sutanu

2016-11-01

Wind is a renewable energy resource that offers several advantages including low pollutant emission and inexpensive construction. Wind turbines operate in conditions dictated by the Atmospheric Boundary Layer (ABL) and that motivates the study of coupling ABL simulations with wind turbine dynamics. The ABL simulations can be used for realistic modeling of the environment which, with the use of fluid-structure interaction, can give realistic predictions of extracted power, rotor loading, and blade structural response. The ABL simulations provide inflow boundary conditions to the wind-turbine simulator which uses arbitrary Lagrangian-Eulerian variational multiscale formulation. In the present work, ABL simulations are performed to examine two different scenarios: (i) A neutral ABL with zero heat-flux and inversion layer at 350m, in which the wind turbine experiences maximum mean shear; (2) A shallow ABL with the surface cooling-rate of -1 K/hr, in which the wind turbine experiences maximum mean velocity at the low-level-jet nose height. We will discuss differences in the unsteady flow between the two different ABL conditions and their impact on the performance of the wind turbine cluster in the coupled ABL-wind turbine simulations.

Recent advances in high-order WENO finite volume methods for compressible multiphase flows

NASA Astrophysics Data System (ADS)

Dumbser, Michael

2013-10-01

We present two new families of better than second order accurate Godunov-type finite volume methods for the solution of nonlinear hyperbolic partial differential equations with nonconservative products. One family is based on a high order Arbitrary-Lagrangian-Eulerian (ALE) formulation on moving meshes, which allows to resolve the material contact wave in a very sharp way when the mesh is moved at the speed of the material interface. The other family of methods is based on a high order Adaptive Mesh Refinement (AMR) strategy, where the mesh can be strongly refined in the vicinity of the material interface. Both classes of schemes have several building blocks in common, in particular: a high order WENO reconstruction operator to obtain high order of accuracy in space; the use of an element-local space-time Galerkin predictor step which evolves the reconstruction polynomials in time and that allows to reach high order of accuracy in time in one single step; the use of a path-conservative approach to treat the nonconservative terms of the PDE. We show applications of both methods to the Baer-Nunziato model for compressible multiphase flows.

Quantum hydrodynamics: capturing a reactive scattering resonance.

PubMed

Derrickson, Sean W; Bittner, Eric R; Kendrick, Brian K

2005-08-01

The hydrodynamic equations of motion associated with the de Broglie-Bohm formulation of quantum mechanics are solved using a meshless method based upon a moving least-squares approach. An arbitrary Lagrangian-Eulerian frame of reference and a regridding algorithm which adds and deletes computational points are used to maintain a uniform and nearly constant interparticle spacing. The methodology also uses averaged fields to maintain unitary time evolution. The numerical instabilities associated with the formation of nodes in the reflected portion of the wave packet are avoided by adding artificial viscosity to the equations of motion. A new and more robust artificial viscosity algorithm is presented which gives accurate scattering results and is capable of capturing quantum resonances. The methodology is applied to a one-dimensional model chemical reaction that is known to exhibit a quantum resonance. The correlation function approach is used to compute the reactive scattering matrix, reaction probability, and time delay as a function of energy. Excellent agreement is obtained between the scattering results based upon the quantum hydrodynamic approach and those based upon standard quantum mechanics. This is the first clear demonstration of the ability of moving grid approaches to accurately and robustly reproduce resonance structures in a scattering system.

A lagrangian-eulerian description of debris transport by a tsunami in the Lisbon waterfront

NASA Astrophysics Data System (ADS)

Conde, Daniel; Canelas, Ricardo; Baptista, Maria Ana; João Telhado, Maria; Ferreira, Rui M. L.

2013-04-01

Several major tsunamis are known to have struck the Portuguese coast over the past millennia (Baptista and Miranda, 2009). The Tagus estuary has great exposure to tsunami occurrences and, being bordered by the largest metropolitan area in the country, is a particularly worrisome location in what concerns safety of populations and economic losses due to disruption of built infrastructures. The last major earthquake and tsunami combination known to have critically affected the Tagus estuary dates back to November 1st 1755. This catastrophe critically damaged Lisbon's infrastructures, led to numerous casualties and priceless heritage losses. The urban tissue of the present city still bears visible the effects of the catastrophe and of the ensuing protection measures. The objective of this work is to simulate the propagation of debris carried by a 1755-like tsunami along the present-day bathimetric and altimetric conditions of Lisbon waterfront. Particular emphasis was directed to the modeling of vehicles since the tsunami is likely to affect areas that are major traffic nodes such as Alcântara, with more than 1500 vehicles in road network of about 3 km. The simulation tool employed is based on a 2DH spatial (eulerian) shallow-flow approach suited to complex and dynamic bottom boundaries. The discretization technique relies on a finite-volume scheme, based on a flux-splitting technique incorporating a reviewed version of the Roe Riemann solver (Canelas et al. 2013). Two formulations were employed to model the advection of debris: a fully coupled continuum approach, where solid bodies are described by the concentration only and an uncoupled material (lagrangian) formulation where solid bodies are tracked between two time-steps once the flow field is determined by the eulerian solver. In the latter case, concentrations are updated after tracking the solid bodies thus correcting the mass and momentum balance to be used for the next time-step. The urban tissue was thoroughly discretized with a mesh finer than street width so that the buildings would act as obstacles and the streets would bind the incoming flow. To simplify the plan-view geometry, it was assumed that buildings would retain its original shape after the earthquake. The results of the eulerian-continuum and of the lagrangian-discrete solutions are presented, compared and discussed. It was found that the patterns of deposition of the eulerian-continuum model can be considerably different to those obtained by the lagrangian-discrete solution if the latter assumes that vehicles have a small equivalent density and if momentum losses due to inter-particle collisions are neglected. Results become more similar if vehicles are considered much denser than water and that the mixture of water and solid bodies loses momentum due to particle collisions. Acknowledgements: Project PTDC/ECM/117660/2010, funded by the Portuguese Foundation for Science and Technology (FCT) has partially supported this work. References Canelas, R.; Murillo, J. & Ferreira, R.M.L. (2013) 2DH modelling of discontinuous flows over mobile beds. Accepted, Journal of Hydraulic Research, December 2012 Baptista M.A. Miranda, J.M. (2009). Revision of the Portuguese catalog of tsunamis. Nat. Hazards Earth Syst. Sci., 9, 25-42.

Unambiguous formalism for higher order Lagrangian field theories

NASA Astrophysics Data System (ADS)

Campos, Cédric M.; de León, Manuel; Martín de Diego, David; Vankerschaver, Joris

2009-11-01

The aim of this paper is to propose an unambiguous intrinsic formalism for higher order field theories which avoids the arbitrariness in the generalization of the conventional description of field theories, and implies the existence of different Cartan forms and Legendre transformations. We propose a differential-geometric setting for the dynamics of a higher order field theory, based on the Skinner and Rusk formalism for mechanics. This approach incorporates aspects of both the Lagrangian and the Hamiltonian description, since the field equations are formulated using the Lagrangian on a higher order jet bundle and the canonical multisymplectic form on its affine dual. As both of these objects are uniquely defined, the Skinner-Rusk approach has the advantage that it does not suffer from the arbitrariness in conventional descriptions. The result is that we obtain a unique and global intrinsic version of the Euler-Lagrange equations for higher order field theories. Several examples illustrate our construction.

Modified Mixed Lagrangian-Eulerian Method Based on Numerical Framework of MT3DMS on Cauchy Boundary.

PubMed

Suk, Heejun

2016-07-01

MT3DMS, a modular three-dimensional multispecies transport model, has long been a popular model in the groundwater field for simulating solute transport in the saturated zone. However, the method of characteristics (MOC), modified MOC (MMOC), and hybrid MOC (HMOC) included in MT3DMS did not treat Cauchy boundary conditions in a straightforward or rigorous manner, from a mathematical point of view. The MOC, MMOC, and HMOC regard the Cauchy boundary as a source condition. For the source, MOC, MMOC, and HMOC calculate the Lagrangian concentration by setting it equal to the cell concentration at an old time level. However, the above calculation is an approximate method because it does not involve backward tracking in MMOC and HMOC or allow performing forward tracking at the source cell in MOC. To circumvent this problem, a new scheme is proposed that avoids direct calculation of the Lagrangian concentration on the Cauchy boundary. The proposed method combines the numerical formulations of two different schemes, the finite element method (FEM) and the Eulerian-Lagrangian method (ELM), into one global matrix equation. This study demonstrates the limitation of all MT3DMS schemes, including MOC, MMOC, HMOC, and a third-order total-variation-diminishing (TVD) scheme under Cauchy boundary conditions. By contrast, the proposed method always shows good agreement with the exact solution, regardless of the flow conditions. Finally, the successful application of the proposed method sheds light on the possible flexibility and capability of the MT3DMS to deal with the mass transport problems of all flow regimes. © 2016, National Ground Water Association.

a Marker-Based Eulerian-Lagrangian Method for Multiphase Flow with Supersonic Combustion Applications

NASA Astrophysics Data System (ADS)

Fan, Xiaofeng; Wang, Jiangfeng

2016-06-01

The atomization of liquid fuel is a kind of intricate dynamic process from continuous phase to discrete phase. Procedures of fuel spray in supersonic flow are modeled with an Eulerian-Lagrangian computational fluid dynamics methodology. The method combines two distinct techniques and develops an integrated numerical simulation method to simulate the atomization processes. The traditional finite volume method based on stationary (Eulerian) Cartesian grid is used to resolve the flow field, and multi-component Navier-Stokes equations are adopted in present work, with accounting for the mass exchange and heat transfer occupied by vaporization process. The marker-based moving (Lagrangian) grid is utilized to depict the behavior of atomized liquid sprays injected into a gaseous environment, and discrete droplet model 13 is adopted. To verify the current approach, the proposed method is applied to simulate processes of liquid atomization in supersonic cross flow. Three classic breakup models, TAB model, wave model and K-H/R-T hybrid model, are discussed. The numerical results are compared with multiple perspectives quantitatively, including spray penetration height and droplet size distribution. In addition, the complex flow field structures induced by the presence of liquid spray are illustrated and discussed. It is validated that the maker-based Eulerian-Lagrangian method is effective and reliable.

Evaluating the performance of the particle finite element method in parallel architectures

NASA Astrophysics Data System (ADS)

Gimenez, Juan M.; Nigro, Norberto M.; Idelsohn, Sergio R.

2014-05-01

This paper presents a high performance implementation for the particle-mesh based method called particle finite element method two (PFEM-2). It consists of a material derivative based formulation of the equations with a hybrid spatial discretization which uses an Eulerian mesh and Lagrangian particles. The main aim of PFEM-2 is to solve transport equations as fast as possible keeping some level of accuracy. The method was found to be competitive with classical Eulerian alternatives for these targets, even in their range of optimal application. To evaluate the goodness of the method with large simulations, it is imperative to use of parallel environments. Parallel strategies for Finite Element Method have been widely studied and many libraries can be used to solve Eulerian stages of PFEM-2. However, Lagrangian stages, such as streamline integration, must be developed considering the parallel strategy selected. The main drawback of PFEM-2 is the large amount of memory needed, which limits its application to large problems with only one computer. Therefore, a distributed-memory implementation is urgently needed. Unlike a shared-memory approach, using domain decomposition the memory is automatically isolated, thus avoiding race conditions; however new issues appear due to data distribution over the processes. Thus, a domain decomposition strategy for both particle and mesh is adopted, which minimizes the communication between processes. Finally, performance analysis running over multicore and multinode architectures are presented. The Courant-Friedrichs-Lewy number used influences the efficiency of the parallelization and, in some cases, a weighted partitioning can be used to improve the speed-up. However the total cputime for cases presented is lower than that obtained when using classical Eulerian strategies.

Comparisons of Lagrangian and Eulerian PDF methods in simulations of non-premixed turbulent jet flames with moderate-to-strong turbulence-chemistry interactions

NASA Astrophysics Data System (ADS)

Jaishree, J.; Haworth, D. C.

2012-06-01

Transported probability density function (PDF) methods have been applied widely and effectively for modelling turbulent reacting flows. In most applications of PDF methods to date, Lagrangian particle Monte Carlo algorithms have been used to solve a modelled PDF transport equation. However, Lagrangian particle PDF methods are computationally intensive and are not readily integrated into conventional Eulerian computational fluid dynamics (CFD) codes. Eulerian field PDF methods have been proposed as an alternative. Here a systematic comparison is performed among three methods for solving the same underlying modelled composition PDF transport equation: a consistent hybrid Lagrangian particle/Eulerian mesh (LPEM) method, a stochastic Eulerian field (SEF) method and a deterministic Eulerian field method with a direct-quadrature-method-of-moments closure (a multi-environment PDF-MEPDF method). The comparisons have been made in simulations of a series of three non-premixed, piloted methane-air turbulent jet flames that exhibit progressively increasing levels of local extinction and turbulence-chemistry interactions: Sandia/TUD flames D, E and F. The three PDF methods have been implemented using the same underlying CFD solver, and results obtained using the three methods have been compared using (to the extent possible) equivalent physical models and numerical parameters. Reasonably converged mean and rms scalar profiles are obtained using 40 particles per cell for the LPEM method or 40 Eulerian fields for the SEF method. Results from these stochastic methods are compared with results obtained using two- and three-environment MEPDF methods. The relative advantages and disadvantages of each method in terms of accuracy and computational requirements are explored and identified. In general, the results obtained from the two stochastic methods (LPEM and SEF) are very similar, and are in closer agreement with experimental measurements than those obtained using the MEPDF method, while MEPDF is the most computationally efficient of the three methods. These and other findings are discussed in detail.

A three-dimensional finite-volume Eulerian-Lagrangian Localized Adjoint Method (ELLAM) for solute-transport modeling

USGS Publications Warehouse

Heberton, C.I.; Russell, T.F.; Konikow, Leonard F.; Hornberger, G.Z.

2000-01-01

This report documents the U.S. Geological Survey Eulerian-Lagrangian Localized Adjoint Method (ELLAM) algorithm that solves an integral form of the solute-transport equation, incorporating an implicit-in-time difference approximation for the dispersive and sink terms. Like the algorithm in the original version of the U.S. Geological Survey MOC3D transport model, ELLAM uses a method of characteristics approach to solve the transport equation on the basis of the velocity field. The ELLAM algorithm, however, is based on an integral formulation of conservation of mass and uses appropriate numerical techniques to obtain global conservation of mass. The implicit procedure eliminates several stability criteria required for an explicit formulation. Consequently, ELLAM allows large transport time increments to be used. ELLAM can produce qualitatively good results using a small number of transport time steps. A description of the ELLAM numerical method, the data-input requirements and output options, and the results of simulator testing and evaluation are presented. The ELLAM algorithm was evaluated for the same set of problems used to test and evaluate Version 1 and Version 2 of MOC3D. These test results indicate that ELLAM offers a viable alternative to the explicit and implicit solvers in MOC3D. Its use is desirable when mass balance is imperative or a fast, qualitative model result is needed. Although accurate solutions can be generated using ELLAM, its efficiency relative to the two previously documented solution algorithms is problem dependent.

Scalable Methods for Eulerian-Lagrangian Simulation Applied to Compressible Multiphase Flows

NASA Astrophysics Data System (ADS)

Zwick, David; Hackl, Jason; Balachandar, S.

2017-11-01

Multiphase flows can be found in countless areas of physics and engineering. Many of these flows can be classified as dispersed two-phase flows, meaning that there are solid particles dispersed in a continuous fluid phase. A common technique for simulating such flow is the Eulerian-Lagrangian method. While useful, this method can suffer from scaling issues on larger problem sizes that are typical of many realistic geometries. Here we present scalable techniques for Eulerian-Lagrangian simulations and apply it to the simulation of a particle bed subjected to expansion waves in a shock tube. The results show that the methods presented here are viable for simulation of larger problems on modern supercomputers. This material is based upon work supported by the National Science Foundation Graduate Research Fellowship under Grant No. DGE-1315138. This work was supported in part by the U.S. Department of Energy under Contract No. DE-NA0002378.

An investigation of turbulent transport in the extreme lower atmosphere

NASA Technical Reports Server (NTRS)

Koper, C. A., Jr.; Sadeh, W. Z.

1975-01-01

A model in which the Lagrangian autocorrelation is expressed by a domain integral over a set of usual Eulerian autocorrelations acquired concurrently at all points within a turbulence box is proposed along with a method for ascertaining the statistical stationarity of turbulent velocity by creating an equivalent ensemble to investigate the flow in the extreme lower atmosphere. Simultaneous measurements of turbulent velocity on a turbulence line along the wake axis were carried out utilizing a longitudinal array of five hot-wire anemometers remotely operated. The stationarity test revealed that the turbulent velocity is approximated as a realization of a weakly self-stationary random process. Based on the Lagrangian autocorrelation it is found that: (1) large diffusion time predominated; (2) ratios of Lagrangian to Eulerian time and spatial scales were smaller than unity; and, (3) short and long diffusion time scales and diffusion spatial scales were constrained within their Eulerian counterparts.

Comparing Lagrangian and Eulerian models for CO2 transport - a step towards Bayesian inverse modeling using WRF/STILT-VPRM

NASA Astrophysics Data System (ADS)

Pillai, D.; Gerbig, C.; Kretschmer, R.; Beck, V.; Karstens, U.; Neininger, B.; Heimann, M.

2012-10-01

We present simulations of atmospheric CO2 concentrations provided by two modeling systems, run at high spatial resolution: the Eulerian-based Weather Research Forecasting (WRF) model and the Lagrangian-based Stochastic Time-Inverted Lagrangian Transport (STILT) model, both of which are coupled to a diagnostic biospheric model, the Vegetation Photosynthesis and Respiration Model (VPRM). The consistency of the simulations is assessed with special attention paid to the details of horizontal as well as vertical transport and mixing of CO2 concentrations in the atmosphere. The dependence of model mismatch (Eulerian vs. Lagrangian) on models' spatial resolution is further investigated. A case study using airborne measurements during which two models showed large deviations from each other is analyzed in detail as an extreme case. Using aircraft observations and pulse release simulations, we identified differences in the representation of details in the interaction between turbulent mixing and advection through wind shear as the main cause of discrepancies between WRF and STILT transport at a spatial resolution such as 2 and 6 km. Based on observations and inter-model comparisons of atmospheric CO2 concentrations, we show that a refinement of the parameterization of turbulent velocity variance and Lagrangian time-scale in STILT is needed to achieve a better match between the Eulerian and the Lagrangian transport at such a high spatial resolution (e.g. 2 and 6 km). Nevertheless, the inter-model differences in simulated CO2 time series for a tall tower observatory at Ochsenkopf in Germany are about a factor of two smaller than the model-data mismatch and about a factor of three smaller than the mismatch between the current global model simulations and the data.

On tide-induced Lagrangian residual current and residual transport: 2. Residual transport with application in south San Francisco Bay, California

USGS Publications Warehouse

Feng, Shizuo; Cheng, Ralph T.; Pangen, Xi

1986-01-01

The transports of solutes and other tracers are fundamental to estuarine processes. The apparent transport mechanisms are convection by tidal current and current-induced shear effect dispersion for processes which take place in a time period of the order of a tidal cycle. However, as emphasis is shifted toward the effects of intertidal processes, the net transport is mainly determined by tide-induced residual circulation and by residual circulation due to other processes. The commonly used intertidal conservation equation takes the form of a convection-dispersion equation in which the convective velocity is the Eulerian residual current, and the dispersion terms are often referred to as the phase effect dispersion or, sometimes, as the “tidal dispersion.” The presence of these dispersion terms is merely the result of a Fickian type hypothesis. Since the actual processes are not Fickian, thus a Fickian hypothesis obscures the physical significance of this equation. Recent research results on residual circulation have suggested that long-term transport phenomena are closely related to the Lagrangian residual current or the Lagrangian residual transport. In this paper a new formulation of an intertidal conservation equation is presented and examined in detail. In a weakly nonlinear tidal estuary the resultant intertidal transport equation also takes the form of a convection-dispersion equation without the ad hoc introduction of phase effect dispersion in a form of dispersion tensor. The convective velocity in the resultant equation is the first-order Lagrangian residual current (the sum of the Eulerian residual current and the Stokes drift). The remaining dispersion terms are important only in higher-order solutions; they are due to shear effect dispersion and turbulent mixing. There exists a dispersion boundary layer adjacent to shoreline boundaries. An order of magnitude estimate of the properties in the dispersion boundary layer is given. The present treatment of intertidal transport processes is illustrated by an analytical solution for an amphidromic system and by a numerical application in South San Francisco Bay, California. The present formulation reveals that the mechanism for long-term transport of solutes is mainly convection due to the Lagrangian residual current in the interior of a tidal estuary. This result also points out the weakness in the tidal dispersion formulation, and explains the large variability of the observed values for tidal dispersion coefficients. Further research on properties of the dispersion boundary layer is needed.

Eulerian and Lagrangian approaches to multidimensional condensation and collection

NASA Astrophysics Data System (ADS)

Li, Xiang-Yu; Brandenburg, A.; Haugen, N. E. L.; Svensson, G.

2017-06-01

Turbulence is argued to play a crucial role in cloud droplet growth. The combined problem of turbulence and cloud droplet growth is numerically challenging. Here an Eulerian scheme based on the Smoluchowski equation is compared with two Lagrangian superparticle (or superdroplet) schemes in the presence of condensation and collection. The growth processes are studied either separately or in combination using either two-dimensional turbulence, a steady flow or just gravitational acceleration without gas flow. Good agreement between the different schemes for the time evolution of the size spectra is observed in the presence of gravity or turbulence. The Lagrangian superparticle schemes are found to be superior over the Eulerian one in terms of computational performance. However, it is shown that the use of interpolation schemes such as the cloud-in-cell algorithm is detrimental in connection with superparticle or superdroplet approaches. Furthermore, the use of symmetric over asymmetric collection schemes is shown to reduce the amount of scatter in the results. For the Eulerian scheme, gravitational collection is rather sensitive to the mass bin resolution, but not so in the case with turbulence.