SEARCH · Engineering Papers

Results for “Formal Methods”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,225 records · Page 68

All sources NASA NTRS DOE OSTI

DOE OSTI2017

The Thermal Proximity Effect: A New Probe of the He ii Reionization History and Quasar Lifetime

Despite decades of effort, the timing and duration of He ii reionization and the properties of the quasars believed to drive it are still not well constrained. We present a new method to study both via the thermal proximity effect—the heating of the intergalactic medium (IGM) around quasars when their radiation doubly ionizes helium. We post-process hydrodynamical simulations with 1D radiative transfer and study how the thermal proximity effect depends on the He ii fraction, x{sub HeII,0}, which prevailed in the IGM before the quasar turned on, and the quasar lifetime t{sub Q}. We find that the amplitude of the temperature boost in the quasar environment depends on x{sub HeII,0}, with a characteristic value of ΔT≃10{sup 4} K for x{sub HeII,0}=1.0, whereas the size of the thermal proximity zone is sensitive to t{sub Q}, with typical sizes of ≃100 cMpc for t{sub Q}=10{sup 8} yr. This temperature boost increases the thermal broadening of H i absorption lines near the quasar. We introduce a new Bayesian statistical method based on measuring the Lyα forest power spectrum as a function of distance from the quasar, and demonstrate that the thermal proximity effect should be easily detectable. For a mock data set of 50 quasars at z≃4, we predict that one can measure x{sub HeII,0} to an (absolute) precision ≈0.04 and t{sub Q} to a precision of ≈0.1 dex. By applying our formalism to existing high-resolution Lyα forest spectra, one should be able to reconstruct the He ii reionization history, providing a global census of hard photons in the high-z universe.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2020

Enforcing constraints for time series prediction in supervised, unsupervised and reinforcement learning

We assume that we are given a time series of data from a dynamical system and our task is to learn the flow map of the dynamical system. We present a collection of results on how to enforce constraints coming from the dynamical system in order to accelerate the training of deep neural networks to represent the flow map of the system as well as increase their predictive ability. In particular, we provide ways to enforce constraints during training for all three major modes of learning, namely supervised, unsupervised and reinforcement learning. In general, the dynamic constraints need to include terms which are analogous to memory terms in model reduction formalisms. Such memory terms act as a restoring force which corrects the errors committed by the learned flow map during prediction. For supervised learning, the constraints are added to the objective function. For the case of unsupervised learning, in particular generative adversarial networks, the constraints are introduced by augmenting the input of the discriminator. Finally, for the case of reinforcement learning and in particular actor-critic methods, the constraints are added to the reward function. In addition, for the reinforcement learning case, we present a novel approach based on homotopy of the action-value function in order to stabilize and accelerate training. We use numerical results for the Lorenz system to illustrate the various constructions.

Stinis, Panagiotis↗

DOE OSTI2023

A Stabilizer Framework for the Contextual Subspace Variational Quantum Eigensolver and the Noncontextual Projection Ansatz

Quantum chemistry is a promising application for noisy intermediate-scale quantum (NISQ) devices. However, quantum computers have thus far not succeeded in providing solutions to problems of real scientific significance, with algorithmic advances being necessary to fully utilize even the modest NISQ machines available today. We discuss a method of ground state energy estimation predicated on a partitioning of the molecular Hamiltonian into two parts: one that is noncontextual and can be solved classically, supplemented by a contextual component that yields quantum corrections obtained via a Variational Quantum Eigensolver (VQE) routine. This approach has been termed Contextual Subspace VQE (CS-VQE); however, there are obstacles to overcome before it can be deployed on NISQ devices. The problem we address here is that of the ansatz, a parametrized quantum state over which we optimize during VQE; it is not initially clear how a splitting of the Hamiltonian should be reflected in the CS-VQE ansätze. We propose a “noncontextual projection” approach that is illuminated by a reformulation of CS-VQE in the stabilizer formalism. This defines an ansatz restriction from the full electronic structure problem to the contextual subspace and facilitates an implementation of CS-VQE that may be deployed on NISQ devices. We validate the noncontextual projection ansatz using a quantum simulator and demonstrate chemically precise ground state energy calculations for a suite of small molecules at a significant reduction in the required qubit count and circuit depth.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2019

EM Program Expectations on Implementation of DOE O 458: Radiation Protection of Public/Environment - 19553

