Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Mathematics and Computing, Physics - Nuclear physics and radiation physics”

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.

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

At least 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Analytical Modeling of Exoplanet Transit Spectroscopy with Dimensional Analysis and Symbolic Regression

Abstract The physical characteristics and atmospheric chemical composition of newly discovered exoplanets are often inferred from their transit spectra, which are obtained from complex numerical models of radiative transfer. Alternatively, simple analytical expressions provide insightful physical intuition into the relevant atmospheric processes. The deep-learning revolution has opened the door for deriving such analytical results directly with a computer algorithm fitting to the data. As a proof of concept, we successfully demonstrate the use of symbolic regression on synthetic data for the transit radii of generic hot-Jupiter exoplanets to derive a corresponding analytical formula. As a preprocessing step, we use dimensional analysis to identify the relevant dimensionless combinations of variables and reduce the number of independent inputs, which improves the performance of the symbolic regression. The dimensional analysis also allowed us to mathematically derive and properly parameterize the most general family of degeneracies among the input atmospheric parameters that affect the characterization of an exoplanet atmosphere through transit spectroscopy.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Positrons for Our Future

JLab Users want beams of polarized (and unpolarized) positrons, with beam qualities and modes of operation similar to those provided in the electron beams at CEBAF today [Gra19]. The creation of positrons with a high degree of spin polarization is however very challenging – limited to beta decay, self-polarization in rings, or via e+/e- pair-production, a process requiring a source of polarized gamma rays. At JLab we have a unique source of polarized gamma rays – the bremsstrahlung radiation produced when the spin polarized electron beams of CEBAF interacts with matter – called Polarized Electrons for Polarized Positrons (or PEPPo) [Abb04]. This project sought to extend the PEPPo concept by exploring a possible layout and parameter space for a polarized positron injector at CEBAF based on this approach.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2021

Building an AI-enhanced modeling framework to address multiscale predictability challenges

Focal Area(s): Build an AI-enhanced modeling framework that integrates the three focus areas in the solicitation. Science Challenge: The 4M-2N complexities (Multiscale, Multiphysics, Multibody, Multidimension, Non-linearity, and Non-Gaussianality) of atmospheric aerosol-cloud-precipitation-turbulence-radiation system poses physical and computational challenges to further advance predictive models; We plan to address the challenges by developing an AI-enhanced modeling framework that facilitates automated calibration and improvement of subgrid parameterizations, enhance data assimilation of measurements to improve initial and boundary conditions used to drive the physical model, and optimally blends data-driven and physics-based forecasting models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Study of Radiative Heat Transfer and Flow Physics from Medium-scale Methanol Pool Fire Simulations [Slides]

Analysis of methanol pool fire conducted as part of validation study for SIERRA/Fuego. Radiation model was effectively calibrated by modifying radiation model parameters for methanol. Computing integrated buoyancy flux, entrainment rate, and turbulent kinetic energy allowed for evaluation of less typical quantities in this validation study. Quantities were compared with experimental data or correlations and generally showed agreement. TKE needed fine mesh to be computed accurately. Predicted flame height less sensitive to variations in mixture fraction than temperature. Mixture fraction is a preferable threshold variable for this application.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Physics-guided machine learning for improved accuracy of the National Solar Radiation Database

The National Solar Radiation Database (NSRDB) provides high-resolution spatiotemporal solar irradiance data for the entire globe. The NSRDB uses a two-step Physical Solar Model (PSM) to compute the effects of clouds and other atmospheric variables on the solar radiation reaching the surface of the Earth. Physical and optical cloud properties are fundamental inputs to the PSM and are derived from the National Oceanic and Atmospheric Administration's Geostationary Operational Environmental Satellites. This paper describes recent improvements to the NSRDB driven by physics-guided machine learning methods for cloud property retrieval. The impacts of these new methods on the NSRDB irradiance data are validated using an extensive set of ground measurement sites, showing significant improvement for all sites. We report on average, the mean absolute percentage error for global horizontal irradiance and direct normal irradiance show reductions of 2.16 and 3.95 percentage points respectively for all daylight conditions, 5.92 and 17.39 percentage points respectively for cloudy conditions, and 9.00 and 22.59 percentage points respectively for gap-filled cloudy conditions. These new methods will help improve the quality and accuracy of the irradiance and cloud data in the NSRDB.

14 SOLAR ENERGY↗
DOE OSTI2021

mystic : software for autonomous discovery and design under uncertainty

