Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “interpolation models”

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 847 records · Page 47
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Improved Prediction of Cold-Air Pools in the Weather Research and Forecasting Model Using a Truly Horizontal Diffusion Scheme for Potential Temperature

The terrain-following vertical coordinate system used by many atmospheric models, including the Weather Research and Forecasting (WRF) Model, is prone to errors in regions of complex terrain. These errors stem, in part, from the calculation of horizontal gradients within the diffusion term of the momentum or scalar evolution equations. In WRF, such gradients can be calculated along coordinate surfaces, or using metric terms that help account for grid skewness. However, neither of these options ensures a truly horizontal gradient calculation, especially if a grid cell is skewed enough that the heights of the neighboring grid points used in the calculation fall outside the vertical range of the cell. In this work, an improved scheme that uses Taylor series approximations to vertically interpolate variables to the level necessary for a truly horizontal gradient calculation is implemented in WRF for the diffusion of potential temperature. The scheme is validated using an atmosphere-at-rest configuration, in which spurious flows develop only as a result of numerical errors and can thus be used as a proxy for model performance. Following validation, the method is applied to the simulation of cold-air pools (CAPs), which occur in regions of complex terrain and are characterized by strong near-surface temperature gradients. Using the truly horizontal scheme, idealized simulations demonstrate reduced numerical mixing in a quiescent CAP, and a realistic case study in the Columbia River basin shows a reduction in positive wind speed bias by up to roughly 20% compared to observations from the Second Wind Forecast Improvement Project.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2010

Regional Model Nesting Within GFS Daily Forecasts Over West Africa

The study uses the RM3, the regional climate model at the Center for Climate Systems Research of Columbia University and the NASA/Goddard Institute for Space Studies (CCSR/GISS). The paper evaluates 30 48-hour RM3 weather forecasts over West Africa during September 2006 made on a 0.5 grid nested within 1 Global Forecast System (GFS) global forecasts. September 2006 was the Special Observing Period #3 of the African Monsoon Multidisciplinary Analysis (AMMA). Archived GFS initial conditions and lateral boundary conditions for the simulations from the US National Weather Service, National Oceanographic and Atmospheric Administration were interpolated four times daily. Results for precipitation forecasts are validated against Tropical Rainfall Measurement Mission (TRMM) satellite estimates and data from the Famine Early Warning System (FEWS), which includes rain gauge measurements, and forecasts of circulation are compared to reanalysis 2. Performance statistics for the precipitation forecasts include bias, root-mean-square errors and spatial correlation coefficients. The nested regional model forecasts are compared to GFS forecasts to gauge whether nesting provides additional realistic information. They are also compared to RM3 simulations driven by reanalysis 2, representing high potential skill forecasts, to gauge the sensitivity of results to lateral boundary conditions. Nested RM3/GFS forecasts generate excessive moisture advection toward West Africa, which in turn causes prodigious amounts of model precipitation. This problem is corrected by empirical adjustments in the preparation of lateral boundary conditions and initial conditions. The resulting modified simulations improve on the GFS precipitation forecasts, achieving time-space correlations with TRMM of 0.77 on the first day and 0.63 on the second day. One realtime RM3/GFS precipitation forecast made at and posted by the African Centre of Meteorological Application for Development (ACMAD) in Niamey, Niger is shown.

Druyan, Leonard M.↗
NASA NTRS1996

Aspirating Seal Development: Analytical Modeling and Seal Test Rig

