Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “time integrator”

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 757 records · Page 42
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Rocks, water, and noble liquids: Unfolding the flavor contents of supernova neutrinos

Measuring core-collapse supernova neutrinos, both from individual supernovae within the Milky Way and from past core collapses throughout the Universe (the diffuse supernova neutrino background, or DSNB), is one of the main goals of current and next generation neutrino experiments. Detecting the heavy-lepton flavor (muon and tau types, collectively v x ) component of the flux is particularly challenging due to small statistics and large backgrounds. Further, while the next galactic neutrino burst will be observed in a plethora of neutrino channels, allowing us to measure a small number of v x events, only upper limits are anticipated for the diffuse v x flux even after decades of data taking with conventional detectors. However, paleo detectors could measure the time-integrated flux of neutrinos from galactic core-collapse supernovae via flavor-blind neutral current interactions. In this work, we show how combining a measurement of the average galactic core-collapse supernova flux with paleo detectors and measurements of the DSNB electron-type neutrino fluxes with the next-generation water Cherenkov detector Hyper-Kamiokande and the liquid noble gas detector DUNE will allow to determine the mean supernova vx flux parameters with precision of order ten percent. Realizing this potential requires both the cosmic supernova rate out to z~1 and the integrated Galactic supernova rate over the last ~1 Gyr to be established at the ~10% level.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Realizability-preserving discontinuous Galerkin method for spectral two-moment radiation transport in special relativity

Here we present a realizability-preserving numerical method for solving a spectral two-moment model to simulate the transport of massless, neutral particles interacting with a steady background material moving with relativistic velocities. The model is obtained as the special relativistic limit of a four-momentum-conservative general relativistic two-moment model. Using a maximum-entropy closure, we solve for the Eulerian-frame energy and momentum. The proposed numerical method is designed to preserve moment realizability, which corresponds to moments defined by a nonnegative phase-space density. The realizability-preserving method is achieved with the following key components: (i) a discontinuous Galerkin phase-space discretization with specially constructed numerical fluxes in the spatial and energy dimensions; (ii) a strong stability-preserving implicit-explicit time-integration method; (iii) a realizability-preserving conserved to primitive moment solver; (iv) a realizability-preserving implicit collision solver; and (v) a realizability-enforcing limiter. Component (iii) is necessitated by the closure procedure, which closes higher order moments nonlinearly in terms of primitive moments. The nonlinear conserved to primitive and the implicit collision solves are formulated as fixed-point problems, which are solved with custom iterative solvers designed to preserve the realizability of each iterate. With a series of numerical tests, we demonstrate the accuracy and robustness of this discontinuous-Galerkin-implicit-explicit method.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Search for Λ − Λ ¯ oscillation in J / ψ → Λ Λ ¯ decay

Using ( 10087 ± 44 ) × 10 6 J / ψ decays collected by the BESIII detector at the BEPCII collider, we search for baryon number violation via Λ − Λ ¯ oscillation in the decay J / ψ → Λ Λ ¯ . No evidence for Λ − Λ ¯ oscillation is observed. The upper limit on the time-integrated probability of Λ − Λ ¯ oscillation is estimated to be 1.4 × 10 − 6 , corresponding to an oscillation parameter less than 2.1 × 10 − 18 GeV at 90% confidence level. Published by the American Physical Society 2025

Ablikim, M.↗
DOE OSTI2020

Temperature distributions and gradients in laser-heated plasmas relevant to magnetized liner inertial fusion

We present two-dimensional temperature measurements of magnetized and unmagnetized plasma experiments performed at Z relevant to the preheat stage in magnetized liner inertial fusion. The deuterium gas fill was doped with a trace amount of argon for spectroscopy purposes, and time-integrated spatially resolved spectra and narrow-band images were collected in both experiments. The spectrum and image data were included in two separate multiobjective analysis methods to extract the electron temperature spatial distribution T e ( r , z ) . The results indicate that the magnetic field increases T e , the axial extent of the laser heating, and the magnitude of the radial temperature gradients. Comparisons with simulations reveal that the simulations overpredict the extent of the laser heating and underpredict the temperature. Temperature gradient scale lengths extracted from the measurements also permit an assessment of the importance of nonlocal heat transport.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Time-resolved investigation of subnanosecond radiation from Al wire hybrid X pinches

