Search NASA⌕ Search

SEARCH · Search NASA

Results for “EXTRAPOLATION”

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 361 records · Page 20
All sourcesNASA NTRSDOE OSTI
NASA NTRS1994

Sonic boom predictions using a solution-adaptive full-potential code

The TranAir full-potential code, which utilizes a solution-adaptive, unstructured Cartesian grid, was used to compute near-field pressure signatures of three generic geometries: (1) a cone-cylinder; (2) a low aspect ratio rectangular wing; and (3) a delta-wing/body. Pressure signatures computed by TranAir are extrapolated to near- and midfield distances using a quasilinear extrapolation technique. Results are compared with wind-tunnel data and, where possible, with extrapolated wind-tunnel data and Whitham-method results. In general, extrapolated TranAir data is in good agreement with experimental and Whitham-method pressure signatures, but peak pressures are slightly underpredicted. It is felt that data computed too close to the model caused some of the discrepancies. The flowfield grid structure was also found to be an important issue in obtaining accurate off-body data. Results were encouraging enough to proceed with the application of TranAir to conceptual configurations in support of the NASA high speed research program (HSRP).

Madson, Michael D.↗
NASA NTRS1996

Sonic Boom Propagation Codes Validated by Flight Test

The sonic boom propagation codes reviewed in this study, SHOCKN and ZEPHYRUS, implement current theory on air absorption using different computational concepts. Review of the codes with a realistic atmosphere model confirm the agreement of propagation results reported by others for idealized propagation conditions. ZEPHYRUS offers greater flexibility in propagation conditions and is thus preferred for practical aircraft analysis. The ZEPHYRUS code was used to propagate sonic boom waveforms measured approximately 1000 feet away from an SR-71 aircraft flying at Mach 1.25 to 5000 feet away. These extrapolated signatures were compared to measurements at 5000 feet. Pressure values of the significant shocks (bow, canopy, inlet and tail) in the waveforms are consistent between extrapolation and measurement. Of particular interest is that four (independent) measurements taken under the aircraft centerline converge to the same extrapolated result despite differences in measurement conditions. Agreement between extrapolated and measured signature duration is prevented by measured duration of the 5000 foot signatures either much longer or shorter than would be expected. The duration anomalies may be due to signature probing not sufficiently parallel to the aircraft flight direction.

Poling, Hugh W.↗
NASA NTRS1994

Accurate Evaluation of Quantum Integrals

Combining an appropriate finite difference method with Richardson's extrapolation results in a simple, highly accurate numerical method for solving a Schr\"{o}dinger's equation. Important results are that error estimates are provided, and that one can extrapolate expectation values rather than the wavefunctions to obtain highly accurate expectation values. We discuss the eigenvalues, the error growth in repeated Richardson's extrapolation, and show that the expectation values calculated on a crude mesh can be extrapolated to obtain expectation values of high accuracy.

Galant, David C.↗
NASA NTRS1995

Accurate Evaluation of Quantum Integrals

Combining an appropriate finite difference method with Richardson's extrapolation results in a simple, highly accurate numerical method for solving a Schrodinger's equation. Important results are that error estimates are provided, and that one can extrapolate expectation values rather than the wavefunctions to obtain highly accurate expectation values. We discuss the eigenvalues, the error growth in repeated Richardson's extrapolation, and show that the expectation values calculated on a crude mesh can be extrapolated to obtain expectation values of high accuracy.

Galant, D. C.↗
NASA NTRS2002

Wind Tunnel Measurements of Shuttle Orbiter Global Heating with Comparisons to Flight

An aerothermodynamic database of global heating images was acquired of the Shuttle Orbiter in the NASA Langley Research Center 20-Inch Mach 6 Air Tunnel. These results were obtained for comparison to the global infrared images of the Orbiter in flight from the infrared sensing aeroheating flight experiment (ISAFE). The most recent ISAFE results from STS-103, consisted of port side images, at hypersonic conditions, of the surface features that result from the strake vortex scrubbing along the side of the vehicle. The wind tunnel results were obtained with the phosphor thermography system, which also provides global information and thus is ideally suited for comparison to the global flight results. The aerothermodynamic database includes both windward and port side heating images of the Orbiter for a range of angles of attack (20 to 40 deg), freestream unit Reynolds number (1 x 10(exp 6))/ft to 8 x 10(exp 6)/ft, body flap deflections (0, 5, and 10 deg), speed brake deflections (0 and 45 deg), as well as with boundary layer trips for forced transition to turbulence heating results. Sample global wind tunnel heat transfer images were extrapolated to flight conditions for comparison to Orbiter flight data. A windward laminar case for an angle of attack of 40 deg was extrapolated to Mach 11.6 flight conditions for comparison to STS-2 flight thermocouple results. A portside wind tunnel image for an angle of attack of 25 deg was extrapolated for Mach 5 flight conditions for comparison to STS-103 global surface temperatures. The comparisons showed excellent qualitative agreement, however the extrapolated wind tunnel results over-predicted the flight surface temperatures on the order of 5% on the windward surface and slightly higher on the portside.