This effort is to develop large diameter (22 - 36 inch) Aspirating Seals for application in aircraft engines. Stein Seal Co. will be fabricating the 36-inch seal(s) for testing. GE's task is to establish a thorough understanding of the operation of Aspirating Seals through analytical modeling and full-scale testing. The two primary objectives of this project are to develop the analytical models of the aspirating seal system, to upgrade using GE's funds, GE's 50-inch seal test rig for testing the Aspirating Seal (back-to-back with a corresponding brush seal), test the aspirating seal(s) for seal closure, tracking and maneuver transients (tilt) at operating pressures and temperatures, and validate the analytical model. The objective of the analytical model development is to evaluate the transient and steady-state dynamic performance characteristics of the seal designed by Stein. The transient dynamic model uses a multi-body system approach: the Stator, Seal face and the rotor are treated as individual bodies with relative degrees of freedom. Initially, the thirty-six springs are represented as a single one trying to keep open the aspirating face. Stops (Contact elements) are provided between the stator and the seal (to compensate the preload in the fully-open position) and between the rotor face and Seal face (to detect rub). The secondary seal is considered as part of the stator. The film's load, damping and stiffness characteristics as functions of pressure and clearance are evaluated using a separate (NASA) code GFACE. Initially, a laminar flow theory is used. Special two-dimensional interpolation routines are written to establish exact film load and damping values at each integration time step. Additionally, other user-routines are written to read-in actual pressure, rpm, stator-growth and rotor growth data and, later, to transfer these as appropriate loads/motions in the system-dynamic model. The transient dynamic model evaluates the various motions, clearances and forces as the seals are subjected to different aircraft maneuvers: Windmilling restart; start-ground idle; ground idle-takeoff; takeoff-burst chop, etc. Results of this model show that the seal closes appropriately and does not ram into the rotor for all of the conditions analyzed. The rig upgrade design for testing Aspirating Seals has been completed. Long lead-time items (forgings, etc.) have been ordered.

Bharat Bagepalli↗
DOE OSTI2025

Data-driven high-dimensional statistical inference with generative models

Crucial to many measurements at the LHC is the use of correlated multi-dimensional information to distinguish rare processes from large backgrounds, which is complicated by the poor modeling of many of the crucial backgrounds in Monte Carlo simulations. In this work, we introduce HI-SIGMA, a method to perform unbinned high-dimensional statistical inference with data-driven background distributions. In contradistinction to many applications of Simulation Based Inference in High Energy Physics, HI-SIGMA relies on generative ML models, rather than classifiers, to learn the signal and background distributions in the high-dimensional space. These ML models allow for interpretable inference while also incorporating model errors and other sources of systematic uncertainties. We showcase this methodology on a simplified version of a di-Higgs measurement in the bbγγ final state, where the di-photon resonance allows for background interpolation from sidebands into the signal region. We demonstrate that HI-SIGMA provides improved sensitivity as compared to standard classifier-based methods, and that systematic uncertainties can be straightforwardly incorporated by extending methods which have been used for histogram based analyses.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

LandCast Mosaic: Reconstructing Global Population Distributions, 1975-2025

LandCast Mosaic (LCM) provides a global, high-resolution gridded population dataset spanning 1975–2025, representing annual, scenario-consistent estimates of daytime, nighttime, and ambient population distributions. LCM builds on the 2025 LandScan Mosaic (LSM) population data by backcasting to earlier years using historical changes in built-surface area derived from the Global Human Settlement Layer (GHSL) and authoritative population counts from international datasets. The workflow scales 2025 building-informed gridded population estimates according to observed changes in built surface, applies linear interpolation for intermediate years, and normalizes estimates to match administrative- and country-level totals. The resulting dataset offers consistent, globally gridded population estimates over fifty years, suitable for temporal analyses of population dynamics, disaster risk modeling, and urban planning applications.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS2017

The Flight Optimization System Weights Estimation Method

FLOPS has been the primary aircraft synthesis software used by the Aeronautics Systems Analysis Branch at NASA Langley Research Center. It was created for rapid conceptual aircraft design and advanced technology impact assessments. FLOPS is a single computer program that includes weights estimation, aerodynamics estimation, engine cycle analysis, propulsion data scaling and interpolation, detailed mission performance analysis, takeoff and landing performance analysis, noise footprint estimation, and cost analysis. It is well known as a baseline and common denominator for aircraft design studies. FLOPS is capable of calibrating a model to known aircraft data, making it useful for new aircraft and modifications to existing aircraft. The weight estimation method in FLOPS is known to be of high fidelity for conventional tube with wing aircraft and a substantial amount of effort went into its development. This report serves as a comprehensive documentation of the FLOPS weight estimation method. The development process is presented with the weight estimation process.

Wells, Douglas P.↗
DOE OSTI2024

Evaluation of Thermal Scattering Law for U 3 Si 2 Using Temperature-Dependent Phonon Density of States

