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.

SourceResearch connectionFrom yearThrough year

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

At least 1,117 records · Page 62

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Spinning black hole scattering at $ \mathcal{O} $(G 3 S 2 )(G 3 S 2 ): Casimir terms, radial action and hidden symmetry

We resolve subtleties in calculating the post-Minksowskian dynamics of binary systems, as a spin expansion, from massive scattering amplitudes of fixed finite spin. In particular, the apparently ambiguous spin Casimir terms can be fully determined from the gradient of the spin-diagonal part of the amplitudes with respect to S 2 = −s(s+1)ħ 2 , using an interpolation between massive amplitudes with different spin representations. From two-loop amplitudes of spin-0 and spin-1 particles minimally coupled to gravity, we extract the spin Casimir terms in the conservative scattering angle between a spinless and a spinning black hole at $ \mathcal{O} $(G 3 S 2 ), finding agreement with known results in the literature. This completes an earlier study [Phys. Rev. Lett. 130 (2023), 021601] that calculated the non-Casimir terms from amplitudes. We also illustrate our methods using a model of spinning bodies in electrodynamics, finding agreement between scattering amplitude predictions and classical predictions in a root-Kerr electromagnetic background up to $ \mathcal{O} $(α 3 S 2 ). For both gravity and electrodynamics, the finite part of the amplitude coincides with the two-body radial action in the aligned spin limit, generalizing the amplitude-action relation beyond the spinless case. Surprisingly, the two-loop amplitude displays a hidden spin-shift symmetry in the probe limit, which was previously observed at one loop. We conjecture that the symmetry holds to all orders in the coupling constant and is a consequence of integrability of Kerr orbits in the probe limit at the first few orders in spin.

Classical Theories of Gravity↗

DOE OSTI2022

Learning two-phase microstructure evolution using neural operators and autoencoder architectures

Abstract Phase-field modeling is an effective but computationally expensive method for capturing the mesoscale morphological and microstructure evolution in materials. Hence, fast and generalizable surrogate models are needed to alleviate the cost of computationally taxing processes such as in optimization and design of materials. The intrinsic discontinuous nature of the physical phenomena incurred by the presence of sharp phase boundaries makes the training of the surrogate model cumbersome. We develop a framework that integrates a convolutional autoencoder architecture with a deep neural operator (DeepONet) to learn the dynamic evolution of a two-phase mixture and accelerate time-to-solution in predicting the microstructure evolution. We utilize the convolutional autoencoder to provide a compact representation of the microstructure data in a low-dimensional latent space. After DeepONet is trained in the latent space, it can be used to replace the high-fidelity phase-field numerical solver in interpolation tasks or to accelerate the numerical solver in extrapolation tasks.

36 MATERIALS SCIENCE↗

NASA NTRS2018

Surface Reconstruction of a Challenging Region with Permanent Shadows on the Moon

As the Earth is orbiting the Sun, there is the region a sunlight can't be reached on the Moon, a satellite of the Earth. Usually, the permanent shadows are located near a pole area and the surface reconstruction in the permanent shadows region is challenging because of the lack of sunlight.In this research, the reconstruction method containing the permanent shadows is suggested. To apply this method, the time-varying shadow region and a permanently shadowed region are needed to be separated. In the time-varying shadow region, the pixel intensity of the image can be composed of the albedo and the reflectance model. Using the light source direction from ISIS3 of USGS, the shape of the region can be extracted from the reflectance model. Though this region is visible in the camera, the change of illumination is quite large in this region, it is hard to reconstruct the shape using stereo-photogrammetry or SfS. So, it is needed of collecting more surface normal information from several pixel level-aligned images. In this process, the albedo information can be initialized and enhanced easily from many images.In the middle of permanent shadows, the high-resolution DEM is generated by LOLA DEM. The LOLA DEM is generated by interpolation between measured data point strips, the actual scene of the permanent shadows. Using the sensor like DIVINER, the Imaginary scene generation via a trained network between the sunlight sensor and other kinds of a sensor is suggested.

Moon, Sunghyun↗

NASA NTRS2019

Soil Moisture Data Assimilation