Throughout the diverse range of science and engineering applications, there is a growing desire to develop computational methods that can reliably predict the behavior of complex systems. Specifically, there is a strategic need for tools that can robustly forecast the behavior of complex physical systems, where data may be high-dimensional, noisy, or sparse, and models of the system may be time-dependent or include uncertainty. We use mystic to build tools that leverage statistical learning, physics-informed learning, and active learning in the efficient generation of reliably predictive surrogates for complex physical systems. mystic is a robust, proven, open-source optimization and uncertainty quantification toolkit with over a decade of use in the design and optimization of neutron instrumentation, solar-powered drones, and gasguns, and in iterative tuning of models for Raman spectroscopy and elastoplastic materials strength. Recent developments have focused on automated learning of statistically robust surrogates under uncertainty, with applications in materials in extreme environments, nanostructures, materials simulations and strength models, and the failure of shielding under particle radiation. In 2020, McKerns demonstrated active learning of optimally robust surrogates with respect to new simulated data for molecular dynamics simulations of materials mixing in warm dense matter, and is currently applying active learning to the automated steering of particle accelerator beams and the optimal design and control of quantum optical sensor instrumentation.

42 ENGINEERING↗
DOE OSTI2022

Transport error estimation using residual Monte Carlo

The residual Monte Carlo (RMC) method is also known in the literature as sequential Monte Carlo and reduced-source Monte Carlo. Given a Monte Carlo method for solving a linear equation and an approximate solution to that system, the residual method enables use of essentially the same Monte Carlo algorithm to directly compute the additive error or “defect” associated with the approximate solution. As the size of the defect decreases relative to the size of the solution, the residual Monte Carlo method becomes increasingly efficient relative to the standard Monte Carlo (SMC) method. Here we present a new RMC algorithm for evaluating the space-angle error in S n radiation transport solutions, and provide computational examples demonstrating that it can be far more efficient than SMC for this purpose. Herein we also describe a particular pitfall that must be avoided if RMC is to be efficient, and explain why the performance of RMC can significantly differ between different transport problems and different quantities of interest for the same problem.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Reference solutions for linear radiation transport: the Hohlraum and Lattice Benchmarks

Radiation transport describes the propagation of energetic particles through space as they interact with a surrounding material medium. In a kinetic description, radiation transport is modeled by a radiation transport equation (RTE) that prescribes the density of the radiation in position-momentum phase space. The purpose of this dataset is to provide highly resolved solutions to two benchmark problems. These two benchmarks do not possess exact solutions; moreover, the construction of a manufactured solution may require a non-physical source that is not desirable, especially if it spoils the physical nature of the solution. Thus the goal of this computational study is to provide a highly resolved reference solution for testing newer, more cost efficient methods that are currently being developed in the research community.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

PAL Experimental Report Form

We constructed the first simultaneous dark- and bright-field X-ray microscope in this experiment, coupling it to simultaneous measurements of the wide-angle Xray diffraction, and the XFEL pulse intensity and energy spectrum on each shot. This experiment established proof of concept for this new advanced measurement technique and demonstrated DFXM at an XFEL for the first time. Within this complex experiment, we also acquired data that advanced our research goals within materials science and mathematics. These included: A preliminary study of how local effects reduce the X-ray damage threshold in diamond; A preliminary study of how X-ray radiation imparts heat into materials based on irreversible phase transitions that occur in boric acid with increasing temperature (based on hydration); A study of X-ray radiation damage from ultrafast thermal cycling and nonthermal melting in bismuth selenide; We also began developing software to automate alignment of compound refractive lenses (CRLs) using stochastic Nelder-Mead optimization methods. This work will expedite and improve the accuracy of alignments for many future experiments using CRLs.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Raytracing

The raytracing (spectrally specific radiative transfer that does not couple back to the hydrodynamic evolution) code was developed because it helps improve our understanding of the physics of explosive events, validates computational physics codes, and estimates optical signals that are cheaper and safer than experiments involving high explosives. The raytracing code can also produce simulations that are useful to analyze what we can understand with sensors, such as determining what we can discriminate with sensors, like the presence of certain metals or the amount of soot. In addition, simulations can provide extra validations to computational physics to help us better understand discrepancies.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
DOE OSTI2024

Classifying photonic topology using the spectral localizer and numerical K -theory

