Search NASA⌕ Search

SEARCH · Search NASA

Results for “High dimensional modeling”

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
DOE OSTI2026

A scalable variational method for estimating the latent infection-rate field of an outbreak

In this paper, we explore whether the infection-rate of a disease can serve as a robust monitoring variable in epidemiological surveillance algorithms. The infection-rate is dependent on population mixing patterns that do not vary erratically day-to-day; in contrast, daily case-counts used in contemporary surveillance algorithms are corrupted by reporting errors. The technical challenge lies in estimating the latent infection-rate from case-counts. Here we devise a Bayesian method to estimate the infection-rate across multiple adjoining areal units, and then use it, via an anomaly detector, to discern a change in epidemiological dynamics. We extend an existing model for estimating the infection-rate in an areal unit by incorporating a Markov random field model, so that we may estimate infection-rates across multiple areal units, while preserving spatial correlations observed in the epidemiological dynamics. To carry out the high-dimensional Bayesian inverse problem, we develop an implementation of mean-field variational inference specific to the infection model and integrate it with the random field model to incorporate correlations across counties. The method is tested on estimating the COVID-19 infection-rates across all 33 counties in New Mexico using data from the summer of 2020, and then employing them to detect the arrival of the Fall 2020 COVID-19 wave. We perform the detection using a temporal algorithm that is applied county-by-county. We also show how the infection-rate field can be used to cluster counties with similar epidemiological dynamics.

60 APPLIED LIFE SCIENCES↗
NASA NTRS1982

A high-resolution model of the planetary boundary layer - Sensitivity tests and comparisons with SESAME-79 data

A one-dimensional, planetary boundary layer (PBL) model is presented and verified using April 10, 1979 SESAME data. The model contains two modules to account for two different regimes of turbulent mixing. Separate parameterizations are made for stable and unstable conditions, with a predictive slab model for surface temperature. Atmospheric variables in the surface layer are calculated with a prognostic model, with moisture included in the coupled surface/PBL modeling. Sensitivity tests are performed for factors such as moisture availability, albedo, surface roughness, and thermal capacity, and a 24 hr simulation is summarized for day and night conditions. The comparison with the SESAME data comprises three hour intervals, using a time-dependent geostrophic wind. Close correlations were found with daytime conditions, but not in nighttime thermal structure, while the turbulence was faithfully predicted. Both geostrophic flow and surface characteristics were shown to have significant effects on the model predictions

Zhang, D.↗
NASA NTRS1986

Polar warming in the middle atmosphere of Mars

During the 1984 Mars opposition, ground-based laser heterodyne spectroscopy was obtained for the nonthermal core emission of the 10.33-micron R(8) and 10.72-micron P(32) lines of C-12(O-16)2 at 23 locations on the Martian disk. It is deduced on the basis of these data that the temperature of the middle Martian atmosphere varies with latitude, and a meridional gradient of 0.4-0.9 K/deg latitude is indicated. The highest temperatures are noted to lie at high latitudes in the winter hemisphere; as in the terrestrial case of seasonal effects at the menopause, this winter polar warming in the Martian middle atmosphere requires departures from radiative equilibrium. Two-dimensional circulation model comparisons with these results indicate that atmospheric dust may enhance this dynamical heating at high winter latitudes.

Deming, D.↗
DOE OSTI2025

Evaluating Chemical Kinetics Predictions for Propane Using 3-D and 0-D Models in a Boosted Spark-Ignited Engine

Propane has been shown to be a promising alternative fuel to reduce emissions while simultaneously achieving high efficiencies in medium- and heavy-duty engines. These high-power density applications require boosted engines which, combined with high compression ratio, can lead to auto-ignition and knock. While three-dimensional (3-D) computational fluid dynamics (CFD) models are often used for resolving the complex fluid flow in engines, these models can become computationally expensive when simulating detailed chemical kinetics. Likewise, zero-dimensional (0-D) models are computationally concise enough for kinetics development, but lack any flow-field information which governs the flame propagation processes in spark ignition (SI) engines. This work presents a comprehensive comparison between 3-D and 0-D closed cycle simulations at knocking conditions in a high compression ratio high stroke-to-bore ratio propane engine. In order to initialize the flow-field for the 3-D closed cycle (intake valve closing, (IVC) to exhaust valve opening, (EVO)) simulation, a motored multi-cycle 3-D model was run using Converge to create a map at IVC, reducing the computational time. The map allowed a non-homogeneous 3-D closed cycle simulation to be satisfactorily validated against experiments, while a homogeneous case using only the turbulence field mapping was also simulated, mimicking 0-D modeling. The 3-D simulations were used to prescribe the initial conditions (e.g., IVC thermodynamics, speciation, burn-rate profile) for a 2-zone 0-D SI engine model in Chemkin Pro for both cases. It was found that 2-zone 0-D modeling underpredicted the knock onset timing, likely due to the lack of thermal stratification in the unburned gas region. Future work will carry multi-zone 0-D modeling to capture the fuel auto-ignition in the unburned region.