DOE Order 458.1, Radiation Protection of the Public and Environment, was issued in 2011. DOE sites were given 18 months to implement. The revised Order promulgates requirements for property clearance actions and independent verification. DOE Order 458.1 allows for preapproved authorized limits for the clearance of property. However, when property with the potential to contain residual radioactive material does not meet preapproved authorized limits, an authorized limits package must be developed. The Order contains a range of requirements for demonstrating with reasonable assurance that dose constraints are met for any actual or likely future use of the property to be released. EM has developed a draft guidance document that provides the elements of an authorized limits package. These elements are used to review all authorized limits packages that are submitted to EM. Expected elements include property description and process knowledge documentation, current and likely future use, the conceptual site model, dose modeling methodologies, specific authorized limits, dose- or concentration-based constraints, collective dose estimate to the potentially exposed population, stakeholder communications, description of the As Low As Reasonably Achievable (ALARA) analysis, implementation plan, evaluation of difficult to assess surfaces, and independent verification plan. The draft guidance document also provides further explanation of the desired elements of an independent verification (IV) plan. An IV plan should describe the planning phase, implementation phase, assessment phase, and decision phase of an activity. The planning phase includes scope and cost, concurrent with cleanup contractor procurement, and review and approval of site radiological survey clearance program and IV program plans. The implementation phase assures sufficient collection of data and adjustment can be made as necessary. During the assessment phase, data and reports are reviewed, and coordination between the clean-up contractor and IV personnel occurs. A decision of whether or not the clearance actions need to be altered occurs during the decision phase. IV certification that the clearance action is complete also occurs during the decision phase. Issuance of the final guidance document is meant to improve the consultation process and minimize iterative interactions with EM. Standardizing and formalizing property clearance and IV actions also helps assure that cleanup actions are transparent and sustainable for future scrutiny. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2015

Oscillator strengths, first-order properties, and nuclear gradients for local ADC(2)

We describe theory and implementation of oscillator strengths, orbital-relaxed first-order properties, and nuclear gradients for the local algebraic diagrammatic construction scheme through second order. The formalism is derived via time-dependent linear response theory based on a second-order unitary coupled cluster model. The implementation presented here is a modification of our previously developed algorithms for Laplace transform based local time-dependent coupled cluster linear response (CC2LR); the local approximations thus are state specific and adaptive. The symmetry of the Jacobian leads to considerable simplifications relative to the local CC2LR method; as a result, a gradient evaluation is about four times less expensive. Test calculations show that in geometry optimizations, usually very similar geometries are obtained as with the local CC2LR method (provided that a second-order method is applicable). As an exemplary application, we performed geometry optimizations on the low-lying singlet states of chlorophyllide a.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2013

DISCOVERY OF AN H{alpha} EMITTING DISK AROUND THE SUPERMASSIVE BLACK HOLE OF M31

Due to its proximity, the mass of the supermassive black hole in the nucleus of the Andromeda galaxy (M31), the most massive black hole in the Local Group of galaxies, has been measured by several methods involving the kinematics of a stellar disk which surrounds it. We report here the discovery of an eccentric H{alpha} emitting disk around the black hole at the center of M31 and show how modeling this disk can provide an independent determination of the mass of the black hole. Our model implies a mass of 5.0{sup +0.8}{sub -1.0} Multiplication-Sign 10{sup 7} M{sub Sun} for the central black hole, consistent with the average of determinations by methods involving stellar dynamics, and compatible (at 1{sigma} level) with measurements obtained from the most detailed models of the stellar disk around the central black hole. This value is also consistent with the M-{sigma} relation. In order to make a comparison, we applied our simulation on the stellar kinematics in the nucleus of M31 and concluded that the parameters obtained for the stellar disk are not formally compatible with the parameters obtained for the H{alpha} emitting disk. This result suggests that the stellar and the H{alpha} emitting disks are intrinsically different from each other. A plausible explanation is that the H{alpha} emission is associated with a gaseous disk. This hypothesis is supported by the detection of traces of weaker nebular lines in the nuclear region of M31. However, we cannot exclude the possibility that the H{alpha} emission is, at least partially, generated by stars.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2018

N-BAR PROBLEMS AS APPROXIMATIONS TO THE BREE PROBLEM