Accurate knowledge of soil moisture at the continental scale is important for improving predictions of weather, agricultural productivity and natural hazards, but observations of soil moisture at such scales are limited to indirect measurements, either obtained through satellite remote sensing or from meteorological networks. Land surface models simulate soil moisture processes, using observation-based meteorological forcing data, and auxiliary information about soil, terrain and vegetation characteristics. Enhanced estimates of soil moisture and other land surface variables, along with their uncertainty, can be obtained by assimilating observations of soil moisture into land surface models. These assimilation results are of direct relevance for the initialization of hydro-meteorological ensemble forecasting systems. The success of the assimilation depends on the choice of the assimilation technique, the nature of the model and the assimilated observations, and, most importantly, the characterization of model and observation error. Systematic differences between satellite-based microwave observations or satellite-retrieved soil moisture and their simulated counterparts require special attention. Other challenges include inferring root-zone soil moisture information from observations that pertain to a shallow surface soil layer, propagating information to unobserved areas and downscaling of coarse information to finer-scale soil moisture estimates. This chapter summarizes state-of-the-art solutions to these issues with conceptual data assimilation examples, using techniques ranging from simplified optimal interpolation to spatial ensemble Kalman filtering. In addition, operational soil moisture assimilation systems are discussed that support numerical weather prediction at ECMWF and provide value-added soil moisture products for the NASA Soil Moisture Active Passive mission.

radar backscatter↗

DOE OSTI2020

Computational modeling of charge hopping dynamics along a disordered one-dimensional wire with energy gradients in quantum environments

This computational study investigates the effects of energy gradients on charge hopping dynamics along a one-dimensional chain of discrete sites coupled to quantum bath, which is modeled at the level of Pauli master equation (PME). This study also assesses the performance of different approximations for the hopping rates. Three different methods for solving the PME, a fourth order Runge–Kutta method, numerical diagonalization of the rate matrix followed by analytic propagation, and kinetic Monte Carlo simulation method, are tested and confirmed to produce virtually identical values of time dependent mean square displacement, diffusion constant, and mobility. Five different rate expressions, exact numerical evaluation of Fermi’s Golden Rule (FGR) rate, stationary phase interpolation (SPI) approximation, semiclassical approximation, classical Marcus rate, and Miller–Abrahams rate, are tested to help understand the effects of approximations in representing quantum environments in the presence of energy gradients. The results based on direct numerical evaluation of FGR rate exhibit transition from diffusive to non-diffusive behavior with the increase in the gradient and show that the charge transport in the quantum bath is more sensitive to the magnitude of the gradient and the disorder than in the classical bath. Among all the four approximations for the hopping rates, the SPI approximation is confirmed to work best overall. A comparison of two different methods to calculate the mobility identifies drift motion of the population distribution as the major source of non-diffusive behavior and provides more reliable information on the contribution of quantum bath.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

Elastic Wave Propagation in Curvilinear Coordinates with Mesh Refinement Interfaces by a Fourth Order Finite Difference Method

In this work, we develop a fourth order accurate finite difference method for the three dimensional elastic wave equation in isotropic media with the piecewise smooth material property. In our model, the material property can be discontinuous at curved interfaces. The governing equations are discretized in second order form on curvilinear meshes by using a fourth order finite difference operator satisfying a summation-by-parts property. The method is energy stable and high order accurate. The highlight is that mesh sizes can be chosen according to the velocity structure of the material so that computational efficiency is improved. At the mesh refinement interfaces with hanging nodes, physical interface conditions are imposed by using ghost points and interpolation. With a fourth order predictor-corrector time integrator, the fully discrete scheme is energy conserving. Numerical experiments are presented to verify the fourth order convergence rate and the energy conserving property.

58 GEOSCIENCES↗

DOE OSTI2024

Data for "Design of Diverse, Functional Mitochondrial Targeting Sequences Across Eukaryotic Organisms Using Variational Autoencoder"

Mitochondria play a key role in energy production and metabolism, making them a promising target for metabolic engineering and disease treatment. However, despite the known influence of passenger proteins on localization efficiency, only a few protein-localization tags have been characterized for mitochondrial targeting. To address this limitation, we leverage a Variational Autoencoder to design novel mitochondrial targeting sequences. In silico analysis reveals that a high fraction of the generated peptides (90.14%) are functional and possess features important for mitochondrial targeting. We characterize artificial peptides in four eukaryotic organisms and, as a proof-of-concept, demonstrate their utility in increasing 3-hydroxypropionic acid titers through pathway compartmentalization and improving 5-aminolevulinate synthase delivery by 1.62-fold and 4.76-fold, respectively. Moreover, we employ latent space interpolation to shed light on the evolutionary origins of dual-targeting sequences. Overall, our work demonstrates the potential of generative artificial intelligence for both fundamental research and practical applications in mitochondrial biology.

AI/ML↗

NASA NTRS2004

Regional Comparison and Assimilation of GOCART and MODIS Aerosol Optical Depth across the Eastern U.S.