The thermal scattering law for uranium silicide, a candidate material for high-density nuclear fuel, was evaluated from the temperature-dependent phonon density of states obtained using the quasi-harmonic approximation. For uranium silicide model systems of various volumes near the equilibrium volume, Hellmann-Feynman forces were calculated by performing ab initio lattice dynamics simulations based on density functional theory. Subsequently, the quasi-harmonic approximation was applied by interpolating the results of the harmonic approximation for each volume system. As a sensitivity test, the phonon density of states was calculated for the structures with temperature-dependent volumes at 298 K and 1000 K, obtained by the quasi-harmonic approximation, and the corresponding thermal neutron scattering cross section was evaluated. The temperature-dependent volume and heat capacity of uranium silicide obtained by the quasi-harmonic approximation showed reasonable agreement with experimental data while lattice thermal expansion showed deviation from experimental observations. Thermal scattering cross sections evaluated from 0 K and temperature-dependent phonon density of states exhibited differences of up to 4.1% at 1000 K. This study examines the impact of temperature on the thermal scattering law, employing ab initio lattice dynamics simulations through the quasi-harmonic approximation.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

Hugoniot properties of porous stainless steel: Insights from molecular dynamics simulations

Here, we present a systematic study of Hugoniot properties of porous 316L stainless steel using both a simple interpolation scheme and direct shock simulations in order to probe pore collapse kinetics as well as final thermodynamic states. Both methods indicate that equilibrated Hugoniot properties depend on pore density only and not on the pore distribution or size. We then create a simple porous equation of state model that is shown to be accurate for a range of validation data. This allows us to extend our simulations to make direct comparison to experimental data that have generally significantly larger system sizes and durations. In addition, our direct shock simulations indicate that the relaxation time after hotspot formation is system size dependent and can reach nanosecond timescales for the largest pores investigated in our study, thereby possibly having a measurable effect on fast dynamic loading experiments.

36 MATERIALS SCIENCE↗
NASA NTRS1989

Model-based quantification of image quality

In 1982, Park and Schowengerdt published an end-to-end analysis of a digital imaging system quantifying three principal degradation components: (1) image blur - blurring caused by the acquisition system, (2) aliasing - caused by insufficient sampling, and (3) reconstruction blur - blurring caused by the imperfect interpolative reconstruction. This analysis, which measures degradation as the square of the radiometric error, includes the sample-scene phase as an explicit random parameter and characterizes the image degradation caused by imperfect acquisition and reconstruction together with the effects of undersampling and random sample-scene phases. In a recent paper Mitchell and Netravelli displayed the visual effects of the above mentioned degradations and presented subjective analysis about their relative importance in determining image quality. The primary aim of the research is to use the analysis of Park and Schowengerdt to correlate their mathematical criteria for measuring image degradations with subjective visual criteria. Insight gained from this research can be exploited in the end-to-end design of optical systems, so that system parameters (transfer functions of the acquisition and display systems) can be designed relative to each other, to obtain the best possible results using quantitative measurements.

Hazra, Rajeeb↗
DOE OSTI2020

Development and Validation of Passive Yaw in the Open-Source WEC-Sim Code

A passive yaw implementation is developed, validated, and explored for the WEC-Sim, an open-source wave energy converter modeling tool that works within MATLAB/Simulink. The Reference Model 5 (RM5) is selected for this investigation, and a WEC-Sim model of the device is modified to allow yaw motion. A boundary element method (BEM) code was used to calculate the excitation force coefficients for a range of wave headings. An algorithm was implemented in WEC-Sim to determine the equivalent wave heading from a body’s instantaneous yaw angle and interpolate the appropriate excitation coefficients to ensure the correct time-domain excitation force. This approach is able to determine excitation force for a body undergoing large yaw displacement. For the mathematically simple case of regular wave excitation, the dynamic equation was integrated numerically and found to closely approximate the results from this implementation in WEC-Sim. A case study is presented for the same device in irregular waves. In this case, computation time is increased by 32x when this interpolation is performed at every time step. To reduce this expense, a threshold yaw displacement can be set to reduce the number of interpolations performed. A threshold of 0.01° was found to increase computation time by only 22x without significantly affecting time domain results. Similar amplitude spectra for yaw force and displacements are observed for all threshold values less than 1°, for which computation time is only increased by 2.2x.

