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 289 records · Page 16
All sourcesNASA NTRSDOE OSTI
NASA NTRS1995

Post-entry and volcanic contaminant abundances of zinc, copper, selenium, germanium and gallium in stratospheric micrometeorites

Some fraction of Zn, Cu, Se, Ga and Ge in chondritic interplanetary dust particles (IDPs) collected in the lower stratosphere between 1981 May and 1984 June has a volcanic origin. I present a method to evaluate the extent of this unavoidable type of stratospheric contamination for individual particles. The mass-normalized abundances for Cu and Ge as a function of mass-normalized stratospheric residence time show their time-integrated stratospheric aerosol abundances. The Zn, Se and Ga abundances show a subdivision into two groups that span approximately two-year periods following the eruptions of the Mount St. Helens (1980 May) and El Chichon (1982 April) volcanoes. Elemental abundances in particles collected at the end of each two-year period indicate low, but not necessarily ambient, volcanic stratospheric abundances. Using this time-integrated baseline, I calculate the straospheric contaminant fractions in nine IDPs and show that Zn, SE and Ga abundances in chondritic IDPs derive in part from stratospheric aerosol contaminants. Post-entry elemental abundances (i.e., the amount that survived atmospheric entry heating of the IDP) show enrichments relative to the CI abundances but in a smaller number of particles than previously suggested.

Rietmeijer, Frans J. M.↗
DOE OSTI2025

Particle signal considerations for isotope ratio analysis with single particle multi-collector inductively coupled plasma mass spectrometry

Particle analysis has benefitted from the advent of single particle inductively coupled plasma-mass spectrometry (spICP-MS) due to its robustness, sensitivity, and high-throughput nature. Previous methods of spICP-MS have typically utilized quadrupole or time-of-flight mass analyzers and therefore employ electron multiplier-based detectors (such as secondary electron multipliers or microchannel plates). However, to obtain precise measurements on elemental or isotopic ratios within individual particles, multi-collector ICP-MS (MC-ICP-MS) can be used. Here, we investigate Ce isotope ratios, specifically 142 Ce/ 140 Ce, by spMC-ICP-MS using an all-Faraday cup collector array. Using 1 μm (diameter) cerium dioxide particles, integration times of the Faraday cup detectors were varied from 50–500 ms. The signal from the cerium isotopes in the particles was used to determine isotope ratios, which closely matched the expected natural isotopic abundances. Due to the signal decay response from the Faraday cups, the signal from particles lasts much longer than the expected 1–2 ms (up to 100 s of ms). To explore this effect on isotope ratio analysis, multiple ratio analysis methods were used to determine how to obtain optimal precision and accuracy. Relative differences were around 2% for methods that calculated isotope ratios from summing the total signal of an individual particle before calculating the ratio (rather than using every data point individually). It was found that summing all data points per particle, or integrating under the signal peak, yielded both accurate and precise isotope ratios within the particle population. Particles were also sampled off a solid substrate via microextraction, and isotope ratios were determined with relative differences of 0.13% to 9%. This demonstrates the ability to use spMC-ICP-MS to obtain isotope ratios on particles, with little to no relative difference in comparison to the expected ratio, even when operating Faraday detectors at fast 50 ms integration times.