Douvry-Rabjeau, Julien [Oakland University, Roches↗
DOE OSTI2025

Hyperplane decision trees as piecewise linear surrogate models for chemical process design

Recent trends in chemical engineering research point towards an increasing reliance on data-driven modeling approaches. Neural networks, for instance, have proven to be accurate when data is plentiful and high-dimensional, but in many cases, they require computationally-intensive training procedures. Here, in this work, we describe hyperplane decision trees (HT) as a highly expressive and low-compute machine learning model architecture. These models are locally linear and have linear decision boundaries, resulting in a piecewise linear model of the data. This property allows them to be converted into mixed-integer linear constraints which can be globally optimized. Our open-source PyTorch implementation of this method is a fast, flexible, and accessible way to build accurate piecewise linear models of data.

Decision trees↗
DOE OSTI2025

Improving the Quasi‐Biennial Oscillation via a Surrogate‐Accelerated Multi‐Objective Optimization

Accurate simulation of the quasi-biennial oscillation (QBO) is challenging due to uncertainties in representing convectively generated gravity waves. We develop an end-to-end uncertainty quantification workflow that calibrates these gravity wave processes in E3SM for a realistic QBO. Central to our approach is a domain knowledge-informed, compressed representation of high-dimensional spatio-temporal wind fields. By employing a parsimonious statistical model that learns the fundamental frequency from complex observations, we extract interpretable and physically meaningful quantities capturing key attributes. Building on this, we train a probabilistic surrogate model that approximates the fundamental characteristics of the QBO as functions of critical physics parameters governing gravity wave generation. Leveraging the Karhunen–Loève decomposition, our surrogate efficiently represents these characteristics as a set of orthogonal features, capturing cross-correlations among multiple physics quantities evaluated at different pressure levels and enabling rapid surrogate-based inference at a fraction of the computational cost of full-scale simulations. Finally, we analyze the inverse problem using a multi-objective approach. Our study reveals a tension between amplitude and period that constrains the QBO representation, precluding a single optimal solution. To navigate this, we quantify the bi-criteria trade-off and generate a set of Pareto optimal parameter values that balance the conflicting objectives. This integrated workflow improves the fidelity of QBO simulations and offers a versatile template for uncertainty quantification in complex geophysical models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Neural refinement of sample weights

Monte Carlo simulations are an essential tool in particle physics data analysis. Events are typically generated alongside weights that redistribute the cross section of the simulated process across the phase space. These weights can be negative, and several post hoc methods have been developed to eliminate or mitigate the negative values. All of these methods share the common strategy of approximating the average weight as a function of phase space. We introduce an alternative approach, which, instead of reweighting to the average, refines the initial weights with a scaling transformation, utilizing a phase space-dependent factor. Since this new refinement method does not need to model the full weight distribution, it can be more accurate. High-dimensional and unbinned phase space is processed using neural networks for the refinement method. In addition to the refinement method, we introduce a new resampling protocol, which can be used in conjunction with any weight transformation to not only preserve the average weight but also the statistical uncertainties of the initial distribution. Using both realistic and synthetic examples, we show that the new neural refinement method is able to match or exceed the accuracy of similar weight transformations and that the new resampling protocol is simpler in implementation than previous methods while exhibiting equivalent statistical properties.

Artificial neural networks↗
NASA NTRS1989

The transient behavior of weldpools with a deformed free surface

Transient heat transfer, fluid flow and phase change in a weldpool with a deformed free surface are modeled using a two-dimensional coordinate transfornmation. The theoretical approach makes it possible to accurately model high Marangoni number flows even when the free surface is not horizontal. The method is used to examine the role of surface deformation and surface tension forces in controlling weld penetration during cathode spot welding of steel and aluminum. The calculated flow fields, temperature distributions and weldpool shapes are consistent with experimental findings, and are in reasonable agreement with previously published modeling efforts.

Thompson, M. E.↗
NASA NTRS1990

Observed effects of landscape variability on convective clouds

A case is presented in which clouds are observed to form first over a mesoscale-size area (100 x 300 km) of harvested wheat in Oklahoma, where the ground temperature is warmer than adjoining areas dominated by growing vegetation. In addition, clouds are suppressed over relatively long bands downwind of small man-made lakes and areas characterized by heavy tree cover. The observed variability of cloud relative to landscape type is compared with that simulated with a one-dimensional boundary-layer model. Clouds form earliest over regions characterized by high, sensible heat flux, and are suppressed over regions characterized by high, latent heat flux during relatively dry atmospheric conditions. This observation has significance in gaining understanding of the feedback mechanisms of land modification on climate, as well as understanding relatively short-range weather forecasting.

Rabin, Robert M.↗
NASA NTRS1995

One-dimensional hybrid satellite track model for the Dynamics Explorer 2 (DE 2) satellite

A one-dimensional hybrid satellite track model has been developed to calculate the high-latitude thermospheric/ionospheric structure below the satellite altitude using Dynamics Explorer 2 (DE 2) satellite measurements and theory. This model is based on Emery et al. (1985) satellite track code but also includes elements of Roble et al. (1987b) global mean thermosphere/ionosphere model. A number of parameterizations and data handling techniques are used to input satellite data from several DE 2 instruments into this model. Profiles of neutral atmospheric densities are determined from the Mass Spectrometer Incoherent Scatter 1990 (MSIS-90) model and measured neutral temperatures. Measured electron precipitation spectra are used in an auroral model to calculate particle impact ionization rates below the satellite. These rates are combined with a solar ionization rate profile and used to solve the O(+) diffusion equation, with the measured electron density as an upper boundary condition. The calculated O(+) density distribution, as well as the ionization profiles, are then used in a photochemical equilibrium model to calculate the electron and molecular ion densities. The electron temperature is also calculated by solving the electron energy equation with an upper boundary condition determined by the DE 2 measurement. The model enables calculations of altitude profiles of conductivity and Joule heating rate along and below the satellite track. In a first application of the new model, a study is made of thermospheric and ionospheric structure below the DE 2 satellite for a single orbit which occurred on October 25, 1981. The field-aligned Poynting flux, which is independently obtained for this orbit, is compared with the model predictions of the height-integrated energy conversion rate. Good quantitative agreement between these two estimates has been reached. In addition, measurements taken at the incoherent scatter radar site at Chatanika (65.1 deg N, 147.4 deg W) during a DE 2 overflight are compared with the model calculations. A good agreement was found in lower thermospheric conductivities and Joule heating rate.

Deng, Wei↗
NASA NTRS1996

Endwall Heat Transfer Measurements in a Transonic Turbine Cascade

Turbine blade endwall heat transfer measurements are given for a range of Reynolds and Mach numbers. Data were obtained for Reynolds numbers based on inlet conditions of 0.5 and 1.0 x 106, for isentropic exit Mach numbers of 1.0 and 1.3, and for freestream turbulence intensities of 0.25% and 7.0%. Tests were conducted in a linear cascade at the NASA Lewis Transonic Turbine Blade Cascade Facility. The test article was a turbine rotor with 136' of turning and an axial chord of 12.7 cm. The large scale allowed for very detailed measurements of both flow field and surface phenomena. The intent of the work is to provide benchmark quality data for computational fluid dynamics (CFD) code and model verification. The flow field in the cascade is highly three-dimensional as a result of thick boundary layers at the test section inlet. Endwall heat transfer data were obtained using a steady-state liquid crystal technique.

Giel, P. W.↗
NASA NTRS1995

One-Dimensional Hybrid Satellite Track Model for the Dynamics Explorer 2 (DE 2) Satellite

A one-dimensional hybrid satellite track model has been developed to calculate the high-latitude thermospheric/ionospheric structure below the satellite altitude using Dynamics Explorer 2 (DE 2) satellite measurements and theory. This model is based on Emery et al. satellite track code but also includes elements of Roble et al. global mean thermosphere/ionosphere model. A number of parameterizations and data handling techniques are used to input satellite data from several DE 2 instruments into this model. Profiles of neutral atmospheric densities are determined from the MSIS-90 model and measured neutral temperatures. Measured electron precipitation spectra are used in an auroral model to calculate particle impact ionization rates below the satellite. These rates are combined with a solar ionization rate profile and used to solve the O(+) diffusion equation, with the measured electron density as an upper boundary condition. The calculated O(+) density distribution, as well as the ionization profiles, are then used in a photochemical equilibrium model to calculate the electron and molecular ion densities. The electron temperature is also calculated by solving the electron energy equation with an upper boundary condition determined by the DE 2 measurement. The model enables calculations of altitude profiles of conductivity and Joule beating rate along and below the satellite track. In a first application of the new model, a study is made of thermospheric and ionospheric structure below the DE 2 satellite for a single orbit which occurred on October 25, 1981. The field-aligned Poynting flux, which is independently obtained for this orbit, is compared with the model predictions of the height-integrated energy conversion rate. Good quantitative agreement between these two estimates has been reached. In addition, measurements taken at the incoherent scatter radar site at Chatanika (65.1 deg N, 147.4 deg W) during a DE 2 overflight are compared with the model calculations. A good agreement was found in lower thermospheric conductivities and Joule heating rate.

Deng, Wei↗
NASA NTRS2000

Monitoring Global Sea Level: Eustatic Variations, Local Variations, and Solid Earth Effects

Project BIFROST (Baseline Inferences for Fennoscandian Rebound Observations. Sea-level and Tectonics) combines networks of continuously operating CPS receivers in Sweden and Finland to measure ongoing crustal deformation due to glacial isostatic adjustment, (CIA). We present an analysis of data collected in the years 1993-1998. We compare the CPS determinations of three-dimensional crustal motion to predictions calculated using the high resolution Fennoscandian deglaciation model recently proposed by Lambeck et al. We find that the the maximum observed uplift rate (approx. 10 mm/ yr) and the maximum predicted uplift rate agree to better than 1 mm/ yr. The patterns of uplift also agree quite well, although differences are discernible. The chi(exp 2) difference between predicted and GPS-observed radial rates is reduced by a factor of 5-6 compared to that for the "null" (no uplift) model, depending on whether a mean difference is first removed. The north components of velocity agree at about the same relative level. whereas the agreement for the east components is worse, a problem possibly related to the lack of bias fixing. We have also compared the values for the observed radial deformation rates to those based on sea-level rates from Baltic tide gauges. The weighted RMS difference between CPS and tide-gauge rates (after removing a mean) is 0.6 mm/ yr, giving an indication of the combined accuracy of the CPS and tide-gauge measurement systems. Spectral analysis of the time series of position estimates yields spectral indices in the range -1 to -2. An EOF analysis indicates, however, that much of this power is correlated among the sites. The correlation appears to be regional and falls off only slightly with distance. Some of this correlated noise is associated with snow accumulation on the antennas or, for those antennas with radomes, on the radomes. This problem has caused us to modify the radomes used several times, leading to one of our more significant sources of uncertainty.

Davis, James L.↗
NASA NTRS2004

Evaluation of GOME Satellite Measurements of Tropospheric NO2 and HCHO using Regional Data from Aircraft Campaigns in the Southeastern United States

We compare tropospheric measurements of nitrogen dioxide (NO2) and formaldehyde (HCHO) from the Global Ozone Monitoring Experiment (GOME) satellite instrument with in situ measurements over eastern Texas and the southeast United States. On aveiage, the GOME and in situ measurements of tropospheric NO2 and HCHO columns are consistent despite pronounced sampling differences. The geometric mean in situ to GOME ratios over the campaign are 1.08 for NO2 and 0.84 for HCHO, with corresponding geometric standard deviations of 1.27 and 1.38. The correlation of the observed column spatial variability between the two NO2 measurement sets is encouraging before (r2 = 0.54, n = 18) and after (r2 = 0.67, n = 18) correcting for a sampling bias. Mean relative vertical profiles of HCHO and NO2 calculated with a global three-dimensional model (GEOS-CHEM) and used in the GOME retrieval are highly consistent with in situ measurements; differences would affect the retrieved NO2 and HCHO columns by a few percent. GOME HCHO columns over eastern Texas include contributions from anthropogenic volatile organic compound (VOC) emissions but are dominated by biogenic VOC emissions at the regional scale in August-September when HCHO columns are within 20% of those over the southeastern United States. In situ measurements show that during summer the lowest 1500 m (the lower mixed layer) contains 75% of the tropospheric NO2 column over Houston and Nashville, and 60% of the HCHO column over Houston. Future validation of space-based measurements of tropospheric NO2 and HCHO columns over polluted regions should include coincident in situ measurements that span the entire satellite footprint, especially in the heterogeneous mixed layer.

Martin, R. V.↗
NASA NTRSDate not supplied

Statistical learning framework for safety and failure analysis of a DNN-based autonomous aircraft system

Deep Neural Networks (DNNs) and Machine Learning technology is increasingly used for safety-critical applications in the Aerospace domain. To ensure safe operations, the DNN and the system must undergo rigorous verification and validation, including advanced statistical analyses. Performance and safety of the DNN and system behavior must not only be analyzed for the nominal case, but under numerous off-nominal and failure cases. In this paper we will describe how our statistical learning framework SYSAI can efficiently perform such analyses using the tool’s unique combination of advanced learning modeling and statistical analysis techniques. SYSAI can effectively explore the high-dimensional state and failure space of the system under test; geometrical shape detection of safety regions and boundaries support explainability of the results to the designer. In this paper, we report experiments and results obtained with a vision-based DNN control system (ACT) that is capable of autonomously steering an aircraft down a runway.

Yuning He↗
NASA NTRSDate not supplied

Instabilities in the Wake of a Pseudorandom Roughness on a Supersonic Flat Plate

A roughness patch defined by a two-dimensional Fourier series was installed in a highly-polished flat plate model immersed in a supersonic quiet flow. The NASA Langley Research Center Supersonic Low Disturbance Tunnel is capable of providing a Mach 3.5 flow with low freestream acoustic noise and low turbulence levels. Measurements made with a hot-wire probe downstream of the roughness element showed the presence of instability modes that did not cause transition in the measurement region at a freestream unit Reynolds number of Re =12.6×10^6/m. The largest amplitude mode appeared to be similar to a symmetric mode. The same roughness patch immersed in a higher freestream acoustic noise environment showed little to no evidence of the instability modes. Instead, disturbances found in the boundary layer matched closely to the freestream spectra measured in the empty tunnel under the same conditions, indicating some evidence of freestream forcing in the boundary layer.

Roughness induced transition↗
NASA NTRS2022

Instabilities in the Wake of a Pseudorandom Roughness on a Supersonic Flat Plate

A roughness patch defined by a two-dimensional Fourier series was installed in a highly-polished flat plate model immersed in a supersonic quiet flow. The NASA Langley Research Center Supersonic Low Disturbance Tunnel is capable of providing a Mach 3.5 flow with low freestream acoustic noise and low turbulence levels. Measurements made with a hot-wire probe downstream of the roughness element showed the presence of instability modes that did not cause transition in the measurement region at a freestream unit Reynolds number of Re =12.6×10^6/m. The largest amplitude mode appeared to be similar to a symmetric mode. The same roughness patch immersed in a higher freestream acoustic noise environment showed little to no evidence of the instability modes. Instead, disturbances found in the boundary layer matched closely to the freestream spectra measured in the empty tunnel under the same conditions, indicating some evidence of freestream forcing in the boundary layer.

Roughness induced transition↗
NASA NTRS2022

Instabilities in the Wake of a Pseudorandom Roughness on a Supersonic Flat Plate

A roughness patch defined by a two-dimensional Fourier series was installed in a highly-polished flat plate model immersed in a supersonic quiet flow. The NASA Langley Research Center Supersonic Low Disturbance Tunnel is capable of providing a Mach 3.5 flow with low freestream acoustic noise and low turbulence levels. Measurements made with a hot-wire probe downstream of the roughness element showed the presence of instability modes that did not cause transition in the measurement region at a freestream unit Reynolds number of Re =12.6×10^6/m. The largest amplitude mode appeared to be similar to a symmetric mode. The same roughness patch immersed in a higher freestream acoustic noise environment showed little to no evidence of the instability modes. Instead, disturbances found in the boundary layer matched closely to the freestream spectra measured in the empty tunnel under the same conditions, indicating some evidence of freestream forcing in the boundary layer.

Roughness induced transition↗