Berry, Scott A.↗
NASA NTRS1999

Loss tolerant speech decoder for telecommunications

A method and device for extrapolating past signal-history data for insertion into missing data segments in order to conceal digital speech frame errors. The extrapolation method uses past-signal history that is stored in a buffer. The method is implemented with a device that utilizes a finite-impulse response (FIR) multi-layer feed-forward artificial neural network that is trained by back-propagation for one-step extrapolation of speech compression algorithm (SCA) parameters. Once a speech connection has been established, the speech compression algorithm device begins sending encoded speech frames. As the speech frames are received, they are decoded and converted back into speech signal voltages. During the normal decoding process, pre-processing of the required SCA parameters will occur and the results stored in the past-history buffer. If a speech frame is detected to be lost or in error, then extrapolation modules are executed and replacement SCA parameters are generated and sent as the parameters required by the SCA. In this way, the information transfer to the SCA is transparent, and the SCA processing continues as usual. The listener will not normally notice that a speech frame has been lost because of the smooth transition between the last-received, lost, and next-received speech frames.

Prieto, Jr., Jaime L.↗
NASA NTRS2011

Bounding Species Distribution Models

Species distribution models are increasing in popularity for mapping suitable habitat for species of management concern. Many investigators now recognize that extrapolations of these models with geographic information systems (GIS) might be sensitive to the environmental bounds of the data used in their development, yet there is no recommended best practice for "clamping" model extrapolations. We relied on two commonly used modeling approaches: classification and regression tree (CART) and maximum entropy (Maxent) models, and we tested a simple alteration of the model extrapolations, bounding extrapolations to the maximum and minimum values of primary environmental predictors, to provide a more realistic map of suitable habitat of hybridized Africanized honey bees in the southwestern United States. Findings suggest that multiple models of bounding, and the most conservative bounding of species distribution models, like those presented here, should probably replace the unbounded or loosely bounded techniques currently used [Current Zoology 57 (5): 642-647, 2011].

Stohlgren, Thomas J.↗
DOE OSTI2015

Accurate ab initio vibrational energies of methyl chloride

Two new nine-dimensional potential energy surfaces (PESs) have been generated using high-level ab initio theory for the two main isotopologues of methyl chloride, CH{sub 3}{sup 35}Cl and CH{sub 3}{sup 37}Cl. The respective PESs, CBS-35{sup HL}, and CBS-37{sup HL}, are based on explicitly correlated coupled cluster calculations with extrapolation to the complete basis set (CBS) limit, and incorporate a range of higher-level (HL) additive energy corrections to account for core-valence electron correlation, higher-order coupled cluster terms, scalar relativistic effects, and diagonal Born-Oppenheimer corrections. Variational calculations of the vibrational energy levels were performed using the computer program TROVE, whose functionality has been extended to handle molecules of the form XY {sub 3}Z. Fully converged energies were obtained by means of a complete vibrational basis set extrapolation. The CBS-35{sup HL} and CBS-37{sup HL} PESs reproduce the fundamental term values with root-mean-square errors of 0.75 and 1.00 cm{sup −1}, respectively. An analysis of the combined effect of the HL corrections and CBS extrapolation on the vibrational wavenumbers indicates that both are needed to compute accurate theoretical results for methyl chloride. We believe that it would be extremely challenging to go beyond the accuracy currently achieved for CH{sub 3}Cl without empirical refinement of the respective PESs.

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

Quantum calculation of inelastic co collisions with H. III. Rate coefficients for RO-vibrational transitions