The Bree solution to the problem of a ratcheting cylinder under constant pressure and cyclic thermal load is a fundamental result in nuclear engineering and widely used as the technical basis for the ASME Boiler and Pressure Vessel Code and other design methods. However, because the loading conditions in the Bree problem are difficult to achieve experimentally there have been relatively few works experimentally examining the problem and extending it to other relevant design situations, for example cladded components. In contrast, 2-bar problems are widely studied experimentally and are relatively easy to setup. These 2-bar problems are thought to be representative of Bree-type geometries, but a formal connection has not been demonstrated. This work formally establishes the connection between the Bree cylinder and an n-bar problem – a coupled bar experiment with, in general, more than two bars linked in parallel. The connection suggests that n-bar experiments using a fairly limited number of bars might be an experimentally-accessible setup that better represents ratcheting phenomenon in actual nuclear pressurized components. Such experiments could test surrogate cladded or multi-material components by using bars of different materials. Finally, this work suggests control schemes that yield optimally efficient n-bar experiments – experiments that best replicates a Bree cylinder with a limited number of bars.

Messner, Mark C.↗

DOE OSTI2019

{ital n}-bar problems as approximations to the Bree problem

The Bree solution to the problem of a ratcheting cylinder under constant pressure and cyclic thermal load is a fundamental result in nuclear engineering and widely used as the technical basis for the ASME Boiler and Pressure Vessel Code and other design methods. However, because the loading conditions in the Bree problem are difficult to achieve experimentally there have been relatively few works experimentally examining the problem and extending it to other relevant design situations, for example cladded components. In contrast, 2-bar problems are widely studied experimentally and are relatively easy to setup. These 2-bar problems are thought to be representative of Bree-type geometries, but a formal connection has not been demonstrated. This work formally establishes the connection between the Bree cylinder and an n-bar problem a coupled bar experiment with, in general, more than two bars linked in parallel. The connection suggests that n-bar experiments using a fairly limited number of bars might be an experimentally-accessible setup that better represents ratcheting phenomenon in actual nuclear pressurized components. Such experiments could test surrogate cladded or multi-material components by using bars of different materials. Finally, this work suggests control schemes that yield optimally efficient n-bar experiments experiments that best replicates a Bree cylinder with a limited number of bars.

Messner, M. C.↗

DOE OSTI2018

Weak and strong orders of linear recurring sequences

The concept of the strong order of linear recurring sequence (LRS) is introduced in this paper. Necessary and sufficient conditions for the existence of the strong LRS order are derived. The strong LRS order is exploited for the formalization of the problem of the extension of a sequence from the available fragment (fragments) of that sequence. The definition of the strong LRS order opens new possibilities for formal sequence analysis whenever the weak LRS order of that sequence exists. Computational experiments with discrete iterative maps are used to illustrate the applicability of the strong LRS order in nonlinear system analysis.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2018

Investigation of deformation effects on the decay properties of {sup 12}C+α Cluster states in {sup 16}O

We have analyzed the elastic scattering angular distributions data of the α+{sup 12}C reaction over a wide energy range (E{sub lab}=28.2 to 35.5 MeV) within the framework of the Optical Model formalism. A double folding (DF) type real potential was used with a phenomenological Woods–Saxon-squared (WS2) type imaginary potential. Good agreement between the calculations and experimental data was obtained. By using the real DF potential we have calculated the properties of the α-cluster states in {sup 16}O by using the Gamow code as well as the α-decay widths by using the WKB method. We implemented a {sup 12}C+α cluster framework for the calculation of the excitation energies and decay widths of {sup 16}O as a function of the orientation of the planar {sup 12}C nucleus with respect to the α-particle. These calculations showed strong sensitivity of the widths and excitation energies to the orientation. Branching ratios were also calculated and though less sensitive to the {sup 12}C orientation, it was found that {sup 12}C{sub gs}+α structure, with the α-particle orbiting the {sup 12}C in its ground state, is dominant. This work demonstrates that deformation, and the orientation, of {sup 12}C plays a crucial role in the understanding of the nature of the α-cluster states in {sup 16}O.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2021

Modeling fission product diffusion in TRISO fuel particles with BISON