Szakas, Sarah E. [Oak Ridge National Laboratory (O↗
DOE OSTI2021

Acquisition and Processing Considerations for Infrared Images of Rotating Turbine Blades

Abstract As designers aim to increase efficiency in gas turbines for aircraft propulsion and power generation, spatially resolved experimental measurements are needed to validate computational models and compare improvement gains of new cooling designs. Infrared (IR) thermography is one such method for obtaining spatially resolved temperature measurements. As technological advances in thermal detectors enable faster integration times, surface temperature measurements of rotating turbine blades become possible to capture including the smallest features. This paper outlines opportunities enabled by the latest IR detector technologies for capturing spatially resolved rotating blade temperatures, while also addressing some of the challenges of implementing IR for turbine rigs such as the one in the Steady Thermal Aero Research Turbine (START) Laboratory. This paper documents critical steps in achieving accurate measurements including calibration, integration times, spatial noise, and motion blur. From these results, recommendations are provided for achieving accurate IR measurements collected in a rotating turbine facility to study film cooling.

Engineering↗
DOE OSTI2020

Improving Solution Accuracy and Convergence for Stochastic Physics Parameterizations with Colored Noise

Stochastic parameterizations are used in numerical weather prediction and climate modeling to help capture the uncertainty in the simulations and improve their statistical properties. Convergence issues can arise when time integration methods originally developed for deterministic differential equations are applied naively to stochastic problems. In previous studies, it has been demonstrated that a correction term, known in stochastic analysis as the Itô correction, can help improve solution accuracy for various deterministic numerical schemes and ensure convergence to the physically relevant solution without substantial computational overhead. The usual formulation of the Itô correction is valid only when the stochasticity is represented by white noise. In this study, a generalized formulation of the Itô correction is derived for noises of any color. The formulation is applied to a test problem described by an advection–diffusion equation forced with a spectrum of fast processes. We present numerical results for cases with both constant and spatially varying advection velocities to show that, for the same time step sizes, the introduction of the generalized Itô correction helps to substantially reduce time integration error and significantly improve the convergence rate of the numerical solutions when the forcing term in the governing equation is rough (fast varying); alternatively, for the same target accuracy, the generalized Itô correction allows for the use of significantly longer time steps and, hence, helps to reduce the computational cost of the numerical simulation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

21 cm Power Spectrum Analysis of North Celestial Pole Observations with the Tianlai Dish Pathfinder Array

The Tianlai Dish Pathfinder Array (TDPA) is a radio interferometer designed to test techniques for 21 cm intensity mapping in the post-reionization Universe as a means of measuring large-scale cosmic structure. Using nine nights of observations targeting the North Celestial Pole field, totaling approximately 107 hr of integration time, we analyze data in the frequency range 700–800 MHz (corresponding to redshift z ∼ 0.9). We do the data format conversion, radio frequency interference flagging, calibration, imaging and point source subtraction, and foreground removal via Singular Value Decomposition. The spherically averaged power spectrum Δ 2 (k) is obtained. Furthermore, this work successfully establishes and validates a comprehensive data analysis framework for the TDPA. We identify key improvements including sky model refinement, increased integration time, and pipeline optimization that will enable future detection of the 21 cm signal through auto-correlation and cross-correlation with optical galaxy surveys.

cosmology: large-scale structure of universe↗
NASA NTRS2004

High-speed Imaging of Global Surface Temperature Distributions on Hypersonic Ballistic-Range Projectiles

The NASA-Ames ballistic range provides a unique capability for aerothermodynamic testing of configurations in hypersonic, real-gas, free-flight environments. The facility can closely simulate conditions at any point along practically any trajectory of interest experienced by a spacecraft entering an atmosphere. Sub-scale models of blunt atmospheric entry vehicles are accelerated by a two-stage light-gas gun to speeds as high as 20 times the speed of sound to fly ballistic trajectories through an 24 m long vacuum-rated test section. The test-section pressure (effective altitude), the launch velocity of the model (flight Mach number), and the test-section working gas (planetary atmosphere) are independently variable. The model travels at hypersonic speeds through a quiescent test gas, creating a strong bow-shock wave and real-gas effects that closely match conditions achieved during actual atmospheric entry. The challenge with ballistic range experiments is to obtain quantitative surface measurements from a model traveling at hypersonic speeds. The models are relatively small (less than 3.8 cm in diameter), which limits the spatial resolution possible with surface mounted sensors. Furthermore, since the model is in flight, surface-mounted sensors require some form of on-board telemetry, which must survive the massive acceleration loads experienced during launch (up to 500,000 gravities). Finally, the model and any on-board instrumentation will be destroyed at the terminal wall of the range. For these reasons, optical measurement techniques are the most practical means of acquiring data. High-speed thermal imaging has been employed in the Ames ballistic range to measure global surface temperature distributions and to visualize the onset of transition to turbulent-flow on the forward regions of hypersonic blunt bodies. Both visible wavelength and infrared high-speed cameras are in use. The visible wavelength cameras are intensified CCD imagers capable of integration times as short as 2 ns. The infrared camera uses an Indium Antimonide (InSb) sensor in the 3 to 5 micron band and is capable of integration times as short as 500 ns. The projectiles are imaged nearly head-on using expendable mirrors offset slightly from the flight path. The proposed paper will discuss the application of high-speed digital imaging systems in the NASA-Ames hypersonic ballistic range, and the challenges encountered when applying these systems. Example images of the thermal radiation from the blunt nose of projectiles flying at nearly 14 times the speed of sound will be given.

Wilder, Michael C.↗
NASA NTRS1993

Adaptive methods for nonlinear structural dynamics and crashworthiness analysis

The objective is to describe three research thrusts in crashworthiness analysis: adaptivity; mixed time integration, or subcycling, in which different timesteps are used for different parts of the mesh in explicit methods; and methods for contact-impact which are highly vectorizable. The techniques are being developed to improve the accuracy of calculations, ease-of-use of crashworthiness programs, and the speed of calculations. The latter is still of importance because crashworthiness calculations are often made with models of 20,000 to 50,000 elements using explicit time integration and require on the order of 20 to 100 hours on current supercomputers. The methodologies are briefly reviewed and then some example calculations employing these methods are described. The methods are also of value to other nonlinear transient computations.

Belytschko, Ted↗
NASA NTRS2015

A Microwave Radiometric Method to Obtain the Average Path Profile of Atmospheric Temperature and Humidity Structure Parameters and Its Application to Optical Propagation System Assessment

The values of the key atmospheric propagation parameters Ct2, Cq2, and Ctq are highly dependent upon the vertical height within the atmosphere thus making it necessary to specify profiles of these values along the atmospheric propagation path. The remote sensing method suggested and described in this work makes use of a rapidly integrating microwave profiling radiometer to capture profiles of temperature and humidity through the atmosphere. The integration times of currently available profiling radiometers are such that they are approaching the temporal intervals over which one can possibly make meaningful assessments of these key atmospheric parameters. Since these parameters are fundamental to all propagation conditions, they can be used to obtain Cn2 profiles for any frequency, including those for an optical propagation path. In this case the important performance parameters of the prevailing isoplanatic angle and Greenwood frequency can be obtained. The integration times are such that Kolmogorov turbulence theory and the Taylor frozen-flow hypothesis must be transcended. Appropriate modifications to these classical approaches are derived from first principles and an expression for the structure functions are obtained. The theory is then applied to an experimental scenario and shows very good results.

Atmospheric turbulence↗
NASA NTRSDate not supplied

An Improved Approach to the Predictability & Reliability of the Onset of Turbulence With Shocks

The construction of numerical schemes for (a) stable and accurate simulation of turbulence with strong shocks, and for (b) obtaining correct propagation speed of discontinuities in the presence of stiff source terms share one important ingredient – minimization of numerical dissipation while maintaining numerical stability. The dual requirements to achieve both numerical stability and minimal numerical dissipation are often conflicting since existing shock capturing schemes were designed mainly to be robust for rapidly developed turbulence-free flows and for shock waves without stiff source term. For the past two decades, Yee and collaborators have focused on an improved understanding of the nonlinear behavior of different high order shock-capturing methods. It was found that even very high order methods without proper nonlinear stability and numerical dissipation control can either numerically smear the onset of turbulence due to excess numerical dissipation, or induce (onset) numerical turbulence that is not physical turbulence due to lack of proper numerical dissipation to improve nonlinear stability for long time integration. Our approach is to combine (I) and (II) below for obtaining the physically correct onset of turbulence with shocks, including problems with stiff source terms: (I) Nonlinear dynamics is utilized to complement the traditional linearized stability theory (Yee & Sweby, Yee et al., Griffiths et al., Lafon & Yee, Yee, Wang et al., Kotov et al. 1990- 2015) in order to (i) Minimize numerically induced false transition to turbulence, (ii) Minimize numerical instability due to long time integration of turbulent flows, (iii) Minimize numerically induced standing wave solutions, and (iv) Minimize wrong propagation of speed of discontinuities due to the presence of stiff source terms. (II) Our recently developed physical preserving (structural preserving) high order methods with improved nonlinear stability & accuracy that are essential in minimizing spurious numerics are used.

HECC↗
DOE OSTI2021

A DG-IMEX Method for Two-moment Neutrino Transport: Nonlinear Solvers for Neutrino–Matter Coupling

Neutrino-matter interactions play an important role in core-collapse supernova (CCSN) explosions as they contribute to both lepton number and/or four-momentum exchange between neutrinos and matter, and thus act as the agent for neutrino-driven explosions. Due to the multiscale nature of neutrino transport in CCSN simulations, an implicit treatment of neutrino-matter interactions is desired, which requires solutions of coupled nonlinear systems in each step of the time integration scheme. In this paper we design and compare nonlinear iterative solvers for implicit systems with energy coupling neutrino-matter interactions commonly used in CCSN simulations. Specifically, we consider electron neutrinos and antineutrinos, which interact with static matter configurations through the Bruenn 85 opacity set. The implicit systems arise from the discretization of a nonrelativistic two-moment model for neutrino transport, which employs the discontinuous Galerkin (DG) method for phase-space discretization and an implicit-explicit (IMEX) time integration scheme. In the context of this DG-IMEX scheme, we propose two approaches to formulate the nonlinear systems — a coupled approach and a nested approach. For each approach, the resulting systems are solved with Anderson-accelerated fixedpoint iteration and Newton’s method. The performance of these four iterative solvers has been compared on relaxation problems with various degree of collisionality, as well as proto-neutron star deleptonization problems with several matter profiles adopted from spherically symmetric CCSN simulations. Here, numerical results suggest that the nested Anderson-accelerated fixed-point solver is more efficient than other tested solvers for solving implicit nonlinear systems with energy coupling neutrino-matter interactions.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1992

Eno-Osher schemes for Euler equations

The combination of the Osher approximate Riemann solver for the Euler equations and various ENO schemes is discussed for one-dimensional flow. The three basic approaches, viz. the ENO scheme using primitive variable reconstruction, either with Cauchy-Kowalewski procedure for time integration or the TVD Runge-Kutta scheme, and the flux-ENO method are tested on different shock tube cases. The shock tube cases were chosen to present a serious challenge to the ENO schemes in order to test their ability to capture flow discontinuities, such as shocks. Also the effect of the ordering of the eigen values, viz. natural or reversed ordering, in the Osher scheme is investigated. The ENO schemes are tested up to fifth order accuracy in space and time. The ENO-Osher scheme using the Cauchy-Kowalewski procedure for time integration is found to be the most accurate and robust compared with the other methods and is also computationally efficient. The tests showed that the ENO schemes perform reasonably well, but have problems in cases where two discontinuities are close together. In that case there are not enough points in the smooth part of the flow to create a non-oscillatory interpolation.

Vandervegt, Jacobus J.↗
DOE OSTI2023

A Vision for Coupling Operation of US Fusion Facilities with HPC Systems and the Implications for Workflows and Data Management

The operation of large US Department of Energy (DOE) research facilities, like the DIII-D National Fusion Facility, results in the collection of complex multi-dimensional scientific datasets, both experimental and model-generated. In the future, it is envisioned that integrated data analysis coupled with large-scale high performance computing (HPC) simulations will be used to improve experimental planning and operation. Practically, massive data sets from these simulations provide the physics basis for generation of both reduced semi-analytic and machine-learning-based models. Storage of both HPC simulation datasets (generated from US DOE leadership computing facilities) and experimental datasets presents significant challenges. In this paper, we present a vision for a DOE-wide data management workflow that integrates US DOE fusion facilities with leadership computing facilities. Data persistence and long-term availability beyond the length of allocated projects is essential, particularly for verification and recalibration of artificial intelligence and machine learning (AI/ML) models. Because these data sets are often generated and shared among hundreds of users across multiple leadership computing facility centers, they would benefit from cross-platform accessibility, persistent identifiers (e.g. DOI, or digital object identifier), and provenance tracking. Here, the ability to handle different data access patterns suggests that a combination of low cost, high latency (e.g. for storing ML training sets) and high cost, low latency systems (e.g. for real-time, integrated machine control feedback) may be needed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1983

Stress and Fracture Analyses Under Elastic-plastic and Creep Conditions: Some Basic Developments and Computational Approaches

A new hybrid-stress finite element algorith, suitable for analyses of large quasi-static deformations of inelastic solids, is presented. Principal variables in the formulation are the nominal stress-rate and spin. A such, a consistent reformulation of the constitutive equation is necessary, and is discussed. The finite element equations give rise to an initial value problem. Time integration has been accomplished by Euler and Runge-Kutta schemes and the superior accuracy of the higher order schemes is noted. In the course of integration of stress in time, it has been demonstrated that classical schemes such as Euler's and Runge-Kutta may lead to strong frame-dependence. As a remedy, modified integration schemes are proposed and the potential of the new schemes for suppressing frame dependence of numerically integrated stress is demonstrated. The topic of the development of valid creep fracture criteria is also addressed.

Reed, K. W.↗
NASA NTRS1976

Fluid dynamic research at NASA-Ames Research Center related to transonic wind tunnel design and testing techniques

Fluid dynamic research with the objective of developing new and improved technology in both test facility concepts and test techniques is being reported. A summary of efforts and results thus far obtained in four areas is presented. The four area are: (1) the use of heavy gases to obtain high Reynolds numbers at transonic speeds: (2) high Reynolds number tests of the C-141A wing configuration; (3) performance and flow quality of the pilot injector driven wind tunnel; and (4) integration time required to extract accurate static and dynamic data from tests in transonic wind tunnels. Some of the principal conclusions relative to each of the four areas are: (1) Initial attempts to apply analytical corrections to test results using gases with gamma other than 1.4 to simulate conditions in air show promise but need significant improvement; (2) for the C-141A configuration, no Reynolds number less than the full scale flight value provides an accurate simulation of the full scale flow; (3) high ratios of tunnel mass flow rate to injection mass flow rate and high flow quality can be obtained in an injector driven transonic wind tunnel; and (4) integration times of 0.5 to 1.0 sec may be required for static force and pressure tests, respectively, at some transonic test conditions in order to obtain the required data accuracy.

Muhlstein, L., Jr.↗
NASA NTRS1993

Variability of (O I) 6300-A emission near Io

Observations of (O I) 6300-A emission near Io have been obtained in 1990, 1991, and 1992 by the National Solar Observatory staff using the solar-stellar spectrography on the McMath-Pierce telescope at Kitt Peak. High-resolution spectra with a resolving power of about 1.2 x 10(exp 5) were obtained with an integration time of 10-15 min each. The viewing aperture for the observations was 5.2 arc sec x 5.2 arc sec centered on Io, with spatial resolution limited within this area by seeing conditions. Observations thus far have been reduced to obtain average brightness values over the aperture which range from approximately 200 to 1000 R for a number of different Io phase angles and Io system III longitudes. The (O I) 6300-A emission brightness exhibits an east-west asymmetry, where the average intensity for Io phase angles in the west (receding ansa) is 1.5 times brighter than in the east (approaching ansa). Similar east-west intensity ratios have also been observed for neutrals near Io and ions in the plasma torus for a number of other optical and ultraviolet emission lines which are excited by electron impact. In addition to the east-west asymmetry, the (O I) 6300-A emission brightness exhibits a strong dependence on the Io system III longitude angle, with a maximum value occurring in the range 200 deg +/- 50 deg. Earlier IUE observations of ultraviolet emission lines of O and S near Io obtained over a number of years have measured the east-west asymmetry, but the long IUE integration times of approximately 7-14 hours masked any detection of system III variability. For the (O I) 6300-A emission, the O(1D) state may be excited by electron impact of atomic oxygen and by electron impact dissociation of SO. The molecule SO may be present at the exobase or may be produced above the exobase as the dominant product of SO2 dissociation by electron impact. Preliminary assessment indicates that production of O(1D) by molecular dissociation may be more important. The (O I) 6300-A emission may therefore provide a remote signature for monitoring (1) the upward transport rates of molecular species in Io's atmosphere, (2) the relative abundance and time-variable dissociation of SO2 and/or SO at the exobase and in the corona of Io, and (3) the spatial distribution of these escaping molecular and atomic species and their ion production rates in the planetary magnetosphere.

Scherb, F.↗
DOE OSTI2022

Effect of Coaxial Electrode Geometry on the Electric Field Enhancement Factor for a High Voltage Vacuum Gap

We present an experimental analysis of the change in the electric field enhancement factor with varying gap size and penetration depth (P.D) of cathode into anode for a coaxial vacuum gap, diagnosed using Fowler–Nordheim analysis and optical imaging via scanning electron microscope (SEM) and time integrated Digital single lens reflex camera (DSLR). Data were collected on the Coaxial Gap Breakdown Machine (240 A, 25 kV, 150 ns, 0.1 Hz). Experiments using five different gap sizes at nine different P.Ds are compared over runs comprising 50 shots for each case. The results show a strong link between enhancement factor and gap size, with P.D and surface topology. For large gap sizes, 150, 330, and 700 μm, the average enhancement factor value increases with increasing P.D. For smaller gap sizes, 50 and 100 μm, the average enhancement factor decreases with P.D. SEM imaging before and after plasma formation for each gap size allows for quantifying surface finish, microprotrusion growth, average blast diameter, and an estimation of the surface area breakdowns occupy. Time integrated DSLR imaging analysis of the gap at each shot allows for a determination of the distribution of breakdowns about the circumference of the gap for each case tested. Here, the Fowler–Nordheim analysis allows for a quantitative analysis of the surface roughness of all gap sizes tested. Results show that for large gap sizes, the gap geometry and increasing area of breakdown is the main cause for increasing average enhancement factor. For small gap sizes, the dominant driving factor for small average enhancement factors—that subsequently decrease with P.D—is significant changes in surface topology due to an increased number of breakdowns.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1988

A modified ATI technique for nowcasting convective rain volumes over areas

This paper explores the applicability of the area-time-integral (ATI) technique for the estimation of the growth portion only of a convective storm (while the rain volume is computed using the entire life history of the event) and for nowcasting the total rain volume of a convective system at the stage of its maximum development. For these purposes, the ATIs were computed from the digital radar data (for 1981-1982) from the North Dakota Cloud Modification Project, using the maximum echo area (ATIA) no less than 25 dBz, the maximum reflectivity, and the maximum echo height as the end of the growth portion of the convective event. Linear regression analysis demonstrated that correlations between total rain volume or the maximum rain volume versus ATIA were the strongest. The uncertainties obtained were comparable to the uncertainties which typically occur in rain volume estimates obtained from radar data employing Z-R conversion followed by space and time integration. This demonstrates that the total rain volume of a storm can be nowcasted at its maximum stage of development.

Makarau, Amos↗