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.
Simplified Solar Modulation Model of Inner Trapped Belt Proton Flux As a Function of Atmospheric Density
No simple algorithm seems to exist for calculating proton fluxes and lifetimes in the Earth's inner, trapped radiation belt throughout the solar cycle. Most models of the inner trapped belt in use depend upon AP8 which only describes the radiation environment at solar maximum and solar minimum in Cycle 20. One exception is NOAAPRO which incorporates flight data from the TIROS/NOAA polar orbiting spacecraft. The present study discloses yet another, simple formulation for approximating proton fluxes at any time in a given solar cycle, in particular between solar maximum and solar minimum. It is derived from AP8 using a regression algorithm technique from nuclear physics. From flux and its time integral fluence, one can then approximate dose rate and its time integral dose. It has already been published in this journal that the absorbed dose rate, D, in the trapped belts exhibits a power law relationship, D = A(rho)(sup -n), where A is a constant, rho is the atmospheric density, and the index n is weakly dependent upon shielding. However, that method does not work for flux and fluence. Instead, we extend this idea by showing that the power law approximation for flux J is actually bivariant in energy E as well as density rho. The resulting relation is J(E,rho)approx.(sum of)A(E(sup n))rho(sup -n), with A itself a power law in E. This provides another method for calculating approximate proton flux and lifetime at any time in the solar cycle. These in turn can be used to predict the associated dose and dose rate.
Chasing Gamma-Ray Signals from Binary Neutron Star Coalescences with the Cherenkov Telescope Array: Prospects and Observing Strategies
The detection of gravitational waves (GWs) from a binary neutron star (BNS) merger by Advanced LIGO and Advanced Virgo (GW170817), together with its electromagnetic counterpart, the short gamma-ray burst GRB 170817A, heralded the birth of multimessenger astronomy. The detection of TeV emission from GRBs motivates follow-up observations with the Cherenkov Telescope Array Observatory (CTAO), which is ideal for detecting such signals due to its unprecedented sensitivity, rapid response, and wide-field survey capabilities. The aim of this work is to evaluate GeV–TeV GW follow-up strategies for CTAO using a multistep simulation pipeline and to estimate the expected rate of joint GW–GRB detections during observing run O5. Using a simulated sample of BNS systems with corresponding GW detections, gamma-ray emission is simulated through phenomenological prescriptions based on the observed population of short GRBs, including off-axis jet scenarios. CTAO observations are simulated to account for instrument response, sky tiling strategies, integration times, and varying observing conditions. Strategies with variable and constant integration times are investigated. We find that, via an optimized follow-up strategy, about 5% of simulated GW-associated short GRBs produce GeV–TeV radiation detectable by CTAO. Detectability is strongly influenced by the jet opening angle and viewing angle, suggesting that even rough estimates of the viewing angle in GW alerts could enhance targeting. This framework motivates future follow-ups of GW-detectable events, including neutron star–black hole mergers, and further supports the development of advanced strategies incorporating galaxy distributions and synergies with future detectors such as the Einstein Telescope.
Solar Modulation of Inner Trapped Belt Radiation Flux as a Function of Atmospheric Density
No simple algorithm seems to exist for calculating proton fluxes and lifetimes in the Earth's inner, trapped radiation belt throughout the solar cycle. Most models of the inner trapped belt in use depend upon AP8 which only describes the radiation environment at solar maximum and solar minimum in Cycle 20. One exception is NOAAPRO which incorporates flight data from the TIROS/NOAA polar orbiting spacecraft. The present study discloses yet another, simple formulation for approximating proton fluxes at any time in a given solar cycle, in particular between solar maximum and solar minimum. It is derived from AP8 using a regression algorithm technique from nuclear physics. From flux and its time integral fluence, one can then approximate dose rate and its time integral dose.
Three-dimensional flux-split Euler schemes involving unstructured dynamic meshes
Improved algorithms for the solution of the 3-D time dependent Euler equations are presented for aerodynamic analysis involving unstructured dynamic meshes. The improvements were developed recently to the spatial and temporal discretizations used by unstructured grid flow solvers. The spatial discretization involves a flux split approach which is naturally dissipative and captures shock waves sharply with at most one grid point within the shock structure. The temporal discretization involves either an explicit time integration scheme using a multistage Runge-Kutta procedure or an implicit time integration scheme using a Gauss-Seidel relaxation procedure which is computationally efficient for either steady or unsteady flow problems. With the implicit Gauss-Seidel procedure, very large time steps may be used for rapid convergence to steady state, and the step size for unsteady cases may be selected for temporal accuracy rather than for numerical stability. Steady flow results are presented for both the NACA 0012 airfoil and the ONERA M6 wing to demonstrate applications of the new Euler solvers. A description of the Euler solvers is presented along with results and comparisons which assess the capability.
From Simple Labels to Time-Use Integrations: Supporting the Spectrum of Qualitative Travel Behavior Data
In transportation research, applications of travel behavior data collection are context-specific and require different types of qualitative inputs. These inputs can be viewed as spanning a spectrum of user burden and data quality, from simple trip labels to complex time-use surveys. However, each currently active smartphone-based travel diary platform appears to only support one type of qualitative input, and the effort required for customization is unclear. In this paper, we characterize the spectrum by defining four canonical use cases: (i) trip labels, (ii) trip questionnaire, (iii) counterfactual trips, and (iv) time-use surveys. We then outline a mechanism for supporting configurable user inputs on the same underlying smartphone-based sensing mechanism and demonstrate that it can support all the use cases without any code changes. We further demonstrate that the flexible data model that underpins this mechanism can enable real-time monitoring and analysis. Finally, we evaluate per-user data collection and engagement metrics for large-scale deployments of three canonical use cases, spanning 10 programs, 435 users, and 251,041 trips, and a maximum duration of 800 days. Future efforts may support additional use cases through an expanded configuration and provide greater insight into user engagement. We hope that these insights enable the research community to look at qualitative inputs through a new lens and experiment with novel use cases to fill in the spectrum.
Measured Indoor Nitrogen Oxides in Homes with a Child with Asthma and Gas or Induction Electric Cooking
The Cooking Energy and Ventilation Impacts on Children's Asthma (CEVICA) study is a randomized control trial investigating the effects of replacing gas with induction electric ranges in the homes of children with asthma in California's San Joaquin Valley. Indoor air quality parameters and respiratory health indicators were measured over 2-week intensive at Baseline and after consecutive 3-month study phases, with half getting electric cooking at the start of Phase 1 and others in Phase 2. Indoor air measurements included time-integrated NO 2 and NO X by passive sampler and time-resolved NO 2 by electrochemical sensors. In the first 40 homes, data were collected during 61 intensives for gas cooking and 55 for induction. Cooking was identified using temperature sensors above the cooktop. Pollutant events were identified from sharp rises in concentrations Time integrated NOX species were lower with electric cooking, with mean differences of 14.2 ppb for NO 2 (95% CI: 9.6, 18.7; p < 0.001), 55.6 ppb for NO X (95% CI: 30.4, 80.9; p < 0.001), and 41.5 ppb for derived NO (95% CI: 19.7, 63.2; p < 0.001). A paired, within home analysis showed larger reductions in Group 2 (gas to induction) than Group 1 that remained electric across Phase 1 and Phase 2 for NO 2 (mean difference in As of 10.5 ppb; 95% CI: 5.6, 15.4; p < 0.001), NO X (32.3 ppb; 95% CI: 10.0, 54.5; p=0.007), and derived NO (21.8 ppb; 95% CI: 3.6, 40.0; p=0.022). Compared to electric cooking, gas had much higher rates of associated NO 2 events, and larger above-baseline NO 2 peaks. These preliminary results are consistent with prior findings that shifting to induction cooking can substantially lower indoor NO 2 compared with gas cooking.
Near-field supersonic flow predictions by an adaptive unstructured tetrahedral grid solver
Applicability of a three-dimensional solution adaptive unstructured tetrahedral Euler flow solver about generic models for near-field sonic boom pressure signature predictions is evaluated. Comparisons of computational and experimental data demonstrates the capability of the method for predicting inviscid solutions useful for high speed calculations about simple 3-D geometries. The approach has promising features and results indicate potential for application to more complex configurations. The mesh generation is based on the advancing front technique, and steady state solutions of the Euler equations are achieved by explicit time integration. Spatial discretization uses the Taylor-Galerkin approach; an alternate time integration, based on the Runge-Kutta method, is also included. The solution-adaptive grid procedure is based on either remeshing or mesh refinement techniques. An alternative geometry-adaptive grid procedure has also been incorporated.
On Orbit Receiver Performance Assessment of the Geoscience Laser Altimeter System (GLAS) on ICESAT
The GLAS instrument on the NASA's ICESat mission has provided over a billion measurements of the Earth surface elevation and atmosphere backscattering at both 532 and 1064-nm wavelengths. The receiver performance has stayed nearly unchanged since ICESat launch in January 2003. The altimeter receiver has achieved a less than 3-cm ranging accuracy when excluding the effects of the laser beam pointing angle determination uncertainties. The receiver can also detect surface echoes through clouds of one-way transmission as low as 5%. The 532-nm atmosphere backscattering receiver can measure aerosol and clouds with cross section as low as 1e-7/m.sr with a 1 second integration time and molecular backscattering from upper atmosphere with a 60 second integration time. The 1064-nm atmosphere backscattering receiver can measure aerosol and clouds with a cross section as low as 4e-6/m.sr. This paper gives a detailed assessment of the GLAS receiver performance based on the in-orbit calibration tests.
A novel integrated time-resolved array avalanche photodiode detection system for nuclear resonant scattering measurements
Here, the nuclear resonant scattering (NRS) experiment requires photon-counting detectors with high time resolution, short dead time, large dynamic range, low noise, and large detection area. An 8-channel avalanche photodiode (APD) array detector system with high integrity, flexibility, and reliability has been developed to adapt to the demands of NRS experiments. The detector system mainly consists of four key parts: (i) an array-APD sensor, (ii) 8-channel integrated fast preamplifiers, (iii) the time-to-digital converter readout electronics, and (iv) a data acquisition system and EPICS support software. Remarkably, the system exhibits a time resolution of better than 500 ps and has a sufficiently low noise level, allowing for the lowest detection energy threshold of 4 keV. The performance of the new array-APD system as well as its real application in nuclear forward scattering (NFS) and nuclear resonant inelastic x-ray scattering (NRIXS) experiments was tested in two synchrotron facilities. With the new system, the NFS signal very close to the prompt electronic scattering signal can be extracted. Thanks to the customized EPICS-areaDetector-based control software, NRIXS spectra can be readily measured with time and energy information of the NRIXS signal stored in the raw data, which is promising for developing NRIXS data analysis in the time domain. The array-APD detector can be deployed for nuclear resonant scattering experiments at various synchrotron radiation facilities.
Accuracy of an unstructured-grid upwind-Euler algorithm for the ONERA M6 wing
Improved algorithms for the solution of the three-dimensional, time-dependent Euler equations are presented for aerodynamic analysis involving unstructured dynamic meshes. The improvements have been developed recently to the spatial and temporal discretizations used by unstructured-grid flow solvers. The spatial discretization involves a flux-split approach that is naturally dissipative and captures shock waves sharply with at most one grid point within the shock structure. The temporal discretization involves either an explicit time-integration scheme using a multistage Runge-Kutta procedure or an implicit time-integration scheme using a Gauss-Seidel relaxation procedure, which is computationally efficient for either steady or unsteady flow problems. With the implicit Gauss-Seidel procedure, very large time steps may be used for rapid convergence to steady state, and the step size for unsteady cases may be selected for temporal accuracy rather than for numerical stability. Steady flow results are presented for both the NACA 0012 airfoil and the Office National d'Etudes et de Recherches Aerospatiales M6 wing to demonstrate applications of the new Euler solvers. The paper presents a description of the Euler solvers along with results and comparisons that assess the capability.
An iterative dynamic chemical stiffness removal method for reacting flow simulations
Abstract An iterative dynamic chemical stiffness removal method (IDCSR) based on quasi-steady-state approximation (QSSA) is proposed. The IDCSR method is built on a previously developed non-iterative method which has proved to work well for small timestep sizes. A novel iterative procedure is designed in IDCSR to enable explicit time integration of stiff chemistry at relatively large timestep sizes relevant to practical reacting flow simulations. The effectiveness of the iterative procedure is first demonstrated with a toy problem and homogeneous auto-ignition with fixed integration step sizes, showing that larger timestep sizes can be allowed for explicit time integration using IDCSR compared with the previous non-iterative method. IDCSR is then compared with existing explicit chemistry solvers for simulations of homogeneous auto-ignition and shows similar or lower computational cost but significantly higher accuracy across a wide range of timestep sizes. IDCSR is further combined with an automatic adaptive time-stepping scheme for simulations of 0-D homogeneous auto-ignition and a 2-D laminar lifted n -dodecane jet flame. For the 0-D auto-ignition simulations, IDCSR is shown to reduce both the error (by 43%–90%) and computational cost (by 6–15 times) compared with existing explicit solvers, while achieving speed-up factors of up to 400 compared with VODE for a wide range of timestep sizes and reaction mechanisms. For the 2-D jet flame simulations, speed-up factors of 15 and 31 for chemistry integration, and 5 and 9 for overall simulation, are achieved by IDCSR compared with CVODE with and without analytic Jacobian, respectively.
A Fundamentally New Coupled Approach to Contact Mechanics via the Dirichlet‐Neumann Schwarz Alternating Method
Contact phenomena are crucial for understanding the behavior of mechanical systems. However, existing computational approaches for simulating mechanical contact often face numerical challenges, such as inaccurate physical predictions, energy conservation errors, and unwanted oscillations. Here, we introduce an alternative technique for simulating dynamic contact based on the non‐overlapping Schwarz alternating method, originally developed for domain decomposition. In multibody contact scenarios, this method treats each body as a separate, non‐overlapping domain and prevents interpenetration using an alternating Dirichlet–Neumann iterative process. This approach has a strong theoretical foundation, eliminates the need for contact constraints, and offers flexibility, making it ideal for multiscale and multiphysics applications. We conducted a numerical comparison between the Schwarz method and traditional methods, such as the Lagrange multiplier and penalty methods, focusing on a benchmark impact problem. Our results indicate that the Schwarz alternating method outperforms traditional methods in several key areas: it provides more accurate predictions for various measurable quantities and demonstrates exceptional energy conservation capabilities. To address unwanted oscillations in contact velocities and forces, we explored various algorithms and stabilization techniques, ultimately opting for the naïve‐stabilized Newmark scheme for its simplicity and effectiveness. Additionally, we validated the efficiency of the Schwarz method in a three‐dimensional impact problem, highlighting its inherent capacity to accommodate different mesh topologies, time‐integration schemes, and time steps for each interacting body.
Probing atomic physics at ultrahigh pressure using laser-driven implosions
Abstract Spectroscopic measurements of dense plasmas at billions of atmospheres provide tests to our fundamental understanding of how matter behaves at extreme conditions. Developing reliable atomic physics models at these conditions, benchmarked by experimental data, is crucial to an improved understanding of radiation transport in both stars and inertial fusion targets. However, detailed spectroscopic measurements at these conditions are rare, and traditional collisional-radiative equilibrium models, based on isolated-atom calculations and ad hoc continuum lowering models, have proved questionable at and beyond solid density. Here we report time-integrated and time-resolved x-ray spectroscopy measurements at several billion atmospheres using laser-driven implosions of Cu-doped targets. We use the imploding shell and its hot core at stagnation to probe the spectral changes of Cu-doped witness layer. These measurements indicate the necessity and viability of modeling dense plasmas with self-consistent methods like density-functional theory, which impact the accuracy of radiation transport simulations used to describe stellar evolution and the design of inertial fusion targets.
Addendum: Unified framework for open quantum dynamics with memory
This Addendum presents a detailed analysis of the discretization error in time-integration and time-derivative that appear in the Nakajima-Zwanzig equation. This was brought to our attention by Makri et al. [arXiv:2410.08239]. Our analysis in the Addendum shows that the relationship derived in our earlier work [Nat. Commun. 15, 8087 (2024)] is valid within the choice of discretization and is not contaminated by the discretization error.
A time-resolved, in-chamber x-ray pinhole imager for Z
We have commissioned a new time-resolved, x-ray imaging diagnostic for the Z facility. The primary intended application is for diagnosing the stagnation behavior of Magnetized Liner Inertial Fusion (MagLIF) and similar targets. We have a variety of imaging systems at Z, both time-integrated and time-resolved, that provide valuable x-ray imaging information, but no system at Z up to this time provides a combined high-resolution imaging with multi-frame time resolution; this new diagnostic, called TRICXI for Time Resolved In-Chamber X-ray Imager, is meant to provide time-resolved spatial imaging with high resolution. The multi-frame camera consists of a microchannel plate camera. A key component to achieving the design goals is to place the instrument inside the Z vacuum chamber within 2 m of the load, which necessitates a considerable amount of x-ray shielding as well as a specially designed, independent vacuum system. Here, a demonstration of the imaging capability for a series of MagLIF shots is presented. Predictions are given for resolution and relative image irradiance to guide experimenters in choosing the desired configuration for their experiments.
VERITAS : A density-functional theory-based multiband kinetic model for understanding x-ray spectroscopy of dense plasmas
X-ray spectroscopy has long been a powerful diagnostic tool for hot, dilute plasmas, providing insights into plasma conditions by measuring line shifts and broadenings of atomic transitions. The technique critically depends on the accuracy of atomic physics models used to interpret spectroscopic measurements for inferring plasma properties such as free-electron density and temperature. Over the past decades, the atomic and plasma physics communities have developed robust atomic physics models to account for various processes in hot, dilute classical plasmas. While these models have been successful in that regime, their applicability becomes uncertain when interpreting x-ray spectroscopy experiments of above-solid-density plasmas. Given that finite-temperature density-functional theory (DFT) offers a more accurate description of dense plasma environments, we present the development of a DFT-based multi-band kinetic model, VERITAS, designed to improve the interpretation of x-ray spectroscopic measurements in high-density plasmas produced by laser-driven spherical implosions. This work details the VERITAS model and its application to both time-integrated and time-resolved x-ray spectra from implosion experiments on OMEGA. The advantages and limitations of the VERITAS model will also be discussed, along with potential directions for advancing x-ray spectroscopy of dense and superdense plasmas.
Advances in phosphor two-color ratio method thermography for full-field surface temperature measurements
We report thermographic phosphors can be employed for optical sensing of surface, gas phase, and bulk material temperatures through different strategies including the time-decay method, time-integrated method, and frequency-domain method. We focus on the time-integrated method, also known as the ratio method, as it can be more practical in many situations. This work advances the ratio method using two machine vision cameras with CMOS detectors for full-field temperature measurements of a solid surface. A phosphor calibration coupon is fabricated using aerosol deposition and employed for in situ determination of the temperature-versus-intensity ratio relationship. Algorithms from digital image correlation are employed to determine the stereoscopic imaging system intrinsic and extrinsic parameters, and accurately register material points on the sample to subpixel locations in each image with 0.07 px or better accuracy. Detector nonlinearity is carefully characterized and corrected. Temperature-dependent, spatial non-uniformity of the full-field intensity ratio—posited to be caused by a blue-shift effect of the bandpass filter for non-collimated light and/or a wavelength-dependent transmission efficiency of the lens—is assessed and treated for cases where a standard flat-field correction fails to correct the non-uniformity. In sum, pixel-wise calibration curves relating the computed intensity ratio to temperature in the range of T = 300–430 K are generated, with an embedded error of less than 3 K. This work offers a full calibration methodology and several improvements on two-color phosphor thermography, opening the door for full-field temperature measurements in dynamic tests with deforming test articles.