Recently, the spectral localizer framework has emerged as an efficient approach for classifying topology in photonic systems featuring local nonlinearities and radiative environments. In nonlinear systems, this framework provides rigorous definitions for concepts such as topological solitons and topological dynamics, where a system’s occupation induces a local change in its topology due to nonlinearity. For systems embedded in radiative environments that do not possess a shared bulk spectral gap, this framework enables the identification of local topology and shows that local topological protection is preserved despite the lack of a common gap. However, as the spectral localizer framework is rooted in the mathematics of C*-algebras, and not vector bundles, understanding and using this framework requires developing intuition for a somewhat different set of underlying concepts than those that appear in traditional approaches for classifying material topology. In this tutorial, we introduce the spectral localizer framework from a ground-up perspective and provide physically motivated arguments for understanding its local topological markers and associated local measure of topological protection. In doing so, we provide numerous examples of the framework’s application to a variety of topological classes, including crystalline and higher-order topology. We then show how Maxwell’s equations can be reformulated to be compatible with the spectral localizer framework, including the possibility of radiative boundary conditions. To aid in this introduction, we also provide a physics-oriented introduction to multi-operator pseudospectral methods and numerical K-theory, two mathematical concepts that form the foundation for the spectral localizer framework. Finally, we provide some mathematically oriented comments on the C*-algebraic origins of this framework, including a discussion of real C*-algebras and graded C*-algebras that are necessary for incorporating physical symmetries. Looking forward, we hope that this tutorial will serve as an approachable starting point for learning the foundations of the spectral localizer framework.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Excited-State and Nonadiabatic Molecular Dynamics Methods with Broad Applicability

Electronically excited states play a pivotal role in the capture, transformation, storage, and emission of radiative energy. Despite much recent progress for small model systems, excited state properties and dynamics in many larger molecular devices and materials have remained elusive for experiment and theory. Nonadiabatic transitions between electronic states turn radiative energy into heat and are particularly difficult to control. This project aims to develop new electronic structure and nonadiabatic molecular dynamics (NAMD) methods to simulate, predict, and ultimately control the flow and transformation of excitation energy in a wide-range of light-driven molecular devices such as dyes, photocatalysts, and fluorescence emitters.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Initial-State Dependent Optimization of Controlled Gate Operations with Quantum Computer

There is no unique way to encode a quantum algorithm into a quantum circuit. With limited qubit counts, connectivity, and coherence times, a quantum circuit optimization is essential to make the best use of near-term quantum devices. We introduce a new circuit optimizer called AQCEL, which aims to remove redundant controlled operations from controlled gates, depending on initial states of the circuit. Especially, the AQCEL can remove unnecessary qubit controls from multi-controlled gates in polynomial computational resources, even when all the relevant qubits are entangled, by identifying zero-amplitude computational basis states using a quantum computer. As a benchmark, the AQCEL is deployed on a quantum algorithm designed to model final state radiation in high energy physics. For this benchmark, we have demonstrated that the AQCEL-optimized circuit can produce equivalent final states with much smaller number of gates. Moreover, when deploying AQCEL with a noisy intermediate scale quantum computer, it efficiently produces a quantum circuit that approximates the original circuit with high fidelity by truncating low-amplitude computational basis states below certain thresholds. Our technique is useful for a wide variety of quantum algorithms, opening up new possibilities to further simplify quantum circuits to be more effective for real devices.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Evaluation of 3D Simulation Capabilities within xRAGE using a Suite of Test Problems

This presentation describes a 3D verification and validation suite for xRAGE. Three verification test problems for hydrodynamics are presented: the Kidder ball problem, the Verney shell problem, and a 5-material compression problem. There is excellent agreement between 2D and 3D XRAGE simulation results for these problems, and between the xRAGE results and the benchmark solutions. In addition, two 3D ICF test problems are presented, based on an OMEGA direct drive capsule experiment, and on a NIF indirect drive capsule experiment. These two ICF test problems are used to evaluate two hydrodynamic methods within xRAGE. In addition to hydrodynamics the simulations also include three temperature plasma physics, thermal conduction and radiation diffusion. When simulating either capsule in 3D the newer unsplit hydrodynamic method in xRAGE produces more vorticity relative to the older default method. For the indirect drive capsule the 3D simulations are in reasonable agreement with the experimental values of ion temperature and neutron production. Unstable 3D vortex rings are present near the DT ice/carbon interface which enhance the confinement of the DT ice.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Computational exploration of DIII-D QH-mode with extended MHD