K-shell x-ray spectra from Al wire hybrid X pinches have been studied using an x-ray streak camera with better than 0.1-ns time resolution together with a Focusing Spectrograph with Spatial Resolution (FSSR) spectrograph. High-intensity radiation with a continuumlike spectrum was observed in the subnanosecond initial phase of the x-ray pulse generated by the hybrid X pinch (HXP). Here, the absence of spectral lines in this phase and the extremely small x-ray source size indicates the importance of radiative processes in the final phase implosion dynamics. Plasma parameters in the following phases of the HXP were determined from analysis of the line intensities. Point-projection radiography together with a slit-step wedge camera and an FSSR spectrograph without time resolution were used to show the number of radiation sources, and to give information on the time-integrated photon energy spectrum.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Field-Dependent Ionic Conductivities from Generalized Fluctuation-Dissipation Relations

We derive a relationship for the electric field dependent ionic conductivity in terms of fluctuations of time integrated microscopic variables. We demonstrate this formalism with molecular dynamics simulations of solutions of differing ionic strength with implicit solvent conditions and molten salts. These calculations are aided by a novel nonequilibrium statistical reweighting scheme that allows for the conductivity to be computed as a continuous function of the applied field. In strong electrolytes, we find the fluctuations of the ionic current are Gaussian, and subsequently, the conductivity is constant with applied field. In weaker electrolytes and molten salts, we find the fluctuations of the ionic current are strongly non-Gaussian, and the conductivity increases with applied field. Furthermore, this nonlinear behavior, known phenomenologically for dilute electrolytes as the Onsager-Wien effect, is general and results from the suppression of ionic correlations at large applied fields, as we elucidate through both dynamic and static correlations within nonequilibrium steady states.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Nonuniform Absorption and Scattered Light in Direct-Drive Implosions Driven by Polarization Smoothing

Laser-direct-drive symmetric implosions on OMEGA illuminate a target with 60 laser beams and are designed to produce spherical implosions. Each beam is smoothed using orthogonal polarizations obtained by passing through distributed polarization rotators (DPR’s). Observations of the scattered light from OMEGA implosions do not show the expected symmetry and have much larger variation than standard predictions. For the first time, we have quantified the scattered-light nonuniformity from individual beams and identified the DPR’s as the source of the enhanced nonuniformity. An instrument was invented that isolated and measured the variation in the intensity and polarization of the light scattered from each OMEGA beam. The asymmetric intensity and polarization measurements are explained when the on-target offsets between the two orthogonal polarizations produced by the DPR’s are modeled using a 3-D cross-beam energy transfer (CBET) code that tracks the polarizations of each beam. The time-integrated nonuniformity in laser absorption and scattered light due to CBET and the DPR polarization offsets during high-performance OMEGA implosions is predicted to be significant and dominated by low spherical harmonic mode numbers. Further, the nonuniformity is predicted to be greatly reduced by replacing the DPR’s with new optics that create smaller offsets.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Hot Spot Evolution Measured by High-Resolution X-Ray Spectroscopy at the National Ignition Facility

Evolution of the hot spot plasma conditions was measured using high-resolution x-ray spectroscopy at the National Ignition Facility. The capsules were filled with DD gas with trace levels of Kr and had either a high-density-carbon (HDC) ablator or a tungsten (W)-doped HDC ablator. Time-resolved measurement of the Kr He β spectra, absolutely calibrated by a simultaneous time-integrated measurement, allows inference of the electron density and temperature through observing Stark broadening and the relative intensities of dielectronic satellites. By matching the calculated hot spot emission using a collisional-radiative code to experimental observations, the hot spot size and areal density are determined. Here these advanced spectroscopy techniques further reveal the effect of W dopant in the ablator on the hot spot parameters for their improved implosion performance.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Measurements of the Branching Fraction, Polarization, and C P Asymmetry for the Decay B 0 → ω ω