Diffusion of fission products in intact TRISO particles depends on particle geometry, fission product source rates, time, temperature, and temperature-dependent diffusion coefficients. Simulating this diffusion process requires models for source rates and diffusion coefficients, plus computation of the temperature field if not prescribed. In addition, simulation quality depends on discretization of the geometry, appropriate time stepping, and the accuracy of the solution method. In this paper, we explore the simulation of fission product diffusion in TRISO fuel particles using the finite element method via the fuel performance code Bison. Recent material model development has occurred in Bison for each material present in tri-structural isotropic (TRISO) fuel particles: the buffer, inner pyrolytic carbon, silicon carbide, and outer pyrolytic carbon layers, as well as the fuel kernel. Also, new mesh generation and fission product release fraction capabilities have been added. Diffusion capabilities are shown to converge to the correct solution via formal verification tests. A large number of code benchmarking problems are also given, with good results, showing that Bison’s computed release fractions closely match those of other software tools. Finally, a significant validation effort is detailed in which fission product release, measured as part of the AGR-1 capsule experiments, is compared to Bison outputs. Bison outputs compare very well to the experimental data and to PARFUME results.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2014

Dynamics in the quantum/classical limit based on selective use of the quantum potential

A classical limit of quantum dynamics can be defined by compensation of the quantum potential in the time-dependent Schrödinger equation. The quantum potential is a non-local quantity, defined in the trajectory-based form of the Schrödinger equation, due to Madelung, de Broglie, and Bohm, which formally generates the quantum-mechanical features in dynamics. Selective inclusion of the quantum potential for the degrees of freedom deemed “quantum,” defines a hybrid quantum/classical dynamics, appropriate for molecular systems comprised of light and heavy nuclei. The wavefunction is associated with all of the nuclei, and the Ehrenfest, or mean-field, averaging of the force acting on the classical degrees of freedom, typical of the mixed quantum/classical methods, is avoided. The hybrid approach is used to examine evolution of light/heavy systems in the harmonic and double-well potentials, using conventional grid-based and approximate quantum-trajectory time propagation. The approximate quantum force is defined on spatial domains, which removes unphysical coupling of the wavefunction fragments corresponding to distinct classical channels or configurations. The quantum potential, associated with the quantum particle, generates forces acting on both quantum and classical particles to describe the backreaction.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2013

Three independent one-dimensional margins for single-fraction frameless stereotactic radiosurgery brain cases using CBCT

Purpose: Setting a proper margin is crucial for not only delivering the required radiation dose to a target volume, but also reducing the unnecessary radiation to the adjacent organs at risk. This study investigated the independent one-dimensional symmetric and asymmetric margins between the clinical target volume (CTV) and the planning target volume (PTV) for linac-based single-fraction frameless stereotactic radiosurgery (SRS).Methods: The authors assumed a Dirac delta function for the systematic error of a specific machine and a Gaussian function for the residual setup errors. Margin formulas were then derived in details to arrive at a suitable CTV-to-PTV margin for single-fraction frameless SRS. Such a margin ensured that the CTV would receive the prescribed dose in 95% of the patients. To validate our margin formalism, the authors retrospectively analyzed nine patients who were previously treated with noncoplanar conformal beams. Cone-beam computed tomography (CBCT) was used in the patient setup. The isocenter shifts between the CBCT and linac were measured for a Varian Trilogy linear accelerator for three months. For each plan, the authors shifted the isocenter of the plan in each direction by ±3 mm simultaneously to simulate the worst setup scenario. Subsequently, the asymptotic behavior of the CTV V{sub 80%} for each patient was studied as the setup error approached the CTV-PTV margin.Results: The authors found that the proper margin for single-fraction frameless SRS cases with brain cancer was about 3 mm for the machine investigated in this study. The isocenter shifts between the CBCT and the linac remained almost constant over a period of three months for this specific machine. This confirmed our assumption that the machine systematic error distribution could be approximated as a delta function. This definition is especially relevant to a single-fraction treatment. The prescribed dose coverage for all the patients investigated was 96.1%± 5.5% with an extreme 3-mm setup error in all three directions simultaneously. It was found that the effect of the setup error on dose coverage was tumor location dependent. It mostly affected the tumors located in the posterior part of the brain, resulting in a minimum coverage of approximately 72%. This was entirely due to the unique geometry of the posterior head.Conclusions: Margin expansion formulas were derived for single-fraction frameless SRS such that the CTV would receive the prescribed dose in 95% of the patients treated for brain cancer. The margins defined in this study are machine-specific and account for nonzero mean systematic error. The margin for single-fraction SRS for a group of machines was also derived in this paper.