This study compares aerosol optical depths (AOD) products from the Moderate Resolution Imaging Spectroradiometer (MODIS) and the Goddard Chemistry Aerosol Radiation and Transport (GOCART) model and their integrated products with ground measurements across the eastern U.S. from March 1, 2000 to December 31, 2001. The Terra MODIS Level-3 (collection 4) AOD at 0.55 pm has better correlation, but consistently overestimates the values of the Aerosol Robotic Network (AERONET) measurements. GOCART has small biases for a 22-month integration, and slight positive biases are appeared for the cold season. These results are also supported by the comparison with the IMPROVE (Interagency Monitoring of Protected Visual Environments) light extinction index. The optimal interpolation improves the daily-scale RMSE from either MODIS or GOCART alone. However, the regional biases in the aerosol products constitute a major constraint to the optimal estimate of AOD.

Matsui, Toshihisa↗

NASA NTRS2016

Development and Evaluation of High-Resolution Climate Simulations Over the Mountainous Northeastern United States

The mountain regions of the northeastern United States are a critical socioeconomic resource for Vermont, New York State, New Hampshire, Maine, and southern Quebec. While global climate models (GCMs) are important tools for climate change risk assessment at regional scales, even the increased spatial resolution of statistically downscaled GCMs (commonly approximately 1/ 8 deg) is not sufficient for hydrologic, ecologic, and land-use modeling of small watersheds within the mountainous Northeast. To address this limitation, an ensemble of topographically downscaled, high-resolution (30"), daily 2-m maximum air temperature; 2-m minimum air temperature; and precipitation simulations are developed for the mountainous Northeast by applying an additional level of downscaling to intermediately downscaled (1/ 8 deg) data using high-resolution topography and station observations. First, observed relationships between 2-m air temperature and elevation and between precipitation and elevation are derived. Then, these relationships are combined with spatial interpolation to enhance the resolution of intermediately downscaled GCM simulations. The resulting topographically downscaled dataset is analyzed for its ability to reproduce station observations. Topographic downscaling adds value to intermediately downscaled maximum and minimum 2-m air temperature at high-elevation stations, as well as moderately improves domain-averaged maximum and minimum 2-m air temperature. Topographic downscaling also improves mean precipitation but not daily probability distributions of precipitation. Overall, the utility of topographic downscaling is dependent on the initial bias of the intermediately downscaled product and the magnitude of the elevation adjustment. As the initial bias or elevation adjustment increases, more value is added to the topographically downscaled product.

climatology↗

DOE OSTI2025

Persistent Classification: Understanding Adversarial Attacks by Studying Decision Boundary Dynamics

ABSTRACT There are a number of hypotheses underlying the existence of adversarial examples for classification problems. These include the high‐dimensionality of the data, the high codimension in the ambient space of the data manifolds of interest, and that the structure of machine learning models may encourage classifiers to develop decision boundaries close to data points. This article proposes a new framework for studying adversarial examples that does not depend directly on the distance to the decision boundary. Similarly to the smoothed classifier literature, we define a (natural or adversarial) data point to be ( γ , σ)‐stable if the probability of the same classification is at least for points sampled in a Gaussian neighborhood of the point with a given standard deviation . We focus on studying the differences between persistence metrics along interpolants of natural and adversarial points. We show that adversarial examples have significantly lower persistence than natural examples for large neural networks in the context of the MNIST and ImageNet datasets. We connect this lack of persistence with decision boundary geometry by measuring angles of interpolants with respect to decision boundaries. Finally, we connect this approach with robustness by developing a manifold alignment gradient metric and demonstrating the increase in robustness that can be achieved when training with the addition of this metric.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

NASA NTRS1995

Development of improved space sampling strategies for ocean chemical properties: Total carbon dioxide and dissolved nitrate

Large-scale ocean observing programs such as the Joint Global Ocean Flux Study (JGOFS) and the World Ocean Circulation Experiment (WOCE) today, must face the problem of designing an adequate sampling strategy. For ocean chemical variables, the goals and observing technologies are quite different from ocean physical variables (temperature, salinity, pressure). We have recently acquired data on the ocean CO2 properties on WOCE cruises P16c and P17c that are sufficiently dense to test for sampling redundancy. We use linear and quadratic interpolation methods on the sampled field to investigate what is the minimum number of samples required to define the deep ocean total inorganic carbon (TCO2) field within the limits of experimental accuracy (+/- 4 micromol/kg). Within the limits of current measurements, these lines were oversampled in the deep ocean. Should the precision of the measurement be improved, then a denser sampling pattern may be desirable in the future. This approach rationalizes the efficient use of resources for field work and for estimating gridded (TCO2) fields needed to constrain geochemical models.

Goyet, Catherine↗

NASA NTRS1981

Optimal interpolation and the Kalman filter

The estimation theory of stochastic-dynamic systems is described and used in a numerical study of optimal interpolation. The general form of data assimilation methods is reviewed. The Kalman-Bucy, KB filter, and optimal interpolation (OI) filters are examined for effectiveness in performance as gain matrices using a one-dimensional form of the shallow-water equations. Control runs in the numerical analyses were performed for a ten-day forecast in concert with the OI method. The effects of optimality, initialization, and assimilation were studied. It was found that correct initialization is necessary in order to localize errors, especially near boundary points. Also, the use of small forecast error growth rates over data-sparse areas was determined to offset inaccurate modeling of correlation functions near boundaries.

Cohn, S.↗

NASA NTRS1992

Stability correction functions for the mean wind speed and temperature in the unstable surface layer

In accordance with a theoretical analysis by Kader and Yaglom, interpolation expressions are proposed for the Monin-Obukhov functions for the wind speed and temperature profiles under unstable conditions in the surface layer of the atmospheric boundary layer. The values of the parameters are adjusted to agree with earlier work under weak instability, as reviewed by Hogstrom. The resulting formulations are valid for considerably larger values of their argument than those that were hitherto available. Their advantage is that they yield concise and closed-form stability correction functions which can be readily applied in profile analysis and model calculations.

Brutsaert, Wilfried↗

DOE OSTI2023

Recent Hydrodynamic Modeling Enhancements in OpenFAST

To improve the hydrodynamic modeling capabilities of the offshore wind design and modeling tool OpenFAST, a restructuring of the HydroDyn hydrodynamics module was undertaken with several new features implemented. The generation of the wave field is now separated from HydroDyn into a new module called SeaState. Unlike previous versions of HydroDyn, which precomputed the wave kinematics at the undisplaced positions of the hydrodynamic nodes, the new SeaState module computes the complete time history of the wave field at the vertices of a user-defined wave grid. During the simulation, the wave kinematics at any point within the grid can be efficiently interpolated, allowing the strip-theory wave loads to be evaluated based on the wave kinematics at the displaced structure position. A phase correction was also implemented for the potential-flow wave-exciting loads to account for large structure displacements. Several wave-stretching methods were implemented for the strip-theory solution along with a wave-load redistribution method that ensures smooth variation of the nodal loads as the nodes of a discrete hydrodynamic mesh enter and exit the water. The load smoothing is included to avoid the excitation of unphysical high-frequency structural vibration during hydroelastic simulations. The MacCamy-Fuchs diffraction correction for the strip-theory solution and the capability to insert a constrained NewWave into extreme stochastic sea states with directional spreading were also added. The present article documents the formulations of the new improvements to HydroDyn with example applications and numerical results.

constrained wave↗

NASA NTRS1992

Gross-fragmentation of meteoroids and bulk density of Geminids from photographic fireball records

The explicit solution of the drag and ablation equations of a single nonfragmenting meteoroid moving in any actual atmosphere was published several years ago. The solution yields the theoretical relation of l, the distance flown by the meteoroid in its trajectory, as a function of time, t, assuming that the height, h, is a known function of l. The photographic records of meteors and fireballs are coded by time marks, using a rotating shutter or a similar device to break the moving image. Time is, thus, the independent variable and for each time mark on a meteoroid trajector, the observed distance along the trajectory, l sub obs, as well as the double- or multiple- station photographs of the same meteoroid. Applying this solution to all available Prairie Network (PN) fireball-records, we recognized that the majority of them gave good solutions with standard deviations somewhat bigger than the intrinsic geometrical precision of the data. We also noticed that, on an average, previous methods of evaluation of the meteoroid velocities (interpolation polynomials, numerical differenciation of the observed l sub obs) used up to only several tens of percent of the intrinsic precision of the PN observational data. When residuals of these solutions, i.e. l sub obs - l sub com, were represented as a function of time for about 75 percent of solutions. The remaining 25 percent of residuals showed systematic changes with time exceeding one standard deviation. We tried to explain these systematic time course of residuals by using different meteoroids first computed theoretically and then analyzed by the same model as the natural PN fireballs were. The conclusion of these model computations is that systematic time changes of residuals in the nonfragmenting model exceeding one standard of deviation are caused by sudden gross fragmentation at one or more trajectory points. Thus, we generalized the explicit solution of the drag and ablation equations of a single nonfragmenting meteoroid by allowing for one or more points, where sudden gross fragmentation can occur. Using this generalized solution, the distances along the meteoroid trajectory can be computed for any choice of input parameters and compared with the observed distances flown by the meteoroid. For the most precise and long fireball trajectories, the least-squares solution can, thus, yield the initial velocities, the ablation coefficients, the dynamical masses, the positions of gross-fragmentation points, and the terminal mass. At a gross-fragmentation point, the ratio of the main mass to all the remaining fragments can be compared with the dynamic mass determined from our gross-fragmentation model and, thus, the meteoroid bulk density can be evaluated. We applied the gross-fragmentation model to sever PN fireballs showing time changes of residuals, and we recognized that, in all these cases, the new computed bulk densities of meteoroids resulted higher in comparison with the meteoroid densities determined with the non-gross-fragmentation model. Other aspects of the study are discussed.

Ceplecha, Zdenek↗

NASA NTRS2012

Correction to "Influence of Dust and Black Carbon on the Snow Albedo in the NASA Goddard Earth Observing System Version 5 Land Surface Model"

The website information describing the forcing meteorological data used for the land surface model (LSM) simulation, which were observed at an Automated Meteorological Station CAWS) at the Sapporo District Meteorological Observatory maintained by the Japan Meteorological Agency (JMA), was missing from the text. The 1-hourly data were obtained from the website of Kisyoutoukeijouhou (Information for available JMA-observed meteorological data in the past) on the website of JMA (in Japanese) (available at: http://www.jma.go.jpijmaimenulreport.html). The measurement height information of 59.5 m for the anemometer at the Sapporo Observatory was also obtained from the website of JMA (in Japanese) (available at: http://www.jma.go.jp/jma/menu/report.html). In addition, the converted 10-m wind speed, based on the AWS/JMA data, was further converted to a 2-m wind speed prior to its use with the land model as a usual treatment of off-line Catchment simulation. Please ignore the ice absorption data on the website mentioned in paragraph [15] which was not used for our calculations (but the data on the website was mostly the same as the estimated ice absorption coefficients by the following method because they partially used the same data by Warren [1984]). We calculated the ice absorption coefficients with the method mentioned in the same paragraph, for which some of the refractive index data by Warren [1984] were used and then interpolated between wavelengths, and also mentioned in paragraph [20] for the visible (VIS) and near-infrared (NIR) ranges. The optical data we used were interpolated between wavelengths as necessary.

Yasunari, Teppei J.↗

DOE OSTI2022

Recent Hydrodynamic Modeling Enhancements in OpenFAST: Preprint

constrained wave↗

DOE OSTI2024

UFNet: Joint U-Net and Fully Connected Neural Network to Bias Correct Precipitation Predictions from

Paper information. Shuang Yu, Indrasis Chakraborty, Gemma J. Anderson, Donald D. Lucas, Yannic Lops, and Daniel Galea. UFNet: Joint U-Net and fully connected neural network to bias correct precipitation predictions from climate models. Artificial Intelligence for the Earth Systems, 2024. Overview. This work develops the UFNet methodology to correct E3SM historical precipitation projection bias. The UFNet deep learning framework consists of a two-part architecture: a U-Net convolutional network to capture the spatiotemporal distribution of precipitation and a fully connected network to capture the distribution of higher-order statistics. The joint network, termed UFNet, can simultaneously improve the spatial structure of the modeled precipitation and capture the distribution of extreme precipitation values. Below we provide guidance for applying UFNet to correct the Energy Exascale Earth System Model (E3SM; Golaz et al. 2019) daily precipitation projection over the contiguous United States (CONUS). Getting started 1. Obtain the historical climate simulation and observation data. The E3SM historical simulation data are available through https://aims2.llnl.gov/search/cmip6/. The CPC unified gauge-based analysis of daily precipitation can be found through https://psl.noaa.gov/data/gridded/data.cpc.globalprecip.html. The ECMWF atmospheric reanalysis of the 20th century (ERA-20C) data are available through https://www.ecmwf.int/en/forecasts/datasets/reanalysis-datasets/era-20c. The spatial resolution of E3SM and observed datasets are both regridded to a common 1° resolution grid using conservative interpolation. The regridded E3SM, CPC and ERA-20C with 1° resolution can be found throught ./data/. 2. Train the fully connected network (DNN) Python train_dnn.py 3. Train the UFNet Python train_ufnet.py 4. Evaluation and compared with the baseline Python evaluation.py

Lucas, DonaldD↗