We present calculated rate coefficients for ro-vibrational transitions of CO in collisions with H atoms for a gas temperature range of 10 K ≤ T ≤ 3000 K, based on the recent three-dimensional ab initio H–CO interaction potential of Song et al. Rate coefficients for ro-vibrational v=1,j=0−30→v{sup ′}=0,j{sup ′} transitions were obtained from scattering cross sections previously computed with the close-coupling (CC) method by Song et al. Combining these with the rate coefficients for vibrational v=1−5→v{sup ′}<v quenching obtained with the infinite-order sudden approximation, we propose a new extrapolation scheme that yields the rate coefficients for ro-vibrational v=2−5,j=0−30→v{sup ′},j{sup ′} de-excitation. Cross sections and rate coefficients for ro-vibrational v=2,j=0−30→v{sup ′}=1,j{sup ′} transitions calculated with the CC method confirm the effectiveness of this extrapolation scheme. Our calculated and extrapolated rates are very different from those that have been adopted in the modeling of many astrophysical environments. The current work provides the most comprehensive and accurate set of ro-vibrational de-excitation rate coefficients for the astrophysical modeling of the H–CO collision system. The application of the previously available and new data sets in astrophysical slab models shows that the line fluxes typically change by 20%–70% in high temperature environments (800 K) with an H/H{sub 2} ratio of 1; larger changes occur for lower temperatures.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Observational Evidence of Magnetic Reconnection for Brightenings and Transition Region Arcades in IRIS Observations

By using a new method of forced-field extrapolation, we study the emerging flux region AR11850 observed by the Interface Region Imaging Spectrograph and Solar Dynamical Observatory . Our results suggest that the bright points (BPs) in this emerging region exhibit responses in lines formed from the upper photosphere to the transition region, which have relatively similar morphologies. They have an oscillation of several minutes according to the Atmospheric Imaging Assembly data at 1600 and 1700 Å. The ratio between the BP intensities measured in 1600 and 1700 Å filtergrams reveals that these BPs are heated differently. Our analysis of the Helioseismic and Magnetic Imager vector magnetic field and the corresponding topology in AR11850 indicates that the BPs are located at the polarity inversion line and most of them are related to magnetic reconnection or cancelation. The heating of the BPs might be different due to different magnetic topology. We find that the heating due to the magnetic cancelation would be stronger than the case of bald patch reconnection. The plasma density rather than the magnetic field strength could play a dominant role in this process. Based on physical conditions in the lower atmosphere, our forced-field extrapolation shows consistent results between the bright arcades visible in slit-jaw image 1400 Å and the extrapolated field lines that pass through the bald patches. It provides reliable observational evidence for testing the mechanism of magnetic reconnection for the BPs and arcades in the emerging flux region, as proposed in simulation studies.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Comparison of Two Coronal Magnetic Field Models to Reconstruct a Sigmoidal Solar Active Region with Coronal Loops

Magnetic field extrapolation is an important tool to study the three-dimensional (3D) solar coronal magnetic field, which is difficult to directly measure. Various analytic models and numerical codes exist, but their results often drastically differ. Thus, a critical comparison of the modeled magnetic field lines with the observed coronal loops is strongly required to establish the credibility of the model. Here we compare two different non-potential extrapolation codes, a nonlinear force-free field code (CESE–MHD–NLFFF) and a non-force-free field (NFFF) code, in modeling a solar active region (AR) that has a sigmoidal configuration just before a major flare erupted from the region. A 2D coronal-loop tracing and fitting method is employed to study the 3D misalignment angles between the extrapolated magnetic field lines and the EUV loops as imaged by SDO /AIA. It is found that the CESE–MHD–NLFFF code with preprocessed magnetogram performs the best, outputting a field that matches the coronal loops in the AR core imaged in AIA 94 Å with a misalignment angle of ∼10°. This suggests that the CESE–MHD–NLFFF code, even without using the information of the coronal loops in constraining the magnetic field, performs as good as some coronal-loop forward-fitting models. For the loops as imaged by AIA 171 Å in the outskirts of the AR, all the codes including the potential field give comparable results of the mean misalignment angle (∼30°). Thus, further improvement of the codes is needed for a better reconstruction of the long loops enveloping the core region.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Measurement of the central exclusive production of charged particle pairs in proton-proton collisions at $\sqrt{s}$ = 200 GeV with the STAR detector at RHIC

