EUPDF: An Eulerian-Based Monte Carlo Probability Density Function (PDF) Solver. User's Manual
NASA Technical Reports Server (NTRS)
Raju, M. S.
1998-01-01
EUPDF is an Eulerian-based Monte Carlo PDF solver developed for application with sprays, combustion, parallel computing and unstructured grids. It is designed to be massively parallel and could easily be coupled with any existing gas-phase flow and spray solvers. The solver accommodates the use of an unstructured mesh with mixed elements of either triangular, quadrilateral, and/or tetrahedral type. The manual provides the user with the coding required to couple the PDF code to any given flow code and a basic understanding of the EUPDF code structure as well as the models involved in the PDF formulation. The source code of EUPDF will be available with the release of the National Combustion Code (NCC) as a complete package.
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.
NASA Astrophysics Data System (ADS)
Druzhinin, O. A.; Elghobashi, S. E.
1999-09-01
In a recent study we showed that the two-fluid (TF) formulation can be used in the direct numerical simulation (DNS) of bubble- (or particle-) laden decaying isotropic turbulence with considerable saving in CPU-time and memory as compared to the trajectory approach employed by many researchers. In the present paper, we develop a Lagrangian-Eulerian mapping (LEM) solver for DNS of bubble-laden turbulent shear flows using TF. The purpose of LEM is to resolve the large gradients of bubble velocity and concentration which result from the absence of the diffusion terms in the equations of bubble-phase motion and the preferential accumulation of bubbles. A standard finite-difference scheme (FDS) fails to resolve these gradients. We examine the performance of the new method in DNS of a bubble-laden Taylor-Green vortex, spatially developing plane mixing layer, and homogeneous shear turbulent flow.
NASA Astrophysics Data System (ADS)
Pathak, Ashish; Raessi, Mehdi
2016-04-01
We present a three-dimensional (3D) and fully Eulerian approach to capturing the interaction between two fluids and moving rigid structures by using the fictitious domain and volume-of-fluid (VOF) methods. The solid bodies can have arbitrarily complex geometry and can pierce the fluid-fluid interface, forming contact lines. The three-phase interfaces are resolved and reconstructed by using a VOF-based methodology. Then, a consistent scheme is employed for transporting mass and momentum, allowing for simulations of three-phase flows of large density ratios. The Eulerian approach significantly simplifies numerical resolution of the kinematics of rigid bodies of complex geometry and with six degrees of freedom. The fluid-structure interaction (FSI) is computed using the fictitious domain method. The methodology was developed in a message passing interface (MPI) parallel framework accelerated with graphics processing units (GPUs). The computationally intensive solution of the pressure Poisson equation is ported to GPUs, while the remaining calculations are performed on CPUs. The performance and accuracy of the methodology are assessed using an array of test cases, focusing individually on the flow solver and the FSI in surface-piercing configurations. Finally, an application of the proposed methodology in simulations of the ocean wave energy converters is presented.
NASA Technical Reports Server (NTRS)
Raju, M. S.
1998-01-01
The success of any solution methodology used in the study of gas-turbine combustor flows depends a great deal on how well it can model the various complex and rate controlling processes associated with the spray's turbulent transport, mixing, chemical kinetics, evaporation, and spreading rates, as well as convective and radiative heat transfer and other phenomena. The phenomena to be modeled, which are controlled by these processes, often strongly interact with each other at different times and locations. In particular, turbulence plays an important role in determining the rates of mass and heat transfer, chemical reactions, and evaporation in many practical combustion devices. The influence of turbulence in a diffusion flame manifests itself in several forms, ranging from the so-called wrinkled, or stretched, flamelets regime to the distributed combustion regime, depending upon how turbulence interacts with various flame scales. Conventional turbulence models have difficulty treating highly nonlinear reaction rates. A solution procedure based on the composition joint probability density function (PDF) approach holds the promise of modeling various important combustion phenomena relevant to practical combustion devices (such as extinction, blowoff limits, and emissions predictions) because it can account for nonlinear chemical reaction rates without making approximations. In an attempt to advance the state-of-the-art in multidimensional numerical methods, we at the NASA Lewis Research Center extended our previous work on the PDF method to unstructured grids, parallel computing, and sprays. EUPDF, which was developed by M.S. Raju of Nyma, Inc., was designed to be massively parallel and could easily be coupled with any existing gas-phase and/or spray solvers. EUPDF can use an unstructured mesh with mixed triangular, quadrilateral, and/or tetrahedral elements. The application of the PDF method showed favorable results when applied to several supersonic
Eulerian action principles for linearized reduced dynamical equations
NASA Astrophysics Data System (ADS)
Brizard, Alain
1994-08-01
New Eulerian action principles for the linearized gyrokinetic Maxwell-Vlasov equations and the linearized kinetic-magnetohydrodynamic (kinetic-MHD) equations are presented. The variational fields for the linearized gyrokinetic Vlasov-Maxwell equations are the perturbed electromagnetic potentials (φ1,A1) and the gyroangle-independent gyrocenter (gy) function Sgy, while the variational fields for the linearized kinetic-MHD equations are the ideal MHD fluid displacement ξ and the gyroangle-independent drift-kinetic (dk) function Sdk (defined as the drift-kinetic limit of Sgy). According to the Lie-transform approach to Vlasov perturbation theory, Sgy generates first-order perturbations in the gyrocenter distribution F1≡{Sgy, F0}gc, where F1 satisfies the linearized gyrokinetic Vlasov equation and {, }gc denotes the unperturbed guiding-center (gc) Poisson bracket. Previous quadratic variational forms were constructed ad hoc from the linearized equations, and required the linearized gyrokinetic (or drift-kinetic) Vlasov equation to be solved a priori (e.g., by integration along an unperturbed guiding-center orbit) through the use of the normal-mode and ballooning-mode representations. The presented action principles ignore these requirements and, thus, apply to more general perturbations.
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.
Objective Eulerian coherent structures
NASA Astrophysics Data System (ADS)
Serra, Mattia; Haller, George
2016-05-01
We define objective Eulerian Coherent Structures (OECSs) in two-dimensional, non-autonomous dynamical systems as the instantaneously most influential material curves. Specifically, OECSs are stationary curves of the averaged instantaneous material stretching-rate or material shearing-rate functionals. From these objective (frame-invariant) variational principles, we obtain explicit differential equations for hyperbolic, elliptic, and parabolic OECSs. As an illustration, we compute OECSs in an unsteady ocean velocity data set. In comparison to structures suggested by other common Eulerian diagnostic tools, we find OECSs to be the correct short-term cores of observed trajectory deformation patterns.
Objective Eulerian coherent structures.
Serra, Mattia; Haller, George
2016-05-01
We define objective Eulerian Coherent Structures (OECSs) in two-dimensional, non-autonomous dynamical systems as the instantaneously most influential material curves. Specifically, OECSs are stationary curves of the averaged instantaneous material stretching-rate or material shearing-rate functionals. From these objective (frame-invariant) variational principles, we obtain explicit differential equations for hyperbolic, elliptic, and parabolic OECSs. As an illustration, we compute OECSs in an unsteady ocean velocity data set. In comparison to structures suggested by other common Eulerian diagnostic tools, we find OECSs to be the correct short-term cores of observed trajectory deformation patterns. PMID:27249950
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.
A multigrid fluid pressure solver handling separating solid boundary conditions.
Chentanez, Nuttapong; Müller-Fischer, Matthias
2012-08-01
We present a multigrid method for solving the linear complementarity problem (LCP) resulting from discretizing the Poisson equation subject to separating solid boundary conditions in an Eulerian liquid simulation’s pressure projection step. The method requires only a few small changes to a multigrid solver for linear systems. Our generalized solver is fast enough to handle 3D liquid simulations with separating boundary conditions in practical domain sizes. Previous methods could only handle relatively small 2D domains in reasonable time, because they used expensive quadratic programming (QP) solvers. We demonstrate our technique in several practical scenarios, including nonaxis-aligned containers and moving solids in which the omission of separating boundary conditions results in disturbing artifacts of liquid sticking to solids. Our measurements show, that the convergence rate of our LCP solver is close to that of a standard multigrid solver. PMID:22411885
NASA Astrophysics Data System (ADS)
Parmentier, Philippe; Winckelmans, Gregoire; Chatelain, Philippe; Hillewaert, Koen
2015-11-01
A hybrid approach, coupling a compressible vortex particle-mesh method (CVPM, also with efficient Poisson solver) and a high order compressible discontinuous Galerkin Eulerian solver, is being developed in order to efficiently simulate flows past bodies; also in the transonic regime. The Eulerian solver is dedicated to capturing the anisotropic flow structures in the near-wall region whereas the CVPM solver is exploited away from the body and in the wake. An overlapping domain decomposition approach is used. The Eulerian solver, which captures the near-body region, also corrects the CVPM solution in that region at every time step. The CVPM solver, which captures the region away from the body and the wake, also provides the outer boundary conditions to the Eulerian solver. Because of the coupling, a boundary element method is also required for consistency. The approach is assessed on typical 2D benchmark cases. Supported by the Fund for Research Training in Industry and Agriculture (F.R.I.A.).
Stanley, Vendall S.; Heroux, Michael A.; Hoekstra, Robert J.; Sala, Marzio
2004-03-01
Amesos is the Direct Sparse Solver Package in Trilinos. The goal of Amesos is to make AX=S as easy as it sounds, at least for direct methods. Amesos provides interfaces to a number of third party sparse direct solvers, including SuperLU, SuperLU MPI, DSCPACK, UMFPACK and KLU. Amesos provides a common object oriented interface to the best sparse direct solvers in the world. A sparse direct solver solves for x in Ax = b. where A is a matrix and x and b are vectors (or multi-vectors). A sparse direct solver flrst factors A into trinagular matrices L and U such that A = LU via gaussian elimination and then solves LU x = b. Switching amongst solvers in Amesos roquires a change to a single parameter. Yet, no solver needs to be linked it, unless it is used. All conversions between the matrices provided by the user and the format required by the underlying solver is performed by Amesos. As new sparse direct solvers are created, they will be incorporated into Amesos, allowing the user to simpty link with the new solver, change a single parameter in the calling sequence, and use the new solver. Amesos allows users to specify whether the matrix has changed. Amesos can be used anywhere that any sparse direct solver is needed.
Energy Science and Technology Software Center (ESTSC)
2004-03-01
Amesos is the Direct Sparse Solver Package in Trilinos. The goal of Amesos is to make AX=S as easy as it sounds, at least for direct methods. Amesos provides interfaces to a number of third party sparse direct solvers, including SuperLU, SuperLU MPI, DSCPACK, UMFPACK and KLU. Amesos provides a common object oriented interface to the best sparse direct solvers in the world. A sparse direct solver solves for x in Ax = b. wheremore » A is a matrix and x and b are vectors (or multi-vectors). A sparse direct solver flrst factors A into trinagular matrices L and U such that A = LU via gaussian elimination and then solves LU x = b. Switching amongst solvers in Amesos roquires a change to a single parameter. Yet, no solver needs to be linked it, unless it is used. All conversions between the matrices provided by the user and the format required by the underlying solver is performed by Amesos. As new sparse direct solvers are created, they will be incorporated into Amesos, allowing the user to simpty link with the new solver, change a single parameter in the calling sequence, and use the new solver. Amesos allows users to specify whether the matrix has changed. Amesos can be used anywhere that any sparse direct solver is needed.« less
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
Comparison of Joint Modeling Approaches Including Eulerian Sliding Interfaces
Lomov, I; Antoun, T; Vorobiev, O
2009-12-16
Accurate representation of discontinuities such as joints and faults is a key ingredient for high fidelity modeling of shock propagation in geologic media. The following study was done to improve treatment of discontinuities (joints) in the Eulerian hydrocode GEODYN (Lomov and Liu 2005). Lagrangian methods with conforming meshes and explicit inclusion of joints in the geologic model are well suited for such an analysis. Unfortunately, current meshing tools are unable to automatically generate adequate hexahedral meshes for large numbers of irregular polyhedra. Another concern is that joint stiffness in such explicit computations requires significantly reduced time steps, with negative implications for both the efficiency and quality of the numerical solution. An alternative approach is to use non-conforming meshes and embed joint information into regular computational elements. However, once slip displacement on the joints become comparable to the zone size, Lagrangian (even non-conforming) meshes could suffer from tangling and decreased time step problems. The use of non-conforming meshes in an Eulerian solver may alleviate these difficulties and provide a viable numerical approach for modeling the effects of faults on the dynamic response of geologic materials. We studied shock propagation in jointed/faulted media using a Lagrangian and two Eulerian approaches. To investigate the accuracy of this joint treatment the GEODYN calculations have been compared with results from the Lagrangian code GEODYN-L which uses an explicit treatment of joints via common plane contact. We explore two approaches to joint treatment in the code, one for joints with finite thickness and the other for tight joints. In all cases the sliding interfaces are tracked explicitly without homogenization or blending the joint and block response into an average response. In general, rock joints will introduce an increase in normal compliance in addition to a reduction in shear strength. In the
Nonlinear Eulerian thermoelasticity for anisotropic crystals
NASA Astrophysics Data System (ADS)
Clayton, J. D.
2013-10-01
A complete continuum thermoelastic theory for large deformation of crystals of arbitrary symmetry is developed. The theory incorporates as a fundamental state variable in the thermodynamic potentials what is termed an Eulerian strain tensor (in material coordinates) constructed from the inverse of the deformation gradient. Thermodynamic identities and relationships among Eulerian and the usual Lagrangian material coefficients are derived, significantly extending previous literature that focused on materials with cubic or hexagonal symmetry and hydrostatic loading conditions. Analytical solutions for homogeneous deformations of ideal cubic crystals are studied over a prescribed range of elastic coefficients; stress states and intrinsic stability measures are compared. For realistic coefficients, Eulerian theory is shown to predict more physically realistic behavior than Lagrangian theory under large compression and shear. Analytical solutions for shock compression of anisotropic single crystals are derived for internal energy functions quartic in Lagrangian or Eulerian strain and linear in entropy; results are analyzed for quartz, sapphire, and diamond. When elastic constants of up to order four are included, both Lagrangian and Eulerian theories are capable of matching Hugoniot data. When only the second-order elastic constant is known, an alternative theory incorporating a mixed Eulerian-Lagrangian strain tensor provides a reasonable approximation of experimental data.
An Eulerian/Lagrangian method for computing blade/vortex impingement
NASA Technical Reports Server (NTRS)
Steinhoff, John; Senge, Heinrich; Yonghu, Wenren
1991-01-01
A combined Eulerian/Lagrangian approach to calculating helicopter rotor flows with concentrated vortices is described. The method computes a general evolving vorticity distribution without any significant numerical diffusion. Concentrated vortices can be accurately propagated over long distances on relatively coarse grids with cores only several grid cells wide. The method is demonstrated for a blade/vortex impingement case in 2D and 3D where a vortex is cut by a rotor blade, and the results are compared to previous 2D calculations involving a fifth-order Navier-Stokes solver on a finer grid.
CASTRO: A NEW COMPRESSIBLE ASTROPHYSICAL SOLVER. II. GRAY RADIATION HYDRODYNAMICS
Zhang, W.; Almgren, A.; Bell, J.; Howell, L.; Burrows, A.
2011-10-01
We describe the development of a flux-limited gray radiation solver for the compressible astrophysics code, CASTRO. CASTRO uses an Eulerian grid with block-structured adaptive mesh refinement based on a nested hierarchy of logically rectangular variable-sized grids with simultaneous refinement in both space and time. The gray radiation solver is based on a mixed-frame formulation of radiation hydrodynamics. In our approach, the system is split into two parts, one part that couples the radiation and fluid in a hyperbolic subsystem, and another parabolic part that evolves radiation diffusion and source-sink terms. The hyperbolic subsystem is solved explicitly with a high-order Godunov scheme, whereas the parabolic part is solved implicitly with a first-order backward Euler method.
NASA Astrophysics Data System (ADS)
Fox, R. O.; Laurent, F.; Massot, M.
2008-03-01
The scope of the present study is Eulerian modeling and simulation of polydisperse liquid sprays undergoing droplet coalescence and evaporation. The fundamental mathematical description is the Williams spray equation governing the joint number density function f(v,u;x,t) of droplet volume and velocity. Eulerian multi-fluid models have already been rigorously derived from this equation in Laurent et al. [F. Laurent, M. Massot, P. Villedieu, Eulerian multi-fluid modeling for the numerical simulation of coalescence in polydisperse dense liquid sprays, J. Comput. Phys. 194 (2004) 505-543]. The first key feature of the paper is the application of direct quadrature method of moments (DQMOM) introduced by Marchisio and Fox [D.L. Marchisio, R.O. Fox, Solution of population balance equations using the direct quadrature method of moments, J. Aerosol Sci. 36 (2005) 43-73] to the Williams spray equation. Both the multi-fluid method and DQMOM yield systems of Eulerian conservation equations with complicated interaction terms representing coalescence. In order to focus on the difficulties associated with treating size-dependent coalescence and to avoid numerical uncertainty issues associated with two-way coupling, only one-way coupling between the droplets and a given gas velocity field is considered. In order to validate and compare these approaches, the chosen configuration is a self-similar 2D axisymmetrical decelerating nozzle with sprays having various size distributions, ranging from smooth ones up to Dirac delta functions. The second key feature of the paper is a thorough comparison of the two approaches for various test-cases to a reference solution obtained through a classical stochastic Lagrangian solver. Both Eulerian models prove to describe adequately spray coalescence and yield a very interesting alternative to the Lagrangian solver. The third key point of the study is a detailed description of the limitations associated with each method, thus giving criteria for
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)
Performance Analysis of GYRO: A Tool Evaluation
Worley, P.; Roth, P.; Candy, J.; Shan, Hongzhang; Mahinthakumar,G.; Sreepathi, S.; Carrington, L.; Kaiser, T.; Snavely, A.; Reed, D.; Zhang, Y.; Huck, K.; Malony, A.; Shende, S.; Moore, S.; Wolf, F.
2005-06-26
The performance of the Eulerian gyrokinetic-Maxwell solver code GYRO is analyzed on five high performance computing systems. First, a manual approach is taken, using custom scripts to analyze the output of embedded wall clock timers, floating point operation counts collected using hardware performance counters, and traces of user and communication events collected using the profiling interface to Message Passing Interface (MPI) libraries. Parts of the analysis are then repeated or extended using a number of sophisticated performance analysis tools: IPM, KOJAK, SvPablo, TAU, and the PMaC modeling tool suite. The paper briefly discusses what has been discovered via this manual analysis process, what performance analyses are inconvenient or infeasible to attempt manually, and to what extent the tools show promise in accelerating or significantly extending the manual performance analyses.
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)
Parallel Multigrid Equation Solver
Energy Science and Technology Software Center (ESTSC)
2001-09-07
Prometheus is a fully parallel multigrid equation solver for matrices that arise in unstructured grid finite element applications. It includes a geometric and an algebraic multigrid method and has solved problems of up to 76 mullion degrees of feedom, problems in linear elasticity on the ASCI blue pacific and ASCI red machines.
Eulerian network modeling of longitudinal dispersion
NASA Astrophysics Data System (ADS)
Mehmani, Yashar; Balhoff, Matthew T.
2015-10-01
A novel Eulerian network model that incorporates "shear dispersion," the stretching of solute due to nonuniform velocity profiles within pore throats is developed. The superposing transport method (STM) is nonlocal in time (i.e., uses information from several previous time steps) and is equivalent to performing network-wide time convolutions of elementary throat response functions. Predicted macroscopic longitudinal dispersion coefficients for disordered sphere packs are in good agreement with published experimental data. We further investigate the impact of mixing assumptions within pores on macroscopic longitudinal dispersion and find the dependence to be weak for disordered sphere packs. Limitations of Eulerian network models as a whole are also discussed, and their inappropriateness for ordered porous media concluded.
Kotulski, Joseph D.; Womble, David E.; Greenberg, David; Driessen, Brian
2004-03-01
PLIRIS is an object-oriented solver built on top of a previous matrix solver used in a number of application codes. Puns solves a linear system directly via LU factorization with partial pivoting. The user provides the linear system in terms of Epetra Objects including a matrix and right-hand-sides. The user can then factor the matrix and perform the forward and back solve at a later time or solve for multiple right-hand-sides at once. This package is used when dense matrices are obtained in the problem formulation. These dense matrices occur whenever boundary element techniques are chosen for the solution procedure. This has been used in electromagnetics for both static and frequency domain problems.
Energy Science and Technology Software Center (ESTSC)
2004-03-01
PLIRIS is an object-oriented solver built on top of a previous matrix solver used in a number of application codes. Puns solves a linear system directly via LU factorization with partial pivoting. The user provides the linear system in terms of Epetra Objects including a matrix and right-hand-sides. The user can then factor the matrix and perform the forward and back solve at a later time or solve for multiple right-hand-sides at once. This packagemore » is used when dense matrices are obtained in the problem formulation. These dense matrices occur whenever boundary element techniques are chosen for the solution procedure. This has been used in electromagnetics for both static and frequency domain problems.« less
Arbitrary Lagrangian Eulerian Adaptive Mesh Refinement
Koniges, A.; Eder, D.; Masters, N.; Fisher, A.; Anderson, R.; Gunney, B.; Wang, P.; Benson, D.; Dixit, P.
2009-09-29
This is a simulation code involving an ALE (arbitrary Lagrangian-Eulerian) hydrocode with AMR (adaptive mesh refinement) and pluggable physics packages for material strength, heat conduction, radiation diffusion, and laser ray tracing developed a LLNL, UCSD, and Berkeley Lab. The code is an extension of the open source SAMRAI (Structured Adaptive Mesh Refinement Application Interface) code/library. The code can be used in laser facilities such as the National Ignition Facility. The code is alsi being applied to slurry flow (landslides).
NASA Technical Reports Server (NTRS)
Dang, A. L.; Navaz, H. K.; Rangel, R. H.
1992-01-01
The liquid thrust chambers performance (LTCP) code is used for parametric studies of flow and combustion in liquid rocket engines. Multiphase flow equations are solved in an Eulerian-Eulerian framework, and multistep finite rate chemistry is incorporated. The discretization scheme is fully implicit and is based on the total variation diminishing (TVD) scheme, which is accurate, robust, very efficient and capable of handling steep gradients and stiff chemistry. Effects of injection velocity and chamber size have been considered, and the effect of group combustion on the evaporation rate has been studied for a dense spray.
Energy Science and Technology Software Center (ESTSC)
2007-03-01
HPCCG is a simple PDE application and preconditioned conjugate gradient solver that solves a linear system on a beam-shaped domain. Although it does not address many performance issues present in real engineering applications, such as load imbalance and preconditioner scalability, it can serve as a first "sanity test" of new processor design choices, inter-connect network design choices and the scalability of a new computer system. Because it is self-contained, easy to compile and easily scaledmore » to 100s or 1000s of porcessors, it can be an attractive study code for computer system designers.« less
Arbitrary Lagrangian Eulerian Adaptive Mesh Refinement
Energy Science and Technology Software Center (ESTSC)
2009-09-29
This is a simulation code involving an ALE (arbitrary Lagrangian-Eulerian) hydrocode with AMR (adaptive mesh refinement) and pluggable physics packages for material strength, heat conduction, radiation diffusion, and laser ray tracing developed a LLNL, UCSD, and Berkeley Lab. The code is an extension of the open source SAMRAI (Structured Adaptive Mesh Refinement Application Interface) code/library. The code can be used in laser facilities such as the National Ignition Facility. The code is alsi being appliedmore » to slurry flow (landslides).« less
Numerical methods for the weakly compressible Generalized Langevin Model in Eulerian reference frame
NASA Astrophysics Data System (ADS)
Azarnykh, Dmitrii; Litvinov, Sergey; Adams, Nikolaus A.
2016-06-01
A well established approach for the computation of turbulent flow without resolving all turbulent flow scales is to solve a filtered or averaged set of equations, and to model non-resolved scales by closures derived from transported probability density functions (PDF) for velocity fluctuations. Effective numerical methods for PDF transport employ the equivalence between the Fokker-Planck equation for the PDF and a Generalized Langevin Model (GLM), and compute the PDF by transporting a set of sampling particles by GLM (Pope (1985) [1]). The natural representation of GLM is a system of stochastic differential equations in a Lagrangian reference frame, typically solved by particle methods. A representation in a Eulerian reference frame, however, has the potential to significantly reduce computational effort and to allow for the seamless integration into a Eulerian-frame numerical flow solver. GLM in a Eulerian frame (GLMEF) formally corresponds to the nonlinear fluctuating hydrodynamic equations derived by Nakamura and Yoshimori (2009) [12]. Unlike the more common Landau-Lifshitz Navier-Stokes (LLNS) equations these equations are derived from the underdamped Langevin equation and are not based on a local equilibrium assumption. Similarly to LLNS equations the numerical solution of GLMEF requires special considerations. In this paper we investigate different numerical approaches to solving GLMEF with respect to the correct representation of stochastic properties of the solution. We find that a discretely conservative staggered finite-difference scheme, adapted from a scheme originally proposed for turbulent incompressible flow, in conjunction with a strongly stable (for non-stochastic PDE) Runge-Kutta method performs better for GLMEF than schemes adopted from those proposed previously for the LLNS. We show that equilibrium stochastic fluctuations are correctly reproduced.
Parallel tridiagonal equation solvers
NASA Technical Reports Server (NTRS)
Stone, H. S.
1974-01-01
Three parallel algorithms were compared for the direct solution of tridiagonal linear systems of equations. The algorithms are suitable for computers such as ILLIAC 4 and CDC STAR. For array computers similar to ILLIAC 4, cyclic odd-even reduction has the least operation count for highly structured sets of equations, and recursive doubling has the least count for relatively unstructured sets of equations. Since the difference in operation counts for these two algorithms is not substantial, their relative running times may be more related to overhead operations, which are not measured in this paper. The third algorithm, based on Buneman's Poisson solver, has more arithmetic operations than the others, and appears to be the least favorable. For pipeline computers similar to CDC STAR, cyclic odd-even reduction appears to be the most preferable algorithm for all cases.
Amesos2 Templated Direct Sparse Solver Package
Energy Science and Technology Software Center (ESTSC)
2011-05-24
Amesos2 is a templated direct sparse solver package. Amesos2 provides interfaces to direct sparse solvers, rather than providing native solver capabilities. Amesos2 is a derivative work of the Trilinos package Amesos.
Eulerian derivation of the Coriolis force
NASA Astrophysics Data System (ADS)
Kageyama, Akira; Hyodo, Mamoru
2006-02-01
In textbooks of geophysical fluid dynamics, the Coriolis force and the centrifugal force in a rotating fluid system are derived by making use of the fluid parcel concept. In contrast to this intuitive derivation to the apparent forces, more rigorous derivation would be useful not only for the pedagogical purpose but also for the applications to other kinds of rotating geophysical systems rather than the fluid. The purpose of this paper is to show a general procedure to derive the transformed equations in the rotating frame of reference based on the local Galilean transformation and rotational coordinate transformation of field quantities. The generality and usefulness of this Eulerian approach is demonstrated in the derivation of apparent forces in rotating fluids as well as the transformed electromagnetic field equation in the rotating system.
Fragmentation statistics from Eulerian hydrocode calculations
Trucano, T.G.; Grady, D.E.; McGlaun, J.M.
1989-01-01
In this paper, we discuss a procedure for computing discrete fragmentation information for terminal ballistics events from the continuum data that emerges from two-dimensional Eulerian hydrocode simulations of these events. The present examples deal with the normal impact of lead cylinders on lead plates at velocities below 1600 m/s. In this regime, the resulting debris is almost exclusively solid lead fragments. We have experimental data of sufficient accuracy to evaluate the extraction of such fragmentation information from code simulations. The problem is interesting because the observed distribution of fragment size would require extreme subgrid resolution in the hydrocode if the statistics were computed directly. Our approach is different. We ask whether or not the general continuum description predicted by the code contains enough information to allow coupling of an analytic fragmentation theory that successfully predicts fragmentation statistics. We believe that this approach is valid, and discuss our current success in matching experimental data. 10 refs., 7 figs., 5 tabs.
NASA Technical Reports Server (NTRS)
Ilin, Andrew V.
2006-01-01
The Magnetic Field Solver computer program calculates the magnetic field generated by a group of collinear, cylindrical axisymmetric electromagnet coils. Given the current flowing in, and the number of turns, axial position, and axial and radial dimensions of each coil, the program calculates matrix coefficients for a finite-difference system of equations that approximates a two-dimensional partial differential equation for the magnetic potential contributed by the coil. The program iteratively solves these finite-difference equations by use of the modified incomplete Cholesky preconditioned-conjugate-gradient method. The total magnetic potential as a function of axial (z) and radial (r) position is then calculated as a sum of the magnetic potentials of the individual coils, using a high-accuracy interpolation scheme. Then the r and z components of the magnetic field as functions of r and z are calculated from the total magnetic potential by use of a high-accuracy finite-difference scheme. Notably, for the finite-difference calculations, the program generates nonuniform two-dimensional computational meshes from nonuniform one-dimensional meshes. Each mesh is generated in such a way as to minimize the numerical error for a benchmark one-dimensional magnetostatic problem.
2D Resistive Magnetohydrodynamics Calculations with an Arbitrary Lagrange Eulerian Code
NASA Astrophysics Data System (ADS)
Rousculp, C. L.; Gianakon, T. A.; Lipnikov, K. N.; Nelson, E. M.
2015-11-01
Single fluid resistive MHD is useful for modeling Z-pinch configurations in cylindrical geometry. One such example is thin walled liners for shock physics or HEDP experiments driven by capacitor banks such as the LANL's PHELIX or Sandia-Z. MHD is also useful for modeling high-explosive-driven flux compression generators (FCGs) and their high-current switches. The resistive MHD in our arbitrary Lagrange Eulerian (ALE) code operates in one and two dimensions in both Cartesian and cylindrical geometry. It is implemented as a time-step split operator, which consists of, ideal MHD connected to the explicit hydro momentum and energy equations and a second order mimetic discretization solver for implicit solution of the magnetic diffusion equation. In a staggered grid scheme, a single-component of cell-centered magnetic flux is conserved in the Lagrangian frame exactly, while magnetic forces are accumulated at the nodes. Total energy is conserved to round off. Total flux is conserved under the ALE relaxation and remap. The diffusion solver consistently computes Ohmic heating. Both Neumann and Dirichlet boundary conditions are available with coupling to external circuit models. Example calculations will be shown.
Sherlock Holmes, Master Problem Solver.
ERIC Educational Resources Information Center
Ballew, Hunter
1994-01-01
Shows the connections between Sherlock Holmes's investigative methods and mathematical problem solving, including observations, characteristics of the problem solver, importance of data, questioning the obvious, learning from experience, learning from errors, and indirect proof. (MKR)
Eulerian Formulation of Spatially Constrained Elastic Rods
NASA Astrophysics Data System (ADS)
Huynen, Alexandre
Slender elastic rods are ubiquitous in nature and technology. For a vast majority of applications, the rod deflection is restricted by an external constraint and a significant part of the elastic body is in contact with a stiff constraining surface. The research work presented in this doctoral dissertation formulates a computational model for the solution of elastic rods constrained inside or around frictionless tube-like surfaces. The segmentation strategy adopted to cope with this complex class of problems consists in sequencing the global problem into, comparatively simpler, elementary problems either in continuous contact with the constraint or contact-free between their extremities. Within the conventional Lagrangian formulation of elastic rods, this approach is however associated with two major drawbacks. First, the boundary conditions specifying the locations of the rod centerline at both extremities of each elementary problem lead to the establishment of isoperimetric constraints, i.e., integral constraints on the unknown length of the rod. Second, the assessment of the unilateral contact condition requires, in principle, the comparison of two curves parametrized by distinct curvilinear coordinates, viz. the rod centerline and the constraint axis. Both conspire to burden the computations associated with the method. To streamline the solution along the elementary problems and rationalize the assessment of the unilateral contact condition, the rod governing equations are reformulated within the Eulerian framework of the constraint. The methodical exploration of both types of elementary problems leads to specific formulations of the rod governing equations that stress the profound connection between the mechanics of the rod and the geometry of the constraint surface. The proposed Eulerian reformulation, which restates the rod local equilibrium in terms of the curvilinear coordinate associated with the constraint axis, describes the rod deformed configuration
CASTRO: A NEW COMPRESSIBLE ASTROPHYSICAL SOLVER. III. MULTIGROUP RADIATION HYDRODYNAMICS
Zhang, W.; Almgren, A.; Bell, J.; Howell, L.; Burrows, A.; Dolence, J.
2013-01-15
We present a formulation for multigroup radiation hydrodynamics that is correct to order O(v/c) using the comoving-frame approach and the flux-limited diffusion approximation. We describe a numerical algorithm for solving the system, implemented in the compressible astrophysics code, CASTRO. CASTRO uses a Eulerian grid with block-structured adaptive mesh refinement based on a nested hierarchy of logically rectangular variable-sized grids with simultaneous refinement in both space and time. In our multigroup radiation solver, the system is split into three parts: one part that couples the radiation and fluid in a hyperbolic subsystem, another part that advects the radiation in frequency space, and a parabolic part that evolves radiation diffusion and source-sink terms. The hyperbolic subsystem and the frequency space advection are solved explicitly with high-order Godunov schemes, whereas the parabolic part is solved implicitly with a first-order backward Euler method. Our multigroup radiation solver works for both neutrino and photon radiation.
CASTRO: A New Compressible Astrophysical Solver. III. Multigroup Radiation Hydrodynamics
NASA Astrophysics Data System (ADS)
Zhang, W.; Howell, L.; Almgren, A.; Burrows, A.; Dolence, J.; Bell, J.
2013-01-01
We present a formulation for multigroup radiation hydrodynamics that is correct to order O(v/c) using the comoving-frame approach and the flux-limited diffusion approximation. We describe a numerical algorithm for solving the system, implemented in the compressible astrophysics code, CASTRO. CASTRO uses a Eulerian grid with block-structured adaptive mesh refinement based on a nested hierarchy of logically rectangular variable-sized grids with simultaneous refinement in both space and time. In our multigroup radiation solver, the system is split into three parts: one part that couples the radiation and fluid in a hyperbolic subsystem, another part that advects the radiation in frequency space, and a parabolic part that evolves radiation diffusion and source-sink terms. The hyperbolic subsystem and the frequency space advection are solved explicitly with high-order Godunov schemes, whereas the parabolic part is solved implicitly with a first-order backward Euler method. Our multigroup radiation solver works for both neutrino and photon radiation.
SPH and Eulerian underwater bubble collapse simulations
Swegle, J.W.; Kipp, M.E.
1998-05-01
SPH (Smoothed Particle Hydrodynamics) is a gridless Lagrangian technique which is appealing as a possible alternative to numerical techniques currently used to analyze high deformation impulsive loading events. Previously, the SPH algorithm has been subjected to detailed testing and analysis to determine the feasibility of using the coupled finite-element/SPH code PRONTO/SPH for the analysis of various types of underwater explosion problems involving fluid-structure and shock-structure interactions. Here, SPH and Eulerian simulations are used to study the details of underwater bubble collapse, particularly the formation of re-entrant jets during collapse, and the loads generated on nearby structures by the jet and the complete collapse of the bubble. Jet formation is shown to be due simply to the asymmetry caused by nearby structures which disrupt the symmetry of the collapse. However, the load generated by the jet is a minor precursor to the major loads which occur at the time of complete collapse of the bubble.
AUTOMATIC DIFFERENTIATION OF AN EULERIAN HYDROCODE
R. HENNINGER; A. CARLE; P. MAUDLIN
2000-11-01
Automatic differentiation (AD) is applied to a two-dimensional Eulerian hydrodynamics computer code (hydrocode) to provide gradients that will be used for design optimization and uncertainty analysis. We examine AD in both the forward and adjoint (reverse) mode using Automatic Differentiation of Fortran (ADIFOR, version 3.0). Setup time, accuracy, and run times are described for three problems. The test set consists of a one-dimensional shock-propagation problem, a two-dimensional metal-jet-formation problem and a two-dimensional shell-collapse problem. Setup time for ADIFOR was approximately one month starting from a simplified, fixed-dimension version of the original code. ADIFOR produced accurate (as compared to finite difference) gradients in both modes for all of the problems. These test problems had 17 independent variables. We find that the forward mode is up to 39% slower and the adjoint mode is at least 11% faster than finding the gradient by means of finite differences. Problems of real interest will certainly have more independent variables. The adjoint mode is thus favored since the computational time increases only slightly for additional independent variables.
Eulerian and Lagrangian statistics in fully developed rotating turbulent flows.
NASA Astrophysics Data System (ADS)
Biferale, Luca; Bonaccorso, Fabio; Mazzitelli, Irene; Lanotte, Alessandra; Perlekar, Prasad; Musacchio, Stefano; Hinsberg, Michel; Toschi, Federico
2015-11-01
We present results concerning both Eulerian and Lagrangian statistics for turbulent under rotation at small and large Rossby numbers. Concerning the Eulerian statistics we discuss the effects of the presence of strong coherent large-scale vortical structures on the small-scale statistics. Concerning Lagrangian properties, we discuss the effects of preferential sampling at changing the inertial properties of the particles also due to the centrifugal and Coriolis forces. Supported by the ERC AdG NewTURB num. 339032.
Coupled Eulerian/Lagrangian Simulation for Overpressure Structural Response
NASA Astrophysics Data System (ADS)
Lloyd, Andrew; Pan, Hua; Miller, David; Cogar, John
2011-06-01
Accurately modeling blast dynamics is critical in the assessment of structures subjected to blast loading. The current industry standard for modeling blast effects in Lagrangian based Finite Element simulations is CONWEP; tabulated pressure data taken directly from blast events. CONWEP is limited, however, and may not always be physically representative of the blast/structural interaction that occurs in the field. Eulerian hydrocodes provide advantages over CONWEP in that they can capture shock front interaction and model blast surface interfaces with fidelity due to the presence of the working fluid. Eulerian codes, however, break down over larger time scales; whereas, Lagrangian modeling allows for discrete finite elements with definable boundary interfaces that can be tracked out to longer time scales. Hence, a hybrid approach that couples the Eulerian blast modeling with Lagrangian system dynamics is necessary. The objective of this paper is to demonstrate improvements in overpressure structural response modeling using a Coupled Eulerian/Lagrangian algorithm implemented in VelodyneTM. Velodyne results using the Coupled Eulerian/Lagrangian algorithm are compared to results from Eulerian hydrocode simulations and Velodyne simulations using the CONWEP algorithm.
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
Self-correcting Multigrid Solver
Jerome L.V. Lewandowski
2004-06-29
A new multigrid algorithm based on the method of self-correction for the solution of elliptic problems is described. The method exploits information contained in the residual to dynamically modify the source term (right-hand side) of the elliptic problem. It is shown that the self-correcting solver is more efficient at damping the short wavelength modes of the algebraic error than its standard equivalent. When used in conjunction with a multigrid method, the resulting solver displays an improved convergence rate with no additional computational work.
Scalable Parallel Algebraic Multigrid Solvers
Bank, R; Lu, S; Tong, C; Vassilevski, P
2005-03-23
The authors propose a parallel algebraic multilevel algorithm (AMG), which has the novel feature that the subproblem residing in each processor is defined over the entire partition domain, although the vast majority of unknowns for each subproblem are associated with the partition owned by the corresponding processor. This feature ensures that a global coarse description of the problem is contained within each of the subproblems. The advantages of this approach are that interprocessor communication is minimized in the solution process while an optimal order of convergence rate is preserved; and the speed of local subproblem solvers can be maximized using the best existing sequential algebraic solvers.
General complex polynomial root solver
NASA Astrophysics Data System (ADS)
Skowron, J.; Gould, A.
2012-12-01
This general complex polynomial root solver, implemented in Fortran and further optimized for binary microlenses, uses a new algorithm to solve polynomial equations and is 1.6-3 times faster than the ZROOTS subroutine that is commercially available from Numerical Recipes, depending on application. The largest improvement, when compared to naive solvers, comes from a fail-safe procedure that permits skipping the majority of the calculations in the great majority of cases, without risking catastrophic failure in the few cases that these are actually required.
NASA Astrophysics Data System (ADS)
Erzincanli, Belkis; Sahin, Mehmet
2013-12-01
An Arbitrary Lagrangian-Eulerian (ALE) formulation based on the unstructured finite volume method is proposed for solving moving boundary problems with large displacements and rotations. The numerical method is based on the side-centered arrangement of the primitive variables that does not require any ad-hoc modifications in order to enhance pressure coupling. The continuity equation is satisfied within each element at machine precision and the summation of the continuity equations can be exactly reduced to the domain boundary, which is important for the global mass conservation. A special attention is given to construct an ALE algorithm obeying the discrete geometric conservation law (DGCL). The mesh deformation algorithm is based on the indirect Radial Basis Function (RBF) algorithm at each time level while avoiding remeshing in order to enhance numerical robustness. For the parallel solution of resulting large-scale algebraic equations in a fully coupled form, a matrix factorization is introduced similar to that of the projection method for the whole system and the parallel algebraic multigrid solver BoomerAMG is used for the scaled discrete Laplacian provided by the HYPRE library which we access through the PETSc library. The present numerical algorithm is initially validated for the decaying Taylor-Green vortex flow, the flow past an oscillating circular cylinder in a channel and the flow induced by an oscillating sphere in a cubic cavity. Then the numerical algorithm is applied to the numerical simulation of flow field around a pair of flapping Drosophila wings in hover flight. The time variation of the Eulerian coherent structures in the near wake is shown along with the aerodynamic loads.
NASA Astrophysics Data System (ADS)
Del Zanna, L.; Zanotti, O.; Bucciantini, N.; Londrillo, P.
2007-10-01
Aims:We present a new numerical code, ECHO, based on a Eulerian conservative high-order scheme for time dependent three-dimensional general relativistic magnetohydrodynamics (GRMHD) and magnetodynamics (GRMD). ECHO is aimed at providing a shock-capturing conservative method able to work at an arbitrary level of formal accuracy (for smooth flows), where the other existing GRMHD and GRMD schemes yield an overall second order at most. Moreover, our goal is to present a general framework based on the 3+1 Eulerian formalism, allowing for different sets of equations and different algorithms and working in a generic space-time metric, so that ECHO may be easily coupled to any solver for Einstein's equations. Methods: Our finite-difference conservative scheme previously developed for special relativistic hydrodynamics and MHD is extended here to the general relativistic case. Various high-order reconstruction methods are implemented and a two-wave approximate Riemann solver is used. The induction equation is treated by adopting the upwind constrained transport (UCT) procedures, appropriate to preserving the divergence-free condition of the magnetic field in shock-capturing methods. The limiting case of magnetodynamics (also known as force-free degenerate electrodynamics) is implemented by simply replacing the fluid velocity with the electromagnetic drift velocity and by neglecting the contribution of matter to the stress tensor. Results: ECHO is particularly accurate, efficient, versatile, and robust. It has been tested against several astrophysical applications, like magnetized accretion onto black holes and constant angular momentum thick disks threaded by toroidal fields. A novel test of the propagation of large-amplitude, circularly polarized Alfvén waves is proposed, and this allows us to prove the spatial and temporal high-order properties of ECHO very accurately. In particular, we show that reconstruction based on a monotonicity-preserving (MP) filter applied to a
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.
EULERIAN DECORRELATION OF FLUCTUATIONS IN THE INTERPLANETARY MAGNETIC FIELD
Matthaeus, W. H.; Osman, K. T.; Dasso, S.; Weygand, J. M.; Kivelson, M. G.
2010-09-20
A method is devised for estimating the two-time correlation function and the associated Eulerian decorrelation timescale in turbulence. With the assumptions of a single decorrelation time and a frozen-in flow approximation for the single-point analysis, the method compares two-point correlation measurements with single-point correlation measurements at the corresponding spatial lag. This method is applied to interplanetary magnetic field measurements from the Advanced Composition Explorer and Wind spacecraft. An average Eulerian decorrelation time of 2.9 hr is found. This measures the total rate of distortion of turbulent fluid elements-including sweeping, nonlinear distortion, and wave propagation.
Constructing the Lagrangian in the Eulerian coordinate for relativistic hydrodynamics
NASA Astrophysics Data System (ADS)
Chiueh, Tzihong
1994-02-01
The Lorentz-covariant Lagrangian for an ideal relativistic flow is constructed in the Eulerian coordinate. In contrast to the Lagrangian of nonrelativistic flows in the Eulerian formulation, for which the continuity equation is required to be externally imposed as a constraint [Mittag, Stephen, and Yourgrau, in Variational Principles in Dynamics .ul and Quantum Theory, edited by W. Yourgrau and X. Mandelstam (Dover, New York, 1968)], this Lorentz-covariant Lagrangian automatically yields the continuity equation as well as the equation of state. In addition, the relativistic generalization of the Bernoulli equation can also be derived from the present formulation.
Eulerian Time-Domain Filtering for Spatial LES
NASA Technical Reports Server (NTRS)
Pruett, C. David
1997-01-01
Eulerian time-domain filtering seems to be appropriate for LES (large eddy simulation) of flows whose large coherent structures convect approximately at a common characteristic velocity; e.g., mixing layers, jets, and wakes. For these flows, we develop an approach to LES based on an explicit second-order digital Butterworth filter, which is applied in,the time domain in an Eulerian context. The approach is validated through a priori and a posteriori analyses of the simulated flow of a heated, subsonic, axisymmetric jet.
Time-domain Raman analytical forward solvers.
Martelli, Fabrizio; Binzoni, Tiziano; Sekar, Sanathana Konugolu Venkata; Farina, Andrea; Cavalieri, Stefano; Pifferi, Antonio
2016-09-01
A set of time-domain analytical forward solvers for Raman signals detected from homogeneous diffusive media is presented. The time-domain solvers have been developed for two geometries: the parallelepiped and the finite cylinder. The potential presence of a background fluorescence emission, contaminating the Raman signal, has also been taken into account. All the solvers have been obtained as solutions of the time dependent diffusion equation. The validation of the solvers has been performed by means of comparisons with the results of "gold standard" Monte Carlo simulations. These forward solvers provide an accurate tool to explore the information content encoded in the time-resolved Raman measurements. PMID:27607645
A generalized gyrokinetic Poisson solver
Lin, Z.; Lee, W.W.
1995-03-01
A generalized gyrokinetic Poisson solver has been developed, which employs local operations in the configuration space to compute the polarization density response. The new technique is based on the actual physical process of gyrophase-averaging. It is useful for nonlocal simulations using general geometry equilibrium. Since it utilizes local operations rather than the global ones such as FFT, the new method is most amenable to massively parallel algorithms.
On unstructured grids and solvers
NASA Technical Reports Server (NTRS)
Barth, T. J.
1990-01-01
The fundamentals and the state-of-the-art technology for unstructured grids and solvers are highlighted. Algorithms and techniques pertinent to mesh generation are discussed. It is shown that grid generation and grid manipulation schemes rely on fast multidimensional searching. Flow solution techniques for the Euler equations, which can be derived from the integral form of the equations are discussed. Sample calculations are also provided.
Numerical Study of Coal Gasification Using Eulerian-Eulerian Multiphase Model
Shi, S.; Guenther, C.; Orsino, S.
2007-09-01
Gasification converts the carbon-containing material into a synthesis gas (syngas) which can be used as a fuel to generate electricity or used as a basic chemical building block for a large number of uses in the petrochemical and refining industries. Based on the mode of conveyance of the fuel and the gasifying medium, gasification can be classified into fixed or moving bed, fluidized bed, and entrained flow reactors. Entrained flow gasifiers normally feature dilute flow with small particle size and can be successfully modeled with the Discrete Phase Method (DPM). For the other types, the Eulerian-Eulerian (E-E) or the so called two-fluid multiphase model is a more appropriate approach. The E-E model treats the solid phase as a distinct interpenetrating granular “fluid” and it is the most general-purposed multi-fluid model. This approach provides transient, three-dimensional, detailed information inside the reactor which would otherwise be unobtainable through experiments due to the large scale, high pressure and/or temperature. In this paper, a transient, three-dimensional model of the Power Systems Development Facility (PSDF) transport gasifier will be presented to illustrate how Computational Fluid Dynamics (CFD) can be used for large-scale complicated geometry with detailed physics and chemistry. In the model, eleven species are included in the gas phase while four pseudo-species are assumed in the solid phase. A total of sixteen reactions, both homogeneous (involving only gas phase species) and heterogeneous (involving species in both gas and solid phases), are used to model the coal gasification chemistry. Computational results have been validated against PSDF experimental data from lignite to bituminous coals under both air and oxygen blown conditions. The PSDF gasifier geometry was meshed with about 70,000, hexahedra-dominated cells. A total of six cases with different coal, feed gas, and/or operation conditions have been performed. The predicted and
Numerical Simulation of Drophila Flight Based on Arbitrary Langrangian-Eulerian Method
NASA Astrophysics Data System (ADS)
Erzincanli, Belkis; Sahin, Mehmet
2012-11-01
A parallel unstructured finite volume algorithm based on Arbitrary Lagrangian Eulerian (ALE) method has been developed in order to investigate the wake structure around a pair of flapping Drosophila wings. The numerical method uses a side-centered arrangement of the primitive variables that does not require any ad-hoc modifications in order to enhance pressure coupling. A radial basis function (RBF) interpolation method is also implemented in order to achieve large mesh deformations. For the parallel solution of resulting large-scale algebraic equations, a matrix factorization is introduced similar to that of the projection method for the whole coupled system and two-cycle of BoomerAMG solver is used for the scaled discrete Laplacian provided by the HYPRE library which we access through the PETSc library. The present numerical algorithm is initially validated for the flow past an oscillating circular cylinder in a channel and the flow induced by an oscillating sphere in a cubic cavity. Then the numerical algorithm is applied to the numerical simulation of flow field around a pair of flapping Drosophila wing in hover flight. The time variation of the near wake structure is shown along with the aerodynamic loads and particle traces. The authors acknowledge financial support from Turkish National Scientific and Technical Research Council (TUBITAK) through project number 111M332. The authors would like to thank Michael Dickinson and Michael Elzinga for providing the experimental data.
QUANTIFYING SUBGRID POLLUTANT VARIABILITY IN EULERIAN AIR QUALITY MODELS
In order to properly assess human risk due to exposure to hazardous air pollutants or air toxics, detailed information is needed on the location and magnitude of ambient air toxic concentrations. Regional scale Eulerian air quality models are typically limited to relatively coar...
Crystal level simulations using Eulerian finite element methods
Becker, R; Barton, N R; Benson, D J
2004-02-06
Over the last several years, significant progress has been made in the use of crystal level material models in simulations of forming operations. However, in Lagrangian finite element approaches simulation capabilities are limited in many cases by mesh distortion associated with deformation heterogeneity. Contexts in which such large distortions arise include: bulk deformation to strains approaching or exceeding unity, especially in highly anisotropic or multiphase materials; shear band formation and intersection of shear bands; and indentation with sharp indenters. Investigators have in the past used Eulerian finite element methods with material response determined from crystal aggregates to study steady state forming processes. However, Eulerian and Arbitrary Lagrangian-Eulerian (ALE) finite element methods have not been widely utilized for simulation of transient deformation processes at the crystal level. The advection schemes used in Eulerian and ALE codes control mesh distortion and allow for simulation of much larger total deformations. We will discuss material state representation issues related to advection and will present results from ALE simulations.
Eulerian-Lagrangian solution of the convection-dispersion equation in natural co-ordinates.
Cheng, R.T.; Casulli, V.; Milford, S.N.
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.-from Authors
Extension of the Time-Spectral Approach to Overset Solvers for Arbitrary Motion
NASA Technical Reports Server (NTRS)
Leffell, Joshua Isaac; Murman, Scott M.; Pulliam, Thomas H.
2012-01-01
demonstrated marked success in reducing the computational costs associated with simulating periodic forced flows, but have yet to be fully applied to overset or Cartesian solvers for arbitrary motion with dynamic hole-cutting. Overset and Cartesian grid methodologies are versatile techniques capable of handling complex geometry configurations in practical engineering applications, and the combination of the Time-Spectral approach with this general capability potentially provides an enabling new design and analysis tool. In an arbitrary moving-body scenario for these approaches, a Lagrangian body moves through a fixed Eulerian mesh and mesh points in the Eulerian mesh interior to the solid body are removed (cut or blanked), leaving a hole in the Eulerian mesh. During the dynamic motion some gridpoints in the domain are blanked and do not have a complete set of time-samples preventing a direct implementation of the Time-Spectral method. Murman[6] demonstrated the Time-Spectral approach for a Cartesian solver with a rigid domain motion, wherein the hole cutting remains constant. Similarly, Custer et al. [15, 16] used the NASA overset OVERFLOW solver and limited the amount of relative motion to ensure static hole-cutting and interpolation. Recently, Mavriplis and Mundis[17] demonstrated a qualitative method for applying the Time-Spectral approach to an unstructured overset solver for arbitrary motion. The goal of the current work is to develop a robust and general method for handling arbitrary motion with the Time-Spectral approach within an overset or Cartesian mesh method, while still approaching the spectral convergence rate of the original Time-Spectral approach. The viscous OVERFLOW solver will be augmented with the new Time-Spectral algorithm and the capability of the method for benchmark problems in rotorcraft and turbomachinery will be demonstrated. This abstract begins with a brief synopsis of the Time-Spectral approach for overset grids and provides details of e current
Finite Element Interface to Linear Solvers
Williams, Alan
2005-03-18
Sparse systems of linear equations arise in many engineering applications, including finite elements, finite volumes, and others. The solution of linear systems is often the most computationally intensive portion of the application. Depending on the complexity of problems addressed by the application, there may be no single solver capable of solving all of the linear systems that arise. This motivates the desire to switch an application from one solver librwy to another, depending on the problem being solved. The interfaces provided by solver libraries differ greatly, making it difficult to switch an application code from one library to another. The amount of library-specific code in an application Can be greatly reduced by having an abstraction layer between solver libraries and the application, putting a common "face" on various solver libraries. One such abstraction layer is the Finite Element Interface to Linear Solvers (EEl), which has seen significant use by finite element applications at Sandia National Laboratories and Lawrence Livermore National Laboratory.
Iannaccone, Francesco; Degroote, Joris; Vierendeels, Jan; Segers, Patrick
2016-01-01
In recent years the role of FSI (fluid-structure interaction) simulations in the analysis of the fluid-mechanics of heart valves is becoming more and more important, being able to capture the interaction between the blood and both the surrounding biological tissues and the valve itself. When setting up an FSI simulation, several choices have to be made to select the most suitable approach for the case of interest: in particular, to simulate flexible leaflet cardiac valves, the type of discretization of the fluid domain is crucial, which can be described with an ALE (Arbitrary Lagrangian-Eulerian) or an Eulerian formulation. The majority of the reported 3D heart valve FSI simulations are performed with the Eulerian formulation, allowing for large deformations of the domains without compromising the quality of the fluid grid. Nevertheless, it is known that the ALE-FSI approach guarantees more accurate results at the interface between the solid and the fluid. The goal of this paper is to describe the same aortic valve model in the two cases, comparing the performances of an ALE-based FSI solution and an Eulerian-based FSI approach. After a first simplified 2D case, the aortic geometry was considered in a full 3D set-up. The model was kept as similar as possible in the two settings, to better compare the simulations’ outcomes. Although for the 2D case the differences were unsubstantial, in our experience the performance of a full 3D ALE-FSI simulation was significantly limited by the technical problems and requirements inherent to the ALE formulation, mainly related to the mesh motion and deformation of the fluid domain. As a secondary outcome of this work, it is important to point out that the choice of the solver also influenced the reliability of the final results. PMID:27128798
Bavo, Alessandra M; Rocatello, Giorgia; Iannaccone, Francesco; Degroote, Joris; Vierendeels, Jan; Segers, Patrick
2016-01-01
In recent years the role of FSI (fluid-structure interaction) simulations in the analysis of the fluid-mechanics of heart valves is becoming more and more important, being able to capture the interaction between the blood and both the surrounding biological tissues and the valve itself. When setting up an FSI simulation, several choices have to be made to select the most suitable approach for the case of interest: in particular, to simulate flexible leaflet cardiac valves, the type of discretization of the fluid domain is crucial, which can be described with an ALE (Arbitrary Lagrangian-Eulerian) or an Eulerian formulation. The majority of the reported 3D heart valve FSI simulations are performed with the Eulerian formulation, allowing for large deformations of the domains without compromising the quality of the fluid grid. Nevertheless, it is known that the ALE-FSI approach guarantees more accurate results at the interface between the solid and the fluid. The goal of this paper is to describe the same aortic valve model in the two cases, comparing the performances of an ALE-based FSI solution and an Eulerian-based FSI approach. After a first simplified 2D case, the aortic geometry was considered in a full 3D set-up. The model was kept as similar as possible in the two settings, to better compare the simulations' outcomes. Although for the 2D case the differences were unsubstantial, in our experience the performance of a full 3D ALE-FSI simulation was significantly limited by the technical problems and requirements inherent to the ALE formulation, mainly related to the mesh motion and deformation of the fluid domain. As a secondary outcome of this work, it is important to point out that the choice of the solver also influenced the reliability of the final results. PMID:27128798
Analysis Tools for CFD Multigrid Solvers
NASA Technical Reports Server (NTRS)
Mineck, Raymond E.; Thomas, James L.; Diskin, Boris
2004-01-01
Analysis tools are needed to guide the development and evaluate the performance of multigrid solvers for the fluid flow equations. Classical analysis tools, such as local mode analysis, often fail to accurately predict performance. Two-grid analysis tools, herein referred to as Idealized Coarse Grid and Idealized Relaxation iterations, have been developed and evaluated within a pilot multigrid solver. These new tools are applicable to general systems of equations and/or discretizations and point to problem areas within an existing multigrid solver. Idealized Relaxation and Idealized Coarse Grid are applied in developing textbook-efficient multigrid solvers for incompressible stagnation flow problems.
Accurate direct Eulerian simulation of dynamic elastic-plastic flow
Kamm, James R; Walter, John W
2009-01-01
The simulation of dynamic, large strain deformation is an important, difficult, and unsolved computational challenge. Existing Eulerian schemes for dynamic material response are plagued by unresolved issues. We present a new scheme for the first-order system of elasto-plasticity equations in the Eulerian frame. This system has an intrinsic constraint on the inverse deformation gradient. Standard Godunov schemes do not satisfy this constraint. The method of Flux Distributions (FD) was devised to discretely enforce such constraints for numerical schemes with cell-centered variables. We describe a Flux Distribution approach that enforces the inverse deformation gradient constraint. As this approach is new and novel, we do not yet have numerical results to validate our claims. This paper is the first installment of our program to develop this new method.
A Dynamically Adaptive Arbitrary Lagrangian-Eulerian Method for Hydrodynamics
Anderson, R W; Pember, R B; Elliott, N S
2004-01-28
A new method that combines staggered grid Arbitrary Lagrangian-Eulerian (ALE) techniques with structured local adaptive mesh refinement (AMR) has been developed for solution of the Euler equations. The novel components of the combined ALE-AMR method hinge upon the integration of traditional AMR techniques with both staggered grid Lagrangian operators as well as elliptic relaxation operators on moving, deforming mesh hierarchies. Numerical examples demonstrate the utility of the method in performing detailed three-dimensional shock-driven instability calculations.
A Dynamically Adaptive Arbitrary Lagrangian-Eulerian Method for Hydrodynamics
Anderson, R W; Pember, R B; Elliott, N S
2002-10-19
A new method that combines staggered grid Arbitrary Lagrangian-Eulerian (ALE) techniques with structured local adaptive mesh refinement (AMR) has been developed for solution of the Euler equations. The novel components of the combined ALE-AMR method hinge upon the integration of traditional AMR techniques with both staggered grid Lagrangian operators as well as elliptic relaxation operators on moving, deforming mesh hierarchies. Numerical examples demonstrate the utility of the method in performing detailed three-dimensional shock-driven instability calculations.
Modified Eulerian-Lagrangian formulation for hydrodynamic modeling
NASA Astrophysics Data System (ADS)
Sorek, Shaul; Borisov, Vyacheslav
2012-04-01
We present the modified Eulerian-Lagrangian (MEL) formulation, based on non-divergent forms of partial differential balance equations, for simulating transport of extensive quantities in a porous medium. Hydrodynamic derivatives are written in terms of modified velocities for particles propagating phase and component quantities along their respective paths. The particles physically interpreted velocities also address the heterogeneity of the matrix and fluid properties. The MEL formulation is also implemented to parabolic Partial Differential Equations (PDE's) as these are shown to be interchangeable with equivalent PDE's having hyperbolic - parabolic characteristics, without violating the same physical concepts. We prove that the MEL schemes provide a convergent and monotone approximation also to PDE's with discontinuous coefficients. An extension to the Peclet number is presented that also accounts for advective dominant PDE's with no reference to the fluid velocity or even when this velocity is not introduced. In Sorek et al. [27], a mathematical analysis for a linear system of coupled PDE's and an example of nonlinear PDE's, proved that the finite difference MEL, unlike an Eulerian scheme, guaranties the absence of spurious oscillations. Currently, we present notions of monotone interpolation associated with the MEL particle tracking procedure and prove the convergence of the MEL schemes to the original balance equation also for discontinuous coefficients on the basis of difference schemes approximating PDE's. We provide numerical examples, also with highly random fields of permeabilities and/or dispersivities, suggesting that the MEL scheme produces resolutions that are more consistent with the physical phenomenon in comparison to the Eulerian and the Eulerian-Lagrangian (EL) schemes.
A domain decomposition scheme for Eulerian shock physics codes
Bell, R.L.; Hertel, E.S. Jr.
1994-08-01
A new algorithm which allows for complex domain decomposition in Eulerian codes was developed at Sandia National Laboratories. This new feature allows a user to customize the zoning for each portion of a calculation and to refine volumes of the computational space of particular interest This option is available in one, two, and three dimensions. The new technique will be described in detail and several examples of the effectiveness of this technique will also be discussed.
An overview of mesoscale material modeling with Eulerian hydrocodes
NASA Astrophysics Data System (ADS)
Olney, K.; Benson, D. J.; Nesterenko, V.
2014-05-01
Eulerian hydrocodes were originally developed for simulating strong shocks in solids and fluids, but their ability to handle arbitrarily large deformations and the formation of new free surfaces makes them attractive for simulating the deformation and failure of materials at the mesoscopic scale. A summary of several numerical techniques that have been developed to address issues that commonly arise for this class of problems is presented.
The hybrid Eulerian Lagrangian numerical scheme tested with Chemistry
NASA Astrophysics Data System (ADS)
Hansen, A. B.; Sørensen, B.; Tarning-Andersen, P.; Christensen, J. H.; Brandt, J.; Kaas, E.
2012-12-01
A newly developed transport scheme, the Hybrid Eulerian Lagrangian (HEL) scheme, has been tested using a module for atmospheric chemistry, including 58 chemical species, and compared to two other traditional advection schemes; a classical pseudospectral Eulerian method the Accurate Space Derivative (ASD) scheme and the bi-cubic semi-Lagrangian (SL) scheme using classical rotation tests. The rotation tests have been designed to test and compare the advection schemes for different spatial and temporal resolutions in different chemical conditions (rural and urban) and for different shapes (cone and slotted cylinder). This gives the advection schemes different challenges with respect to relatively slow or fast chemistry and smooth or sharp gradients. In every test, error measures have been calculated and used for ranking of the advection schemes with respect to performance, i.e. lowest overall errors for all chemical species. The results presented show that the new transport scheme, HEL, by far outperforms both the Eulerian and semi-Lagrangian schemes with very low error estimates compared to the two other schemes.
The hybrid Eulerian Lagrangian numerical scheme tested with Chemistry
NASA Astrophysics Data System (ADS)
Hansen, A. B.; Sørensen, B.; Tarning-Andersen, P.; Christensen, J. H.; Brandt, J.; Kaas, E.
2012-11-01
A newly developed advection scheme, the Hybrid Eulerian Lagrangian (HEL) scheme, has been tested, including a module for atmospheric chemistry, including 58 chemical species, and compared to two other traditional advection schemes; a classical pseudospectral Eulerian method the Accurate Space Derivative (ASD) scheme and the bi-cubic semi-Lagrangian (SL) scheme using classical rotation tests. The rotation tests have been designed to test and compare the advection schemes for different spatial and temporal resolutions in different chemical conditions (rural and urban) and for different shapes (cone and slotted cylinder) giving the advection schemes different challenges with respect to relatively slow or fast chemistry and smooth or sharp gradients, respectively. In every test, error measures have been calculated and used for ranking of the advection schemes with respect to performance, i.e. lowest overall errors for all chemical species. Furthermore, the HEL and SL schemes have been compared in a shallow water model, demonstrating the performance in a more realistic non-linear deformation flow. The results in this paper show that the new advection scheme, HEL, by far outperforms both the Eulerian and semi-Lagrangian schemes with very low error estimates compared to the two other schemes. Although no analytic solution can be obtained for the performance in the non-linear shallow water model flow, the tracer distribution appears realistic as compared to LMCSL when a mixing between local parcel concentrations is introduced in HEL.
MACSYMA's symbolic ordinary differential equation solver
NASA Technical Reports Server (NTRS)
Golden, J. P.
1977-01-01
The MACSYMA's symbolic ordinary differential equation solver ODE2 is described. The code for this routine is delineated, which is of interest because it is written in top-level MACSYMA language, and may serve as a good example of programming in that language. Other symbolic ordinary differential equation solvers are mentioned.
KLU2 Direct Linear Solver Package
Energy Science and Technology Software Center (ESTSC)
2012-01-04
KLU2 is a direct sparse solver for solving unsymmetric linear systems. It is related to the existing KLU solver, (in Amesos package and also as a stand-alone package from University of Florida) but provides template support for scalar and ordinal types. It uses a left looking LU factorization method.
Improving Resource-Unaware SAT Solvers
NASA Astrophysics Data System (ADS)
Hölldobler, Steffen; Manthey, Norbert; Saptawijaya, Ari
The paper discusses cache utilization in state-of-the-art SAT solvers. The aim of the study is to show how a resource-unaware SAT solver can be improved by utilizing the cache sensibly. The analysis is performed on a CDCL-based SAT solver using a subset of the industrial SAT Competition 2009 benchmark. For the analysis, the total cycles, the resource stall cycles, the L2 cache hits and the L2 cache misses are traced using sample based profiling. Based on the analysis, several techniques - some of which have not been used in SAT solvers so far - are proposed resulting in a combined speedup up to 83% without affecting the search path of the solver. The average speedup on the benchmark is 60%. The new techniques are also applied to MiniSAT2.0 improving its runtime by 20% on average.
Belos Block Linear Solvers Package
Energy Science and Technology Software Center (ESTSC)
2004-03-01
Belos is an extensible and interoperable framework for large-scale, iterative methods for solving systems of linear equations with multiple right-hand sides. The motivation for this framework is to provide a generic interface to a collection of algorithms for solving large-scale linear systems. Belos is interoperable because both the matrix and vectors are considered to be opaque objects--only knowledge of the matrix and vectors via elementary operations is necessary. An implementation of Balos is accomplished viamore » the use of interfaces. One of the goals of Belos is to allow the user flexibility in specifying the data representation for the matrix and vectors and so leverage any existing software investment. The algorithms that will be included in package are Krylov-based linear solvers, like Block GMRES (Generalized Minimal RESidual) and Block CG (Conjugate-Gradient).« less
On two-phase flow solvers in irregular domains with contact line
NASA Astrophysics Data System (ADS)
Lepilliez, Mathieu; Popescu, Elena Roxana; Gibou, Frederic; Tanguy, Sébastien
2016-09-01
We present numerical methods that enable the direct numerical simulation of two-phase flows in irregular domains. A method is presented to account for surface tension effects in a mesh cell containing a triple line between the liquid, gas and solid phases. Our numerical method is based on the level-set method to capture the liquid-gas interface and on the single-phase Navier-Stokes solver in irregular domain proposed in [35] to impose the solid boundary in an Eulerian framework. We also present a strategy for the implicit treatment of the viscous term and how to impose both a Neumann boundary condition and a jump condition when solving for the pressure field. Special care is given on how to take into account the contact angle, the no-slip boundary condition for the velocity field and the volume forces. Finally, we present numerical results in two and three spatial dimensions evaluating our simulations with several benchmarks.
ALPS - A LINEAR PROGRAM SOLVER
NASA Technical Reports Server (NTRS)
Viterna, L. A.
1994-01-01
Linear programming is a widely-used engineering and management tool. Scheduling, resource allocation, and production planning are all well-known applications of linear programs (LP's). Most LP's are too large to be solved by hand, so over the decades many computer codes for solving LP's have been developed. ALPS, A Linear Program Solver, is a full-featured LP analysis program. ALPS can solve plain linear programs as well as more complicated mixed integer and pure integer programs. ALPS also contains an efficient solution technique for pure binary (0-1 integer) programs. One of the many weaknesses of LP solvers is the lack of interaction with the user. ALPS is a menu-driven program with no special commands or keywords to learn. In addition, ALPS contains a full-screen editor to enter and maintain the LP formulation. These formulations can be written to and read from plain ASCII files for portability. For those less experienced in LP formulation, ALPS contains a problem "parser" which checks the formulation for errors. ALPS creates fully formatted, readable reports that can be sent to a printer or output file. ALPS is written entirely in IBM's APL2/PC product, Version 1.01. The APL2 workspace containing all the ALPS code can be run on any APL2/PC system (AT or 386). On a 32-bit system, this configuration can take advantage of all extended memory. The user can also examine and modify the ALPS code. The APL2 workspace has also been "packed" to be run on any DOS system (without APL2) as a stand-alone "EXE" file, but has limited memory capacity on a 640K system. A numeric coprocessor (80X87) is optional but recommended. The standard distribution medium for ALPS is a 5.25 inch 360K MS-DOS format diskette. IBM, IBM PC and IBM APL2 are registered trademarks of International Business Machines Corporation. MS-DOS is a registered trademark of Microsoft Corporation.
GARDNER, P.R.
2006-04-01
Sudoku, also known as Number Place, is a logic-based placement puzzle. The aim of the puzzle is to enter a numerical digit from 1 through 9 in each cell of a 9 x 9 grid made up of 3 x 3 subgrids (called ''regions''), starting with various digits given in some cells (the ''givens''). Each row, column, and region must contain only one instance of each numeral. Completing the puzzle requires patience and logical ability. Although first published in a U.S. puzzle magazine in 1979, Sudoku initially caught on in Japan in 1986 and attained international popularity in 2005. Last fall, after noticing Sudoku puzzles in some newspapers and magazines, I attempted a few just to see how hard they were. Of course, the difficulties varied considerably. ''Obviously'' one could use Trial and Error but all the advice was to ''Use Logic''. Thinking to flex, and strengthen, those powers, I began to tackle the puzzles systematically. That is, when I discovered a new tactical rule, I would write it down, eventually generating a list of ten or so, with some having overlap. They served pretty well except for the more difficult puzzles, but even then I managed to develop an additional three rules that covered all of them until I hit the Oregonian puzzle shown. With all of my rules, I could not seem to solve that puzzle. Initially putting my failure down to rapid mental fatigue (being unable to hold a sufficient quantity of information in my mind at one time), I decided to write a program to implement my rules and see what I had failed to notice earlier. The solver, too, failed. That is, my rules were insufficient to solve that particular puzzle. I happened across a book written by a fellow who constructs such puzzles and who claimed that, sometimes, the only tactic left was trial and error. With a trial and error routine implemented, my solver successfully completed the Oregonian puzzle, and has successfully solved every puzzle submitted to it since.
SIERRA framework version 4 : solver services.
Williams, Alan B.
2005-02-01
Several SIERRA applications make use of third-party libraries to solve systems of linear and nonlinear equations, and to solve eigenproblems. The classes and interfaces in the SIERRA framework that provide linear system assembly services and access to solver libraries are collectively referred to as solver services. This paper provides an overview of SIERRA's solver services including the design goals that drove the development, and relationships and interactions among the various classes. The process of assembling and manipulating linear systems will be described, as well as access to solution methods and other operations.
A scalable 2-D parallel sparse solver
Kothari, S.C.; Mitra, S.
1995-12-01
Scalability beyond a small number of processors, typically 32 or less, is known to be a problem for existing parallel general sparse (PGS) direct solvers. This paper presents a parallel general sparse PGS direct solver for general sparse linear systems on distributed memory machines. The algorithm is based on the well-known sequential sparse algorithm Y12M. To achieve efficient parallelization, a 2-D scattered decomposition of the sparse matrix is used. The proposed algorithm is more scalable than existing parallel sparse direct solvers. Its scalability is evaluated on a 256 processor nCUBE2s machine using Boeing/Harwell benchmark matrices.
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.
An Overview of Mesoscale Material Modeling with Eulerian Hydrocodes
NASA Astrophysics Data System (ADS)
Benson, David
2013-06-01
Eulerian hydrocodes were originally developed for simulating strong shocks in solids and fluids, but their ability to handle arbitrarily large deformations and the formation of new free surfaces makes them attractive for simulating the deformation and failure of materials at the mesoscopic scale. A summary of some of the numerical techniques that have been developed to address common issues for this class of problems is presented with the shock compression of powders used as a model problem. Achieving the correct packing density with the correct statistical distribution of particle sizes and shapes is, in itself, a challenging problem. However, since Eulerian codes permit multiple materials within each element, or cell, the material interfaces do not have to follow the mesh lines. The use of digital image processing to map the pixels of micrographs to the Eulerian mesh has proven to be a popular and useful means of creating accurate models of complex microstructures. Micro CT scans have been used to extend this approach to three dimensions for several classes of materials. The interaction between the particles is of considerable interest. During shock compression, individual particles may melt and form jets, and the voids between them collapse. Dynamic interface ordering has become a necessity, and many codes now have a suite of options for handling multi-material mechanics. True contact algorithms are now replacing multi-material approximations in some cases. At the mesoscale, material properties often vary spatially due to sub-scale effects. Using a large number of material species to represent the variations is usually unattractive. Directly specifying the properties point-wise as history variables has not proven successful because the limiters in the transport algorithms quickly smooth out the variations. Circumventing the limiter problem is shown to be relatively simple with the use of a reference configuration and the transport of the initial coordinates
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.
MESA: A 3-D Eulerian hydrocode for penetration mechanics studies
Mandell, D.A.; Holian, K.S.; Henninger, R.
1991-01-01
We describe an explicit, finite-difference hydrocode, called MESA, and compare calculations to metal and ceramic plate impacts with spall and to Taylor cylinder tests. The MESA code was developed with support from DARPA, the Army and the Marine Corps for use in armor/anti-armor problems primarily, but the code has been used for a number of other applications. MESA includes 2-D and 3-D Eulerian hydrodynamics, a number of material strength and fracture models, and a programmed burn high explosives model. 15 refs., 4 figs.
NASA Technical Reports Server (NTRS)
Ferencz, Donald C.; Viterna, Larry A.
1991-01-01
ALPS is a computer program which can be used to solve general linear program (optimization) problems. ALPS was designed for those who have minimal linear programming (LP) knowledge and features a menu-driven scheme to guide the user through the process of creating and solving LP formulations. Once created, the problems can be edited and stored in standard DOS ASCII files to provide portability to various word processors or even other linear programming packages. Unlike many math-oriented LP solvers, ALPS contains an LP parser that reads through the LP formulation and reports several types of errors to the user. ALPS provides a large amount of solution data which is often useful in problem solving. In addition to pure linear programs, ALPS can solve for integer, mixed integer, and binary type problems. Pure linear programs are solved with the revised simplex method. Integer or mixed integer programs are solved initially with the revised simplex, and the completed using the branch-and-bound technique. Binary programs are solved with the method of implicit enumeration. This manual describes how to use ALPS to create, edit, and solve linear programming problems. Instructions for installing ALPS on a PC compatible computer are included in the appendices along with a general introduction to linear programming. A programmers guide is also included for assistance in modifying and maintaining the program.
Euler solvers for transonic applications
NASA Technical Reports Server (NTRS)
Vanleer, Bram
1989-01-01
The 1980s may well be called the Euler era of applied aerodynamics. Computer codes based on discrete approximations of the Euler equations are now routinely used to obtain solutions of transonic flow problems in which the effects of entropy and vorticity production are significant. Such codes can even predict separation from a sharp edge, owing to the inclusion of artificial dissipation, intended to lend numerical stability to the calculation but at the same time enforcing the Kutta condition. One effect not correctly predictable by Euler codes is the separation from a smooth surface, and neither is viscous drag; for these some form of the Navier-Stokes equation is needed. It, therefore, comes as no surprise to observe that the Navier-Stokes has already begun before Euler solutions were fully exploited. Moreover, most numerical developments for the Euler equations are now constrained by the requirement that the techniques introduced, notably artificial dissipation, must not interfere with the new physics added when going from an Euler to a full Navier-Stokes approximation. In order to appreciate the contributions of Euler solvers to the understanding of transonic aerodynamics, it is useful to review the components of these computational tools. Space discretization, time- or pseudo-time marching and boundary procedures, the essential constituents are discussed. The subject of grid generation and grid adaptation to the solution are touched upon only where relevant. A list of unanswered questions and an outlook for the future are covered.
Parallelizing alternating direction implicit solver on GPUs
Technology Transfer Automated Retrieval System (TEKTRAN)
We present a parallel Alternating Direction Implicit (ADI) solver on GPUs. Our implementation significantly improves existing implementations in two aspects. First, we address the scalability issue of existing Parallel Cyclic Reduction (PCR) implementations by eliminating their hardware resource con...
Optimization of solver for gas flow modeling
NASA Astrophysics Data System (ADS)
Savichkin, D.; Dodulad, O.; Kloss, Yu
2014-05-01
The main purpose of the work is optimization of the solver for rarefied gas flow modeling based on the Boltzmann equation. Optimization method is based on SIMD extensions for ×86 processors. Computational code is profiled and manually optimized with SSE instructions. Heat flow, shock waves and Knudsen pump are modeled with optimized solver. Dependencies of computational time from mesh sizes and CPU capabilities are provided.
A parallel PCG solver for MODFLOW.
Dong, Yanhui; Li, Guomin
2009-01-01
In order to simulate large-scale ground water flow problems more efficiently with MODFLOW, the OpenMP programming paradigm was used to parallelize the preconditioned conjugate-gradient (PCG) solver with in this study. Incremental parallelization, the significant advantage supported by OpenMP on a shared-memory computer, made the solver transit to a parallel program smoothly one block of code at a time. The parallel PCG solver, suitable for both MODFLOW-2000 and MODFLOW-2005, is verified using an 8-processor computer. Both the impact of compilers and different model domain sizes were considered in the numerical experiments. Based on the timing results, execution times using the parallel PCG solver are typically about 1.40 to 5.31 times faster than those using the serial one. In addition, the simulation results are the exact same as the original PCG solver, because the majority of serial codes were not changed. It is worth noting that this parallelizing approach reduces cost in terms of software maintenance because only a single source PCG solver code needs to be maintained in the MODFLOW source tree. PMID:19563427
Finite Element Interface to Linear Solvers
Energy Science and Technology Software Center (ESTSC)
2005-03-18
Sparse systems of linear equations arise in many engineering applications, including finite elements, finite volumes, and others. The solution of linear systems is often the most computationally intensive portion of the application. Depending on the complexity of problems addressed by the application, there may be no single solver capable of solving all of the linear systems that arise. This motivates the desire to switch an application from one solver librwy to another, depending on themore » problem being solved. The interfaces provided by solver libraries differ greatly, making it difficult to switch an application code from one library to another. The amount of library-specific code in an application Can be greatly reduced by having an abstraction layer between solver libraries and the application, putting a common "face" on various solver libraries. One such abstraction layer is the Finite Element Interface to Linear Solvers (EEl), which has seen significant use by finite element applications at Sandia National Laboratories and Lawrence Livermore National Laboratory.« less
An Eulerian approach for computing the finite time Lyapunov exponent
NASA Astrophysics Data System (ADS)
Leung, Shingyu
2011-05-01
We propose efficient Eulerian methods for approximating the finite-time Lyapunov exponent (FTLE). The idea is to compute the related flow map using the Level Set Method and the Liouville equation. There are several advantages of the proposed approach. Unlike the usual Lagrangian-type computations, the resulting method requires the velocity field defined only at discrete locations. No interpolation of the velocity field is needed. Also, the method automatically stops a particle trajectory in the case when the ray hits the boundary of the computational domain. The computational complexity of the algorithm is O(Δ x-( d+1) ) with d the dimension of the physical space. Since there are the same number of mesh points in the x- t space, the computational complexity of the proposed Eulerian approach is optimal in the sense that each grid point is visited for only O(1) time. We also extend the algorithm to compute the FTLE on a co-dimension one manifold. The resulting algorithm does not require computation on any local coordinate system and is simple to implement even for an evolving manifold.
Eulerian hydrocode modeling of a dynamic tensile extrusion experiment (u)
Burkett, Michael W; Clancy, Sean P
2009-01-01
Eulerian hydrocode simulations utilizing the Mechanical Threshold Stress flow stress model were performed to provide insight into a dynamic extrusion experiment. The dynamic extrusion response of copper (three different grain sizes) and tantalum spheres were simulated with MESA, an explicit, 2-D Eulerian continuum mechanics hydrocode and compared with experimental data. The experimental data consisted of high-speed images of the extrusion process, recovered extruded samples, and post test metallography. The hydrocode was developed to predict large-strain and high-strain-rate loading problems. Some of the features of the features of MESA include a high-order advection algorithm, a material interface tracking scheme and a van Leer monotonic advection-limiting. The Mechanical Threshold Stress (MTS) model was utilized to evolve the flow stress as a function of strain, strain rate and temperature for copper and tantalum. Plastic strains exceeding 300% were predicted in the extrusion of copper at 400 m/s, while plastic strains exceeding 800% were predicted for Ta. Quantitative comparisons between the predicted and measured deformation topologies and extrusion rate were made. Additionally, predictions of the texture evolution (based upon the deformation rate history and the rigid body rotations experienced by the copper during the extrusion process) were compared with the orientation imaging microscopy measurements. Finally, comparisons between the calculated and measured influence of the initial texture on the dynamic extrusion response of tantalum was performed.
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.
Advanced three-dimensional Eulerian hydrodynamic algorithm development
Rider, W.J.; Kothe, D.B.; Mosso, S.
1998-11-01
This is the final report of a three-year, Laboratory Directed Research and Development (LDRD) project at the Los Alamos National Laboratory (LANL). The purpose of this project is to investigate, implement, and evaluate algorithms that have high potential for improving the robustness, fidelity and accuracy of three-dimensional Eulerian hydrodynamic simulations. Eulerian computations are necessary to simulate a number of important physical phenomena ranging from the molding process for metal parts to nuclear weapons safety issues to astrophysical phenomena such as that associated with a Type 2 supernovae. A number of algorithmic issues were explored in the course of this research including interface/volume tracking, surface physics integration, high resolution integration techniques, multilevel iterative methods, multimaterial hydrodynamics and coupling radiation with hydrodynamics. This project combines core strengths of several Laboratory divisions. The project has high institutional benefit given the renewed emphasis on numerical simulations in Science-Based Stockpile Stewardship and the Accelerated Strategic Computing Initiative and LANL`s tactical goals related to high performance computing and simulation.
PSPIKE: A Parallel Hybrid Sparse Linear System Solver
NASA Astrophysics Data System (ADS)
Manguoglu, Murat; Sameh, Ahmed H.; Schenk, Olaf
The availability of large-scale computing platforms comprised of tens of thousands of multicore processors motivates the need for the next generation of highly scalable sparse linear system solvers. These solvers must optimize parallel performance, processor (serial) performance, as well as memory requirements, while being robust across broad classes of applications and systems. In this paper, we present a new parallel solver that combines the desirable characteristics of direct methods (robustness) and effective iterative solvers (low computational cost), while alleviating their drawbacks (memory requirements, lack of robustness). Our proposed hybrid solver is based on the general sparse solver PARDISO, and the “Spike” family of hybrid solvers. The resulting algorithm, called PSPIKE, is as robust as direct solvers, more reliable than classical preconditioned Krylov subspace methods, and much more scalable than direct sparse solvers. We support our performance and parallel scalability claims using detailed experimental studies and comparison with direct solvers, as well as classical preconditioned Krylov methods.
Heberton, C.I.; Russell, T.F.; Konikow, L.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.
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.
Relativistic perturbations in ΛCDM: Eulerian & Lagrangian approaches
NASA Astrophysics Data System (ADS)
Villa, Eleonora; Rampf, Cornelius
2016-01-01
We study the relativistic dynamics of a pressure-less and irrotational fluid of dark matter (CDM) with a cosmological constant (Λ), up to second order in cosmological perturbation theory. In our analysis we also account for vector and tensor perturbations and include primordial non-Gaussianity. We consider three gauges: the synchronous-comoving gauge, the Poisson gauge and the total matter gauge, where the first is the unique relativistic Lagrangian frame of reference, and the latters are convenient gauge choices for Eulerian frames. Our starting point is the metric and fluid variables in the Poisson gauge up to second order. We then perform the gauge transformations to the synchronous-comoving gauge and subsequently to the total matter gauge. Our expressions for the metrics, densities, velocities, and the gauge generators are novel and coincide with known results in the limit of a vanishing cosmological constant.
The Eulerian time correlation function from direct simulation data
NASA Astrophysics Data System (ADS)
Rubinstein, Robert; He, Guowei
2001-11-01
The Eulerian time correlation function in homogeneous isotropic turbulence is obtained from direct numerical simulation. We develop curvefits of this function using a producure suggested by Boon and Yip (Molecular Hydrodynamics), which develops a continued fraction expansion of the Laplace transform of the time correlation function. Results of different two-pole expressions are compared with the results of the simulations. Good agreement using one such expression is obtained. The curvefit is developed both for the DNS dataset and for the time correlation function computed from LES. The dynamic meaning of the time correlation function in turbulence is compared to the role of the time correlation function in molecular hydrodynamics, where it is associated with the hydrodynamic modes of the fluid.
Full Eulerian lattice Boltzmann model for conjugate heat transfer
NASA Astrophysics Data System (ADS)
Hu, Yang; Li, Decai; Shu, Shi; Niu, Xiaodong
2015-12-01
In this paper a full Eulerian lattice Boltzmann model is proposed for conjugate heat transfer. A unified governing equation with a source term for the temperature field is derived. By introducing the source term, we prove that the continuity of temperature and its normal flux at the interface is satisfied automatically. The curved interface is assumed to be zigzag lines. All physical quantities are recorded and updated on a Cartesian grid. As a result, any complicated treatment near the interface is avoided, which makes the proposed model suitable to simulate the conjugate heat transfer with complex interfaces efficiently. The present conjugate interface treatment is validated by several steady and unsteady numerical tests, including pure heat conduction, forced convection, and natural convection problems. Both flat and curved interfaces are also involved. The obtained results show good agreement with the analytical and/or finite volume results.
Vortex and strain skeletons in Eulerian and Lagrangian frames.
Sahner, Jan; Weinkauf, Tino; Teuber, Nathalie; Hege, Hans-Christian
2007-01-01
We present an approach to analyze mixing in flow fields by extracting vortex and strain features as extremal structures of derived scalar quantities that satisfy a duality property: they indicate vortical as well as high-strain (saddletype) regions. Specifically, we consider the Okubo-Weiss criterion and the recently introduced MZ-criterion. While the first is derived from a purely Eulerian framework, the latter is based on Lagrangian considerations. In both cases high values indicate vortex activity whereas low values indicate regions of high strain. By considering the extremal features of those quantities, we define the notions of a vortex and a strain skeleton in a hierarchical manner: the collection of maximal 0D, 1D and 2D structures assemble the vortex skeleton; the minimal structures identify the strain skeleton. We extract those features using scalar field topology and apply our method to a number of steady and unsteady 3D flow fields. PMID:17622681
Noninvasive Free Flap Monitoring Using Eulerian Video Magnification
Liu, Yuan Fang; Vuong, Christopher; Walker, Paul Charles; Peterson, Nathaniel Ray; Inman, Jared Christian; Filho, Pedro Alcantara Andrade; Lee, Steve Choon-Sung
2016-01-01
Eulerian Video Magnification (EVM) can enhance subtle changes in videos to reveal what was once invisible to the naked eye. In this proof of concept study, we investigated using EVM as a novel form of free flap monitoring. Free flaps with skin paddles were filmed in the operating room with manipulation of their pedicles. In a representative 77-year-old female who received a latissimus dorsi-serratus-rib composite free flap, EVM was able to detect blockage of arterial or venous supply instantaneously, providing a visible representation through degree of color change in videos. EVM has the potential to serve as a powerful free flap monitoring tool with the benefit of being noninvasive, sensitive, easy-to-use, and nearly cost-free. PMID:27092284
Hamiltonian magnetohydrodynamics: Lagrangian, Eulerian, and dynamically accessible stability—Theory
Andreussi, T.; Morrison, P. J.; Pegoraro, F.
2013-09-15
Stability conditions of magnetized plasma flows are obtained by exploiting the Hamiltonian structure of the magnetohydrodynamics (MHD) equations and, in particular, by using three kinds of energy principles. First, the Lagrangian variable energy principle is described and sufficient stability conditions are presented. Next, plasma flows are described in terms of Eulerian variables and the noncanonical Hamiltonian formulation of MHD is exploited. For symmetric equilibria, the energy-Casimir principle is expanded to second order and sufficient conditions for stability to symmetric perturbation are obtained. Then, dynamically accessible variations, i.e., variations that explicitly preserve invariants of the system, are introduced and the respective energy principle is considered. General criteria for stability are obtained, along with comparisons between the three different approaches.
Adaptive reconnection-based arbitrary Lagrangian Eulerian method
Bo, Wurigen; Shashkov, Mikhail
2015-07-21
We present a new adaptive Arbitrary Lagrangian Eulerian (ALE) method. This method is based on the reconnection-based ALE (ReALE) methodology of Refs. [35], [34] and [6]. The main elements in a standard ReALE method are: an explicit Lagrangian phase on an arbitrary polygonal (in 2D) mesh in which the solution and positions of grid nodes are updated; a rezoning phase in which a new grid is defined by changing the connectivity (using Voronoi tessellation) but not the number of cells; and a remapping phase in which the Lagrangian solution is transferred onto the new grid. Furthermore, in the standard ReALE method, the rezoned mesh is smoothed by using one or several steps toward centroidal Voronoi tessellation, but it is not adapted to the solution in any way.
Adaptive reconnection-based arbitrary Lagrangian Eulerian method
Bo, Wurigen; Shashkov, Mikhail
2015-07-21
We present a new adaptive Arbitrary Lagrangian Eulerian (ALE) method. This method is based on the reconnection-based ALE (ReALE) methodology of Refs. [35], [34] and [6]. The main elements in a standard ReALE method are: an explicit Lagrangian phase on an arbitrary polygonal (in 2D) mesh in which the solution and positions of grid nodes are updated; a rezoning phase in which a new grid is defined by changing the connectivity (using Voronoi tessellation) but not the number of cells; and a remapping phase in which the Lagrangian solution is transferred onto the new grid. Furthermore, in the standard ReALEmore » method, the rezoned mesh is smoothed by using one or several steps toward centroidal Voronoi tessellation, but it is not adapted to the solution in any way.« less
An implicit-explicit Eulerian Godunov scheme for compressible flow
Collins, J.P.; Colella, P.; Glaz, H.M.
1995-02-01
A hybrid implicit-explicit scheme is developed for Eulerian hydrodynamics. The hybridization is a continuous switch and operates on each characteristic field separately. The explicit scheme is a version of the second-order Geodunov scheme; the implicit method is only first-order accurate in time but leads to a block tridiagonal matrix inversion for efficiency and is unconditionally stable for the case of linear advection. The methodology is described for the cases of linear advection, for nonlinear scalar problems, and for gas dynamics. An important element of our work is the use of a modified Engquist-Osher flux function in place of the Godunov flux. Several numerical results are presented to demonstrate the properties of the method, especially stable numerical shocks at very high CFL numbers and second-order accurate steady states. 24 refs., 6 figs., 2 tabs.
Full Eulerian lattice Boltzmann model for conjugate heat transfer.
Hu, Yang; Li, Decai; Shu, Shi; Niu, Xiaodong
2015-12-01
In this paper a full Eulerian lattice Boltzmann model is proposed for conjugate heat transfer. A unified governing equation with a source term for the temperature field is derived. By introducing the source term, we prove that the continuity of temperature and its normal flux at the interface is satisfied automatically. The curved interface is assumed to be zigzag lines. All physical quantities are recorded and updated on a Cartesian grid. As a result, any complicated treatment near the interface is avoided, which makes the proposed model suitable to simulate the conjugate heat transfer with complex interfaces efficiently. The present conjugate interface treatment is validated by several steady and unsteady numerical tests, including pure heat conduction, forced convection, and natural convection problems. Both flat and curved interfaces are also involved. The obtained results show good agreement with the analytical and/or finite volume results. PMID:26764851
Techniques to derive geometries for image-based Eulerian computations
Dillard, Seth; Buchholz, James; Vigmostad, Sarah; Kim, Hyunggun; Udaykumar, H.S.
2014-01-01
Purpose The performance of three frequently used level set-based segmentation methods is examined for the purpose of defining features and boundary conditions for image-based Eulerian fluid and solid mechanics models. The focus of the evaluation is to identify an approach that produces the best geometric representation from a computational fluid/solid modeling point of view. In particular, extraction of geometries from a wide variety of imaging modalities and noise intensities, to supply to an immersed boundary approach, is targeted. Design/methodology/approach Two- and three-dimensional images, acquired from optical, X-ray CT, and ultrasound imaging modalities, are segmented with active contours, k-means, and adaptive clustering methods. Segmentation contours are converted to level sets and smoothed as necessary for use in fluid/solid simulations. Results produced by the three approaches are compared visually and with contrast ratio, signal-to-noise ratio, and contrast-to-noise ratio measures. Findings While the active contours method possesses built-in smoothing and regularization and produces continuous contours, the clustering methods (k-means and adaptive clustering) produce discrete (pixelated) contours that require smoothing using speckle-reducing anisotropic diffusion (SRAD). Thus, for images with high contrast and low to moderate noise, active contours are generally preferable. However, adaptive clustering is found to be far superior to the other two methods for images possessing high levels of noise and global intensity variations, due to its more sophisticated use of local pixel/voxel intensity statistics. Originality/value It is often difficult to know a priori which segmentation will perform best for a given image type, particularly when geometric modeling is the ultimate goal. This work offers insight to the algorithm selection process, as well as outlining a practical framework for generating useful geometric surfaces in an Eulerian setting. PMID
Development and deployment of constitutive softening routines in Eulerian hydrocodes.
Fuller, Timothy Jesse; Dewers, Thomas A.; Swan, Matthew Scot
2013-03-01
The state of the art in failure modeling enables assessment of crack nucleation, propagation, and progression to fragmentation due to high velocity impact. Vulnerability assessments suggest a need to track material behavior through failure, to the point of fragmentation and beyond. This eld of research is particularly challenging for structures made of porous quasi-brittle materials, such as ceramics used in modern armor systems, due to the complex material response when loading exceeds the quasi-brittle material's elastic limit. Further complications arise when incorporating the quasi-brittle material response in multi-material Eulerian hydrocode simulations. In this report, recent e orts in coupling a ceramic materials response in the post-failure regime with an Eulerian hydro code are described. Material behavior is modeled by the Kayenta material model [2]and Alegra as the host nite element code [14]. Kayenta, a three invariant phenomenological plasticity model originally developed for modeling the stress response of geologic materials, has in recent years been used with some success in the modeling of ceramic and other quasi-brittle materials to high velocity impact. Due to the granular nature of ceramic materials, Kayenta allows for signi cant pressures to develop due to dilatant plastic ow, even in shear dominated loading where traditional equations of state predict little or no pressure response. When a material's ability to carry further load is compromised, Kayenta allows the material's strength and sti ness to progressively degrade through the evolution of damage to the point of material failure. As material dilatation and damage progress, accommodations are made within Alegra to treat in a consistent manner the evolving state.
New iterative solvers for the NAG Libraries
Salvini, S.; Shaw, G.
1996-12-31
The purpose of this paper is to introduce the work which has been carried out at NAG Ltd to update the iterative solvers for sparse systems of linear equations, both symmetric and unsymmetric, in the NAG Fortran 77 Library. Our current plans to extend this work and include it in our other numerical libraries in our range are also briefly mentioned. We have added to the Library the new Chapter F11, entirely dedicated to sparse linear algebra. At Mark 17, the F11 Chapter includes sparse iterative solvers, preconditioners, utilities and black-box routines for sparse symmetric (both positive-definite and indefinite) linear systems. Mark 18 will add solvers, preconditioners, utilities and black-boxes for sparse unsymmetric systems: the development of these has already been completed.
Using SPARK as a Solver for Modelica
Wetter, Michael; Wetter, Michael; Haves, Philip; Moshier, Michael A.; Sowell, Edward F.
2008-06-30
Modelica is an object-oriented acausal modeling language that is well positioned to become a de-facto standard for expressing models of complex physical systems. To simulate a model expressed in Modelica, it needs to be translated into executable code. For generating run-time efficient code, such a translation needs to employ algebraic formula manipulations. As the SPARK solver has been shown to be competitive for generating such code but currently cannot be used with the Modelica language, we report in this paper how SPARK's symbolic and numerical algorithms can be implemented in OpenModelica, an open-source implementation of a Modelica modeling and simulation environment. We also report benchmark results that show that for our air flow network simulation benchmark, the SPARK solver is competitive with Dymola, which is believed to provide the best solver for Modelica.
Multigrid in energy preconditioner for Krylov solvers
Slaybaugh, R.N.; Evans, T.M.; Davidson, G.G.; Wilson, P.P.H.
2013-06-01
We have added a new multigrid in energy (MGE) preconditioner to the Denovo discrete-ordinates radiation transport code. This preconditioner takes advantage of a new multilevel parallel decomposition. A multigroup Krylov subspace iterative solver that is decomposed in energy as well as space-angle forms the backbone of the transport solves in Denovo. The space-angle-energy decomposition facilitates scaling to hundreds of thousands of cores. The multigrid in energy preconditioner scales well in the energy dimension and significantly reduces the number of Krylov iterations required for convergence. This preconditioner is well-suited for use with advanced eigenvalue solvers such as Rayleigh Quotient Iteration and Arnoldi.
ODE System Solver W. Krylov Iteration & Rootfinding
Hindmarsh, Alan C.
1991-09-09
LSODKR is a new initial value ODE solver for stiff and nonstiff systems. It is a variant of the LSODPK and LSODE solvers, intended mainly for large stiff systems. The main differences between LSODKR and LSODE are the following: (a) for stiff systems, LSODKR uses a corrector iteration composed of Newton iteration and one of four preconditioned Krylov subspace iteration methods. The user must supply routines for the preconditioning operations, (b) Within the corrector iteration, LSODKR does automatic switching between functional (fixpoint) iteration and modified Newton iteration, (c) LSODKR includes the ability to find roots of given functions of the solution during the integration.
ODE System Solver W. Krylov Iteration & Rootfinding
Energy Science and Technology Software Center (ESTSC)
1991-09-09
LSODKR is a new initial value ODE solver for stiff and nonstiff systems. It is a variant of the LSODPK and LSODE solvers, intended mainly for large stiff systems. The main differences between LSODKR and LSODE are the following: (a) for stiff systems, LSODKR uses a corrector iteration composed of Newton iteration and one of four preconditioned Krylov subspace iteration methods. The user must supply routines for the preconditioning operations, (b) Within the corrector iteration,more » LSODKR does automatic switching between functional (fixpoint) iteration and modified Newton iteration, (c) LSODKR includes the ability to find roots of given functions of the solution during the integration.« less
Steady potential solver for unsteady aerodynamic analyses
NASA Technical Reports Server (NTRS)
Hoyniak, Dan
1994-01-01
Development of a steady flow solver for use with LINFLO was the objective of this report. The solver must be compatible with LINFLO, be composed of composite mesh, and have transonic capability. The approaches used were: (1) steady flow potential equations written in nonconservative form; (2) Newton's Method; (3) implicit, least-squares, interpolation method to obtain finite difference equations; and (4) matrix inversion routines from LINFLO. This report was given during the NASA LeRC Workshop on Forced Response in Turbomachinery in August of 1993.
Wave Speeds, Riemann Solvers and Artificial Viscosity
Rider, W.J.
1999-07-18
A common perspective on the numerical solution of the equation Euler equations for shock physics is examined. The common viewpoint is based upon the selection of nonlinear wavespeeds upon which the dissipation (implicit or explicit) is founded. This perspective shows commonality between Riemann solver based method (i.e. Godunov-type) and artificial viscosity (i.e. von Neumann-Richtmyer). As an example we derive an improved nonlinear viscous stabilization of a Richtmyer-Lax-Wendroff method. Additionally, we will define a form of classical artificial viscosity based upon the HLL Riemann solver.
Ames, Thomas L.; Farnsworth, Grant V.; Ketcheson, David Isaac; Robinson, Allen Conrad
2009-09-01
The modeling of solids is most naturally placed within a Lagrangian framework because it requires constitutive models which depend on knowledge of the original material orientations and subsequent deformations. Detailed kinematic information is needed to ensure material frame indifference which is captured through the deformation gradient F. Such information can be tracked easily in a Lagrangian code. Unfortunately, not all problems can be easily modeled using Lagrangian concepts due to severe distortions in the underlying motion. Either a Lagrangian/Eulerian or a pure Eulerian modeling framework must be introduced. We discuss and contrast several Lagrangian/Eulerian approaches for keeping track of the details of material kinematics.
Wang, C.Y.; Zeuch, W.R.
1982-01-01
This paper describes an arbitrary Lagrangian-Eulerian method for analyzing fluid-structure interactions in fast-reactor containment with complex internal structures. The fluid transient can be calculated either implicitly or explicitly, using a finite-difference mesh with vertices that may be moved with the fluid (Lagrangian), held fixed (Eulerian), or moved in any other prescribed manner (hybrid Lagrangian Eulerian). The structural response is computed explicitly by two nonlinear, elastic-plastic finite-element modules formulated in corotational coordinates. Interaction between fluid and structure is accounted for by enforcing the interface boundary conditions. The method has convincing advantages in treating complicated phenomena such as flow through perforated structures, large material distortions, flow around corners and irregularities, and highly contorted fluid boundaries. Several sample problems are given to illustrate the effectiveness of this arbitrary Lagrangian-Eulerian method.
EULERIAN-LAGRANGIAN LOCALIZED ADJOINT METHOD FOR THE ADVECTION-DIFFUSION EQUATION
Many numerical methods use characteristic analysis to accommodate the advective component of transport. uch characteristic methods include Eulerian-Lagrangian methods (ELM), modified method of characteristics (MMOC), and operator splitting methods. eneralization of characteristic...
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.
AN EULERIAN-LAGRANGIAN LOCALIZED ADJOINT METHOD FOR THE ADVECTION-DIFFUSION EQUATION
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...
Leung Shingyu; Qian Jianliang
2010-11-20
We propose the backward phase flow method to implement the Fourier-Bros-Iagolnitzer (FBI)-transform-based Eulerian Gaussian beam method for solving the Schroedinger equation in the semi-classical regime. The idea of Eulerian Gaussian beams has been first proposed in . In this paper we aim at two crucial computational issues of the Eulerian Gaussian beam method: how to carry out long-time beam propagation and how to compute beam ingredients rapidly in phase space. By virtue of the FBI transform, we address the first issue by introducing the reinitialization strategy into the Eulerian Gaussian beam framework. Essentially we reinitialize beam propagation by applying the FBI transform to wavefields at intermediate time steps when the beams become too wide. To address the second issue, inspired by the original phase flow method, we propose the backward phase flow method which allows us to compute beam ingredients rapidly. Numerical examples demonstrate the efficiency and accuracy of the proposed algorithms.
Eulerian simulations of collisional effects on electrostatic plasma waves
Pezzi, Oreste; Valentini, Francesco; Perrone, Denise; Veltri, Pierluigi
2013-09-15
The problem of collisions in a plasma is a wide subject with a huge historical literature. In fact, the description of realistic plasmas is a tough problem to attack, both from the theoretical and the numerical point of view. In this paper, a Eulerian time-splitting algorithm for the study of the propagation of electrostatic waves in collisional plasmas is presented. Collisions are modeled through one-dimensional operators of the Fokker-Planck type, both in linear and nonlinear forms. The accuracy of the numerical code is discussed by comparing the numerical results to the analytical predictions obtained in some limit cases when trying to evaluate the effects of collisions in the phenomenon of wave plasma echo and collisional dissipation of Bernstein-Greene-Kruskal waves. Particular attention is devoted to the study of the nonlinear Dougherty collisional operator, recently used to describe the collisional dissipation of electron plasma waves in a pure electron plasma column [M. W. Anderson and T. M. O'Neil, Phys. Plasmas 14, 112110 (2007)]. Finally, for the study of collisional plasmas, a recipe to set the simulation parameters in order to prevent the filamentation problem can be provided, by exploiting the property of velocity diffusion operators to smooth out small velocity scales.
Complete Eulerian-mean tracer equation for coarse resolution OGCMs
NASA Astrophysics Data System (ADS)
Dubovikov, M. S.; Canuto, V. M.
2006-06-01
McDougall and McIntosh showed that the adiabatic mesoscale mixing is represented incompletely in the tracer Eulerian-averaged equation (EAE) of coarse resolution OGCMs. We show that completing EAE requires an adequate decomposition of the mesoscale tracer flux which is achieved by means of transforming mesoscale fields to isopycnal coordinates (IC) where mesoscale dynamics has the simplest form. The transformation results in splitting Fτ into two components and : the former is determined by buoyancy mesoscale dynamics only and has a trivial kinematic dependence on the mean tracer field, the latter is determined by mesoscale tracer dynamics. Thus, the problem of modelling (parameterizing) Fτ in ZC is divided in two stages which can be termed kinematic and dynamic. The kinematic stage consists in adequate decomposing Fτ, and the result is expressed in terms of mesoscale fields. The dynamic stage consists in applying a specific dynamic mesoscale model to parameterize the components of Fτ. In this article, we show that some components of Fτ are missing in ZC-OGCMs tracer equation and that their contribution is of the same order of magnitude as the mesoscale contribution itself. We also show that Fτ has components across mean isopycnals and that their existence is consistent with the adiabatic approximation which requires vanishing all fluxes across isopycnal surfaces. As for practical results, we derive the complete equation for the large scale tracer in ZC-OGCMs and present the parameterization of the terms which have been missing thus far.
Lagrangian, Eulerian, and Dynamically Accessible Stability of MHD flows
NASA Astrophysics Data System (ADS)
Andreussi, Tommaso; Morrison, Philip; Pegoraro, Francesco
2012-10-01
Stability conditions of magnetized plasma flows are obtained by exploiting the Hamiltonian structure of the magnetohydrodynamics (MHD) equations and, in particular, by using three kinds of energy principles. First, the Lagrangian energy principle of Ref. [1] is introduced and sufficient stability conditions are presented. Next, plasma flows are described in terms of Eulerian variables and the noncanonical Hamiltonian formulation of MHD [2] is exploited. For symmetric equilibria, the energy-Casimir principle of Ref. [3] is expanded to second order and sufficient conditions for stability to symmetric perturbation are obtained. Then, dynamically accessible variations, i.e. variations that explicitly preserve the invariants of the system, are introduced and the respective energy principle is considered. As in Ref. [4], general criteria for stability are obtained. A comparison between the three different approaches is finally presented. [4pt] [1] E.A. Frieman and M. Rotenberg, Rev. Mod. Phys., 32 898 (1960).[0pt] [2] P.J. Morrison, J.M. Greene, Phys. Rev. Lett., 45 790 (1980).[0pt] [3] T. Andreussi, P.J. Morrison, F. Pegoraro, Phys. Plasmas, 19 052102 (2012).[0pt] [4] E. Hameiri, Phys. Plasmas, 10 2643 (2003).
Eulerian simulations of collisional effects on electrostatic plasma waves
NASA Astrophysics Data System (ADS)
Pezzi, Oreste; Valentini, Francesco; Perrone, Denise; Veltri, Pierluigi
2013-09-01
The problem of collisions in a plasma is a wide subject with a huge historical literature. In fact, the description of realistic plasmas is a tough problem to attack, both from the theoretical and the numerical point of view. In this paper, a Eulerian time-splitting algorithm for the study of the propagation of electrostatic waves in collisional plasmas is presented. Collisions are modeled through one-dimensional operators of the Fokker-Planck type, both in linear and nonlinear forms. The accuracy of the numerical code is discussed by comparing the numerical results to the analytical predictions obtained in some limit cases when trying to evaluate the effects of collisions in the phenomenon of wave plasma echo and collisional dissipation of Bernstein-Greene-Kruskal waves. Particular attention is devoted to the study of the nonlinear Dougherty collisional operator, recently used to describe the collisional dissipation of electron plasma waves in a pure electron plasma column [M. W. Anderson and T. M. O'Neil, Phys. Plasmas 14, 112110 (2007)]. Finally, for the study of collisional plasmas, a recipe to set the simulation parameters in order to prevent the filamentation problem can be provided, by exploiting the property of velocity diffusion operators to smooth out small velocity scales.
Characterization of Pacific Ocean Surface Temperatures Using Eulerian Motion Magnification
NASA Astrophysics Data System (ADS)
Rojo Hernandez, J. D.; Mesa, O. J.
2014-12-01
The Eulerian Motion Magnification Method was used in order to identify the spatial-temporal patterns in the variability of sea-surface temperatures (SST) in the Pacific Ocean. This method, developed by a research team at MIT, consists in jointly applying spatial and temporal filters to a sequence of images with a known playback speed, and then amplifying the intensity of a signal associated with a certain frequency, so that periodic phenomena can be easily displayed. Magnifying the SST in the frequency band of 2-7 years - which corresponds to ENSO- various processes can be clearly observed, such as the dynamics of temperature variability in the Pacific Ocean associated with the occurrence of warm and cool episodes of the differentiated ocean warming type (Central-Pacific El Nino and Eastern-Pacific El Nino), the possible interaction between tropical and extra-tropical waves that may enhance or diminish the possible ENSO events, and it displays that the ocean heating and/or cooling patterns can be represented as Kelvin and Rosby wave propagation at inter-annual scale.
Lagrangian and Eulerian description of bed-load particle kinematics
NASA Astrophysics Data System (ADS)
Ballio, Francesco; Sadabadi, Seyed Abbas Hosseini; Pokrajac, Dubravka; Radice, Alessio
2016-04-01
The motion of bed-load sediment particles transported by a flow can be analyzed within a Lagrangian or an Eulerian framework. In the former case, we consider the particles as individual objects in motion and we study their kinematic properties. The latter approach is instead referred to suitably chosen control volumes. Quantities describing sediment motion in the two frameworks are different, and the relationships among the two approaches are not straightforward. In this work, we intend to discuss the kinematic properties of sediment transport: first, a set of quantities is univocally defined; then, relationships among different representations are explored. Proof-of-concept results presented in the study are from a recent experiment involving weak bed-load sediment transport, where the moving particles were released over a fixed rough bed. The bulk flow velocity was 1.4 times the critical value for incipient particle motion, and particles were mostly moving by rolling and sliding, with limited saltation. The particle motion was filmed from the top and the measurements were conducted by image-based methods, obtaining extensive samples of virtually-instantaneous quantities.
Ising-like phase transition of an n-component Eulerian face-cubic model
NASA Astrophysics Data System (ADS)
Ding, Chengxiang; Guo, Wenan; Deng, Youjin
2013-11-01
By means of Monte Carlo simulations and a finite-size scaling analysis, we find a critical line of an n-component Eulerian face-cubic model on the square lattice and the simple cubic lattice in the region v>1, where v is the bond weight. The phase transition belongs to the Ising universality class independent of n. The critical properties of the phase transition can also be captured by the percolation of the complement of the Eulerian graph.
Lagrangian and Eulerian Statistics of Vorticity Dynamics in Turbulent Stratified Shear Flows
NASA Astrophysics Data System (ADS)
Jacobitz, Frank; Schneider, Kai; Farge, Marie
2015-11-01
The Lagrangian and Eulerian time-rate of change statistics of vorticity in homogeneous turbulence with shear and stable stratification are studied. Direct numerical simulations are performed, in which the Richardson number is varied from Ri=0, corresponding to unstratified shear flow, to Ri=1, corresponding to strongly stratified shear flow. The probability density functions (pdfs) of both Lagrangian and Eulerian time-rates of change show a strong influence on the Richardson number. The Lagrangian time-rate of change pdf has a stretched-exponential shape due to the vortex stretching and tilting term in the equation for fluctuating vorticity. The shape of the Eulerian time-rate of change pdf was also observed to be stretched-exponential and the extreme values for the Eulerian time-rate of change are larger than those observed for the Lagrangian counterpart due to the nonlinear term in the vorticity equation. The Lagrangian and Eulerian acceleration pdfs are mainly determined by the pressure-gradient and nonlinear terms in the Navier-Stokes equation, respectively. The Lagrangian time-rate of change pdf of fluctuating density does not show a stretched exponential shape, while its Eulerian counterpart does due to the nonlinear term in the in the density advection-diffusion equation.
Aleph Field Solver Challenge Problem Results Summary.
Hooper, Russell; Moore, Stan Gerald
2015-01-01
Aleph models continuum electrostatic and steady and transient thermal fields using a finite-element method. Much work has gone into expanding the core solver capability to support enriched mod- eling consisting of multiple interacting fields, special boundary conditions and two-way interfacial coupling with particles modeled using Aleph's complementary particle-in-cell capability. This report provides quantitative evidence for correct implementation of Aleph's field solver via order- of-convergence assessments on a collection of problems of increasing complexity. It is intended to provide Aleph with a pedigree and to establish a basis for confidence in results for more challeng- ing problems important to Sandia's mission that Aleph was specifically designed to address.
Verifying a Local Generic Solver in Coq
NASA Astrophysics Data System (ADS)
Hofmann, Martin; Karbyshev, Aleksandr; Seidl, Helmut
Fixpoint engines are the core components of program analysis tools and compilers. If these tools are to be trusted, special attention should be paid also to the correctness of such solvers. In this paper we consider the local generic fixpoint solver RLD which can be applied to constraint systems {x}sqsupseteq fx,{x}in V, over some lattice {D} where the right-hand sides f x are given as arbitrary functions implemented in some specification language. The verification of this algorithm is challenging, because it uses higher-order functions and relies on side effects to track variable dependences as they are encountered dynamically during fixpoint iterations. Here, we present a correctness proof of this algorithm which has been formalized by means of the interactive proof assistant Coq.
Domain decomposition for the SPN solver MINOS
Jamelot, Erell; Baudron, Anne-Marie; Lautard, Jean-Jacques
2012-07-01
In this article we present a domain decomposition method for the mixed SPN equations, discretized with Raviart-Thomas-Nedelec finite elements. This domain decomposition is based on the iterative Schwarz algorithm with Robin interface conditions to handle communications. After having described this method, we give details on how to optimize the convergence. Finally, we give some numerical results computed in a realistic 3D domain. The computations are done with the MINOS solver of the APOLLO3 (R) code. (authors)
A perspective on unstructured grid flow solvers
NASA Technical Reports Server (NTRS)
Venkatakrishnan, V.
1995-01-01
This survey paper assesses the status of compressible Euler and Navier-Stokes solvers on unstructured grids. Different spatial and temporal discretization options for steady and unsteady flows are discussed. The integration of these components into an overall framework to solve practical problems is addressed. Issues such as grid adaptation, higher order methods, hybrid discretizations and parallel computing are briefly discussed. Finally, some outstanding issues and future research directions are presented.
Differential geometry based solvation model I: Eulerian formulation.
Chen, Zhan; Baker, Nathan A; Wei, G W
2010-11-01
This paper presents a differential geometry based model for the analysis and computation of the equilibrium property of solvation. Differential geometry theory of surfaces is utilized to define and construct smooth interfaces with good stability and differentiability for use in characterizing the solvent-solute boundaries and in generating continuous dielectric functions across the computational domain. A total free energy functional is constructed to couple polar and nonpolar contributions to the salvation process. Geometric measure theory is employed to rigorously convert a Lagrangian formulation of the surface energy into an Eulerian formulation so as to bring all energy terms into an equal footing. By minimizing the total free energy functional, we derive coupled generalized Poisson-Boltzmann equation (GPBE) and generalized geometric flow equation (GGFE) for the electrostatic potential and the construction of realistic solvent-solute boundaries, respectively. By solving the coupled GPBE and GGFE, we obtain the electrostatic potential, the solvent-solute boundary profile, and the smooth dielectric function, and thereby improve the accuracy and stability of implicit solvation calculations. We also design efficient second order numerical schemes for the solution of the GPBE and GGFE. Matrix resulted from the discretization of the GPBE is accelerated with appropriate preconditioners. An alternative direct implicit (ADI) scheme is designed to improve the stability of solving the GGFE. Two iterative approaches are designed to solve the coupled system of nonlinear partial differential equations. Extensive numerical experiments are designed to validate the present theoretical model, test computational methods, and optimize numerical algorithms. Example solvation analysis of both small compounds and proteins are carried out to further demonstrate the accuracy, stability, efficiency and robustness of the present new model and numerical approaches. Comparison is given to
MILAMIN 2 - Fast MATLAB FEM solver
NASA Astrophysics Data System (ADS)
Dabrowski, Marcin; Krotkiewski, Marcin; Schmid, Daniel W.
2013-04-01
MILAMIN is a free and efficient MATLAB-based two-dimensional FEM solver utilizing unstructured meshes [Dabrowski et al., G-cubed (2008)]. The code consists of steady-state thermal diffusion and incompressible Stokes flow solvers implemented in approximately 200 lines of native MATLAB code. The brevity makes the code easily customizable. An important quality of MILAMIN is speed - it can handle millions of nodes within minutes on one CPU core of a standard desktop computer, and is faster than many commercial solutions. The new MILAMIN 2 allows three-dimensional modeling. It is designed as a set of functional modules that can be used as building blocks for efficient FEM simulations using MATLAB. The utilities are largely implemented as native MATLAB functions. For performance critical parts we use MUTILS - a suite of compiled MEX functions optimized for shared memory multi-core computers. The most important features of MILAMIN 2 are: 1. Modular approach to defining, tracking, and discretizing the geometry of the model 2. Interfaces to external mesh generators (e.g., Triangle, Fade2d, T3D) and mesh utilities (e.g., element type conversion, fast point location, boundary extraction) 3. Efficient computation of the stiffness matrix for a wide range of element types, anisotropic materials and three-dimensional problems 4. Fast global matrix assembly using a dedicated MEX function 5. Automatic integration rules 6. Flexible prescription (spatial, temporal, and field functions) and efficient application of Dirichlet, Neuman, and periodic boundary conditions 7. Treatment of transient and non-linear problems 8. Various iterative and multi-level solution strategies 9. Post-processing tools (e.g., numerical integration) 10. Visualization primitives using MATLAB, and VTK export functions We provide a large number of examples that show how to implement a custom FEM solver using the MILAMIN 2 framework. The examples are MATLAB scripts of increasing complexity that address a given
On Lagrangian and Eulerian Acceleration in Rotating and Sheared Homogeneous Turbulence
NASA Astrophysics Data System (ADS)
Jacobitz, Frank; Schneider, Kai; Bos, Wouter; Farge, Marie
2013-11-01
The Lagrangian and Eulerian acceleration properties of turbulence are of importance for problems ranging from fundamental theoretical considerations to modeling of dispersion processes. The acceleration statistics of rotating and sheared homogeneous turbulence are studied here using direct numerical simulations. The study focusses in particular on the influence of the Coriolis to shear rate ratio and also on the scale dependence of the statistics. The probability density functions (pdfs) of both Lagrangian and Eulerian acceleration show a strong and similar influence on the rotation ratio. The flatness further quantifies this influence and yields values close to three for strong rotation. For moderate and vanishing rotation, the flatness of the Eulerian acceleration is larger than that of the Lagrangian acceleration, contrary to previous results for isotropic turbulence. A wavelet-based scale-dependent analysis shows that the flatness of both Eulerian and Lagrangian acceleration increases as scale decreases. For strong rotation, the Eulerian acceleration is more intermittent than the Lagrangian acceleration, while the opposite result is obtained for moderate rotation.
NASA Astrophysics Data System (ADS)
Speetjens, M. F. M.; Lauret, M.; Nijmeijer, H.; Anderson, P. D.
2013-05-01
Transport of passive tracers may be described through the spatio-temporal evolution of Eulerian concentration distributions or via the geometrical composition of the Lagrangian flow structure. The present study seeks to deepen insight into the connection between the Eulerian and Lagrangian perspectives by investigating the role of Lagrangian coherent structures (LCSs) in the Eulerian concentration distributions in time-periodic and spatially-periodic mixing flows. Representation of the Eulerian transport by the mapping method, describing realistic transport problems by distribution matrices, admits a generic analysis based on matrix and graph theory. This reveals that LCSs-and the transport barriers that separate them-leave a distinct “footprint” in the eigenmode spectrum of the distribution matrix and, by proxy, of the underlying Eulerian transport operator. Transport barriers impart a block-diagonal structure upon the mapping matrix, where each block matrix A corresponds with a given LCS. Its kind is reflected in the spectrum of A; higher-order periodicity yields a distinct permutation within A. The composition of the distribution matrix versus the Lagrangian flow structure thus predicted is demonstrated by way of examples. These findings increase fundamental understanding of transport phenomena and have great practical potential for e.g. flow and mixing control.
Approximate Riemann solvers for the Godunov SPH (GSPH)
NASA Astrophysics Data System (ADS)
Puri, Kunal; Ramachandran, Prabhu
2014-08-01
The Godunov Smoothed Particle Hydrodynamics (GSPH) method is coupled with non-iterative, approximate Riemann solvers for solutions to the compressible Euler equations. The use of approximate solvers avoids the expensive solution of the non-linear Riemann problem for every interacting particle pair, as required by GSPH. In addition, we establish an equivalence between the dissipative terms of GSPH and the signal based SPH artificial viscosity, under the restriction of a class of approximate Riemann solvers. This equivalence is used to explain the anomalous “wall heating” experienced by GSPH and we provide some suggestions to overcome it. Numerical tests in one and two dimensions are used to validate the proposed Riemann solvers. A general SPH pairing instability is observed for two-dimensional problems when using unequal mass particles. In general, Ducowicz Roe's and HLLC approximate Riemann solvers are found to be suitable replacements for the iterative Riemann solver in the original GSPH scheme.
New Multigrid Solver Advances in TOPS
Falgout, R D; Brannick, J; Brezina, M; Manteuffel, T; McCormick, S
2005-06-27
In this paper, we highlight new multigrid solver advances in the Terascale Optimal PDE Simulations (TOPS) project in the Scientific Discovery Through Advanced Computing (SciDAC) program. We discuss two new algebraic multigrid (AMG) developments in TOPS: the adaptive smoothed aggregation method ({alpha}SA) and a coarse-grid selection algorithm based on compatible relaxation (CR). The {alpha}SA method is showing promising results in initial studies for Quantum Chromodynamics (QCD) applications. The CR method has the potential to greatly improve the applicability of AMG.
DPS--a computerised diagnostic problem solver.
Bartos, P; Gyárfas, F; Popper, M
1982-01-01
The paper contains a short description of the DPS system which is a computerized diagnostic problem solver. The system is under development of the Research Institute of Medical Bionics in Bratislava, Czechoslovakia. Its underlying philosophy yields from viewing the diagnostic process as process of cognitive problem solving. The implementation of the system is based on the methods of Artificial Intelligence and utilisation of production systems and frame theory should be noted in this context. Finally a list of program modules and their characterisation is presented. PMID:6811229
Updates to the NEQAIR Radiation Solver
NASA Technical Reports Server (NTRS)
Cruden, Brett A.; Brandis, Aaron M.
2014-01-01
The NEQAIR code is one of the original heritage solvers for radiative heating prediction in aerothermal environments, and is still used today for mission design purposes. This paper discusses the implementation of the first major revision to the NEQAIR code in the last five years, NEQAIR v14.0. The most notable features of NEQAIR v14.0 are the parallelization of the radiation computation, reducing runtimes by about 30×, and the inclusion of mid-wave CO2 infrared radiation.
Input-output-controlled nonlinear equation solvers
NASA Technical Reports Server (NTRS)
Padovan, Joseph
1988-01-01
To upgrade the efficiency and stability of the successive substitution (SS) and Newton-Raphson (NR) schemes, the concept of input-output-controlled solvers (IOCS) is introduced. By employing the formal properties of the constrained version of the SS and NR schemes, the IOCS algorithm can handle indefiniteness of the system Jacobian, can maintain iterate monotonicity, and provide for separate control of load incrementation and iterate excursions, as well as having other features. To illustrate the algorithmic properties, the results for several benchmark examples are presented. These define the associated numerical efficiency and stability of the IOCS.
NASA Astrophysics Data System (ADS)
Wang, Y.; Shu, C.; Yang, L. M.
2016-02-01
A boundary condition-enforced-immersed boundary-lattice Boltzmann flux solver is proposed in this work for effective simulation of thermal flows with Neumann boundary conditions. In this method, two auxiliary layers of Lagrangian points are introduced and respectively placed inside and outside of the solid body, on which the temperature corrections (related to the heat source) are set as unknowns. To effectively consider the fluid-boundary interaction, these unknowns are expressed as algebraic summations of the temperature correction on Eulerian points, which are in turn obtained from biased distributions of unknown temperature corrections on the immersed boundary. By enforcing the temperature gradient at the solid boundary being equal to that approximated by the corrected temperature field, a set of algebraic equations are formed and solved to obtain all the unknowns simultaneously. They are then distributed biasedly to the inner region of the auxiliary layer so that the diffusion from the smooth delta function can be reduced substantially. In addition, the solutions of the flow and temperature fields are obtained by the thermal lattice Boltzmann flux solver with the second order of accuracy. The proposed method is well validated through its applications to simulate several benchmarks of natural, forced and mixed convection problems. It has been demonstrated that the present solver has about 1.724 order of accuracy and the error between the present result and theoretical value for the temperature gradient on the solid surface is in the order of 10-13, which indicates that the proposed method is able to satisfy the Neumann boundary condition accurately.
Comparison of electromagnetic solvers for antennas mounted on vehicles
NASA Astrophysics Data System (ADS)
Mocker, M. S. L.; Hipp, S.; Spinnler, F.; Tazi, H.; Eibert, T. F.
2015-11-01
An electromagnetic solver comparison for various use cases of antennas mounted on vehicles is presented. For this purpose, several modeling approaches, called transient, frequency and integral solver, including the features fast resonant method and autoregressive filter, offered by CST MWS, are investigated. The solvers and methods are compared for a roof antenna itself, a simplified vehicle, a roof including a panorama window and a combination of antenna and vehicle. With these examples, the influence of different materials, data formats and parameters such as size and complexity are investigated. Also, the necessary configurations for the mesh and the solvers are described.
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.
Simulation on gasification of forestry residues in fluidized beds by Eulerian-Lagrangian approach.
Xie, Jun; Zhong, Wenqi; Jin, Baosheng; Shao, Yingjuan; Liu, Hao
2012-10-01
A comprehensive three-dimensional numerical model is developed to simulate forestry residues gasification in a fluidized bed reactor using Eulerian-Lagrangian approach. The complex granular flow behaviors and chemical reaction characteristics are addressed simultaneously. The model uses an Eulerian method for fluid phase and a discrete particle method for solid phase, which takes particle contact force into account. Heterogeneous and homogenous reaction rates are solved on the Eulerian grid. The numerical model is employed to study the gasification performance in a lab-scale pine gasifier. A series of simulations have been performed with some critical parameters including temperature, equivalence ratio and steam to biomass ratio. The model predicts product gas composition and carbon conversion efficiency in good agreement with experimental data. The formation and development of flow regimes, profiles of particle species, and distributions of gas compositions inside the reactor are also discussed. PMID:22858466
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.
Algebraic Multiscale Solver for Elastic Geomechanical Deformation
NASA Astrophysics Data System (ADS)
Castelletto, N.; Hajibeygi, H.; Tchelepi, H.
2015-12-01
Predicting the geomechanical response of geological formations to thermal, pressure, and mechanical loading is important in many engineering applications. The mathematical formulation that describes deformation of a reservoir coupled with flow and transport entails heterogeneous coefficients with a wide range of length scales. Such detailed heterogeneous descriptions of reservoir properties impose severe computational challenges for the study of realistic-scale (km) reservoirs. To deal with these challenges, we developed an Algebraic Multiscale Solver for ELastic geomechanical deformation (EL-AMS). Constructed on finite element fine-scale system, EL-AMS imposes a coarse-scale grid, which is a non-overlapping decomposition of the domain. Then, local (coarse) basis functions for the displacement vector are introduced. These basis functions honor the elastic properties of the local domains subject to the imposed local boundary conditions. The basis form the Restriction and Prolongation operators. These operators allow for the construction of accurate coarse-scale systems for the displacement. While the multiscale system is efficient for resolving low-frequency errors, coupling it with a fine-scale smoother, e.g., ILU(0), leads to an efficient iterative solver. Numerical results for several test cases illustrate that EL-AMS is quite efficient and applicable to simulate elastic deformation of large-scale heterogeneous reservoirs.
Using the scalable nonlinear equations solvers package
Gropp, W.D.; McInnes, L.C.; Smith, B.F.
1995-02-01
SNES (Scalable Nonlinear Equations Solvers) is a software package for the numerical solution of large-scale systems of nonlinear equations on both uniprocessors and parallel architectures. SNES also contains a component for the solution of unconstrained minimization problems, called SUMS (Scalable Unconstrained Minimization Solvers). Newton-like methods, which are known for their efficiency and robustness, constitute the core of the package. As part of the multilevel PETSc library, SNES incorporates many features and options from other parts of PETSc. In keeping with the spirit of the PETSc library, the nonlinear solution routines are data-structure-neutral, making them flexible and easily extensible. This users guide contains a detailed description of uniprocessor usage of SNES, with some added comments regarding multiprocessor usage. At this time the parallel version is undergoing refinement and extension, as we work toward a common interface for the uniprocessor and parallel cases. Thus, forthcoming versions of the software will contain additional features, and changes to parallel interface may result at any time. The new parallel version will employ the MPI (Message Passing Interface) standard for interprocessor communication. Since most of these details will be hidden, users will need to perform only minimal message-passing programming.
On code verification of RANS solvers
NASA Astrophysics Data System (ADS)
Eça, L.; Klaij, C. M.; Vaz, G.; Hoekstra, M.; Pereira, F. S.
2016-04-01
This article discusses Code Verification of Reynolds-Averaged Navier Stokes (RANS) solvers that rely on face based finite volume discretizations for volumes of arbitrary shape. The study includes test cases with known analytical solutions (generated with the method of manufactured solutions) corresponding to laminar and turbulent flow, with the latter using eddy-viscosity turbulence models. The procedure to perform Code Verification based on grid refinement studies is discussed and the requirements for its correct application are illustrated in a simple one-dimensional problem. It is shown that geometrically similar grids are recommended for proper Code Verification and so the data should not have scatter making the use of least square fits unnecessary. Results show that it may be advantageous to determine the extrapolated error to cell size/time step zero instead of assuming that it is zero, especially when it is hard to determine the asymptotic order of grid convergence. In the RANS examples, several of the features of the ReFRESCO solver are checked including the effects of the available turbulence models in the convergence properties of the code. It is shown that it is required to account for non-orthogonality effects in the discretization of the diffusion terms and that the turbulence quantities transport equations can deteriorate the order of grid convergence of mean flow quantities.
A coupled PFEM-Eulerian approach for the solution of porous FSI problems
NASA Astrophysics Data System (ADS)
Larese, A.; Rossi, R.; Oñate, E.; Idelsohn, S. R.
2012-12-01
This paper aims to present a coupled solution strategy for the problem of seepage through a rockfill dam taking into account the free-surface flow within the solid as well as in its vicinity. A combination of a Lagrangian model for the structural behavior and an Eulerian approach for the fluid is used. The particle finite element method is adopted for the evaluation of the structural response, whereas an Eulerian fixed-mesh approach is employed for the fluid. The free surface is tracked by the use of a level set technique. The numerical results are validated with experiments on scale models rockfill dams.
Development and Evaluation of a Hybrid Eulerian-Lagrangian Modeling Approach
NASA Astrophysics Data System (ADS)
Choi, Y.; Czader, B.; Percell, P.; Byun, D.
2014-12-01
A hybrid Lagrangian-Eulerian modeling tool has been developed using the CMAQ code. It is a moving nest domain that resides in a base Eulerian model with the same vertical CMAQ structure and physical and chemical interactions. It may follow the trajectory defined by the mean mixed-layer wind or any other user defined trajectory. It was developed as a computationally efficient 3-D grid sub-model for the purpose of evaluations of the source-receptor relationship. The correctness of the algorithms and the overall performance was evaluated against CMAQ simulation results.
Equations of nonlinear dynamics of elastic shells in cylindrical Eulerian coordinates
NASA Astrophysics Data System (ADS)
Zubov, L. M.
2016-05-01
The equations of dynamics of elastic shells subjected to large deformations are formulated. The Eulerian coordinates on a circular cylinder and time are accepted as independent variables, and one of the unknown functions is the distance from a point of the shell surface to the cylinder axis. The equations of dynamics of nonlinearly elastic shells in the Eulerian coordinates are convenient for exact formulation of the problem on the interaction of strongly deformable shells with moving fluids and gases. The equations obtained can be used for dynamic calculations of fluids and gases flowings in pipelines, blood vessels, hoses, and other nonlinearly deformable thin-walled tubular elements of constructions.
Modeling of confined turbulent fluid-particle flows using Eulerian and Lagrangian schemes
NASA Technical Reports Server (NTRS)
Adeniji-Fashola, A.; Chen, C. P.
1990-01-01
Two important aspects of fluid-particulate interaction in dilute gas-particle turbulent flows (the turbulent particle dispersion and the turbulence modulation effects) are addressed, using the Eulerian and Lagrangian modeling approaches to describe the particulate phase. Gradient-diffusion approximations are employed in the Eulerian formulation, while a stochastic procedure is utilized to simulate turbulent dispersion in the Lagrangina formulation. The k-epsilon turbulence model is used to characterize the time and length scales of the continuous phase turbulence. Models proposed for both schemes are used to predict turbulent fully-developed gas-solid vertical pipe flow with reasonable accuracy.
Eulerian Mapping Closure Approach for Probability Density Function of Concentration in Shear Flows
NASA Technical Reports Server (NTRS)
He, Guowei; Bushnell, Dennis M. (Technical Monitor)
2002-01-01
The Eulerian mapping closure approach is developed for uncertainty propagation in computational fluid mechanics. The approach is used to study the Probability Density Function (PDF) for the concentration of species advected by a random shear flow. An analytical argument shows that fluctuation of the concentration field at one point in space is non-Gaussian and exhibits stretched exponential form. An Eulerian mapping approach provides an appropriate approximation to both convection and diffusion terms and leads to a closed mapping equation. The results obtained describe the evolution of the initial Gaussian field, which is in agreement with direct numerical simulations.
jShyLU Scalable Hybrid Preconditioner and Solver
Energy Science and Technology Software Center (ESTSC)
2012-09-11
ShyLU is numerical software to solve sparse linear systems of equations. ShyLU uses a hybrid direct-iterative Schur complement method, and may be used either as a preconditioner or as a solver. ShyLU is parallel and optimized for a single compute Solver node. ShyLU will be a package in the Trilinos software framework.
Experiences with linear solvers for oil reservoir simulation problems
Joubert, W.; Janardhan, R.; Biswas, D.; Carey, G.
1996-12-31
This talk will focus on practical experiences with iterative linear solver algorithms used in conjunction with Amoco Production Company`s Falcon oil reservoir simulation code. The goal of this study is to determine the best linear solver algorithms for these types of problems. The results of numerical experiments will be presented.
Shape reanalysis and sensitivities utilizing preconditioned iterative boundary solvers
NASA Technical Reports Server (NTRS)
Guru Prasad, K.; Kane, J. H.
1992-01-01
The computational advantages associated with the utilization of preconditined iterative equation solvers are quantified for the reanalysis of perturbed shapes using continuum structural boundary element analysis (BEA). Both single- and multi-zone three-dimensional problems are examined. Significant reductions in computer time are obtained by making use of previously computed solution vectors and preconditioners in subsequent analyses. The effectiveness of this technique is demonstrated for the computation of shape response sensitivities required in shape optimization. Computer times and accuracies achieved using the preconditioned iterative solvers are compared with those obtained via direct solvers and implicit differentiation of the boundary integral equations. It is concluded that this approach employing preconditioned iterative equation solvers in reanalysis and sensitivity analysis can be competitive with if not superior to those involving direct solvers.
A real-time impurity solver for DMFT
NASA Astrophysics Data System (ADS)
Kim, Hyungwon; Aron, Camille; Han, Jong E.; Kotliar, Gabriel
Dynamical mean-field theory (DMFT) offers a non-perturbative approach to problems with strongly correlated electrons. The method heavily relies on the ability to numerically solve an auxiliary Anderson-type impurity problem. While powerful Matsubara-frequency solvers have been developed over the past two decades to tackle equilibrium situations, the status of real-time impurity solvers that could compete with Matsubara-frequency solvers and be readily generalizable to non-equilibrium situations is still premature. We present a real-time solver which is based on a quantum Master equation description of the dissipative dynamics of the impurity and its exact diagonalization. As a benchmark, we illustrate the strengths of our solver in the context of the equilibrium Mott-insulator transition of the one-band Hubbard model and compare it with iterative perturbation theory (IPT) method. Finally, we discuss its direct application to a nonequilibrium situation.
Linear iterative solvers for implicit ODE methods
NASA Technical Reports Server (NTRS)
Saylor, Paul E.; Skeel, Robert D.
1990-01-01
The numerical solution of stiff initial value problems, which lead to the problem of solving large systems of mildly nonlinear equations are considered. For many problems derived from engineering and science, a solution is possible only with methods derived from iterative linear equation solvers. A common approach to solving the nonlinear equations is to employ an approximate solution obtained from an explicit method. The error is examined to determine how it is distributed among the stiff and non-stiff components, which bears on the choice of an iterative method. The conclusion is that error is (roughly) uniformly distributed, a fact that suggests the Chebyshev method (and the accompanying Manteuffel adaptive parameter algorithm). This method is described, also commenting on Richardson's method and its advantages for large problems. Richardson's method and the Chebyshev method with the Mantueffel algorithm are applied to the solution of the nonlinear equations by Newton's method.
Scalable Adaptive Multilevel Solvers for Multiphysics Problems
Xu, Jinchao
2014-12-01
In this project, we investigated adaptive, parallel, and multilevel methods for numerical modeling of various real-world applications, including Magnetohydrodynamics (MHD), complex fluids, Electromagnetism, Navier-Stokes equations, and reservoir simulation. First, we have designed improved mathematical models and numerical discretizaitons for viscoelastic fluids and MHD. Second, we have derived new a posteriori error estimators and extended the applicability of adaptivity to various problems. Third, we have developed multilevel solvers for solving scalar partial differential equations (PDEs) as well as coupled systems of PDEs, especially on unstructured grids. Moreover, we have integrated the study between adaptive method and multilevel methods, and made significant efforts and advances in adaptive multilevel methods of the multi-physics problems.
Optimising a parallel conjugate gradient solver
Field, M.R.
1996-12-31
This work arises from the introduction of a parallel iterative solver to a large structural analysis finite element code. The code is called FEX and it was developed at Hitachi`s Mechanical Engineering Laboratory. The FEX package can deal with a large range of structural analysis problems using a large number of finite element techniques. FEX can solve either stress or thermal analysis problems of a range of different types from plane stress to a full three-dimensional model. These problems can consist of a number of different materials which can be modelled by a range of material models. The structure being modelled can have the load applied at either a point or a surface, or by a pressure, a centrifugal force or just gravity. Alternatively a thermal load can be applied with a given initial temperature. The displacement of the structure can be constrained by having a fixed boundary or by prescribing the displacement at a boundary.
General purpose nonlinear system solver based on Newton-Krylov method.
Energy Science and Technology Software Center (ESTSC)
2013-12-01
KINSOL is part of a software family called SUNDIALS: SUite of Nonlinear and Differential/Algebraic equation Solvers [1]. KINSOL is a general-purpose nonlinear system solver based on Newton-Krylov and fixed-point solver technologies [2].
Eulerian simulation of the perforation of aluminum plates by nondeforming projectiles
Silling, S.A.
1992-03-01
A new algorithm for the treatment of sliding interfaces between solids with or without friction in an Eulerian wavecode is described. The algorithm has been implemented in the two-dimensional version of the CTH code. The code was used to simulate penetration and perforation of aluminum plates by rigid, conical-nosed tungsten projectiles. Comparison with experimental data is provided.
In this paper, results of Eulerian grid and Lagrangian photochemical model simulations of emissions from a major elevated point source are presented. eries of simulations with grid sizes varying from 30 km to 2 km were performed with the Urban Airshed Model, a photochemical grid ...
A Monte Carlo technique for the simulation of two-particle relative diffusion is developed using the Eulerian statistical properties of the turbulence. Numerical results are presented which clarify the role of interparticle velocity correlations and reveal the importance of an Eu...
Eulerian-Lagrangian localized adjoint methods for reactive transport in groundwater
Ewing, R.E.; Wang, Hong
1996-12-31
In this paper, we present Eulerian-Lagrangian localized adjoint methods (ELLAM) to solve convection-diffusion-reaction equations governing contaminant transport in groundwater flowing through an adsorbing porous medium. These ELLAM schemes can treat various combinations of boundary conditions and conserve mass. Numerical results are presented to demonstrate the strong potential of ELLAM schemes.
The "Preface to the Special Edition on Model Evaluation: Evaluation of Urban and Regional Eulerian Air Quality Models" is a brief introduction to the papers included in a special issue of Atmospheric Environment. The Preface provides a background for the papers, which have thei...
TREATMENT FOR LAGRANGIAN TRANSPORT AND DIFFUSION OF SUBGRID PLUMES IN A EULERIAN GRID FRAMEWORK
Since Eulerian grid models cannot explicitly treat processes existing on spatial scales smaller than the grid cell size, it is imperative to be able to develop innovative approaches to simulate relevant subgrid scale features within the modeling grid framework. otable class of su...
NASA Astrophysics Data System (ADS)
Diggs, Angela; Balachandar, S.
2016-05-01
The present work addresses numerical methods required to compute particle volume fraction or number density. Local volume fraction of the lth particle, αl, is the quantity of foremost importance in calculating the gas-mediated particle-particle interaction effect in multiphase flows. A general multiphase flow with a distribution of Lagrangian particles inside a fluid flow discretized on an Eulerian grid is considered. Particle volume fraction is needed both as a Lagrangian quantity associated with each particle and also as an Eulerian quantity associated with the grid cell for Eulerian-Lagrangian simulations. In Grid-Based (GB) methods the particle volume fraction is first obtained within each grid cell as an Eulerian quantity and then the local particle volume fraction associated with any Lagrangian particle can be obtained from interpolation. The second class of methods presented are Particle-Based (PB) methods, where particle volume fraction will first be obtained at each particle as a Lagrangian quantity, which then can be projected onto the Eulerian grid. Traditionally, the GB methods are used in multiphase flow, but sub-grid resolution can be obtained through use of the PB methods. By evaluating the total error, and its discretization, bias and statistical error components, the performance of the different PB methods is compared against several common GB methods of calculating volume fraction. The standard von Neumann error analysis technique has been adapted for evaluation of rate of convergence of the different methods. The discussion and error analysis presented focus on the volume fraction calculation, but the methods can be extended to obtain field representations of other Lagrangian quantities, such as particle velocity and temperature.
Comparison of open-source linear programming solvers.
Gearhart, Jared Lee; Adair, Kristin Lynn; Durfee, Justin D.; Jones, Katherine A.; Martin, Nathaniel; Detry, Richard Joseph
2013-10-01
When developing linear programming models, issues such as budget limitations, customer requirements, or licensing may preclude the use of commercial linear programming solvers. In such cases, one option is to use an open-source linear programming solver. A survey of linear programming tools was conducted to identify potential open-source solvers. From this survey, four open-source solvers were tested using a collection of linear programming test problems and the results were compared to IBM ILOG CPLEX Optimizer (CPLEX) [1], an industry standard. The solvers considered were: COIN-OR Linear Programming (CLP) [2], [3], GNU Linear Programming Kit (GLPK) [4], lp_solve [5] and Modular In-core Nonlinear Optimization System (MINOS) [6]. As no open-source solver outperforms CPLEX, this study demonstrates the power of commercial linear programming software. CLP was found to be the top performing open-source solver considered in terms of capability and speed. GLPK also performed well but cannot match the speed of CLP or CPLEX. lp_solve and MINOS were considerably slower and encountered issues when solving several test problems.
A non-conforming 3D spherical harmonic transport solver
Van Criekingen, S.
2006-07-01
A new 3D transport solver for the time-independent Boltzmann transport equation has been developed. This solver is based on the second-order even-parity form of the transport equation. The angular discretization is performed through the expansion of the angular neutron flux in spherical harmonics (PN method). The novelty of this solver is the use of non-conforming finite elements for the spatial discretization. Such elements lead to a discontinuous flux approximation. This interface continuity requirement relaxation property is shared with mixed-dual formulations such as the ones based on Raviart-Thomas finite elements. Encouraging numerical results are presented. (authors)
Multi-GPU kinetic solvers using MPI and CUDA
NASA Astrophysics Data System (ADS)
Zabelok, Sergey; Arslanbekov, Robert; Kolobov, Vladimir
2014-12-01
This paper describes recent progress towards porting a Unified Flow Solver (UFS) to heterogeneous parallel computing. The main challenge of porting UFS to graphics processing units (GPUs) comes from the dynamically adapted mesh, which causes irregular data access. We describe the implementation of CUDA kernels for three modules in UFS: the direct Boltzmann solver using discrete velocity method (DVM), the DSMC module, and the Lattice Boltzmann Method (LBM) solver, all using octree Cartesian mesh with adaptive Mesh Refinement (AMR). Double digit speedup on single GPU and good scaling for multi-GPU has been demonstrated.
Elliptic Solvers for Adaptive Mesh Refinement Grids
Quinlan, D.J.; Dendy, J.E., Jr.; Shapira, Y.
1999-06-03
We are developing multigrid methods that will efficiently solve elliptic problems with anisotropic and discontinuous coefficients on adaptive grids. The final product will be a library that provides for the simplified solution of such problems. This library will directly benefit the efforts of other Laboratory groups. The focus of this work is research on serial and parallel elliptic algorithms and the inclusion of our black-box multigrid techniques into this new setting. The approach applies the Los Alamos object-oriented class libraries that greatly simplify the development of serial and parallel adaptive mesh refinement applications. In the final year of this LDRD, we focused on putting the software together; in particular we completed the final AMR++ library, we wrote tutorials and manuals, and we built example applications. We implemented the Fast Adaptive Composite Grid method as the principal elliptic solver. We presented results at the Overset Grid Conference and other more AMR specific conferences. We worked on optimization of serial and parallel performance and published several papers on the details of this work. Performance remains an important issue and is the subject of continuing research work.
Advanced Multigrid Solvers for Fluid Dynamics
NASA Technical Reports Server (NTRS)
Brandt, Achi
1999-01-01
The main objective of this project has been to support the development of multigrid techniques in computational fluid dynamics that can achieve "textbook multigrid efficiency" (TME), which is several orders of magnitude faster than current industrial CFD solvers. Toward that goal we have assembled a detailed table which lists every foreseen kind of computational difficulty for achieving it, together with the possible ways for resolving the difficulty, their current state of development, and references. We have developed several codes to test and demonstrate, in the framework of simple model problems, several approaches for overcoming the most important of the listed difficulties that had not been resolved before. In particular, TME has been demonstrated for incompressible flows on one hand, and for near-sonic flows on the other hand. General approaches were advanced for the relaxation of stagnation points and boundary conditions under various situations. Also, new algebraic multigrid techniques were formed for treating unstructured grid formulations. More details on all these are given below.
NASA Technical Reports Server (NTRS)
Raju, Manthena S.
1998-01-01
Sprays occur in a wide variety of industrial and power applications and in the processing of materials. A liquid spray is a phase flow with a gas as the continuous phase and a liquid as the dispersed phase (in the form of droplets or ligaments). Interactions between the two phases, which are coupled through exchanges of mass, momentum, and energy, can occur in different ways at different times and locations involving various thermal, mass, and fluid dynamic factors. An understanding of the flow, combustion, and thermal properties of a rapidly vaporizing spray requires careful modeling of the rate-controlling processes associated with the spray's turbulent transport, mixing, chemical kinetics, evaporation, and spreading rates, as well as other phenomena. In an attempt to advance the state-of-the-art in multidimensional numerical methods, we at the NASA Lewis Research Center extended our previous work on sprays to unstructured grids and parallel computing. LSPRAY, which was developed by M.S. Raju of Nyma, Inc., is designed to be massively parallel and could easily be coupled with any existing gas-phase flow and/or Monte Carlo probability density function (PDF) solver. The LSPRAY solver accommodates the use of an unstructured mesh with mixed triangular, quadrilateral, and/or tetrahedral elements in the gas-phase solvers. It is used specifically for fuel sprays within gas turbine combustors, but it has many other uses. The spray model used in LSPRAY provided favorable results when applied to stratified-charge rotary combustion (Wankel) engines and several other confined and unconfined spray flames. The source code will be available with the National Combustion Code (NCC) as a complete package.
Handling Vacuum Regions in a Hybrid Plasma Solver
NASA Astrophysics Data System (ADS)
Holmström, M.
2013-04-01
In a hybrid plasma solver (particle ions, fluid mass-less electrons) regions of vacuum, or very low charge density, can cause problems since the evaluation of the electric field involves division by charge density. This causes large electric fields in low density regions that can lead to numerical instabilities. Here we propose a self consistent handling of vacuum regions for hybrid solvers. Vacuum regions can be considered having infinite resistivity, and in this limit Faraday's law approaches a magnetic diffusion equation. We describe an algorithm that solves such a diffusion equation in regions with charge density below a threshold value. We also present an implementation of this algorithm in a hybrid plasma solver, and an application to the interaction between the Moon and the solar wind. We also discuss the implementation of hyperresistivity for smoothing the electric field in a PIC solver.
Parallel iterative solvers and preconditioners using approximate hierarchical methods
Grama, A.; Kumar, V.; Sameh, A.
1996-12-31
In this paper, we report results of the performance, convergence, and accuracy of a parallel GMRES solver for Boundary Element Methods. The solver uses a hierarchical approximate matrix-vector product based on a hybrid Barnes-Hut / Fast Multipole Method. We study the impact of various accuracy parameters on the convergence and show that with minimal loss in accuracy, our solver yields significant speedups. We demonstrate the excellent parallel efficiency and scalability of our solver. The combined speedups from approximation and parallelism represent an improvement of several orders in solution time. We also develop fast and paralellizable preconditioners for this problem. We report on the performance of an inner-outer scheme and a preconditioner based on truncated Green`s function. Experimental results on a 256 processor Cray T3D are presented.
Simulation of the dispersion of nuclear contamination using an adaptive Eulerian grid model.
Lagzi, I; Kármán, D; Turányi, T; Tomlin, A S; Haszpra, L
2004-01-01
Application of an Eulerian model using layered adaptive unstructured grids coupled to a meso-scale meteorological model is presented for modelling the dispersion of nuclear contamination following the accidental release from a single but strong source to the atmosphere. The model automatically places a finer resolution grid, adaptively in time, in regions were high spatial numerical error is expected. The high-resolution grid region follows the movement of the contaminated air over time. Using this method, grid resolutions of the order of 6 km can be achieved in a computationally effective way. The concept is illustrated by the simulation of hypothetical nuclear accidents at the Paks NPP, in Central Hungary. The paper demonstrates that the adaptive model can achieve accuracy comparable to that of a high-resolution Eulerian model using significantly less grid points and computer simulation time. PMID:15149762
NASA Astrophysics Data System (ADS)
Rubin, M. B.; Nadler, B.
2016-03-01
An Eulerian formulation has been developed for the constitutive response of a group of materials that includes anisotropic elastic and viscoelastic solids and viscous fluids. The material is considered to be a composite of an elastic solid and a viscous fluid. Evolution equations are proposed for a triad of vectors m i that represent the stretches and orientations of material line elements in the solid component. Evolution equations for an orthonormal triad of vectors s i are also proposed to characterize anisotropy of the fluid component. In particular, for an elastic solid it is shown that the material response is totally characterized by the functional form of the strain energy and by the current values of m i , which are measurable in the current state of the material. Moreover, it is shown that the proposed Eulerian formulation removes unphysical arbitrariness of the choice of the reference configuration in the standard formulation of constitutive equations for anisotropic elastic solids.
Eulerian CFD modeling and X-ray validation of non-evaporating diesel spray
NASA Astrophysics Data System (ADS)
Xue, Qingluan; Som, Sibendu; Quan, Shaoping; Pomraning, Eric; Senecal, P. K.
2013-11-01
This work implemented an Eulerian single-phase approach by Vallet et al. into CFD software (Convergent) for diesel spray simulations. This Eulerian approach considers liquid and gas phase as a complex mixture of a single flow with a highly variable density to describe the near nozzle dense sprays. The mean density is obtained form the Favre-averaged liquid mass fraction. Liquid mass fraction is transported with a model for the turbulent liquid diffusion flux into the gas. A mean gradient-based model is employed for the diffusion flux in this study. A non-evaporating diesel spray was measured using x-ray radiography at Argonne National Laboratory. The quantitative and time-resolved data of liquid penetration and mass distribution in the dense spray region are used to validate this approach. The different turbulence models are also used for the simulations. The comparison between the simulated results and experimental data and the turbulence model effect are discussed.
Large eddy simulation of Rayleigh-Taylor instability using the arbitrary Lagrangian-Eulerian method
Darlington, R
1999-12-01
This research addresses the application of a large eddy simulation (LES) to Arbitrary Lagrangian Eulerian (ALE) simulations of Rayleigh-Taylor instability. First, ALE simulations of simplified Rayleigh-Taylor instability are studied. The advantages of ALE over Eulerian simulations are shown. Next, the behavior of the LES is examined in a more complicated ALE simulation of Rayleigh-Taylor instability. The effects of eddy viscosity and stochastic backscatter are examined. The LES is also coupled with ALE to increase grid resolution in areas where it is needed. Finally, the methods studied above are applied to two sets of experimental simulations. In these simulations, ALE allows the mesh to follow expanding experimental targets, while LES can be used to mimic the effect of unresolved instability modes.
Cosmological dynamics: from the Eulerian to the Lagrangian frame. Part I. Newtonian approximation
Villa, Eleonora; Maino, Davide; Matarrese, Sabino E-mail: sabino.matarrese@pd.infn.it
2014-06-01
We analyse the non-linear gravitational dynamics of a pressure-less fluid in the Newtonian limit of General Relativity in both the Eulerian and Lagrangian pictures. Starting from the Newtonian metric in the Poisson gauge, we transform to the synchronous and comoving gauge and obtain the Lagrangian metric within the Newtonian approximation. Our approach is fully non-perturbative, which implies that if our quantities are expanded according to the rules of standard perturbation theory, all terms are exactly recovered at any order in perturbation theory, only provided they are Newtonian. We explicitly show this result up to second order and in both gauges. Our transformation clarifies the meaning of the change of spatial and time coordinates from the Eulerian to the Lagrangian frame in the Newtonian approximation.
A three-dimensional Eulerian method for the numerical simulation of high-velocity impact problems
NASA Astrophysics Data System (ADS)
Wu, Shi-Yu; Liu, Kai-Xin; Chen, Qian-Yi
2014-03-01
In the present paper, a three-dimensional (3D) Eulerian technique for the 3D numerical simulation of high-velocity impact problems is proposed. In the Eulerian framework, a complete 3D conservation element and solution element scheme for conservative hyperbolic governing equations with source terms is given. A modified ghost fluid method is proposed for the treatment of the boundary conditions. Numerical simulations of the Taylor bar problem and the ricochet phenomenon of a sphere impacting a plate target at an angle of 60° are carried out. The numerical results are in good agreement with the corresponding experimental observations. It is proved that our computational technique is feasible for analyzing 3D high-velocity impact problems.
Littoral transport in the surf zone elucidated by an Eulerian sediment tracer.
Duane, D.B.; James, W.R.
1980-01-01
An Eulerian, or time integration, sand tracer experiment was designed and carried out in the surf zone near Pt. Mugu, California on April 19, 1972. Data indicate that conditions of stationarity and finite boundaries required for proper application of Eulerian tracer theory exist for short time periods in the surf zone. Grain counts suggest time required for tracer sand to attain equilibrium concentration is on the order of 30-60 minutes. Grain counts also indicate transport (discharge) was strongly dependent upon grain size, with the maximum rate occurring in the size 2.5-2.75 phi, decreasing to both finer and coarser sizes. The measured instantaneous transport was at the annual rate of 2.4 x 106 m3/yr.- Authors
Comparison of Direct Eulerian Godunov and Lagrange Plus Remap, Artificial Viscosity Schemes
Pember, R B; Anderson, R W
2001-03-30
The authors compare two algorithms for solving the equations of unsteady inviscid compressible flow in an Eulerian frame: a staggered grid, Lagrange plus remap artificial viscosity scheme and a cell-centered, direct Eulerian higher-order Godunov scheme. They use the two methods to compute solutions to a number of one- and two-dimensional problems. The results show the accuracy of the two schemes to be generally equivalent. In a 1984 survey paper by Woodward and Colella, the Lagrange plus remap approach did not compare favorably with the higher-order Godunov methodology. They examine, therefore, how certain features of the staggered grid scheme considered here contribute to its improved accuracy. The critical features are shown to be the use of a monotonic artificial viscosity in the Lagrange step and, in the remap step, the use of a corner transport upwind scheme with van Leer limiters in conjunction with separate advection of internal and kinetic energies.
Simulations of needle insertion by using a Eulerian hydrocode FEM and the experimental validations.
Kataoka, Hiroyuki; Noda, Shigeho; Yokota, Hideo; Takagi, Shu; Himeno, Ryutaro; Okazawa, Shigenobu
2008-01-01
In this paper, simulations for needle insertion were performed by using a novel Eulerian hydrocode FEM, which was adaptive for large deformation and tissue fracture. We also performed experiments for the same needle insertion with silicon rubbers and needles, which had conical tips of different angles in order to investigate the accuracy of the simulations. The resistance forces in the simulations accurately followed those in the experiments until the conical portion of the needle was inside the rubbers, and the validation of the Eulerian hydrocode was revealed. However, the present simulation showed that after the conical portion was inside the tissue, the simulated resistance forces became lower than the experimental ones. The proportional increase of the friction forces and the roughly flatness of the tip force along the time were simulated. It was predicted that the tightening force along the needle side was underestimated. PMID:18982649
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.
A Lagrangian-Eulerian finite element method with adaptive gridding for advection-dispersion problems
Ijiri, Y.; Karasaki, K.
1994-02-01
In the present paper, a Lagrangian-Eulerian finite element method with adaptive gridding for solving advection-dispersion equations is described. The code creates new grid points in the vicinity of sharp fronts at every time step in order to reduce numerical dispersion. The code yields quite accurate solutions for a wide range of mesh Peclet numbers and for mesh Courant numbers well in excess of 1.
Unit physics performance of a mix model in Eulerian fluid computations
Vold, Erik; Douglass, Rod
2011-01-25
In this report, we evaluate the performance of a K-L drag-buoyancy mix model, described in a reference study by Dimonte-Tipton [1] hereafter denoted as [D-T]. The model was implemented in an Eulerian multi-material AMR code, and the results are discussed here for a series of unit physics tests. The tests were chosen to calibrate the model coefficients against empirical data, principally from RT (Rayleigh-Taylor) and RM (Richtmyer-Meshkov) experiments, and the present results are compared to experiments and to results reported in [D-T]. Results show the Eulerian implementation of the mix model agrees well with expectations for test problems in which there is no convective flow of the mass averaged fluid, i.e., in RT mix or in the decay of homogeneous isotropic turbulence (HIT). In RM shock-driven mix, the mix layer moves through the Eulerian computational grid, and there are differences with the previous results computed in a Lagrange frame [D-T]. The differences are attributed to the mass averaged fluid motion and examined in detail. Shock and re-shock mix are not well matched simultaneously. Results are also presented and discussed regarding model sensitivity to coefficient values and to initial conditions (IC), grid convergence, and the generation of atomically mixed volume fractions.
A new Eulerian method to estimate "spicy" Agulhas leakage in climate models
NASA Astrophysics Data System (ADS)
Putrasahan, D. A.; Beal, Lisa M.; Kirtman, Ben P.; Cheng, Yu
2015-06-01
The Agulhas leakage brings warm and saline water from the Indian Ocean into the South Atlantic, forming part of the returning branch of the global thermohaline circulation, an important component of climate. We develop a new method using the passive tracer spice, to estimate Agulhas leakage transport in an ocean eddy-resolving climate model that properly captures Agulhas Retroflection and leakage. We identify Agulhas leakage waters as positive spice anomalies (≥ +0.1) and are able to trace leakage to depths greater than 1000 m. Spice-based Eulerian Agulhas leakage captures spikes in Lagrangian leakage transport which coincide with the passage of Agulhas rings, yielding a statistically significant correlation (0.47). On interannual times scales, differences between Eulerian and Lagrangian transports are reduced and correlation increases to 0.77. We obtain a mean spice-based Eulerian Agulhas leakage of 17 ± 4 sverdrup (Sv) and Lagrangian Agulhas leakage of 18 ± 2 Sv, which are within the range of observational estimates.
Eulerian Simulation of Acoustic Waves Over Long Range in Realistic Environments
NASA Astrophysics Data System (ADS)
Chitta, Subhashini; Steinhoff, John
2015-11-01
In this paper, we describe a new method for computation of long-range acoustics. The approach is a hybrid of near and far-field methods, and is unique in its Eulerian treatment of the far-field propagation. The near-field generated by any existing method to project an acoustic solution onto a spherical surface that surrounds a source. The acoustic field on this source surface is then extended to an arbitrarily large distance in an inhomogeneous far-field. This would normally require an Eulerian solution of the wave equation. However, conventional Eulerian methods have prohibitive grid requirements. This problem is overcome by using a new method, ``Wave Confinement'' (WC) that propagates wave-identifying phase fronts as nonlinear solitary waves that live on grid indefinitely. This involves modification of wave equation by the addition of a nonlinear term without changing the basic conservation properties of the equation. These solitary waves can then be used to ``carry'' the essential integrals of the acoustic wave. For example, arrival time, centroid position and other properties that are invariant as the wave passes a grid point. Because of this property the grid can be made as coarse as necessary, consistent with overall accuracy to resolve atmospheric/ground variations. This work is being funded by the U.S. Army under a Small Business Innovation Research (SBIR) program (contract number: # W911W6-12-C-0036). The authors would like to thank Dr. Frank Caradonna and Dr. Ben W. Sim for this support.
NASA Astrophysics Data System (ADS)
Fox, Rodney O.; Vie, Aymeric; Laurent, Frederique; Chalons, Christophe; Massot, Marc
2012-11-01
Numerous applications involve a disperse phase carried by a gaseous flow. To simulate such flows, one can resort to a number density function (NDF) governed a kinetic equation. Traditionally, Lagrangian Monte-Carlo methods are used to solve for the NDF, but are expensive as the number of numerical particles needed must be large to control statistical errors. Moreover, such methods are not well adapted to high-performance computing because of the intrinsic inhomogeneity of the NDF. To overcome these issues, Eulerian methods can be used to solve for the moments of the NDF resulting in an unclosed Eulerian system of hyperbolic conservation laws. To obtain closure, in this work a multivariate bi-Gaussian quadrature is used, which can account for particle trajectory crossing (PTC) over a large range of Stokes numbers. This closure uses up to four quadrature points in 2-D velocity phase space to capture large-scale PTC, and an anisotropic Gaussian distribution around each quadrature point to model small-scale PTC. Simulations of 2-D particle-laden isotropic turbulence at different Stokes numbers are employed to validate the Eulerian models against results from the Lagrangian approach. Good agreement is found for the number density fields over the entire range of Stokes numbers tested. Research carried out at the Center for Turbulence Research 2012 Summer Program.
An Eulerian-Lagrangian Form for the Euler Equations in Sobolev Spaces
NASA Astrophysics Data System (ADS)
Pooley, Benjamin C.; Robinson, James C.
2016-06-01
In 2000 Constantin showed that the incompressible Euler equations can be written in an "Eulerian-Lagrangian" form which involves the back-to-labels map (the inverse of the trajectory map for each fixed time). In the same paper a local existence result is proved in certain Hölder spaces {C^{1,μ}} . We review the Eulerian-Lagrangian formulation of the equations and prove that given initial data in H s for {n ≥ 2} and {s > n/2+1} , a unique local-in-time solution exists on the n-torus that is continuous into H s and C 1 into H s-1. These solutions automatically have C 1 trajectories. The proof here is direct and does not appeal to results already known about the classical formulation. Moreover, these solutions are regular enough that the classical and Eulerian-Lagrangian formulations are equivalent, therefore what we present amounts to an alternative approach to some of the standard theory.
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.
NASA Astrophysics Data System (ADS)
Abalos, M.; Ploeger, F.; Konopka, P.; Randel, W. J.; Serrano, E.
2013-11-01
We aim to reconcile the recently published, apparently contrasting results regarding the relative importance of tropical upwelling versus horizontal transport for the seasonality of ozone above the tropical tropopause. Different analysis methods in the literature (Lagrangian versus Eulerian, and isentropic versus pressure vertical coordinates) yield different perspectives of ozone transport, and the results must be carefully compared in equivalent terms to avoid misinterpretation. By examining the Lagrangian calculations in the Eulerian formulation, we show here that the results are in fact consistent with each other and with a common understanding of the ozone transport processes near and above the tropical tropopause. We further emphasize that the complementary approaches are suited for answering two different scientific questions: (1) what drives the observed seasonal cycle in ozone at a particular level above the tropical tropopause? and (2) how important is horizontal transport from mid-latitudes for ozone concentrations in the tropical lower stratosphere? Regarding the first question, the analysis of the transformed Eulerian mean (TEM) ozone budget shows that the annual cycle in tropical upwelling is the main forcing of the ozone seasonality at altitudes with large vertical gradients in the tropical lower stratosphere. To answer the second question a Lagrangian framework must be used, and the results show that a large fraction (~50%) of the ozone molecules ascending through the tropical lower stratosphere is of extra-tropical origin and has been in-mixed from mid-latitudes.
NASA Astrophysics Data System (ADS)
Abalos, M.; Ploeger, F.; Konopka, P.; Randel, W. J.; Serrano, E.
2013-07-01
We aim to reconcile the recently published, apparently contrasting results regarding the relative importance of tropical upwelling versus horizontal transport for the seasonality of ozone above the tropical tropopause. Different analysis methods in the literature (Lagrangian versus Eulerian, and isentropic versus pressure vertical coordinates) yield different perspectives of ozone transport, and the results must be carefully compared in equivalent terms to avoid misinterpretation. By examining the Lagrangian calculations in the Eulerian formulation, we show here that the results are in fact consistent with each other and with a common understanding of the ozone transport processes near and above the tropical tropopause. We further emphasize that the complementary approaches are suited for answering two different scientific questions: (1) what drives the observed seasonal cycle in ozone at a particular level above the tropical tropopause? and (2) how important is horizontal transport from mid-latitudes for ozone concentrations in the tropical lower stratosphere? Regarding the first question, the analysis of the Transformed Eulerian Mean (TEM) ozone budget shows that the annual cycle in tropical upwelling is the main forcing of the ozone seasonality at altitudes with large vertical gradients in the tropical lower stratosphere. To answer the second question a Lagrangian framework must be used, and the results show that a large fraction (∼50%) of the ozone molecules ascending through the tropical lower stratosphere is of extra-tropical origin and has been in-mixed from mid-latitudes.
Eulerian laser Doppler vibrometry: Online blade damage identification on a multi-blade test rotor
NASA Astrophysics Data System (ADS)
Oberholster, A. J.; Heyns, P. S.
2011-01-01
Laser Doppler vibrometry enables the telemetry-free measurement of online turbomachinery blade vibration. Specifically, the Eulerian or fixed reference frame implementation of laser vibrometry provides a practical solution to the condition monitoring of rotating blades. The short data samples that are characteristic of this measurement approach do however negate the use of traditional frequency domain signal processing techniques. It is therefore necessary to employ techniques such as time domain analysis and non-harmonic Fourier analysis to obtain useful information from the blade vibration signatures. The latter analysis technique allows the calculation of phase angle trends which can be used as indicators of blade health deterioration, as has been shown in previous work for a single-blade rotor. This article presents the results from tests conducted on a five-blade axial-flow test rotor at different rotor speeds and measurement positions. With the aid of artificial neural networks, it is demonstrated that the parameters obtained from non-harmonic Fourier analysis and time domain signal processing on Eulerian laser Doppler vibrometry signals can successfully be used to identify and quantify blade damage from among healthy blades. It is also shown that the natural frequencies of individual blades can be approximated from the Eulerian signatures recorded during rotor run-up and run-down.
NASA Astrophysics Data System (ADS)
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.
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. PMID:26871161
An advanced implicit solver for MHD
NASA Astrophysics Data System (ADS)
Udrea, Bogdan
A new implicit algorithm has been developed for the solution of the time-dependent, viscous and resistive single fluid magnetohydrodynamic (MHD) equations. The algorithm is based on an approximate Riemann solver for the hyperbolic fluxes and central differencing applied on a staggered grid for the parabolic fluxes. The algorithm employs a locally aligned coordinate system that allows the solution to the Riemann problems to be solved in a natural direction, normal to cell interfaces. The result is an original scheme that is robust and reduces the complexity of the flux formulas. The evaluation of the parabolic fluxes is also implemented using a locally aligned coordinate system, this time on the staggered grid. The implicit formulation employed by WARP3 is a two level scheme that was applied for the first time to the single fluid MHD model. The flux Jacobians that appear in the implicit scheme are evaluated numerically. The linear system that results from the implicit discretization is solved using a robust symmetric Gauss-Seidel method. The code has an explicit mode capability so that implementation and test of new algorithms or new physics can be performed in this simpler mode. Last but not least the code was designed and written to run on parallel computers so that complex, high resolution runs can be per formed in hours rather than days. The code has been benchmarked against analytical and experimental gas dynamics and MHD results. The benchmarks consisted of one-dimensional Riemann problems and diffusion dominated problems, two-dimensional supersonic flow over a wedge, axisymmetric magnetoplasmadynamic (MPD) thruster simulation and three-dimensional supersonic flow over intersecting wedges and spheromak stability simulation. The code has been proven to be robust and the results of the simulations showed excellent agreement with analytical and experimental results. Parallel performance studies showed that the code performs as expected when run on parallel
Benchmarking transport solvers for fracture flow problems
NASA Astrophysics Data System (ADS)
Olkiewicz, Piotr; Dabrowski, Marcin
2015-04-01
Fracture flow may dominate in rocks with low porosity and it can accompany both industrial and natural processes. Typical examples of such processes are natural flows in crystalline rocks and industrial flows in geothermal systems or hydraulic fracturing. Fracture flow provides an important mechanism for transporting mass and energy. For example, geothermal energy is primarily transported by the flow of the heated water or steam rather than by the thermal diffusion. The geometry of the fracture network and the distribution of the mean apertures of individual fractures are the key parameters with regard to the fracture network transmissivity. Transport in fractures can occur through the combination of advection and diffusion processes like in the case of dissolved chemical components. The local distribution of the fracture aperture may play an important role for both flow and transport processes. In this work, we benchmark various numerical solvers for flow and transport processes in a single fracture in 2D and 3D. Fracture aperture distributions are generated by a number of synthetic methods. We examine a single-phase flow of an incompressible viscous Newtonian fluid in the low Reynolds number limit. Periodic boundary conditions are used and a pressure difference is imposed in the background. The velocity field is primarly found using the Stokes equations. We systematically compare the obtained velocity field to the results obtained by solving the Reynolds equation. This allows us to examine the impact of the aperture distribution on the permeability of the medium and the local velocity distribution for two different mathematical descriptions of the fracture flow. Furthermore, we analyse the impact of aperture distribution on the front characteristics such as the standard deviation and the fractal dimension for systems in 2D and 3D.
A robust multilevel simultaneous eigenvalue solver
NASA Technical Reports Server (NTRS)
Costiner, Sorin; Taasan, Shlomo
1993-01-01
Multilevel (ML) algorithms for eigenvalue problems are often faced with several types of difficulties such as: the mixing of approximated eigenvectors by the solution process, the approximation of incomplete clusters of eigenvectors, the poor representation of solution on coarse levels, and the existence of close or equal eigenvalues. Algorithms that do not treat appropriately these difficulties usually fail, or their performance degrades when facing them. These issues motivated the development of a robust adaptive ML algorithm which treats these difficulties, for the calculation of a few eigenvectors and their corresponding eigenvalues. The main techniques used in the new algorithm include: the adaptive completion and separation of the relevant clusters on different levels, the simultaneous treatment of solutions within each cluster, and the robustness tests which monitor the algorithm's efficiency and convergence. The eigenvectors' separation efficiency is based on a new ML projection technique generalizing the Rayleigh Ritz projection, combined with a technique, the backrotations. These separation techniques, when combined with an FMG formulation, in many cases lead to algorithms of O(qN) complexity, for q eigenvectors of size N on the finest level. Previously developed ML algorithms are less focused on the mentioned difficulties. Moreover, algorithms which employ fine level separation techniques are of O(q(sub 2)N) complexity and usually do not overcome all these difficulties. Computational examples are presented where Schrodinger type eigenvalue problems in 2-D and 3-D, having equal and closely clustered eigenvalues, are solved with the efficiency of the Poisson multigrid solver. A second order approximation is obtained in O(qN) work, where the total computational work is equivalent to only a few fine level relaxations per eigenvector.
A Comparative Study of Randomized Constraint Solvers for Random-Symbolic Testing
NASA Technical Reports Server (NTRS)
Takaki, Mitsuo; Cavalcanti, Diego; Gheyi, Rohit; Iyoda, Juliano; dAmorim, Marcelo; Prudencio, Ricardo
2009-01-01
The complexity of constraints is a major obstacle for constraint-based software verification. Automatic constraint solvers are fundamentally incomplete: input constraints often build on some undecidable theory or some theory the solver does not support. This paper proposes and evaluates several randomized solvers to address this issue. We compare the effectiveness of a symbolic solver (CVC3), a random solver, three hybrid solvers (i.e., mix of random and symbolic), and two heuristic search solvers. We evaluate the solvers on two benchmarks: one consisting of manually generated constraints and another generated with a concolic execution of 8 subjects. In addition to fully decidable constraints, the benchmarks include constraints with non-linear integer arithmetic, integer modulo and division, bitwise arithmetic, and floating-point arithmetic. As expected symbolic solving (in particular, CVC3) subsumes the other solvers for the concolic execution of subjects that only generate decidable constraints. For the remaining subjects the solvers are complementary.
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.
Quantitative analysis of numerical solvers for oscillatory biomolecular system models
Quo, Chang F; Wang, May D
2008-01-01
Background This article provides guidelines for selecting optimal numerical solvers for biomolecular system models. Because various parameters of the same system could have drastically different ranges from 10-15 to 1010, the ODEs can be stiff and ill-conditioned, resulting in non-unique, non-existing, or non-reproducible modeling solutions. Previous studies have not examined in depth how to best select numerical solvers for biomolecular system models, which makes it difficult to experimentally validate the modeling results. To address this problem, we have chosen one of the well-known stiff initial value problems with limit cycle behavior as a test-bed system model. Solving this model, we have illustrated that different answers may result from different numerical solvers. We use MATLAB numerical solvers because they are optimized and widely used by the modeling community. We have also conducted a systematic study of numerical solver performances by using qualitative and quantitative measures such as convergence, accuracy, and computational cost (i.e. in terms of function evaluation, partial derivative, LU decomposition, and "take-off" points). The results show that the modeling solutions can be drastically different using different numerical solvers. Thus, it is important to intelligently select numerical solvers when solving biomolecular system models. Results The classic Belousov-Zhabotinskii (BZ) reaction is described by the Oregonator model and is used as a case study. We report two guidelines in selecting optimal numerical solver(s) for stiff, complex oscillatory systems: (i) for problems with unknown parameters, ode45 is the optimal choice regardless of the relative error tolerance; (ii) for known stiff problems, both ode113 and ode15s are good choices under strict relative tolerance conditions. Conclusions For any given biomolecular model, by building a library of numerical solvers with quantitative performance assessment metric, we show that it is possible
Euler/Navier-Stokes Solvers Applied to Ducted Fan Configurations
NASA Technical Reports Server (NTRS)
Keith, Theo G., Jr.; Srivastava, Rakesh
1997-01-01
Due to noise considerations, ultra high bypass ducted fans have become a more viable design. These ducted fans typically consist of a rotor stage containing a wide chord fan and a stator stage. One of the concerns for this design is the classical flutter that keeps occurring in various unducted fan blade designs. These flutter are catastrophic and are to be avoided in the flight envelope of the engine. Some numerical investigations by Williams, Cho and Dalton, have suggested that a duct around a propeller makes it more unstable. This needs to be further investigated. In order to design an engine to safely perform a set of desired tasks, accurate information of the stresses on the blade during the entire cycle of blade motion is required. This requirement in turn demands that accurate knowledge of steady and unsteady blade loading be available. Aerodynamic solvers based on unsteady three-dimensional analysis will provide accurate and fast solutions and are best suited for aeroelastic analysis. The Euler solvers capture significant physics of the flowfield and are reasonably fast. An aerodynamic solver Ref. based on Euler equations had been developed under a separate grant from NASA Lewis in the past. Under the current grant, this solver has been modified to calculate the aeroelastic characteristics of unducted and ducted rotors. Even though, the aeroelastic solver based on three-dimensional Euler equations is computationally efficient, it is still very expensive to investigate the effects of multiple stages on the aeroelastic characteristics. In order to investigate the effects of multiple stages, a two-dimensional multi stage aeroelastic solver was also developed under this task, in collaboration with Dr. T. S. R. Reddy of the University of Toledo. Both of these solvers were applied to several test cases and validated against experimental data, where available.
EUPDF-II: An Eulerian Joint Scalar Monte Carlo PDF Module : User's Manual
NASA Technical Reports Server (NTRS)
Raju, M. S.; Liu, Nan-Suey (Technical Monitor)
2004-01-01
EUPDF-II provides the solution for the species and temperature fields based on an evolution equation for PDF (Probability Density Function) and it is developed mainly for application with sprays, combustion, parallel computing, and unstructured grids. It is designed to be massively parallel and could easily be coupled with any existing gas-phase CFD and spray solvers. The solver accommodates the use of an unstructured mesh with mixed elements of either triangular, quadrilateral, and/or tetrahedral type. The manual provides the user with an understanding of the various models involved in the PDF formulation, its code structure and solution algorithm, and various other issues related to parallelization and its coupling with other solvers. The source code of EUPDF-II will be available with National Combustion Code (NCC) as a complete package.
Performance Models for the Spike Banded Linear System Solver
Manguoglu, Murat; Saied, Faisal; Sameh, Ahmed; Grama, Ananth
2011-01-01
With availability of large-scale parallel platforms comprised of tens-of-thousands of processors and beyond, there is significant impetus for the development of scalable parallel sparse linear system solvers and preconditioners. An integral part of this design process is the development of performance models capable of predicting performance and providing accurate cost models for the solvers and preconditioners. There has been some work in the past on characterizing performance of the iterative solvers themselves. In this paper, we investigate the problem of characterizing performance and scalability of banded preconditioners. Recent work has demonstrated the superior convergence properties and robustness of banded preconditioners,more » compared to state-of-the-art ILU family of preconditioners as well as algebraic multigrid preconditioners. Furthermore, when used in conjunction with efficient banded solvers, banded preconditioners are capable of significantly faster time-to-solution. Our banded solver, the Truncated Spike algorithm is specifically designed for parallel performance and tolerance to deep memory hierarchies. Its regular structure is also highly amenable to accurate performance characterization. Using these characteristics, we derive the following results in this paper: (i) we develop parallel formulations of the Truncated Spike solver, (ii) we develop a highly accurate pseudo-analytical parallel performance model for our solver, (iii) we show excellent predication capabilities of our model – based on which we argue the high scalability of our solver. Our pseudo-analytical performance model is based on analytical performance characterization of each phase of our solver. These analytical models are then parameterized using actual runtime information on target platforms. An important consequence of our performance models is that they reveal underlying performance bottlenecks in both serial and parallel formulations. All of our results are validated
Performance of Nonlinear Finite-Difference Poisson-Boltzmann Solvers.
Cai, Qin; Hsieh, Meng-Juei; Wang, Jun; Luo, Ray
2010-01-12
We implemented and optimized seven finite-difference solvers for the full nonlinear Poisson-Boltzmann equation in biomolecular applications, including four relaxation methods, one conjugate gradient method, and two inexact Newton methods. The performance of the seven solvers was extensively evaluated with a large number of nucleic acids and proteins. Worth noting is the inexact Newton method in our analysis. We investigated the role of linear solvers in its performance by incorporating the incomplete Cholesky conjugate gradient and the geometric multigrid into its inner linear loop. We tailored and optimized both linear solvers for faster convergence rate. In addition, we explored strategies to optimize the successive over-relaxation method to reduce its convergence failures without too much sacrifice in its convergence rate. Specifically we attempted to adaptively change the relaxation parameter and to utilize the damping strategy from the inexact Newton method to improve the successive over-relaxation method. Our analysis shows that the nonlinear methods accompanied with a functional-assisted strategy, such as the conjugate gradient method and the inexact Newton method, can guarantee convergence in the tested molecules. Especially the inexact Newton method exhibits impressive performance when it is combined with highly efficient linear solvers that are tailored for its special requirement. PMID:24723843
The novel high-performance 3-D MT inverse solver
NASA Astrophysics Data System (ADS)
Kruglyakov, Mikhail; Geraskin, Alexey; Kuvshinov, Alexey
2016-04-01
We present novel, robust, scalable, and fast 3-D magnetotelluric (MT) inverse solver. The solver is written in multi-language paradigm to make it as efficient, readable and maintainable as possible. Separation of concerns and single responsibility concepts go through implementation of the solver. As a forward modelling engine a modern scalable solver extrEMe, based on contracting integral equation approach, is used. Iterative gradient-type (quasi-Newton) optimization scheme is invoked to search for (regularized) inverse problem solution, and adjoint source approach is used to calculate efficiently the gradient of the misfit. The inverse solver is able to deal with highly detailed and contrasting models, allows for working (separately or jointly) with any type of MT responses, and supports massive parallelization. Moreover, different parallelization strategies implemented in the code allow optimal usage of available computational resources for a given problem statement. To parameterize an inverse domain the so-called mask parameterization is implemented, which means that one can merge any subset of forward modelling cells in order to account for (usually) irregular distribution of observation sites. We report results of 3-D numerical experiments aimed at analysing the robustness, performance and scalability of the code. In particular, our computational experiments carried out at different platforms ranging from modern laptops to HPC Piz Daint (6th supercomputer in the world) demonstrate practically linear scalability of the code up to thousands of nodes.
Adaptive kinetic-fluid solvers for heterogeneous computing architectures
NASA Astrophysics Data System (ADS)
Zabelok, Sergey; Arslanbekov, Robert; Kolobov, Vladimir
2015-12-01
We show feasibility and benefits of porting an adaptive multi-scale kinetic-fluid code to CPU-GPU systems. Challenges are due to the irregular data access for adaptive Cartesian mesh, vast difference of computational cost between kinetic and fluid cells, and desire to evenly load all CPUs and GPUs during grid adaptation and algorithm refinement. Our Unified Flow Solver (UFS) combines Adaptive Mesh Refinement (AMR) with automatic cell-by-cell selection of kinetic or fluid solvers based on continuum breakdown criteria. Using GPUs enables hybrid simulations of mixed rarefied-continuum flows with a million of Boltzmann cells each having a 24 × 24 × 24 velocity mesh. We describe the implementation of CUDA kernels for three modules in UFS: the direct Boltzmann solver using the discrete velocity method (DVM), the Direct Simulation Monte Carlo (DSMC) solver, and a mesoscopic solver based on the Lattice Boltzmann Method (LBM), all using adaptive Cartesian mesh. Double digit speedups on single GPU and good scaling for multi-GPUs have been demonstrated.
General Equation Set Solver for Compressible and Incompressible Turbomachinery Flows
NASA Technical Reports Server (NTRS)
Sondak, Douglas L.; Dorney, Daniel J.
2002-01-01
Turbomachines for propulsion applications operate with many different working fluids and flow conditions. The flow may be incompressible, such as in the liquid hydrogen pump in a rocket engine, or supersonic, such as in the turbine which may drive the hydrogen pump. Separate codes have traditionally been used for incompressible and compressible flow solvers. The General Equation Set (GES) method can be used to solve both incompressible and compressible flows, and it is not restricted to perfect gases, as are many compressible-flow turbomachinery solvers. An unsteady GES turbomachinery flow solver has been developed and applied to both air and water flows through turbines. It has been shown to be an excellent alternative to maintaining two separate codes.
A multiple right hand side iterative solver for history matching
Killough, J.E.; Sharma, Y.; Dupuy, A.; Bissell, R.; Wallis, J.
1995-12-31
History matching of oil and gas reservoirs can be accelerated by directly calculating the gradients of observed quantities (e.g., well pressure) with respect to the adjustable reserve parameters (e.g., permeability). This leads to a set of linear equations which add a significant overhead to the full simulation run without gradients. Direct Gauss elimination solvers can be used to address this problem by performing the factorization of the matrix only once and then reusing the factor matrix for the solution of the multiple right hand sides. This is a limited technique, however. Experience has shown that problems with greater than few thousand cells may not be practical for direct solvers because of computation time and memory limitations. This paper discusses the implementation of a multiple right hand side iterative linear equation solver (MRHS) for a system of adjoint equations to significantly enhance the performance of a gradient simulator.
Gpu Implementation of a Viscous Flow Solver on Unstructured Grids
NASA Astrophysics Data System (ADS)
Xu, Tianhao; Chen, Long
2016-06-01
Graphics processing units have gained popularities in scientific computing over past several years due to their outstanding parallel computing capability. Computational fluid dynamics applications involve large amounts of calculations, therefore a latest GPU card is preferable of which the peak computing performance and memory bandwidth are much better than a contemporary high-end CPU. We herein focus on the detailed implementation of our GPU targeting Reynolds-averaged Navier-Stokes equations solver based on finite-volume method. The solver employs a vertex-centered scheme on unstructured grids for the sake of being capable of handling complex topologies. Multiple optimizations are carried out to improve the memory accessing performance and kernel utilization. Both steady and unsteady flow simulation cases are carried out using explicit Runge-Kutta scheme. The solver with GPU acceleration in this paper is demonstrated to have competitive advantages over the CPU targeting one.
Two Solvers for Tractable Temporal Constraints with Preferences
NASA Technical Reports Server (NTRS)
Rossi, F.; Khatib,L.; Morris, P.; Morris, R.; Clancy, Daniel (Technical Monitor)
2002-01-01
A number of reasoning problems involving the manipulation of temporal information can naturally be viewed as implicitly inducing an ordering of potential local decisions involving time on the basis of preferences. Soft temporal constraints problems allow to describe in a natural way scenarios where events happen over time and preferences are associated to event distances and durations. In general, solving soft temporal problems require exponential time in the worst case, but there are interesting subclasses of problems which are polynomially solvable. We describe two solvers based on two different approaches for solving the same tractable subclass. For each solver we present the theoretical results it stands on, a description of the algorithm and some experimental results. The random generator used to build the problems on which tests are performed is also described. Finally, we compare the two solvers highlighting the tradeoff between performance and representational power.
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.
Stability analysis of Eulerian-Lagrangian methods for the one-dimensional shallow-water equations
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.
An Arbitrary Lagrangian-Eulerian Discretization of MHD on 3D Unstructured Grids
Rieben, R N; White, D A; Wallin, B K; Solberg, J M
2006-06-12
We present an arbitrary Lagrangian-Eulerian (ALE) discretization of the equations of resistive magnetohydrodynamics (MHD) on unstructured hexahedral grids. The method is formulated using an operator-split approach with three distinct phases: electromagnetic diffusion, Lagrangian motion, and Eulerian advection. The resistive magnetic dynamo equation is discretized using a compatible mixed finite element method with a 2nd order accurate implicit time differencing scheme which preserves the divergence-free nature of the magnetic field. At each discrete time step, electromagnetic force and heat terms are calculated and coupled to the hydrodynamic equations to compute the Lagrangian motion of the conducting materials. By virtue of the compatible discretization method used, the invariants of Lagrangian MHD motion are preserved in a discrete sense. When the Lagrangian motion of the mesh causes significant distortion, that distortion is corrected with a relaxation of the mesh, followed by a 2nd order monotonic remap of the electromagnetic state variables. The remap is equivalent to Eulerian advection of the magnetic flux density with a fictitious mesh relaxation velocity. The magnetic advection is performed using a novel variant of constrained transport (CT) that is valid for unstructured hexahedral grids with arbitrary mesh velocities. The advection method maintains the divergence free nature of the magnetic field and is second order accurate in regions where the solution is sufficiently smooth. For regions in which the magnetic field is discontinuous (e.g. MHD shocks) the method is limited using a novel variant of algebraic flux correction (AFC) which is local extremum diminishing (LED) and divergence preserving. Finally, we verify each stage of the discretization via a set of numerical experiments.
Lagrangian and arbitrary Lagrangian Eulerian simulations of complex roll-forming processes
NASA Astrophysics Data System (ADS)
Crutzen, Yanick; Boman, Romain; Papeleux, Luc; Ponthot, Jean-Philippe
2016-04-01
The Arbitrary Lagrangian Eulerian (ALE) formalism is a breakthrough technique in the numerical simulation of the continuous-type roll-forming process. In contrast to the classical Lagrangian approach, the ALE formalism can compute the hopefully stationary state for the entire mill length with definitely effortless set-up tasks thanks to a nearly-stationary mesh. In this paper, advantages of ALE and Lagrangian formalisms are extensively discussed for simulating such continuous-type processes. Through a highly complex industrial application, the ease of use of ALE modelling is illustrated with the in-house code METAFOR. ALE and Lagrangian results are in good agreement with each other.
A High-order Eulerian-Lagrangian Finite Element Method for Coupled Electro-mechanical Systems
NASA Astrophysics Data System (ADS)
Brandstetter, Gerd
The main focus of this work is on the development of a high-order Eulerian-Lagrangian finite element method for the simulation of electro-mechanical systems. The coupled problem is solved by a staggered scheme, where the mechanical motion is discretized by standard Lagrangian finite elements, and the electrical field is solved on a fixed Eulerian grid with embedded boundary conditions. Traditional Lagrangian-Lagrangian or arbitrary Lagrangian-Eulerian (ALE) methods encounter deficiencies, for example, when dealing with mesh distortion due to large deformations, or topology changes due to contacting bodies. The presented Eulerian-Lagrangian approach addresses these issues in a natural way. Within this context we develop a high-order immersed boundary discontinuous-Galerkin (IB-DG) method, which is shown to be necessary for (i) the accurate representation of the electrical gradient along nonlinear boundary features such as singular corners, and (ii) to achieve full convergence during the iterative global solution. We develop an implicit scheme based on the mid-point rule, as well as an explicit scheme based on the centered-difference method, with the incorporation of energy conserving, frictionless contact algorithms for an elastic-to-rigid-surface contact. The performance of the proposed method is assessed for several benchmark tests: the electro-static force vector around a singular corner, the quasi-static pull-in of an electro-mechanically actuated switch, the excitation of a carbon nanotube at resonance, and the cyclic impact simulation of a micro-electro-mechanical resonant-switch. We report improved accuracy for the high-order method as compared to low-order methods, and linear convergence in the iterative solution of the staggered scheme. Additionally, we investigate a Newton-Krylov shooting scheme in order to directly find cyclic steady states of electro-mechanical devices excited at resonance-- as opposed to a naive time-stepping from zero initial
Arbitrary Lagrangian-Eulerian approach in reduced order modeling of a flow with a moving boundary
NASA Astrophysics Data System (ADS)
Stankiewicz, W.; Roszak, R.; Morzyński, M.
2013-06-01
Flow-induced deflections of aircraft structures result in oscillations that might turn into such a dangerous phenomena like flutter or buffeting. In this paper the design of an aeroelastic system consisting of Reduced Order Model (ROM) of the flow with a moving boundary is presented. The model is based on Galerkin projection of governing equation onto space spanned by modes obtained from high-fidelity computations. The motion of the boundary and mesh is defined in Arbitrary Lagrangian-Eulerian (ALE) approach and results in additional convective term in Galerkin system. The developed system is demonstrated on the example of a flow around an oscillating wing.
Modeling erosion in a centrifugal pump in an Eulerian-Lagrangian frame using OpenFOAM®
NASA Astrophysics Data System (ADS)
Lopez, Alejandro; Stickland, Matthew; Dempster, William
2015-07-01
Erosion induced by solid particle impingement is a very commonwear mechanism in turbomachinery and Computational Fluid Dynamics is one of the most widely used tools for its prediction. In this article, erosion is modeled in one of the channels of a centrifugal pump using OpenFOAM®,which is an Open Source CFD package. A review of some of the most commonly used erosion models is carried out in an Eulerian-Lagrangian frame along with a comparative study of the erosion rates obtained with each model. Results yielded some disparities between models due to the different factors taken into consideration. The mesh is then deformed to obtain the resulting eroded geometry.
Fast Euler solver for transonic airfoils. I - Theory. II - Applications
NASA Technical Reports Server (NTRS)
Dadone, Andrea; Moretti, Gino
1988-01-01
Equations written in terms of generalized Riemann variables are presently integrated by inverting six bidiagonal matrices and two tridiagonal matrices, using an implicit Euler solver that is based on the lambda-formulation. The solution is found on a C-grid whose boundaries are very close to the airfoil. The fast solver is then applied to the computation of several flowfields on a NACA 0012 airfoil at various Mach number and alpha values, yielding results that are primarily concerned with transonic flows. The effects of grid fineness and boundary distances are analyzed; the code is found to be robust and accurate, as well as fast.
Numerical System Solver Developed for the National Cycle Program
NASA Technical Reports Server (NTRS)
Binder, Michael P.
1999-01-01
As part of the National Cycle Program (NCP), a powerful new numerical solver has been developed to support the simulation of aeropropulsion systems. This software uses a hierarchical object-oriented design. It can provide steady-state and time-dependent solutions to nonlinear and even discontinuous problems typically encountered when aircraft and spacecraft propulsion systems are simulated. It also can handle constrained solutions, in which one or more factors may limit the behavior of the engine system. Timedependent simulation capabilities include adaptive time-stepping and synchronization with digital control elements. The NCP solver is playing an important role in making the NCP a flexible, powerful, and reliable simulation package.
NASA Astrophysics Data System (ADS)
Cerminara, Matteo; Esposti Ongaro, Tomaso; Carlo Berselli, Luigi
2014-05-01
We have developed a compressible multiphase flow model to simulate the three-dimensional dynamics of turbulent volcanic ash plumes. The model describes the eruptive mixture as a polydisperse fluid, composed of different types of gases and particles, treated as interpenetrating Eulerian phases. Solid phases represent the discrete ash classes into which the total granulometric spectrum is discretized, and can differ by size and density. The model is designed to quickly and accurately resolve important physical phenomena in the dynamics of volcanic ash plumes. In particular, it can simulate turbulent mixing (driving atmospheric entrainment and controlling the heat transfer), thermal expansion (controlling the plume buoyancy), the interaction between solid particles and volcanic gas (including kinetic non-equilibrium effects) and the effects of compressibility (over-pressured eruptions and infrasonic measurements). The model is based on the turbulent dispersed multiphase flow theory for dilute flows (volume concentration <0.001, implying that averaged inter-particle distance is larger than 10 diameters) where particle collisions are neglected. Moreover, in order to speed up the code without losing accuracy, we make the hypothesis of fine particles (Stokes number <0.2 , i.e., volcanic ash particles finer then a millimeter), so that we are able to consider non-equilibrium effects only at the first order. We adopt LES formalism (which is preferable in transient regimes) for compressible flows to model the non-linear coupling between turbulent scales and the effect of sub-grid turbulence on the large-scale dynamics. A three-dimensional numerical code has been developed basing on the OpenFOAM computational framework, a CFD open source parallel software package. Numerical benchmarks demonstrate that the model is able to capture important non-equilibrium phenomena in gas-particle mixtures, such as particle clustering and ejection from large-eddy turbulent structures, as well
NASA Astrophysics Data System (ADS)
Ergun, Sule
In the case of a postulated loss of coolant accident (LOCA) in a nuclear reactor, an accurate prediction of clad temperature is needed to determine the safety margins. The large break LOCA analyses can be divided in to three time periods. These periods are blowdown, refill and reflood. During the blowdown and reflood phases of the LOCA, when the local void fraction is greater than 80% and the wall is at a temperature above minimum film boiling temperature (Tmin), heat is transferred from the fuel rod to a continuous vapor flow with dispersed droplets. The high void fraction mixture of droplets and vapor provide cooling to prevent the clad temperature from exceeding the safety limit. The heat transfer process for high void fraction mixture is called dispersed flow film boiling (DFFB). This thesis has been modeled DFFB in the reflood phase of a LOCA in a pressurized water reactor (PWR) rod bundle. In this study, the modifications and modification requirements for the COBRA-TF code to obtain a five field Eulerian - Eulerian modeling for two-phase DFFB is described. COBRA-TF is a best estimate code developed for the rod bundle analysis and has four fields, namely, vapor, entrained drop and continuous liquid film. COBRA-TF has a detailed reflood package which takes effect of spacer grids on heat transfer into account. This study has a detailed description of code's solution scheme and the models used for dispersed flow film boiling. The dispersed flow film boiling heat transfer model of the COBRA-TF code has been modified by adding a small droplet field to the code as the fifth field. The effect of smaller, thermally more active droplets on heat, mass and momentum transfer during DFFB has been modeled. Since the large drop break up due to spacer grids is one of the reasons for small droplet generation, the spacer grid models of the COBRA-TF have been revised and modified. In addition to small droplet generation, the spacer grid rewet is an important aspect of heat
Li, Xiantao; Wöhlbier, John G; Jin, Shi; Booske, John H
2004-01-01
We provide methods of computing multivalued solutions to the Euler-Poisson system and test them in the context of a klystron amplifier. An Eulerian formulation capable of computing multivalued solutions is derived from a kinetic description of the Euler-Poisson system and a moment closure. The system of the moment equations may be closed due to the special structure of the solution in phase space. The Eulerian moment equations are computed for a velocity modulated electron beam, which has been shown by prior Lagrangian theories to break in a finite time and form multivalued solutions. The results of the Eulerian moment equations are compared to direct computation of the kinetic equations and a Lagrangian method also developed in the paper. We use the Lagrangian formulation for the explicit computation of wave breaking time and location for typical velocity modulation boundary conditions. PMID:15324179
Lee, T.; Leaf, G. K.; Mathematics and Computer Science; The City Univ. of New York
2009-04-01
We propose an Eulerian description of the bounce-back boundary condition based on the high-order implicit time-marching schemes to improve the accuracy of lattice Boltzmann simulation in the vicinity of curved boundary. The Eulerian description requires only one grid spacing between fluid nodes when second-order accuracy in time and space is desired, although high-order accurate boundary conditions can be constructed on more grid-point support. The Eulerian description also provides an analytical framework for several different interpolation-based boundary conditions. For instance, the semi-Lagrangian, linear interpolation boundary condition is found to be a first-order upwind discretization that changes the time-marching schemes from implicit to explicit as the distance between the fluid boundary node and the solid boundary increases.
Analytical Model for Pair Dispersion in Gaussian Models of Eulerian Turbulence
NASA Astrophysics Data System (ADS)
Eyink, Gregory; Benveniste, Damien
2012-11-01
Synthetic models of Eulerian turbulence are often used as computational shortcuts for studying Lagrangian properties of turbulence (e.g. Elliott & Majda, 1996). These models have been criticized by Thomson & Devenish (2005), who argued on physical grounds that their sweeping effects are very different from true turbulence. We give analytical results for Eulerian turbulence modeled by Gaussian fields. Our starting point is an exact integrodifferential equation for the particle pair separation distribution obtained from Gaussian integration-by-parts. When velocity correlation times are short, a Markovian approximation leads to a Richardson-type diffusion model. We obtain a time-dependent pair diffusivity tensor of the form Kij (r , t) =Sij (r) τ (r , t) where Sij (r) is the structure-function tensor and τ (r , t) is an effective correlation time of velocity increments. Crucially, this is found to be the minimum value of three times: the intrinsic turnover time τeddy (r) at separation r, the overall evolution time t , and the sweeping time r /v0 with v0 the rms velocity. We thus verify the main argument of Thomson & Devenish (2005), but we predict scaling laws for pair dispersion different from theirs for zero-mean velocity ensembles.
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.
NASA Astrophysics Data System (ADS)
Smith, Scott C.; Houser, Janet L.; Centrella, Joan M.
1996-02-01
We have carried out three-dimensional numerical simulations of the dynamical bar instability in a rotating star and the resulting gravitational radiation using both an Eulerian code written in cylindrical coordinates and a smooth particle hydrodynamics (SPH) code. The star is modeled initially as a polytrope with index n = 3/2 and Trot|W| ≍ 0.30, where Trot is the rotational kinetic energy and |W| is the gravitational potential energy. In both codes the gravitational field is purely Newtonian, and the gravitational radiation is calculated in the quadrupole approximation. We have run three simulations with the Eulerian code, varying the number of angular zones and the treatment of the boundary between the star and the vacuum. Using the SPH code we did seven runs, varying the number of particles, the artificial viscosity, and the type of initial model. We compare the growth rate and rotation speed of the bar, the mass, and angular momentum distributions, and the gravitational radiation quantities. We highlight the successes and difficulties of both methods and make suggestions for future improvements.
Candy, J.; Waltz, R.E.
2006-03-15
Equations which describe the evolution of volume-averaged gyrokinetic entropy are derived and added to GYRO [J. Candy and R.E. Waltz, J. Comput. Phys. 186, 545 (2003)], a Eulerian gyrokinetic turbulence simulation code. In particular, the creation of entropy through spatial upwind dissipation (there is zero velocity-space dissipation in GYRO) and the reduction of entropy via the production of fluctuations are monitored in detail. This new diagnostic has yielded several key confirmations of the validity of the GYRO simulations. First, fluctuations balance dissipation in the ensemble-averaged sense, thus demonstrating that turbulent GYRO simulations achieve a true statistical steady state. Second, at the standard spatial grid size, neither entropy nor energy flux is significantly changed by a 16-fold increase (from 32 to 512 grid points per cell) in the number of grid points in the two-dimensional velocity space. Third, the measured flux is invariant to an eightfold increase in the upwind dissipation coefficients. A notable conclusion is that the lack of change in entropy with grid refinement refutes the familiar but incorrect notion that Eulerian gyrokinetic codes miss important velocity-space structure. The issues of density and energy conservation and their relation to negligible second-order effects such as the parallel nonlinearity are also discussed.
Reuge, N; Cadoret, L.; Pannala, Sreekanth; Syamlal, M; Coufort, C; Caussat, B
2008-01-01
Computational fluid dynamic (CFD) models must be thoroughly validated before they can be used with confidence for designing fluidized bed reactors. In this study, validation data were collected from a fluidized bed of (Geldart's group B) alumina particles operated at different gas velocities involving two fluidization hydrodynamic regimes (bubbling and slugging). The bed expansion, height of bed fluctuations, and frequency of fluctuations were measured from a videos of the fluidized bed. The Eulerian-Eulerian two fluid model MFIX was then used to simulate the experiments. Two different models for the particle stresses - Schaeffer (Syamlal et al., (1993), Schaeffer (1987)) and Princeton (Srivastava and Sundaresan (2003)) models - and different values of the restitution coefficient and internal angle of friction were evaluated. 3-D simulations are required for getting quantitative and qualitative agreement with experimental data. The results from the Princeton model are in better agreement with data than from the Schaeffer model. Both free-slip and Johnson-Jackson boundary conditions give nearly identical results. An increase in e from 0.8 to 1 leads to larger bed expansions and lower heights of fluctuations in the bubbling regime whereas it leads to unchanged bed expansion and to a massive reduction in the height of fluctuations in the slugging regime. The angle of internal friction (φ) in the range 10 -40 does not affect the bed expansion, but its reduction significantly reduces the height of fluctuations.
Coordinate Projection-based Solver for ODE with Invariants
Serban, Radu
2008-04-08
CPODES is a general purpose (serial and parallel) solver for systems of ordinary differential equation (ODE) with invariants. It implements a coordinate projection approach using different types of projection (orthogonal or oblique) and one of several methods for the decompositon of the Jacobian of the invariant equations.
Navier-Stokes Solvers and Generalizations for Reacting Flow Problems
Elman, Howard C
2013-01-27
This is an overview of our accomplishments during the final term of this grant (1 September 2008 -- 30 June 2012). These fall mainly into three categories: fast algorithms for linear eigenvalue problems; solution algorithms and modeling methods for partial differential equations with uncertain coefficients; and preconditioning methods and solvers for models of computational fluid dynamics (CFD).
Assessment of linear finite-difference Poisson-Boltzmann solvers.
Wang, Jun; Luo, Ray
2010-06-01
CPU time and memory usage are two vital issues that any numerical solvers for the Poisson-Boltzmann equation have to face in biomolecular applications. In this study, we systematically analyzed the CPU time and memory usage of five commonly used finite-difference solvers with a large and diversified set of biomolecular structures. Our comparative analysis shows that modified incomplete Cholesky conjugate gradient and geometric multigrid are the most efficient in the diversified test set. For the two efficient solvers, our test shows that their CPU times increase approximately linearly with the numbers of grids. Their CPU times also increase almost linearly with the negative logarithm of the convergence criterion at very similar rate. Our comparison further shows that geometric multigrid performs better in the large set of tested biomolecules. However, modified incomplete Cholesky conjugate gradient is superior to geometric multigrid in molecular dynamics simulations of tested molecules. We also investigated other significant components in numerical solutions of the Poisson-Boltzmann equation. It turns out that the time-limiting step is the free boundary condition setup for the linear systems for the selected proteins if the electrostatic focusing is not used. Thus, development of future numerical solvers for the Poisson-Boltzmann equation should balance all aspects of the numerical procedures in realistic biomolecular applications. PMID:20063271
Assessment of Linear Finite-Difference Poisson-Boltzmann Solvers
Wang, Jun; Luo, Ray
2009-01-01
CPU time and memory usage are two vital issues that any numerical solvers for the Poisson-Boltzmann equation have to face in biomolecular applications. In this study we systematically analyzed the CPU time and memory usage of five commonly used finite-difference solvers with a large and diversified set of biomolecular structures. Our comparative analysis shows that modified incomplete Cholesky conjugate gradient and geometric multigrid are the most efficient in the diversified test set. For the two efficient solvers, our test shows that their CPU times increase approximately linearly with the numbers of grids. Their CPU times also increase almost linearly with the negative logarithm of the convergence criterion at very similar rate. Our comparison further shows that geometric multigrid performs better in the large set of tested biomolecules. However, modified incomplete Cholesky conjugate gradient is superior to geometric multigrid in molecular dynamics simulations of tested molecules. We also investigated other significant components in numerical solutions of the Poisson-Boltzmann equation. It turns out that the time-limiting step is the free boundary condition setup for the linear systems for the selected proteins if the electrostatic focusing is not used. Thus, development of future numerical solvers for the Poisson-Boltzmann equation should balance all aspects of the numerical procedures in realistic biomolecular applications. PMID:20063271
Advances in computational fluid dynamics solvers for modern computing environments
NASA Astrophysics Data System (ADS)
Hertenstein, Daniel; Humphrey, John R.; Paolini, Aaron L.; Kelmelis, Eric J.
2013-05-01
EM Photonics has been investigating the application of massively multicore processors to a key problem area: Computational Fluid Dynamics (CFD). While the capabilities of CFD solvers have continually increased and improved to support features such as moving bodies and adjoint-based mesh adaptation, the software architecture has often lagged behind. This has led to poor scaling as core counts reach the tens of thousands. In the modern High Performance Computing (HPC) world, clusters with hundreds of thousands of cores are becoming the standard. In addition, accelerator devices such as NVIDIA GPUs and Intel Xeon Phi are being installed in many new systems. It is important for CFD solvers to take advantage of the new hardware as the computations involved are well suited for the massively multicore architecture. In our work, we demonstrate that new features in NVIDIA GPUs are able to empower existing CFD solvers by example using AVUS, a CFD solver developed by the Air Force Research Labratory (AFRL) and the Volcanic Ash Advisory Center (VAAC). The effort has resulted in increased performance and scalability without sacrificing accuracy. There are many well-known codes in the CFD space that can benefit from this work, such as FUN3D, OVERFLOW, and TetrUSS. Such codes are widely used in the commercial, government, and defense sectors.
Intellectual Abilities That Discriminate Good and Poor Problem Solvers.
ERIC Educational Resources Information Center
Meyer, Ruth Ann
1981-01-01
This study compared good and poor fourth-grade problem solvers on a battery of 19 "reference" tests for verbal, induction, numerical, word fluency, memory, perceptual speed, and simple visualization abilities. Results suggest verbal, numerical, and especially induction abilities are important to successful mathematical problem solving. (MP)
NASA Astrophysics Data System (ADS)
Gotovac, Hrvoje; Srzic, Veljko
2014-05-01
Contaminant transport in natural aquifers is a complex, multiscale process that is frequently studied using different Eulerian, Lagrangian and hybrid numerical methods. Conservative solute transport is typically modeled using the advection-dispersion equation (ADE). Despite the large number of available numerical methods that have been developed to solve it, the accurate numerical solution of the ADE still presents formidable challenges. In particular, current numerical solutions of multidimensional advection-dominated transport in non-uniform velocity fields are affected by one or all of the following problems: numerical dispersion that introduces artificial mixing and dilution, grid orientation effects, unresolved spatial and temporal scales and unphysical numerical oscillations (e.g., Herrera et al, 2009; Bosso et al., 2012). In this work we will present Eulerian Lagrangian Adaptive Fup Collocation Method (ELAFCM) based on Fup basis functions and collocation approach for spatial approximation and explicit stabilized Runge-Kutta-Chebyshev temporal integration (public domain routine SERK2) which is especially well suited for stiff parabolic problems. Spatial adaptive strategy is based on Fup basis functions which are closely related to the wavelets and splines so that they are also compactly supported basis functions; they exactly describe algebraic polynomials and enable a multiresolution adaptive analysis (MRA). MRA is here performed via Fup Collocation Transform (FCT) so that at each time step concentration solution is decomposed using only a few significant Fup basis functions on adaptive collocation grid with appropriate scales (frequencies) and locations, a desired level of accuracy and a near minimum computational cost. FCT adds more collocations points and higher resolution levels only in sensitive zones with sharp concentration gradients, fronts and/or narrow transition zones. According to the our recent achievements there is no need for solving the large
Multiscale Universal Interface: A concurrent framework for coupling heterogeneous solvers
Tang, Yu-Hang; Kudo, Shuhei; Bian, Xin; Li, Zhen; Karniadakis, George Em
2015-09-15
Graphical abstract: - Abstract: Concurrently coupled numerical simulations using heterogeneous solvers are powerful tools for modeling multiscale phenomena. However, major modifications to existing codes are often required to enable such simulations, posing significant difficulties in practice. In this paper we present a C++ library, i.e. the Multiscale Universal Interface (MUI), which is capable of facilitating the coupling effort for a wide range of multiscale simulations. The library adopts a header-only form with minimal external dependency and hence can be easily dropped into existing codes. A data sampler concept is introduced, combined with a hybrid dynamic/static typing mechanism, to create an easily customizable framework for solver-independent data interpretation. The library integrates MPI MPMD support and an asynchronous communication protocol to handle inter-solver information exchange irrespective of the solvers' own MPI awareness. Template metaprogramming is heavily employed to simultaneously improve runtime performance and code flexibility. We validated the library by solving three different multiscale problems, which also serve to demonstrate the flexibility of the framework in handling heterogeneous models and solvers. In the first example, a Couette flow was simulated using two concurrently coupled Smoothed Particle Hydrodynamics (SPH) simulations of different spatial resolutions. In the second example, we coupled the deterministic SPH method with the stochastic Dissipative Particle Dynamics (DPD) method to study the effect of surface grafting on the hydrodynamics properties on the surface. In the third example, we consider conjugate heat transfer between a solid domain and a fluid domain by coupling the particle-based energy-conserving DPD (eDPD) method with the Finite Element Method (FEM)
Migration of vectorized iterative solvers to distributed memory architectures
Pommerell, C.; Ruehl, R.
1994-12-31
Both necessity and opportunity motivate the use of high-performance computers for iterative linear solvers. Necessity results from the size of the problems being solved-smaller problems are often better handled by direct methods. Opportunity arises from the formulation of the iterative methods in terms of simple linear algebra operations, even if this {open_quote}natural{close_quotes} parallelism is not easy to exploit in irregularly structured sparse matrices and with good preconditioners. As a result, high-performance implementations of iterative solvers have attracted a lot of interest in recent years. Most efforts are geared to vectorize or parallelize the dominating operation-structured or unstructured sparse matrix-vector multiplication, or to increase locality and parallelism by reformulating the algorithm-reducing global synchronization in inner products or local data exchange in preconditioners. Target architectures for iterative solvers currently include mostly vector supercomputers and architectures with one or few optimized (e.g., super-scalar and/or super-pipelined RISC) processors and hierarchical memory systems. More recently, parallel computers with physically distributed memory and a better price/performance ratio have been offered by vendors as a very interesting alternative to vector supercomputers. However, programming comfort on such distributed memory parallel processors (DMPPs) still lags behind. Here the authors are concerned with iterative solvers and their changing computing environment. In particular, they are considering migration from traditional vector supercomputers to DMPPs. Application requirements force one to use flexible and portable libraries. They want to extend the portability of iterative solvers rather than reimplementing everything for each new machine, or even for each new architecture.