61 RADIATION PROTECTION AND DOSIMETRY↗

DOE OSTI2015

SU-E-J-89: Comparative Analysis of MIM and Velocity’s Image Deformation Algorithm Using Simulated KV-CBCT Images for Quality Assurance

Purpose: Deformable image registration (DIR) is used routinely in the clinic without a formalized quality assurance (QA) process. Using simulated deformations to digitally deform images in a known way and comparing to DIR algorithm predictions is a powerful technique for DIR QA. This technique must also simulate realistic image noise and artifacts, especially between modalities. This study developed an algorithm to create simulated daily kV cone-beam computed-tomography (CBCT) images from CT images for DIR QA between these modalities. Methods: A Catphan and physical head-and-neck phantom, with known deformations, were used. CT and kV-CBCT images of the Catphan were utilized to characterize the changes in Hounsfield units, noise, and image cupping that occur between these imaging modalities. The algorithm then imprinted these changes onto a CT image of the deformed head-and-neck phantom, thereby creating a simulated-CBCT image. CT and kV-CBCT images of the undeformed and deformed head-and-neck phantom were also acquired. The Velocity and MIM DIR algorithms were applied between the undeformed CT image and each of the deformed CT, CBCT, and simulated-CBCT images to obtain predicted deformations. The error between the known and predicted deformations was used as a metric to evaluate the quality of the simulated-CBCT image. Ideally, the simulated-CBCT image registration would produce the same accuracy as the deformed CBCT image registration. Results: For Velocity, the mean error was 1.4 mm for the CT-CT registration, 1.7 mm for the CT-CBCT registration, and 1.4 mm for the CT-simulated-CBCT registration. These same numbers were 1.5, 4.5, and 5.9 mm, respectively, for MIM. Conclusion: All cases produced similar accuracy for Velocity. MIM produced similar values of accuracy for CT-CT registration, but was not as accurate for CT-CBCT registrations. The MIM simulated-CBCT registration followed this same trend, but overestimated MIM DIR errors relative to the CT-CBCT registration.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2012

Precision Measurements of Lifetime and $\bf CP$ Violation Phase in the $\bf B^{0}_{s}$ Meson with the DØ Detector

The Standard Model of Partice Physics (SM) is probably the most successful theory, regarding to his predictions. The SM prediction for $CP$ violation is not enough to explain the overwhelming asymmetry among the matter and anti-matter abundance. Measuring some process where $CP$ violation is different to the one predicted by the SM would be a clear signal for Physics Beyond the Standard Model. The SM prediction for the $CP$ violation phase, $\phi_{s}$, in the $B^{0}_{s}$ meson is practically equal to zero for the current experiments. This means that measuring a deviation from zero in $\phi_{s}$ could be an indication for Physics Beyond the SM. On the other hand, the approximation based on the ``heavy quark symmetry'' let approximated calculations of the fundamental quantities of those hadrons containing a heavy quark, $c,b,t$. These calculations are expressed as expansions on inverse powers of the heavy quark mass in such hadron. This formalism is called `` Heavy Quark Effective Theory'' (HQET), and has been successful predicting some properties of the heavy hadrons. The HQET prediction for the lifetime ratio the $B^{0}_{d}$ and $B^{0}_{s}$ is practically equal to one. So, measuring with good precision the $B^{0}_{s}$ lifetime is also a way to test an approximation based on the SM. In this thesis it is detailed presented the method to measure the $\phi_{s}$ and the lifetime ratio of the $B^{0}_{d}$ and $B^{0}_{s}$, among other quantities, with the DØ located in the Fermi National Accelerator Laboratory, in the United States.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2016

SU-F-T-127: Charged Particle Transport in Condensed Media

