Results for “Discontinuous diffusion coefficients” — Page 2

NASA NTRS1980

Gas phase hydrogen permeation in alpha titanium and carbon steels

Commercially pure titanium and heats of Armco ingot iron and steels containing from 0.008-1.23 w/oC were annealed or normalized and machined into hollow cylinders. Coefficients of diffusion for alpha-Ti and alpha-Fe were determined by the lag-time technique. Steady state permeation experiments yield first power pressure dependence for alpha-Ti and Sievert's law square root dependence for Armco iron and carbon steels. As in the case of diffusion, permeation data confirm that alpha-titanium is subject to at least partial phase boundary reaction control while the steels are purely diffusion controlled. The permeation rate in steels also decreases as the carbon content increases. As a consequence of Sievert's law, the computed hydrogen solubility decreases as the carbon content increases. This decreases in explained in terms of hydrogen trapping at carbide interfaces. Oxidizing and nitriding the surfaces of alpha-titanium membranes result in a decrease in the permeation rate for such treatment on the gas inlet surfaces but resulted in a slight increase in the rate for such treatment on the gas outlet surfaces. This is explained in terms of a discontinuous TiH2 layer.

Johnson, D. L.↗

DOE OSTI2020

A quadratic programming flux correction method for high-order DG discretizations of S transport

In this work, we present a new flux-fixup approach for arbitrarily high-order discontinuous Galerkin (DG) discretizations of the S N transport equation, and we demonstrate the compatibility of this approach with the Variable Eddington Factor (VEF) method. The new fixup approach is sweep-compatible: during a transport sweep (block Gauss-Seidel iteration in which the scattering source is lagged), a local quadratic programming (QP) problem is solved in each spatial element to ensure that the solution satisfies certain physical constraints, including local particle balance. In this paper, we describe two choices of physical constraints, resulting in two variants of the method: QP Zero (QPZ) and QP Maximum Principle (QPMP). In QPZ, the finite element coefficients of the solution are constrained to be nonnegative. In QPMP, they are constrained to adhere to an approximate discrete maximum principle. There are two primary takeaways in this paper. First, when the positive Bernstein basis is used for DG discretization, the QPMP method eliminates negativities, preserves high-order accuracy for smooth problems, and significantly dampens unphysical oscillations in the solution. The latter feature – the dampening of unphysical oscillations – is an improvement upon standard, simpler fixup approaches such as the approach described in (denoted as the “zero and rescale” (ZR) method in this paper). This improvement comes at a moderate computational cost, but it is not prohibitive. Our results show that, even in an unrealistic worst-case scenario where 83% of the spatial elements require a fixup, the computational cost of performing a transport sweep with fixup is only ~31% greater than performing one without fixup. The second takeaway is that the VEF method can be used to accelerate the convergence of transport sweeps even when a fixup is applied. When optically thick regions are present, transport sweeps converge slowly, regardless of whether a fixup is applied, and acceleration is needed. However, attempting to apply standard diffusion synthetic acceleration (DSA) to fixed-up transport sweeps results in divergence for optically thick problems. Our results show that the same is not true for VEF. 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.

97 MATHEMATICS AND COMPUTING↗