We present a comprehensive study of B 0 → ω ω decays using 772 × 10 6 B B ¯ pairs collected with the Belle detector at the KEKB e + e − collider. This process is a suppressed charmless decay into two vector mesons and can exhibit interesting polarization and C P violation. The decay is observed for the first time with a significance of 7.9 standard deviations. We measure a branching fraction B = ( 1.53 ± 0.29 ± 0.17 ) × 10 − 6 , a fraction of longitudinal polarization f L = 0.87 ± 0.13 ± 0.13 , and a time-integrated C P asymmetry A C P = − 0.44 ± 0.43 ± 0.11 , where the first uncertainties listed are statistical and the second are systematic. This is the first observation of B 0 → ω ω and the first measurements of f L and A C P for this decay. Published by the American Physical Society 2024

Physics↗
DOE OSTI2025

High-resolution in-situ characterization of laser powder bed fusion via transmission X-ray microscopy at X-ray free electron lasers

In this work, we describe the instrumentation used to perform the first operando transmission X-ray microscopy (TXM) and simultaneous X-ray diffraction of laser melting simulating laser powder bed fusion on the XCS instrument at the Linac Coherent Light Source (LCLS) X-ray free-electron laser (XFEL). Our TXM with 40× magnification in the X-ray regime at 11 keV gave spatial resolutions down to 940 nm per line pair, with effective pixel sizes down to 206 nm, image integration times of <100 fs, and frame rates tunable between 2.1 and 119 ns for two probe frames (0.48 GHz to 8.4 MHz). Images were recorded on Zyla and Icarus (UXI) detectors to trade off between spatial resolution and time dynamics. A 1 kW CW IR laser was coupled into the interaction point to conduct pump–probe studies of laser melting and solidification dynamics. Our temporal and spatial resolution with attenuation-based contrast exceeds that currently possible with synchrotron-based high-speed radiography. This system was sensitive to feature velocities of 10–12000 m s −1 but we did not observe any motion in this range in the laser melting of Al6061 alloy. Shockwaves were not observed and hot cracking proceeded at velocities below the detection limits. Pore accumulation was observed between successive shots, indicating that bubble escape mechanisms were not active. With proper experimental design, the spatial resolution, contrast and field of view could be further improved or modified. The increased brightness and narrower bandwidth of the XFEL allowed for this imaging technique and it lays the groundwork for a wide range of operando techniques to study additive manufacturing.

47 OTHER INSTRUMENTATION↗
DOE OSTI2021

Voltage Loss Comparison in CdSe/CdTe Solar Cells and Polycrystalline CdSeTe Heterostructures

Cd(Se)Te solar cells have considerable headroom to increase voltage. Voltage losses occur due to incomplete absorption above bandgap E_g and band tail absorption below E_g; such losses are quantified using radiative voltage. The largest voltage losses are attributed to nonradiative recombination, which is quantified via carrier lifetime and radiative efficiency. We compare radiative voltage, radiative efficiency, and carrier lifetime for Cu-doped CdSe/CdTe solar cells and for undoped polycrystalline CdSeTe heterostructures passivated with Al_2O_3. Using external quantum efficiency spectrum and a CdSeTe absorption spectrum obtained from absolute photoluminescence (PL), we show that the radiative voltage is greater than 1.1 V. Time-resolved PL experiments and modeling show that a major part of voltage losses can be attributed to recombination in the absorber bulk. The front interface recombination makes a larger impact within the first few nanoseconds after pulsed excitation, and the comparison of time-gated and time-integrated PL can be used to assess relative contributions of front interface and bulk recombination rates.

14 SOLAR ENERGY↗
DOE OSTI2020

Indoor air quality in California homes with code-required mechanical ventilation

Data were collected in 70 detached houses built in 2011-2017 in compliance with the mechanical ventilation requirements of California's building energy efficiency standards. Each home was monitored for a 1-week period with windows closed and the central mechanical ventilation system operating. Pollutant measurements included time-resolved fine particulate matter (PM 2.5 ) indoors and outdoors and formaldehyde and carbon dioxide (CO 2 ) indoors. Time-integrated measurements were made for formaldehyde, NO 2 , and nitrogen oxides (NO X ) indoors and outdoors. Operation of the cooktop, range hood, and other exhaust fans was continuously recorded during the monitoring period. Onetime diagnostic measurements included mechanical airflows and envelope and duct system air leakage. All homes met or were very close to meeting the ventilation requirements. On average, the dwelling unit ventilation fan moved 50% more airflow than the minimum requirement. Pollutant concentrations were similar to or lower than those reported in a 2006-2007 study of California new homes built in 2002-2005. Mean and median indoor concentrations were lower by 44% and 38% for formaldehyde and 44% and 54% for PM 2.5 . Ventilation fans were operating in only 26% of homes when first visited, and the control switches in many homes did not have informative labels as required by building standards.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2020