Purpose: Provide quality interaction cross sections for charged particle Monte Carlo (MC) track structure codes and evaluate low energy electron transport in condensed media. Methods: MC methods and codes are often used to model or simulate charged particle radiation transport in matter. Detailed (or event-by-event) track structure simulations are of special interest for the modeling of the physical and chemical stages of radiation action with matter and the initial radiation damage to biological systems. They require reliable interaction cross sections of all radiation qualities considered (e.g., electrons, protons, alpha particles, light and heavy ions) with the target material under consideration, mainly in the condensed phase, including liquids and solids. Interaction cross sections are calculated using the dielectric formalism, a mixture of first principles, theoretical modeling and experimental information for bulk and surface transport and implemented into MC track structure codes. Secondary electron emission yields from amorphous solid water (ASW) and thin metal foils (copper and gold) after fast proton impact are simulated and measured experimentally. Results: After considering different transport models for bulk and surface transport and a careful modeling of the experimental geometry, simulated secondary electron emission yields follow the trends of experimental data well. Furthermore, yields for electron emissions below 50 eV are sensitive to differential elastic scattering cross sections. Conclusion: Low-energy electron transport in condensed media is still a challenge for detailed track structure simulation codes. Interaction cross-sections, transport models, and target geometry need to be considered adequately. The research was funded in part by NSF Major Research Instrumentation Program 2009, Award Number: 0923270 and by NIH/NCI R01 CA93351.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2014

WE-E-BRE-11: New Method to Simulate DNA Damage Using Ionization Cross-Sections and a Geometrical Nucleosome Model

Purpose: To obtain probability distributions of various DNA damage types as a function of the incident electron kinetic energy. Methods: Using Geant4-DNA electron ionization cross-sections, we calculated path length distributions for electrons of energies between 10 eV and 1 MeV, defined as the length between two subsequent interactions. These path lengths were then convolved with probability distributions for the creation of same-strand damage, opposite-strand damage, clustered damage, isolated damage, and same DNA strand target damage. These probability distributions of DNA damage were obtained by a Monte Carlo routine calculating probabilities of interaction in DNA targets inside a nucleosome geometrical model. Results represent the probability of a secondary electron, initially created inside a DNA strand target, of undergoing its next interaction: (1) in the opposite strand (DSB), (2) in the same strand (SSB+), (3) in either the opposite or same-strand (clustered), (4) in the same DNA target (multiple-hit) or (5) more than 10 base pairs away (isolated). Results: Electrons with kinetic energy between 50 and 250 eV have a maximal probability of creating DSB, SSB+, clustered damage and multiple-hits in the same target The probabilities for these damage patterns have values of 2.5%, 4.3%, 6.7% and 5.4%, respectively. Isolated damage is most probable between 700 eV to 900 eV with a probability of 0.2%. Conclusion: We obtained DNA damage probability distributions as a function of electron incident energy. We showed that electrons with kinetic energies between 50 and 250 eV have the highest probability of producing complex forms of DNA damage (DSB, SSB+). We also showed that a double ionization within the same DNA target is the most frequent outcome occurring 5% of the time. It is expected that electron slowing down spectra can be convolved with our formalism to calculate source specific DNA damage patterns. Research grants from governments of Canada and Quebec. PP acknowledges partial support by the CREATE Medical Physics Research Training Network grant of the Natural Sciences and Engineering Research Council (Grant number: 432290)

74 ATOMIC AND MOLECULAR PHYSICS↗

DOE OSTI2016

Asymptotic analysis of discrete schemes for non-equilibrium radiation diffusion

Motivated by providing well-behaved fully discrete schemes in practice, this paper extends the asymptotic analysis on time integration methods for non-equilibrium radiation diffusion in [2] to space discretizations. Therein studies were carried out on a two-temperature model with Larsen's flux-limited diffusion operator, both the implicitly balanced (IB) and linearly implicit (LI) methods were shown asymptotic-preserving. In this paper, we focus on asymptotic analysis for space discrete schemes in dimensions one and two. First, in construction of the schemes, in contrast to traditional first-order approximations, asymmetric second-order accurate spatial approximations are devised for flux-limiters on boundary, and discrete schemes with second-order accuracy on global spatial domain are acquired consequently. Then by employing formal asymptotic analysis, the first-order asymptotic-preserving property for these schemes and furthermore for the fully discrete schemes is shown. Finally, with the help of manufactured solutions, numerical tests are performed, which demonstrate quantitatively the fully discrete schemes with IB time evolution indeed have the accuracy and asymptotic convergence as theory predicts, hence are well qualified for both non-equilibrium and equilibrium radiation diffusion. - Highlights: • Provide AP fully discrete schemes for non-equilibrium radiation diffusion. • Propose second order accurate schemes by asymmetric approach for boundary flux-limiter. • Show first order AP property of spatially and fully discrete schemes with IB evolution. • Devise subtle artificial solutions; verify accuracy and AP property quantitatively. • Ideas can be generalized to 3-dimensional problems and higher order implicit schemes.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