50 EE - Wind and Water Power Program - Water (EE-4↗
NASA NTRS2019

Scaling Laws for Impact Craters in Water

Tsunami waves from asteroid impacts into water are of concern from asteroids in the 200 m to 1 km diameter range because this spans the range from asteroids that will likely hit the surface and not airburst, but also be small enough that global climate effects will hopefully be minimal. Current estimates of impact tsunamis depend on either hydrocode simulations or on semi-analytical models. Unfortunately there is significant disagreement between these methods. One of the main reasons for the disparity is that the semi-analytical models such as Holsapple (1993) rely on experimental impacts into deep water. However for asteroids in the 200 - 1000 m range even the deep ocean basins can appear as shallow water impacts where the crater formed in the water reaches all the way to the sea floor. Another reason for the disparity arises from the linear interpolation of data across many orders of magnitude difference between Froude number (ratio of kinetic to gravitational energy) used in the laboratory experiments and what would be seen in an asteroid impact. The Gault & Sonett (1982) experiment shot millimeter sized glass spheres into water at 1 to 6 km/s and the Olevson (1969) experiment dropped millimeter sized water drops at a few meters per second. The goal of this work was to fill in the gaps and conduct experiments and simulations at the correct Froude numbers of interest, and in both deep and shallow water, to help resolve the disparity and extend the semi-analytical models.

Robertson, Darrel↗
DOE OSTI2026

Integrated Neutronics Modeling for Inertial Fusion Energy Systems: Development and Application to LD-FIRST

Lawrence Livermore National Laboratory (LLNL) is proposing a new Laser Driven Fusion Integration Research and Science Test Facility (LD-FIRST) with the goal of providing an experimental testbed for future Inertial Fusion Energy (IFE) systems. However, IFE systems require detailed and accurate multiphysics modeling to quantify material damage, thermal loading, and tritium breeding within complex chamber environments. This article presents the first step in an integrated multiphysics framework that couples meshed CAD-based geometry within Monte Carlo neutronic simulations to enable high-fidelity analysis of IFE chamber concepts, with future coupling to external codes. The neutronics workflow utilizes OpenMC and its third-party capability to use CAD-based geometries through DAGMC and tally on unstructured meshes with Libmesh to evaluate neutron transport behavior, geometric fidelity, and material performance under reactor-relevant conditions. The use of tailored tallies on unstructured meshes in this framework allows direct transfer without interpolating to CFD simulation tools. Two IFE chambers were evaluated, both conceived by LLNL: HYLIFE-II and Laser IFE (LIFE). This work produced high-fidelity conformal surface and volumetric meshes of the HYLIFE-II and LIFE chambers with mapped spatial insight into material damage, thermal loading, and tritium breeding. The HYLIFE-II model was built utilizing available resources and used as a test case to verify that the neutronics framework can handle complex geometries. The LIFE chamber CAD was provided by LLNL and was the main focus of this work. This work analyzes multiple ternary alloy breeding materials for the LIFE chamber, across different 6 Li enrichments to produce data relevant to the LD-FIRST project. This work also investigates the level of model fidelity for the LIFE chamber, and results show that inclusion of detailed first wall and coolant structures increased the predicted tritium breeding ratio (TBR) by ~30%, highlighting the sensitivity of tritium breeding and the need for a high-fidelity simulation framework for IFE chambers. These developments provide a scalable toolset for the design and optimization of next-generation IFE chambers, forming a solid foundation for future coupled multiphysics analysis.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Simulation of adiabatic quantum computing for molecular ground states

Quantum computation promises to provide substantial speedups in many practical applications with a particularly exciting one being the simulation of quantum many-body systems. Adiabatic state preparation (ASP) is one way that quantum computers could recreate and simulate the ground state of a physical system. In this paper, we explore a novel approach for classically simulating the time dynamics of ASP with high accuracy and with only modest computational resources via an adaptive sampling configuration interaction scheme for truncating the Hilbert space to only the most important determinants. We verify that this truncation introduces negligible error and use this new approach to simulate ASP for sets of small molecular systems and Hubbard models. Furthermore, we examine two approaches to speeding up ASP when performed on quantum hardware: (i) using the complete active space configuration interaction (CASCI) wave function instead of the Hartree–Fock initial state and (ii) a nonlinear interpolation between the initial and target Hamiltonians. We find that starting with a CASCI wave function with a limited active space yields substantial speedups for many of the systems examined, while nonlinear interpolation does not. In additional, we observe interesting trends in the minimum gap location (based on the initial state) as well as how state preparation time can depend on certain molecular properties, such as the number of valence electrons. Importantly, we find that the required state preparation times do not show an immediate exponential wall that would preclude an efficient run of ASP on actual hardware.

Kremenetski, Vladimir↗
DOE OSTI2020

Disruption Prediction and Avoidance in High Beta Long Pulse KSTAR Plasmas

DoE Grant DE-SC0016409 (Disruption Prediction and Avoidance in High Beta Long Pulse KSTAR Plasmas) supported a collaboration among Columbia University, the Princeton Plasma Physics Laboratory (PPPL), and the Plasma Science and Fusion Center (PSFC) at MIT to conduct research on the KSTAR tokamak at the National Fusion Research Institute (NFRI) in Korea to improve predictive capability for disruptions in tokamaks. This initiative included upgrades to kinetic diagnostics at KSTAR whose measurements form the basis for the disruption prediction models. The major thrust of the PPPL and PSFC work was implementing a capability at KSTAR to allow partially-polarized background light seen by the KSTAR Motional Stark Effect (MSE) diagnostic to be compensated by ‘wavelength interpolation’ rather than the traditional ‘beam modulation’ technique. This objective was completed on time and on budget, and the MSE Background Polychromator is currently functioning at KSTAR. Primary responsibility for maintaining and operating the MSE-BP was transferred to Nova Photonics, Inc. in 2019.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1985

Experiments on the Assimilation of Patches of Data Using a Barotropic Model

Experiments were performed to examine the impact of adopting various procedures in the assimilation of patches of data during the progress of a numerical forecast. A barotropic model described by Takacs and Balgovind, 1983 is used for the experiments. The model is global and its discretization is similar to the GLA Fourth Order GCM. The effect of orography is included in the model. The experiments are of an identical twin kind. The nature run starts from global fields of wind (u and v) and height (h) at 300 mb. They are real data of January, 1979 OOZ which were interpolated from the adjacent sigma levels of the GLA GCM and then balanced by the bounded derivative method. From this initial state a 6 day forecast is performed. At intervals of 3 hours this information is saved over sectors of 60 width in longitude, from the south pole to the north pole, to simulate data coverage and observational frequency of a polar orbiting satellite.

Semazzi, F. H. M.↗
NASA NTRS1988

Nonlinear modeling and estimation of slew induced structural deformations

A model of the nonlinear dynamics of a deformable maneuvering multibody system is described, whereby elastic deformation are modeled by restoring forces and dissipative forces at point mass appendages. This model is brought into bilinear form. Estimation of deformations occasioned by rapid slewing maneuvers is carried out by a filter based on a globally equivalent linear model of the bilinear dynamics, and is shown to be an improvement over the extended Kalman filter. To further alleviate the computational burden, the estimated deformation state is propagated between observations by a low dimensional operator spline interpolator of bilinear system Volterra series, which is easily implemented.

Dwyer, T. A. W., III↗
NASA NTRS1991

Bilinear modeling and state estimation of slew-induced deformations of a flexible multibody system

A model of the nonlinear dynamics of a deformable maneuvering multibody system is described, whereby elastic deformations are modeled by restoring forces and dissipative forces at point mass appendages. This model is brought into bilinear form. Estimation of deformations occasioned by rapid slewing maneuvers is carried out by a filter based on a globally equivalent linear model of the bilinear dynamics, and is shown to be an improvement over the extended Kalman filter. To further alleviate the computational burden, the estimated deformation state is propagated between observations by a low dimensional operator spline interpolator of bilinear system Volterra series, which is easily implemented.

Dwyer, Thomas A. W., III↗
NASA NTRS1994

Computation of helicopter rotor acoustics in forward flight

This paper presents a new method for computing acoustic signals from helicopter rotors in forward flight. The aerodynamic and acoustic solutions in the near field are computed with a finite-difference solver for the Euler equations. A nonrotating cylindrical Kirchhoff surface is then placed around the entire rotor system. This Kirchhoff surface moves subsonically with the rotor in forward flight. The finite-difference solution is interpolated onto this cylindrical surface at each time step and a Kirchhoff integration is used to carry the acoustic signal to the far field. Computed values for high-speed impulsive noise show excellent agreement with model-rotor and flight-test experimental data. Results from the new method offer high accuracy with reasonable computer resource requirements.

Strawn, Roger↗