Indoor air quality in new and renovated low-income apartments with mechanical ventilation and natural gas cooking in California

This work presents pollutant concentrations and performance data for code-required mechanical ventilation equipment in 23 low-income apartments at 4 properties constructed or renovated 2013-2017. All apartments had natural gas cooking burners. Occupants pledged to not use windows for ventilation during the study but several did. Measured airflows of range hoods and bathroom exhaust fans were lower than product specifications. Only eight apartments operationally met all ventilation code requirements. Pollutants measured over one week in each apartment included time-resolved fine particulate matter (PM 2.5 ), nitrogen dioxide (NO 2 ), formaldehyde and carbon dioxide (CO 2 ) and time-integrated formaldehyde, NO 2 and nitrogen oxides (NO X ). Compared to a recent study of California houses with code-compliant ventilation, apartments were smaller, had fewer occupants, higher densities, and higher mechanical ventilation rates. Mean PM 2.5 , formaldehyde, NO 2 , and CO 2 were 7.7 µg/m 3 , 14.1 ppb, 18.8 ppb, and 741 ppm in apartments; these are 4% lower, 25% lower, 165% higher, and 18% higher compared to houses with similar cooking frequency. Four apartments had weekly PM 2.5 above the California annual outdoor standard of 12 µg/m 3 and also discrete days above the World Health Organization 24-h guideline of 25 µg/m 3 . Two apartments had weekly NO 2 above the California annual outdoor standard of 30 ppb.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2020

A Holistic Algorithmic Approach to Improving Accuracy, Robustness, and Computational Efficiency for Atmospheric Dynamics

Atmospheric weather and climate models must perform simulations very quickly to be useful. Therefore, modelers have traditionally focused on reducing computations as much as possible. However, in our new era of increasingly compute-capable hardware, data movement is now the prohibiting expense. This study examines the computational benefits of a new algorithmic approach to modeling atmospheric dynamics on scales relevant to weather and climate simulation. Rather than minimizing computations, this new approach considers the larger problem more holistically, including spatial accuracy, temporal accuracy, robustness (i.e., oscillations), on-node efficiency, and internode data transfers together at once. Numerical experiments demonstrate how computations can be strategically increased to simultaneously address each of these constraints while reducing data movement to adapt to modern accelerated hardware. The new algorithm can achieve at times up to 80% peak floating point throughput in single precision on the Nvidia Tesla V100 GPU, where the traditional approach is shown to only achieve single-digit floating point efficiency. Further, the new algorithm is twice as fast as a standard Runge--Kutta time integrator, and high-order accuracy with Weighted Essentially Non-Oscillatory (WENO) limiting came at less than 30% additional runtime cost on a GPU, thus increasing the accuracy per degree of freedom.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Elastic Wave Propagation in Curvilinear Coordinates with Mesh Refinement Interfaces by a Fourth Order Finite Difference Method

In this work, we develop a fourth order accurate finite difference method for the three dimensional elastic wave equation in isotropic media with the piecewise smooth material property. In our model, the material property can be discontinuous at curved interfaces. The governing equations are discretized in second order form on curvilinear meshes by using a fourth order finite difference operator satisfying a summation-by-parts property. The method is energy stable and high order accurate. The highlight is that mesh sizes can be chosen according to the velocity structure of the material so that computational efficiency is improved. At the mesh refinement interfaces with hanging nodes, physical interface conditions are imposed by using ghost points and interpolation. With a fourth order predictor-corrector time integrator, the fully discrete scheme is energy conserving. Numerical experiments are presented to verify the fourth order convergence rate and the energy conserving property.

