skip to main content
OSTI.GOV title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: A 3D finite element ALE method using an approximate Riemann solution

Abstract

Arbitrary Lagrangian–Eulerian finite volume methods that solve a multidimensional Riemann-like problem at the cell center in a staggered grid hydrodynamic (SGH) arrangement have been proposed. This research proposes a new 3D finite element arbitrary Lagrangian–Eulerian SGH method that incorporates a multidimensional Riemann-like problem. Here, two different Riemann jump relations are investigated. A new limiting method that greatly improves the accuracy of the SGH method on isentropic flows is investigated. A remap method that improves upon a well-known mesh relaxation and remapping technique in order to ensure total energy conservation during the remap is also presented. Numerical details and test problem results are presented.

Authors:
 [1];  [1]
  1. Los Alamos National Lab. (LANL), Los Alamos, NM (United States)
Publication Date:
Research Org.:
Los Alamos National Lab. (LANL), Los Alamos, NM (United States)
Sponsoring Org.:
USDOE
OSTI Identifier:
1304832
Report Number(s):
LA-UR-16-20312
Journal ID: ISSN 0271-2091
Grant/Contract Number:
AC52-06NA25396
Resource Type:
Journal Article: Accepted Manuscript
Journal Name:
International Journal for Numerical Methods in Fluids
Additional Journal Information:
Journal Name: International Journal for Numerical Methods in Fluids; Journal ID: ISSN 0271-2091
Publisher:
Wiley
Country of Publication:
United States
Language:
English
Subject:
97 MATHEMATICS AND COMPUTING

Citation Formats

Chiravalle, V. P., and Morgan, N. R. A 3D finite element ALE method using an approximate Riemann solution. United States: N. p., 2016. Web. doi:10.1002/fld.4284.
Chiravalle, V. P., & Morgan, N. R. A 3D finite element ALE method using an approximate Riemann solution. United States. doi:10.1002/fld.4284.
Chiravalle, V. P., and Morgan, N. R. 2016. "A 3D finite element ALE method using an approximate Riemann solution". United States. doi:10.1002/fld.4284. https://www.osti.gov/servlets/purl/1304832.
@article{osti_1304832,
title = {A 3D finite element ALE method using an approximate Riemann solution},
author = {Chiravalle, V. P. and Morgan, N. R.},
abstractNote = {Arbitrary Lagrangian–Eulerian finite volume methods that solve a multidimensional Riemann-like problem at the cell center in a staggered grid hydrodynamic (SGH) arrangement have been proposed. This research proposes a new 3D finite element arbitrary Lagrangian–Eulerian SGH method that incorporates a multidimensional Riemann-like problem. Here, two different Riemann jump relations are investigated. A new limiting method that greatly improves the accuracy of the SGH method on isentropic flows is investigated. A remap method that improves upon a well-known mesh relaxation and remapping technique in order to ensure total energy conservation during the remap is also presented. Numerical details and test problem results are presented.},
doi = {10.1002/fld.4284},
journal = {International Journal for Numerical Methods in Fluids},
number = ,
volume = ,
place = {United States},
year = 2016,
month = 8
}

Journal Article:
Free Publicly Available Full Text
Publisher's Version of Record

Save / Share:
  • Coulomb Sturmians are used to calculate the ground state of H{sub 2}{sup +}. The required multicenter integrals are calculated to very high accuracy in momentum space. Differential cross sections describing electron scattering from the ground state of H{sub 2}{sup +} are then calculated using the finite element method. The solution of the resulting linear algebraic problem is accomplished using the iterative-variational method.
  • This paper presents a technique for solving inverse problems in magnetostatic nondestructive evaluation (NDE), using finite element models. In-line inspection of ferromagnetic gas pipelines containing pipe-wall defects, is chosen as the candidate NDE process. The signal inversion technique consists of iteratively solving the forward problem by updating the finite element mesh, rather than the material properties of the finite elements. Preliminary simulation results obtained using a 2D finite element model are presented.
  • The problem of magnetohydrodynamic flow through channels has become important because of several engineering applications such as design of nuclear reactor cooling systems, electromagnetic pumps, MHD flowmeters, MHD generators, blood flow measurements, etc. A numerical code based on the Finite Element Method (FEM) was developed to solve the two-dimensional, steady-state magnetohydrodynamic (MHD) flow in a rectangular channel. In order to apply the FEM, the Galerkin Weak Formulation was used. In this analysis, in contrast with the previous works, the thickness of the duct wall is taken into account and the results are compared to those obtained in the limit casemore » when the thickness is much smaller than a characteristic dimension of the duct. In this case, convergence behavior of several iterative methods, for high Hartmann numbers, was also investigated.« less
  • We introduce an approximate Roe type Riemann solver for the numerical simulation of two-phase fluid flows composed of liquid droplets suspended in gas. We compute a Roe linearization of some well-conditioned approximate Rankine-Hugoniot relations in nonconservation form. The computed solutions are found to be in good agreement with the exact solution in one dimension slab geometry. We extend this solver to two-dimensional geometries using a fininte volume formulation. 24 refs., 15 figs., 2 tabs.
  • The authors briefly describe the solution of the zeroth and first moment equations of particle transport with Minerbo`s maximum entropy closure using a Riemann solver. This solver properly accounts for the direction and speed of propagation of information in this highly nonlinear hyperbolic system. While a direct iteration on the Eddington factor can become unstable and produce unphysical solutions, time-stepping to steady state with a Riemann solver has shown no such unphysical behavior. This establishes that the Minerbo closure is a useful two-moment transport approximation even in strong absorbers. It also establishes the utility of examining the hyperbolic structure ofmore » nonlinear closures and of Riemann solvers in radiation transport applications.« less