This work will enable improved understanding of the scientific basis for DIII-D Quiescent H-mode (QH-mode) by leveraging recent computational developments of the NIMROD code. Experimentally, it is established that the QH-mode tokamak operation regime is free of edge-localized modes (ELMs) with good energy confinement. ELMs produce an impulsive heat-flux on the divertor that is unacceptable under burning-plasma conditions. Beyond DIII-D in future high-power-density tokamaks, accommodation of a high heat flux at the divertor represents a key challenge to overcome. Proposed solutions involve operation at high divertor density to achieve divertor detachment, operation with a liquid lithium wall and/or introduction of impurities to enhance radiative losses. Additionally, the first wall material is likely to be a high-Z material such as Tungsten in future devices as materials such as graphite tiles restrict the lifetime of the first wall due to a high erosion rate. While the atomic species composition of the plasma is known to impact the edge-pedestal structure and stability, the specifics of the underlying physics basis for this interaction is not firm. Studying this interaction is the goal of this proposal without focusing specifically on the fine details of the first wall and divertor solution. It is the goal is to lay the ground work to establish the compatibility of the divertor solutions and wall material with pedestal solutions that eliminate ELMs by leveraging present DIII-D experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Lagrangian particle model for 3D simulation of pellets and SPI fragments in tokamaks

A 3D numerical model for the ablation of pellets and shattered pellet injection (SPI) fragments in tokamaks in the plasma disruption mitigation and fueling parameter space has been developed based on the Lagrangian particle code [R. Samulyak, X. Wang, H.-S. Chen, Lagrangian Particle Method for Compressible Fluid Dynamics, J. Comput. Phys., 362 (2018), 1-19]. The pellet code implements the low magnetic Reynolds number MHD equations, kinetic models for the electronic heating, a pellet surface ablation model, an equation of state that supports multiple ionization states, radiation, and a model for grad-B drift of the ablated material across the magnetic field. The Lagrangian particle algorithm is highly adaptive, capable of simulating a large number of fragments in 3D while eliminating numerical difficulties of dealing with the tokamak background plasma. The code has achieved good agreement with theory for spherically symmetric ablation flows. Axisymmetric simulations of neon and deuterium pellets in magnetic fields ranging from 1 to 6 Tesla have been compared with previous simulations using the FronTier code, and very good agreement has also been obtained. Furthermore, the main physics contribution of the paper is a detailed study of the influence of 3D effects, in particular grad-B drift, on pellet ablation rates and properties of ablation clouds. Smaller reductions of ablation rates in magnetic fields compared to axially symmetric simulations have been demonstrated because the ablated material is not confined to narrowing channels in the presence of grad-B drift. Contribution of various factors in the grad-B drift model has also been quantified.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Stochastic gravitational wave background from global cosmic strings

Global cosmic strings are generically predicted in particle physics beyond the Standard Model, e.g., a post-inflationary global U(1) symmetry breaking which may associate with axion-like dark matter. Here, we demonstrate that although subdominant to Goldstone emission, gravitational waves (GWs) radiated from global strings can be observable with current/future GW detectors. The frequency spectrum of such GWs is also shown to be a powerful tool to probe the Hubble expansion rate of the Universe at times prior to the Big Bang nucleosynthesis where the standard cosmology has yet to be tested.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Heat transfer in directly-irradiated high-temperature solid–gas flows laden with polydisperse particles

Heat transfer in directly-irradiated high-temperature solid–gas flows laden with polydisperse particles is investigated using a novel transient three-dimensional computational fluid dynamics model. The model couples particle–gas hydrodynamics of solid–gas flows laden with polydisperse particles, radiative heat transfer in non-grey absorbing, emitting and anisotropically-scattering multi-component participating media, conduction heat transfer in the gas phase, and interfacial convection heat transfer. The multiphase particle-in-cell method is used to predict high-fidelity solid–gas flow characteristics, such as the local discrete particle size distribution, with increased computational efficiency by combining the advantages of both Eulerian and Lagrangian methods. The multi-component radiative transfer model is implemented using an advanced collision-based Monte Carlo ray-tracing method. The number of the prescribed discrete particle components is found to be the key parameter affecting the computational accuracy and efficiency, which primarily depends on the size distribution of the particles. For the model particle–gas flow featuring free-falling Gamma-distributed ceramic particles exposed to concentrated solar irradiation, the particle volume fraction, radiative, fluid flow and thermal characteristics appear to converge with the increasing number of the discrete particle components. Five particle components are sufficient to obtain physically meaningful results. A further increase in the number of the particle components only slightly increases the accuracy of the numerical predictions at the expense of a rapidly increasing computational time. For five particle components, the particle vertical velocity at the receiver exit for particles with the diameter of 43.4 μm is 57% of that for the particles with the diameter of 202.8 μm. The temperatures of these two particle components increase from the initial ambient values by factors of 2 and 1.2, respectively, during the simulation time. In conclusion, the model developed allows for increased fidelity of particle–gas flow simulations with significant radiative effects.

97 MATHEMATICS AND COMPUTING↗