58 GEOSCIENCES↗
DOE OSTI2022

Performance Analysis of Speculative Parallel Adaptive Local Timestepping for Conservation Laws

Stable simulation of conservation laws, such as those used to model fluid dynamics and plasma physics applications, requires the satisfaction of the so-called Courant-Friedrichs-Lewy condition. By allowing regions of the mesh to advance with different timesteps that locally satisfy this stability constraint, significant work reduction can be attained when compared to a time integration scheme using a single timestep size. However, parallelizing this algorithm presents considerable difficulty. Since the stability condition depends on the state of the system, dependencies become dynamic and potentially non-local. In this article, we present an adaptive local timestepping algorithm using an optimistic (Timewarp-based) parallel discrete event simulation. We introduce waiting heuristics to limit misspeculation and a semi-static load balancing scheme to eliminate load imbalance as parts of the mesh require finer or coarser timesteps. Last, we outline an interface for separating the physics of the specific conservation law from the temporal integration allowing for productive adoption of our proposed algorithm. We present a misspeculation study for three conservation laws, demonstrating both the productivity of the local timestepping API, for which 74% of the lines of code are reused across different conservation laws, and the robustness of the waiting heuristics—at most 1.5% of element updates are rolled back. Our performance studies demonstrate up to a 2.8× speedup versus a baseline unoptimized local timestepping approach, a 4x improvement in per-node throughput compared to an MPI parallelization of synchronous timestepping, and scalability up to 3,072 cores on NERSC’s Cori Haswell partition.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

FIERRO V.X

FIERRO is a modern C++ code intended to simulate quasi-static solid mechanics problems and transient, compressible material dynamic problems with Lagrangian methods, which have meshes with constant mass elements that move with the material. FIERRO is designed to aid a) material model research that has historically been done using commercial implicit and explicit finite element codes, b) numerical methods research, and c) computer science research. FIERRO supports user developed material models that adhere to several industry standard formats by using a C++ to Fortran interface to couple the model to the numerical solvers. FIERRO is built on the ELEMENTS library that supports a diverse suite of element types, including high-order elements, and quadrature rules. The mesh class within the ELEMENTS library is designed for efficient calculations on unstructured meshes and to minimize memory usage. FIERRO is designed to readily accommodate a range of numerical methods including continuous finite element, finite volume, and discontinuous Galerkin methods. FIERRO is designed to support explicit and implicit time integration methods. FIERRO only supports a single material in an element. No physical data exists within the code.

Morgan, Nathaniel↗
DOE OSTI2022

Phase Field Dislocation Dynamics (PFDD) version 2.x

This disclosure is for version 2.x of a mesoscale model called Phase Field Dislocation Dynamics (PFDD). PFDD is used for investigating deformation in nanoscale (grain sizes of ~300 nm and less) materials, such as metals and alloys. This approach models the motion and interaction of individual defects, namely dislocations, in the material using scalar-valued phase field variables, also called order parameters. The system is evolved through energy minimization thus the model calculates the total energy density in terms of the phase field variables. The energy minimization is completed using the Ginzburg-Landau equation, and is implemented with explicit time integration. The total system energy can be comprised of several terms, including the strain energy (which describes dislocation-dislocation interactions), the energy due to an applied stress (dislocation interactions with the applied stress), and a core/lattice (perfect dislocations) or generalized stacking fault (partial dislocations) energy (described the dislocation core structure). The latter term in particular may vary based on the crystal structure being modeled and is typically informed using lower length scale (e.g., atomistic) approaches, although no such (atomistic) calculations are completed within the PFDD algorithm. This basic formulation was previously reviewed by Los Alamos National Laboratory and released under license number C17113. This previously reviewed version we will henceforth refer to as PFDD v1.0. PFDD v1.0 consisted of 2 codes (one parallel and one serial) plus input files, all written in the C language. This new disclosure is addressing the next versions of the PFDD, versions 2.x. There have been several enhancements of PFDD v1.0, which are described here and included in the attached code, which we will refer to as PFDD v2.0. There are also several new features described here that are either planned or already in process and are expected to be subsequent releases, i.e., v2.1, v2.2, ...v2.x.

Hunter, Abigail↗