-
Arbitrary Order Energy and Enstrophy Conserving Finite Element Methods for 2D Incompressible Fluid Dynamics and Drift-Reduced Magnetohydrodynamics
Authors:
Milan Holec,
Ben Zhu,
Ilon Joseph,
Christopher J. Vogl,
Ben S. Southworth,
Alejandro Campos,
Andris M. Dimits,
Will E. Pazner
Abstract:
Maintaining conservation laws in the fully discrete setting is critical for accurate long-time behavior of numerical simulations and requires accounting for discrete conservation properties in both space and time. This paper derives arbitrary order finite element exterior calculus spatial discretizations for the two-dimensional (2D) Navier-Stokes and drift-reduced magnetohydrodynamic equations tha…
▽ More
Maintaining conservation laws in the fully discrete setting is critical for accurate long-time behavior of numerical simulations and requires accounting for discrete conservation properties in both space and time. This paper derives arbitrary order finite element exterior calculus spatial discretizations for the two-dimensional (2D) Navier-Stokes and drift-reduced magnetohydrodynamic equations that conserve both energy and enstrophy to machine precision when coupled with generally symplectic time-integration methods. Both continuous and discontinuous-Galerkin (DG) weak formulations can ensure conservation, but only generally symplectic time integration methods, such as the implicit midpoint method, permit exact conservation in time. Moreover, the symplectic implicit midpoint method yields an order of magnitude speedup over explicit schemes. The methods are implemented using the MFEM library and the solutions are verified for an extensive suite of 2D neutral fluid turbulence test problems. Numerical solutions are verified via comparison to a semi-analytic linear eigensolver as well as to the finite difference Global Drift Ballooning (GDB) code. However, it is found that turbulent simulations that conserve both energy and enstrophy tend to have too much power at high wavenumber and that this part of the spectrum should be controlled by reintroducing artificial dissipation. The DG formulation allows upwinding of the advection operator which dissipates enstrophy while still maintaining conservation of energy. Coupling upwinded DG with implicit symplectic integration appears to offer the best compromise of allowing mid-range wavenumbers to reach the appropriate amplitude while still controlling the high-wavenumber part of the spectrum.
△ Less
Submitted 25 February, 2022;
originally announced February 2022.
-
A New Scheme for Solving High-Order DG Discretizations of Thermal Radiative Transfer using the Variable Eddington Factor Method
Authors:
Ben C. Yee,
Samuel S. Olivier,
Ben S. Southworth,
Milan Holec,
Terry S. Haut
Abstract:
We present a new approach for solving high-order thermal radiative transfer (TRT) using the Variable Eddington Factor (VEF) method (also known as quasidiffusion). Our approach leverages the VEF equations, which consist of the first and second moments of the $S_N$ transport equation, to more efficiently compute the TRT solution for each time step. The scheme consists of two loops - an outer loop to…
▽ More
We present a new approach for solving high-order thermal radiative transfer (TRT) using the Variable Eddington Factor (VEF) method (also known as quasidiffusion). Our approach leverages the VEF equations, which consist of the first and second moments of the $S_N$ transport equation, to more efficiently compute the TRT solution for each time step. The scheme consists of two loops - an outer loop to converge the Eddington tensor and an inner loop to converge the iteration between the temperature equation and the VEF system. By converging the outer iteration, one obtains the fully implicit TRT solution for the given time step with a relatively low number of transport sweeps. However, one could choose to perform exactly one outer iteration (and therefore exactly one sweep) per time step, resulting in a semi-implicit scheme that is both highly efficient and robust. Our results indicate that the error between the one-sweep and fully implicit variants of our scheme may be small enough for consideration in many problems of interest.
△ Less
Submitted 15 April, 2021;
originally announced April 2021.
-
Diffusion synthetic acceleration for heterogeneous domains, compatible with voids
Authors:
Ben S. Southworth,
Milan Holec,
Terry S. Haut
Abstract:
A standard approach to solving the S$_N$ transport equations is to use source iteration with diffusion synthetic acceleration (DSA). Although this approach is widely used and effective on many problems, there remain some practical issues with DSA preconditioning, particularly on highly heterogeneous domains. For large-scale parallel simulation, it is critical that both (i) preconditioned source it…
▽ More
A standard approach to solving the S$_N$ transport equations is to use source iteration with diffusion synthetic acceleration (DSA). Although this approach is widely used and effective on many problems, there remain some practical issues with DSA preconditioning, particularly on highly heterogeneous domains. For large-scale parallel simulation, it is critical that both (i) preconditioned source iteration converges rapidly, and (ii) the action of the DSA preconditioner can be applied using fast, scalable solvers, such as algebraic multigrid (AMG). For heterogeneous domains, these two interests can be at odds. In particular, there exist DSA diffusion discretizations that can be solved rapidly using AMG, but they do not always yield robust/fast convergence of the larger source iteration. Conversely, there exist robust DSA discretizations where source iteration converges rapidly on difficult heterogeneous problems, but fast parallel solvers like AMG tend to struggle applying the action of such operators. Moreover, very few current methods for the solution of deterministic transport are compatible with voids. This paper develops a new heterogeneous DSA preconditioner based on only preconditioning the optically thick subdomains. The resulting method proves robust on a variety of heterogeneous transport problems, including a linearized hohlraum mesh related to inertial confinement fusion. Moreover, the action of the preconditioner is easily computed using $\mathcal{O}(1)$ AMG iterations, {convergence of the transport iteration typically requires $2-5\times$ less iterations than current state-of-the-art ``full DSA,'' and the proposed method is} trivially compatible with voids. On the hohlraum problem, rapid convergence is obtained by preconditioning less than 3\% of the mesh elements with $5-10$ AMG iterations.
△ Less
Submitted 20 July, 2020; v1 submitted 24 January, 2020;
originally announced January 2020.
-
A Quadratic Programming Flux Correction Method for High-Order DG Discretizations of SN Transport
Authors:
Ben C. Yee,
Samuel S. Olivier,
Terry S. Haut,
Milan Holec,
Vladimir Z. Tomov,
Peter G. Maginot
Abstract:
We present a new flux-fixup approach for arbitrarily high-order discontinuous Galerkin discretizations of the SN transport equation. This approach is sweep-compatible: as the transport sweep is performed, a local quadratic programming (QP) problem is solved in each spatial cell to ensure that the solution satisfies certain physical constraints, including local particle balance. The constraints can…
▽ More
We present a new flux-fixup approach for arbitrarily high-order discontinuous Galerkin discretizations of the SN transport equation. This approach is sweep-compatible: as the transport sweep is performed, a local quadratic programming (QP) problem is solved in each spatial cell to ensure that the solution satisfies certain physical constraints, including local particle balance. The constraints can be chosen in two ways, leading to two variants of the method: QP Zero (QPZ) and QP Maximum Principle (QPMP). The coefficients of the solution are constrained to be nonnegative in QPZ, and they are constrained by an approximate discrete maximum principle in QPMP. There are two primary takeaways in this paper. First, it is shown that the QPMP method, when used with the positive Bernstein basis, eliminates negativities, preserves high-order accuracy for smooth problems, and significantly dampens unphysical oscillations in the solution. The second takeaway is that the Variable Eddington Factor (VEF) method can be used to accelerate the convergence of source iteration with fixup for problems with optically thick regions. Regardless of whether a fixup is applied, source iteration converges slowly when optically thick regions are present and acceleration is needed. When VEF is combined with fixed-up transport sweeps, the result is a scheme that produces a nonnegative solution, converges independently of the mean free path, and, in the case of the QPMP fixup, adheres to an approximate discrete maximum principle.
△ Less
Submitted 26 June, 2020; v1 submitted 7 October, 2019;
originally announced October 2019.
-
AWBS kinetic modeling of electrons with nonlocal Ohms law in plasmas relevant to inertial confinement fusion
Authors:
M. Holec,
P. Loiseau,
A. Debayle,
J. P. Brodrick,
D. Del Sorbo,
C. P. Ridgers,
V. Tikhonchuk,
J. -L. Feugeas,
Ph. Nicolai,
B. Dubroca,
R. J. Kingham
Abstract:
The interaction of lasers with plasmas very often leads to nonlocal transport conditions, where the classical hydrodynamic model fails to describe important microscopic physics related to highly mobile particles. In this study we analyze and further propose a modification of the Albritton- Williams-Bernstein-Swartz collision operator Phys. Rev. Lett 57, 1887 (1986) for the nonlocal electron transp…
▽ More
The interaction of lasers with plasmas very often leads to nonlocal transport conditions, where the classical hydrodynamic model fails to describe important microscopic physics related to highly mobile particles. In this study we analyze and further propose a modification of the Albritton- Williams-Bernstein-Swartz collision operator Phys. Rev. Lett 57, 1887 (1986) for the nonlocal electron transport under conditions relevant to ICF. The electron distribution function provided by this modification exhibits some very desirable properties when compared to the full Fokker- Planck operator in the local diffusive regime, and also performs very well when benchmarked against Vlasov-Fokker-Planck and collisional PIC codes in the nonlocal transport regime, where we find that the effect of the electric field via the nonlocal Ohms law is an essential ingredient in order to capture the electron kinetics properly.
△ Less
Submitted 22 December, 2018;
originally announced January 2019.