Design-order, non-conformal low-Mach fluid algorithms using a hybrid CVFEM/DG approach

doi:10.1016/j.jcp.2018.01.007

Title: Design-order, non-conformal low-Mach fluid algorithms using a hybrid CVFEM/DG approach

Abstract

Here, a hybrid, design-order sliding mesh algorithm, which uses a control volume finite element method (CVFEM), in conjunction with a discontinuous Galerkin (DG) approach at non-conformal interfaces, is outlined in the context of a low-Mach fluid dynamics equation set. This novel hybrid DG approach is also demonstrated to be compatible with a classic edge-based vertex centered (EBVC) scheme. For the CVFEM, element polynomial, P, promotion is used to extend the low-order P = 1 CVFEM method to higher-order, i.e., P = 2. An equal-order low-Mach pressure-stabilized methodology, with emphasis on the non-conformal interface boundary condition, is presented. A fully implicit matrix solver approach that accounts for the full stencil connectivity across the non-conformal interface is employed. A complete suite of formal verification studies using the method of manufactured solutions (MMS) is performed to verify the order of accuracy of the underlying methodology. The chosen suite of analytical verification cases range from a simple steady diffusion system to a traveling viscous vortex across mixed-order non-conformal interfaces. Results from all verification studies demonstrate either second- or third-order spatial accuracy and, for transient solutions, second-order temporal accuracy.

Authors:

Domino, Stefan P. ^[1]

Sandia National Lab. (SNL-NM), Albuquerque, NM (United States)

Publication Date:: 2018-01-12

Research Org.:: Sandia National Lab. (SNL-NM), Albuquerque, NM (United States)

Sponsoring Org.:: USDOE Office of Science (SC)

OSTI Identifier:: 1465806

Report Number(s):: SAND-2017-10354J
Journal ID: ISSN 0021-9991; 666503; TRN: US1902560

Grant/Contract Number:: AC04-94AL85000

Resource Type:: Journal Article: Accepted Manuscript

Journal Name:: Journal of Computational Physics

Additional Journal Information:: Journal Volume: 359; Journal Issue: C; Journal ID: ISSN 0021-9991

Publisher:: Elsevier

Country of Publication:: United States

Language:: English

Subject:: 71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSICS; Control volume finite element; Discontinuous Galerkin; Higher-order; Sliding mesh; Non-conformal

Citation Formats


                    Domino, Stefan P. Design-order, non-conformal low-Mach fluid algorithms using a hybrid CVFEM/DG approach.  United States: N. p., 2018. 
        Web.  doi:10.1016/j.jcp.2018.01.007.

Copy to clipboard


                    Domino, Stefan P. Design-order, non-conformal low-Mach fluid algorithms using a hybrid CVFEM/DG approach.  United States.  doi:10.1016/j.jcp.2018.01.007.

Copy to clipboard


                    Domino, Stefan P. Fri .  
        "Design-order, non-conformal low-Mach fluid algorithms using a hybrid CVFEM/DG approach".  United States.  doi:10.1016/j.jcp.2018.01.007.  https://www.osti.gov/servlets/purl/1465806.

Copy to clipboard


                    
@article{osti_1465806,

  title        = {Design-order, non-conformal low-Mach fluid algorithms using a hybrid CVFEM/DG approach},

  author       = {Domino, Stefan P.},

  abstractNote = {Here, a hybrid, design-order sliding mesh algorithm, which uses a control volume finite element method (CVFEM), in conjunction with a discontinuous Galerkin (DG) approach at non-conformal interfaces, is outlined in the context of a low-Mach fluid dynamics equation set. This novel hybrid DG approach is also demonstrated to be compatible with a classic edge-based vertex centered (EBVC) scheme. For the CVFEM, element polynomial, P, promotion is used to extend the low-order P = 1 CVFEM method to higher-order, i.e., P = 2. An equal-order low-Mach pressure-stabilized methodology, with emphasis on the non-conformal interface boundary condition, is presented. A fully implicit matrix solver approach that accounts for the full stencil connectivity across the non-conformal interface is employed. A complete suite of formal verification studies using the method of manufactured solutions (MMS) is performed to verify the order of accuracy of the underlying methodology. The chosen suite of analytical verification cases range from a simple steady diffusion system to a traveling viscous vortex across mixed-order non-conformal interfaces. Results from all verification studies demonstrate either second- or third-order spatial accuracy and, for transient solutions, second-order temporal accuracy.},

  doi          = {10.1016/j.jcp.2018.01.007},

  journal      = {Journal of Computational Physics},

  issn         = {0021-9991},

  number       = C,

  volume       = 359,

  place        = {United States},

  year         = {2018},

  month        = {1}

}

Copy to clipboard

Journal Article:

Free Publicly Available Full Text

Accepted Manuscript (DOE)

Publisher's Version of Record

DOI: 10.1016/j.jcp.2018.01.007

Other availability

Search WorldCat to find libraries that may hold this journal

Citation Metrics:

Cited by: 4 works

Citation information provided by
Web of Science

Figures / Tables:

Figure 1: Overview of the conceptual differences between non-conformal and overset methodologies. In the sliding mesh technology, (a), a set of surfaces are defined that form a nonconformal interface while in the overset method, a background mesh is intersected with an overset mesh to remove inactive elements, (b); overlap regionsmore »

All figures and tables (19 total)

Save / Share:

Export Metadata

Save to My Library

Figures / Tables found in this record:

Figures/Tables have been extracted from DOE-funded journal article accepted manuscripts.

Similar records in OSTI.GOV collections:

Deploy production sliding mesh capability with linear solver benchmarking.

Technical Report Domino, Stefan P. ; Thomas, Stephen ; Barone, Matthew F. ; ...

Wind applications require the ability to simulate rotating blades. To support this use-case, a novel design-order sliding mesh algorithm has been developed and deployed. The hybrid method combines the control volume finite element methodology (CVFEM) with concepts found within a discontinuous Galerkin (DG) finite element method (FEM) to manage a sliding mesh. The method has been demonstrated to be design-order for the tested polynomial basis (P=1 and P=2) and has been deployed to provide production simulation capability for a Vestas V27 (225 kW) wind turbine. Other stationary and canonical rotating ow simulations are also presented. As the majority of wind-energymore »« less
DOI: 10.2172/1423525

Full Text Available
High-order fully implicit solver for all-speed fluid dynamics: AUSM ride from nearly incompressible variable-density flows to shock dynamics

Journal Article Nourgaliev, R. ; Greene, P. ; Weston, B. ; ... - Shock Waves

We introduce a novel Newton–Krylov (NK)-based fully implicit algorithm for solving fluid flows in a wide range of flow conditions—from variable density nearly incompressible to supersonic shock dynamics. The key enabling feature of our all-speed solver is the ability to efficiently solve conservation laws by choosing a set of independent variables that produce a well-conditioned Jacobian matrix for the linear iterations of the global nonlinear iterative solver. In particular, instead of choosing to discretize the conservative variables (density, momentum, total energy), which is traditionally used in Eulerian high-speed compressible fluid dynamics, we demonstrate superior performance by discretizing the primitive variables—pressure–velocity–temperaturemore »« less
DOI: 10.1007/s00193-018-0871-8

Full Text Available
A discontinuous Galerkin conservative level set scheme for interface capturing in multiphase flows

Journal Article Owkes, Mark, E-mail: mfc86@cornell.edu ; Desjardins, Olivier - Journal of Computational Physics

The accurate conservative level set (ACLS) method of Desjardins et al. [O. Desjardins, V. Moureau, H. Pitsch, An accurate conservative level set/ghost fluid method for simulating turbulent atomization, J. Comput. Phys. 227 (18) (2008) 8395–8416] is extended by using a discontinuous Galerkin (DG) discretization. DG allows for the scheme to have an arbitrarily high order of accuracy with the smallest possible computational stencil resulting in an accurate method with good parallel scaling. This work includes a DG implementation of the level set transport equation, which moves the level set with the flow field velocity, and a DG implementation of themore »« less
DOI: 10.1016/J.JCP.2013.04.036
An assessment of atypical mesh topologies for low-Mach large-eddy simulation

Journal Article Domino, Stefan P. ; Sakievich, Philip ; Barone, Matthew F. - Computers and Fluids

An implicit, low-dissipation, low-Mach, variable density control volume finite element formulation is used to explore foundational understanding of numerical accuracy for large-eddy simulation applications on hybrid meshes. Detailed simulation comparisons are made between low-order hexahedral, tetrahedral, pyramid, and wedge/prism topologies against a third-order, unstructured hexahedral topology. Using smooth analytical and manufactured low-Mach solutions, design-order convergence is established for the hexahedral, tetrahedral, pyramid, and wedge element topologies using a new open boundary condition based on energy-stable methodologies previously deployed within a finite-difference context. A wide range of simulations demonstrate that low-order hexahedral- and wedge-based element topologies behave nearly identically in bothmore »« less
Cited by 1
DOI: 10.1016/j.compfluid.2018.12.002

Full Text Available
A Lagrangian discontinuous Galerkin hydrodynamic method

Journal Article Liu, Xiaodong ; Morgan, Nathaniel Ray ; Burton, Donald E. - Computers and Fluids

Here, we present a new Lagrangian discontinuous Galerkin (DG) hydrodynamic method for solving the two-dimensional gas dynamic equations on unstructured hybrid meshes. The physical conservation laws for the momentum and total energy are discretized using a DG method based on linear Taylor expansions. Three different approaches are investigated for calculating the density variation over the element. The first approach evolves a Taylor expansion of the specific volume field. The second approach follows certain finite element methods and uses the strong mass conservation to calculate the density field at a location inside the element or on the element surface. The thirdmore »« less
Cited by 4
DOI: 10.1016/j.compfluid.2017.12.007

Full Text Available

Similar Records