SEARCH · Search NASA

Results for “time integration”

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

SourceResearch connectionFrom yearThrough year

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

At least 199 records · Page 11

All sources NASA NTRS DOE OSTI

NASA NTRS1990

An interactive adaptive remeshing algorithm for the two-dimensional Euler equations

An interactive adaptive remeshing algorithm utilizing a frontal grid generator and a variety of time integration schemes for the two-dimensional Euler equations on unstructured meshes is presented. Several device dependent interactive graphics interfaces have been developed along with a device independent DI-3000 interface which can be employed on any computer that has the supporting software including the Cray-2 supercomputers Voyager and Navier. The time integration methods available include: an explicit four stage Runge-Kutta and a fully implicit LU decomposition. A cell-centered finite volume upwind scheme utilizing Roe's approximate Riemann solver is developed. To obtain higher order accurate results a monotone linear reconstruction procedure proposed by Barth is utilized. Results for flow over a transonic circular arc and flow through a supersonic nozzle are examined.

Slack, David C.↗

NASA NTRS1999

Entropy Splitting and Numerical Dissipation

A rigorous stability estimate for arbitrary order of accuracy of spatial central difference schemes for initial-boundary value problems of nonlinear symmetrizable systems of hyperbolic conservation laws was established recently by Olsson and Oliger (1994) and Olsson (1995) and was applied to the two-dimensional compressible Euler equations for a perfect gas by Gerritsen and Olsson (1996) and Gerritsen (1996). The basic building block in developing the stability estimate is a generalized energy approach based on a special splitting of the flux derivative via a convex entropy function and certain homogeneous properties. Due to some of the unique properties of the compressible Euler equations for a perfect gas, the splitting resulted in the sum of a conservative portion and a non-conservative portion of the flux derivative. hereafter referred to as the "Entropy Splitting." There are several potential desirable attributes and side benefits of the entropy splitting for the compressible Euler equations that were not fully explored in Gerritsen and Olsson. The paper has several objectives. The first is to investigate the choice of the arbitrary parameter that determines the amount of splitting and its dependence on the type of physics of current interest to computational fluid dynamics. The second is to investigate in what manner the splitting affects the nonlinear stability of the central schemes for long time integrations of unsteady flows such as in nonlinear aeroacoustics and turbulence dynamics. If numerical dissipation indeed is needed to stabilize the central scheme, can the splitting help minimize the numerical dissipation compared to its un-split cousin? Extensive numerical study on the vortex preservation capability of the splitting in conjunction with central schemes for long time integrations will be presented. The third is to study the effect of the non-conservative proportion of splitting in obtaining the correct shock location for high speed complex shock-turbulence interactions. The fourth is to determine if this method can be extended to other physical equations of state and other evolutionary equation sets. If numerical dissipation is needed, the Yee, Sandham, and Djomehri (1999) numerical dissipation is employed. The Yee et al. schemes fit in the Olsson and Oliger framework.

Yee, H. C.↗

NASA NTRS1974

A simulation of synthetic aperture radar imaging of ocean waves

A simulation of radar imaging of ocean waves with synthetic aperture techniques is presented. The modelling is simplistic from the oceanographic and electromagnetic viewpoint in order to minimize the computational problems, yet reveal some of the physical problems associated with the imaging of moving ocean waves. The model assumes: (1) The radar illuminates a one-dimensional, one harmonic ocean wave. (2) The scattering is assumed to be governed by geometrical optics. (3) The radar is assumed to be down-looking, with Doppler processing (range processing is suppressed due to the one-dimensional nature of the problem). (4) The beamwidth of the antenna (or integration time) is assumed to be sufficiently narrow to restrict the specular points of the peaks and troughs of the wave. The results show that conventional processing of the image gives familiar results if the ocean waves are stationary. When the ocean wave dispersion relationship is satisfied, the image is smeared due to the motion of the specular points over the integration time. In effect, the image of the ocean is transferred to the near field of the synthetic aperture.

Swift, C. T.↗

NASA NTRS1972

On the possibility of Earth reentry simulation of shallow angle Jupiter entry

Possible Earth reentry simulation of shallow-angle (3 to 30 deg) Jupiter entry was investigated in terms of four parameters of the bow shock layer ahead of a blunt vehicle: peak (equilibrium) temperature, peak pressure, peak inward radiative flux, and time-integrated radiative flux. It is shown that simulation ranging from fair to good can be achieved. Reentry is easiest (lowest Earth reentry speed) at steep Earth reentry, in the Earth entry speed range of 15 to 22 km/sec, for both the Jupiter nominal and cool atmospheres. Increased Earth speed is required, generally, for increasing Jupiter entry angle, and for increasing temperature, radiative flux, time-integrated flux and pressure, in that order.

Mcronald, A. D.↗

NASA NTRS1978

DACOM - A rapid, high sensitivity airborne carbon monoxide monitor

A rapid, high sensitivity differential absorption carbon monoxide monitor (DACOM) has been developed to measure CO from an aircraft. The differential absorption measurement is accomplished by rapidly current tuning the wavelength of a PbSSe diode laser across a single 4.7 micron CO absorption line. As the laser wavelength sweeps periodically across the absorption line, the laser beam is amplitude modulated and is then synchronously detected. With this technique a minimum detectable differential transmission of 0.02% is obtained for a 1 second integration time. This performance is equivalent to a CO sensitivity of 5 ppb for a 10 meter absorption path. The techniques used to achieve this high sensitivity will be discussed. In the aircraft installation a 10 meter absorption path is located in the airstream along the fuselage of a C-54 aircraft. During flight tests of DACOM, a sensitivity of 20 ppb for a 1-second integration time has been achieved.

Sachse, G. W.↗

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 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↗

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.↗

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.↗

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↗

DOE OSTI2025

Optimizing Charge-coupled Device Readout Enabled by the Floating-gate Amplifier

Multiple-Amplifier Sensing (MAS) charge-coupled devices (CCDs) have recently been shown to be promising silicon detectors that meet noise sensitivity requirements for next generation Stage-5 spectroscopic surveys and potentially, future space-based imaging of extremely faint objects on missions such as the Habitable Worlds Observatory. Building upon the capability of the Skipper CCD to achieve deeply sub-electron noise floors, MAS CCDs utilize multiple floating-gate amplifiers along the serial register to increase the readout speed by a factor of the number of output nodes compared to a Skipper CCD. We introduce and experimentally demonstrate on a 16-channel prototype device new readout techniques that exploit the MAS CCD’s floating-gate amplifiers to optimize the correlated double sampling by resetting once per line instead of once per pixel. With this new mode, we find an optimal filter to subtract the noise from the signal during read out. We also take advantage of the MAS CCD’s structure to tune the read time by independently changing integration times for the signal and reference level. Together with optimal weighted averaging of the 16 outputs, these approaches enable us to reach a sub-electron noise of 0.9 e − rms pix −1 at 19 μs pix −1 for a single charge measurement per pixel—simultaneously giving a 30% faster readout time and 10% lower read noise compared to performance previously evaluated without these techniques.

Lin, Kenneth W. (ORCID:0000000189672281)↗