We report on the measurement of the Central Exclusive Production of charged particle pairs h+h- (h = π, K, p) with the STAR detector at RHIC in proton-proton collisions at $\sqrt{s}$ = 200 GeV. The charged particle pairs produced in the reaction pp → p' + h + h - + p' are reconstructed from the tracks in the central detector and identified using the specific energy loss and the time of flight method, while the forward-scattered protons are measured in the Roman Pot system. Exclusivity of the event is guaranteed by requiring the transverse momentum balance of all four final-state particles. Differential cross sections are measured as functions of observables related to the central hadronic final state and to the forward-scattered protons. They are measured in a fiducial region corresponding to the acceptance of the STAR detector and determined by the central particles’ transverse momenta and pseudorapidities as well as by the forward-scattered protons’ momenta. This fiducial region roughly corresponds to the square of the four-momentum transfers at the proton vertices in the range 0.04 GeV 2 < -t 1 , -t 2 < 0.2 GeV 2 , invariant masses of the charged particle pairs up to a few GeV and pseudorapidities of the centrally-produced hadrons in the range |η| < 0.7. The measured cross sections are compared to phenomenological predictions based on the Double Pomeron Exchange (DPE) model. Structures observed in the mass spectra of π + π - and K + K - pairs are consistent with the DPE model, while angular distributions of pions suggest a dominant spin-0 contribution to π + π - production. For π + π - production, the fiducial cross section is extrapolated to the Lorentz-invariant region, which allows decomposition of the invariant mass spectrum into continuum and resonant contributions. The extrapolated cross section is well described by the continuum production and at least three resonances, the f 0 (980), f 2 (1270) and f 0 (1500), with a possible small contribution from the f 0 (1370). Fits to the extrapolated differential cross section as a function of t 1 and t 2 enable extraction of the exponential slope parameters in several bins of the invariant mass of π + π - pairs. These parameters are sensitive to the size of the interaction region.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Accurate Effective Stress Measures: Predicting Creep Life for 3D Stresses Using 2D and 1D Creep Rupture Simulations and Data

Operating structural components experience complex loading conditions resulting in 3D stress states. Current design practice estimates multiaxial creep rupture life by mapping a general state of stress to a uniaxial creep rupture correlation using effective stress measures. The data supporting the development of effective stress measures are nearly always only uniaxial and biaxial, as 3D creep rupture tests are not widely available. This limitation means current effective stress measures must extrapolate from 2D to 3D stress states, potentially introducing extrapolation error. In this work, we use a physics-based, crystal plasticity finite element model to simulate uniaxial, biaxial, and triaxial creep rupture. Here, we use the virtual dataset to assess the accuracy of current and novel effective stress measures in extrapolating from 2D to 3D stresses and also explore how the predictive accuracy of the effective stress measures might change if experimental 3D rupture data was available. We confirm these conclusions, based on simulation data, against multiaxial creep rupture experimental data for several materials, drawn from the literature. The results of the virtual experiments show that calibrating effective stress measures using triaxial test data would significantly improve accuracy and that some effective stress measures are more accurate than others, particularly for highly triaxial stress states. Results obtained using experimental data confirm the numerical findings and suggest that a unified effective stress measure should include an explicit dependence on the first stress invariant, the maximum tensile principal stress, and the von Mises stress.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Not so HOT Triangulations

Here, we propose primal–dual mesh optimization algorithms that overcome shortcomings of the standard algorithm while retaining some of its desirable features. “Hodge-Optimized Triangulations” defines the “HOT energy” as a bound on the discretization error of the diagonalized Delaunay Hodge star operator. HOT energy is a natural choice for an objective function, but unstable for both mathematical and algorithmic reasons: it has minima for collapsed edges, and its extrapolation to non-regular triangulations is inaccurate and has unbounded minima. We propose a different extrapolation with a stronger theoretical foundation, and avoid extrapolation by recalculating the objective just beyond the flip threshold. We propose new objectives, based on normalizations of the HOT energy, with barriers to edge collapses and other undesirable configurations. We propose mesh improvement algorithms coupling these. When HOT optimization nearly collapses an edge, we actually collapse the edge. Otherwise, we use the barrier objective to update positions and weights and remove vertices. By combining discrete connectivity changes with continuous optimization, we more fully explore the space of possible meshes and obtain higher quality solutions.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Data reduction considerations for the burning velocity of spherical constant volume flames of R32 (CH 2 F 2 ) with air

Here, the present work explores data reduction techniques for the measurement of the laminar burning velocities of R32(CH 2 F 2 )-air mixtures using a constant volume combustion device, in which the pressure-time history is the only measured parameter. To allow clear assessment of the accuracy of the data reduction methods, the pressure-time histories used for analysis are synthetically generated via a detailed numerical simulation employing full kinetics and with and without an optically-thin radiation model. Various data reduction models are employed, including a two-zone model and two multi-zone models, and these are compared with the results from the burning velocity obtained from the output of the numerical simulation. The data reduction schemes are shown to be accurate if the same radiation model is employed in the data reduction as was used in the flame simulation to generate the pressure trace used for post-processing. If the incorrect radiation model is employed, however, the errors can be quite large. The effects of stretch, radiation, and different data post-processing methodologies are explored and the errors quantified. Stretch is shown to be important for the early stages and the selected data range that is used for extrapolation has a significant effect on the extrapolated burning velocity. However, with an appropriate choice of data considered for extrapolation, the prediction of the unstretched burning velocity can be quite accurate.

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

An empirical calibration of the serpentine-water oxygen isotope fractionation at T = 25–100°C

The ability to precisely constrain the physical–chemical conditions of serpentinization, such as temperature and fluid source, is limited by the accuracy of calibrations for oxygen isotope fractionation between serpentine and water – i.e., 1000 lnα( serpentine-water ) – which disagree by up to 20‰ when extrapolated to T<200°C. In this study, we present a new empirical calibration of 1000 lnα( serpentine-water ) aiming to improve applications of oxygen isotope geochemistry to very low-T serpentinization (T<100°C). We used the high-spatial resolution capabilities of Secondary Ion Mass Spectrometry (SIMS) to analyze oxygen isotope ratios in mineral pairs of calcite + serpentine, quartz + serpentine and talc + serpentine co-crystallized at scales≤50μm in six serpentinite samples from the Samail ophiolite (Oman). SIMS analysis shows that the mineral pairs are relatively homogenous in oxygen isotope ratios with variability in δ 18 O values ≤2‰ (2s). Clumped isotope thermometry and petrological constraints indicate crystallization temperatures from 25 to 100°C for the investigated samples. These independent constraints on temperature allowed us to derive 1000 lnα( serpentine-water ) by combining Δ 18 O( quartz-serpentine ), Δ 18 O( talc-serpentine ) and Δ 18 O( calcite-serpentine ) measured by SIMS in the investigated samples with Δ 18 O( quartz-water ), Δ 18 O( talc-water ) and Δ 18 O( calcite-water ) calibrations available in the literature. Our empirical calibration of 1000 lnα( serpentine-water ) = 1.04±0.20 (2SE) × 10 6 /T 2 (T in K), from T = 25 to 100°C, is within uncertainty of former high-temperature empirical calibrations extrapolated to T<100°C and matches experimental calibrations when extrapolated to T=250°C. In conclusion, the new serpentine-water calibration enables more accurate reconstructions of fluid-rock interactions occurring during low-temperature serpentinization processes in various tectonic settings.

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

Using dynamic mode decomposition to predict the dynamics of a two-time non-equilibrium Green’s function

Computing the numerical solution of the Kadanoff–Baym equations, a set of nonlinear integral differential equations satisfied by the two-time Green's functions derived from many-body perturbation theory for a quantum many-body system away from equilibrium, is a challenging task. Recently, we have successfully applied dynamic mode decomposition (DMD) to construct a data driven reduced order model that can be used to extrapolate the time-diagonal of a two-time Green's function from numerical solutions of the KBE within a small time window. In this paper, we extend the previous work and use DMD to predict off-diagonal elements of the two-time Green's function. We partition the two-time Green's function into a number of one-time functions along the diagonal and subdiagonals of the two-time window as well as in horizontal and vertical directions. We use DMD to construct separate reduced order models to predict the dynamics of these one-time functions in a two-step procedure. We extrapolate along diagonal and several subdiagonals within a subdiagonal band of a two-time window in the first step. In the second step, we use DMD to extrapolate the Green's function outside of the sub-diagonal band. In conclusion, we demonstrate the efficiency and accuracy of this approach by applying it to a two-band Hubbard model problem.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

CFD-DEM solution verification: Fixed-bed studies

Fixed (frozen) particle configurations are used to study CFD-DEM convergence behavior in the limit of infinite grid resolution. A filtering scheme based on solving a diffusion equation, i.e., a Gaussian kernel, is used to interpolate particle information to the continuum grid, free of cell size limitations. A regression method is used for Richardson extrapolation of the discrete solutions to the grid-free solution. The regression method is shown to predict accurate extrapolation curves in the presence of the observed oscillatory convergence behavior. The extrapolation curve is used to complete the solution verification process by quantifying the numerical uncertainty present in the calculations. Several numerical and model parameters are studied including the spatial discretization method, the kernel width and the mean drag law. Here, a manageable level of discretization error is found with relative discretization errors generally <4%.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