Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “interpolation”

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 1,117 records · Page 62
All sourcesNASA NTRSDOE OSTI
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↗
DOE OSTI2024

The development of thermodynamically consistent and physics-informed equation-of-state model through machine learning

Ab initio molecular dynamics (AIMD) simulations have become an important tool used in the construction of equations of state (EOS) tables for warm dense matter. Due to computational costs, only a limited number of system state conditions can be simulated, and the remaining EOS surface must be interpolated for use in radiation-hydrodynamic simulations of experiments. In this work, we develop a thermodynamically consistent EOS model that utilizes a physics-informed machine learning approach to implicitly learn the underlying Helmholtz free-energy from AIMD generated energies and pressures. The model, referred to as PIML-EOS, was trained and tested on warm dense polystyrene producing a fit within a 1% relative error for both energy and pressure and is shown to satisfy both the Maxwell and Gibbs–Duhem relations. In addition, we provide a path toward obtaining thermodynamic quantities, such as the total entropy and chemical potential (containing both ionic and electronic contributions), which are not available from current AIMD simulations.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2025

A self-consistent electrostatic electrified shallow water model for charged liquid surface dynamics

The dynamics of an electrified liquid surface are investigated using a shallow water model that is self-consistently coupled with an electrostatic solver. To account for the spatiotemporal variation of a curved liquid surface, the electric field on curved surfaces is calculated by solving a two-dimensional electrostatic equation using the weighted least squares (WLSQ) interpolation method. First, the WLSQ implementation is verified with analytical theory obtained from a Laplace solution assuming a half-cylinder liquid surface profile. Then, the coupled electrified shallow water model is used to study the instability of liquid surface perturbation as a function of the potential drop between an electrode and conducting liquid, surface tension, and gravity. We present a linear dispersion theory of liquid surface instability for long-wavelength perturbations, including the effects of liquid viscosity. Furthermore, the effects of multiple sinusoidal surface perturbations on the electrified liquid surface instability are investigated.

Capillary waves↗
DOE OSTI2013

A basis in an invariant subspace of analytic functions

The existence problem for a basis in a differentiation-invariant subspace of analytic functions defined in a bounded convex domain in the complex plane is investigated. Conditions are found for the solvability of a certain special interpolation problem in the space of entire functions of exponential type with conjugate diagrams lying in a fixed convex domain. These underlie sufficient conditions for the existence of a basis in the invariant subspace. This basis consists of linear combinations of eigenfunctions and associated functions of the differentiation operator, whose exponents are combined into relatively small clusters. Necessary conditions for the existence of a basis are also found. Under a natural constraint on the number of points in the groups, these coincide with the sufficient conditions. That is, a criterion is found under this constraint that a basis constructed from relatively small clusters exists in an invariant subspace of analytic functions in a bounded convex domain in the complex plane. Bibliography: 25 titles.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Deep Learning for Multigroup Cross-Section Representation in Two-Step Core Calculations

Here we investigate using deep learning, a type of machine-learning algorithm employing multiple layers of artificial neurons, for the mathematical representation of multigroup cross sections for use in the Griffin reactor multiphysics code for two-step deterministic neutronics calculations. A three-dimensional fuel element typical of a high-temperature gas reactor as well as a two-dimensional sodium-cooled fast reactor lattice are modeled using the Serpent Monte Carlo code, and multigroup macroscopic cross sections are generated for various state parameters to produce a training data set and a separate validation data set. A fully connected, feedforward neural network is trained using the open-source PyTorch machine-learning framework, and its accuracy is compared against the standard piecewise linear interpolation model. Additionally, we provide in this work a generic technique for propagating the cross-section model errors up to the k eff using sensitivity coefficients with the first-order uncertainty propagation rule. Quantifying the eigenvalue error due to the cross-section regression errors is especially practical for appropriately selecting the mathematical representation of the cross sections. We demonstrate that the artificial neural network model produces lower errors and therefore enables better accuracy relative to the piecewise linear model when the cross sections exhibit nonlinear dependencies; especially when a coarse grid is employed, where the errors can be halved by the artificial neural network. However, for linearly dependent multigroup cross sections as found for the sodium-cooled fast reactor case, a simpler linear regression outperforms deeper networks.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2026

ReMU: regional minimal updating for model-based derivative-free optimization

Derivative-free optimization (DFO) problems are optimization problems where derivative information is unavailable or extremely difficult to obtain. Model-based DFO solvers have been applied extensively in scientific computing. Powell's NEWUOA (2004) [Powell, The NEWUOA software for unconstrained optimization without derivatives, in Large-Scale Nonlinear Optimization, Nonconvex Optimization and its Applications Vol. 83, G. Di Pillo and M. Roma, eds., Springer, 2006, pp. 255–297] and Wild's POUNDerS (2014) [Wild, Solving derivative-free nonlinear least squares problems with POUNDERS, in Advances and Trends in Optimization with Engineering Applications, T. Terlaky, M.F. Anjos, and S. Ahmed, eds., SIAM, 2017, pp. 529–540] explore the numerical power of the minimal norm Hessian (MNH) model for DFO and contributed to the open discussion on building better models with fewer data to achieve faster numerical convergence. Another decade later, we propose the regional minimal updating (ReMU) models, and extend the previous models into a broader class, including the H 2 norm models [Xie and Yuan, Least H 2 norm updating of quadratic interpolation models for derivative-free trust-region algorithms, IMA J. Numer. Anal. 46 (2025), pp. 21–50]. This paper shows motivation behind ReMU models, computational details, theoretical and numerical results on particular extreme points and the barycentre of ReMU's weight coefficient region, and the associated KKT matrix error and distance. Novel metrics, such as the truncated Newton step error, are proposed to numerically understand the new models' properties. A new algorithmic strategy, based on iteratively adjusting the ReMU model type, is also proposed, and shows numerical advantages by combining and switching between the barycentric model and the classic least Frobenius norm model in an online fashion.

derivative-free trust-region methods↗
DOE OSTI2013

70 DA WHITE DWARFS IDENTIFIED IN LAMOST PILOT SURVEY

We present a spectroscopically identified catalog of 70 DA white dwarfs (WDs) from the LAMOST pilot survey. Thirty-five are found to be new identifications after cross-correlation with the Eisenstein et al. and Villanova catalogs. The effective temperature and gravity of these WDs are estimated by Balmer lines fitting. Most of them are hot WDs. The cooling times and masses of these WDs are estimated by interpolation in theoretical evolution tracks. The peak of the mass distribution is found to be {approx}0.6 M {sub Sun }, which is consistent with prior work in the literature. The distances of these WDs are estimated using the method of synthetic spectral distances. All of these WDs are found to be in the Galactic disk from our analysis of space motions. Our sample supports the expectation that WDs with high mass are concentrated near the plane of the Galactic disk.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

Relation between the standard perturbation theory and regularized multi-point propagator method

We investigate the relation between the regularized multi-propagator method, called 'Reg PT', and the standard perturbation theory. Reg PT is one of the most successful models to describe nonlinear evolution of dark matter fluctuations. However, Reg PT is a mathematically unproven interpolation formula between the large-scale solution calculated by the standard perturbation theory and the limiting solution in the small scale calculated by the multi-point propagator method. In this paper, we give an alternative explanation for Reg PT in the context of the standard perturbation theory, showing that Reg PT does not ever have more effective information on nonlinear matter evolution than the standard perturbation theory. In other words, the solutions of the standard perturbation theory reproduce the results of N-body simulations better than those of Reg PT, especially at the high-k region. This fact means that the standard perturbation theory at the two-loop level is still one of the best predictions of the nonlinear power spectrum to date. Nevertheless, the standard perturbation theory has not been preferred because of the divergent behavior of the solution at small scales. To solve this problem, we also propose a modified standard perturbation theory which avoids the divergence.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

AB INITIO EQUATION OF STATE FOR HYDROGEN-HELIUM MIXTURES WITH RECALIBRATION OF THE GIANT-PLANET MASS-RADIUS RELATION

Using density functional molecular dynamics simulations, we determine the equation of state (EOS) for hydrogen-helium mixtures spanning density-temperature conditions typical of giant-planet interiors, {approx}0.2-9 g cm{sup -3} and 1000-80,000 K for a typical helium mass fraction of 0.245. In addition to computing internal energy and pressure, we determine the entropy using an ab initio thermodynamic integration technique. A comprehensive EOS table with 391 density-temperature points is constructed and the results are presented in the form of a two-dimensional free energy fit for interpolation. Deviations between our ab initio EOS and the semi-analytical EOS model by Saumon and Chabrier are analyzed in detail, and we use the results for initial revision of the inferred thermal state of giant planets with known values for mass and radius. Changes are most pronounced for planets in the Jupiter mass range and below. We present a revision to the mass-radius relationship that makes the hottest exoplanets increase in radius by {approx}0.2 Jupiter radii at fixed entropy and for masses greater than {approx}0.5 Jupiter mass. This change is large enough to have possible implications for some discrepant ''inflated giant exoplanets''.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

Isochrones for old (>5 Gyr) stars and stellar populations. I. Models for –2.4 ≤ [Fe/H] ≤+0.6, 0.25 ≤ Y ≤ 0.33, and –0.4 ≤ [α/Fe] ≤+0.4

Canonical grids of stellar evolutionary sequences have been computed for the helium mass-fraction abundances Y = 0.25, 0.29, and 0.33, and for iron abundances that vary from –2.4 to +0.4 (in 0.2 dex increments) when [α/Fe] =+0.4, or for the ranges –2.0 ≤ [Fe/H] ≤+0.6, –1.8 ≤ [Fe/H] ≤+0.6 when [α/Fe] =0.0 and –0.4, respectively. The grids, which consist of tracks for masses from 0.12 M{sub ⊙} to 1.1-1.5 M{sub ⊙} (depending on the metallicity) are based on up-to-date physics, including the gravitational settling of helium (but not metals diffusion). Interpolation software is provided to generate isochrones for arbitrary ages between ≈5 and 15 Gyr and any values of Y, [α/Fe], and [Fe/H] within the aformentioned ranges. Comparisons of isochrones with published color-magnitude diagrams (CMDs) for the open clusters M67 ([Fe/H] ≈0.0) and NGC 6791 ([Fe/H] ≈0.3) and for four of the metal-poor globular clusters (47 Tuc, M3, M5, and M92) indicate that the models for the observed metallicities do a reasonably good job of reproducing the locations and slopes of the cluster main sequences and giant branches. The same conclusion is reached from a consideration of plots of nearby subdwarfs that have accurate Hipparcos parallaxes and metallicities in the range –2.0 ≲ [Fe/H] ≲ –1.0 on various CMDs and on the (log T {sub eff}, M{sub V} ) diagram. A relatively hot temperature scale similar to that derived in recent calibrations of the infrared flux method is favored by both the isochrones and the adopted color transformations, which are based on the latest MARCS model atmospheres.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

Optimized blind gamma-ray pulsar searches at fixed computing budget

The sensitivity of blind gamma-ray pulsar searches in multiple years worth of photon data, as from the Fermi LAT, is primarily limited by the finite computational resources available. Addressing this 'needle in a haystack' problem, here we present methods for optimizing blind searches to achieve the highest sensitivity at fixed computing cost. For both coherent and semicoherent methods, we consider their statistical properties and study their search sensitivity under computational constraints. The results validate a multistage strategy, where the first stage scans the entire parameter space using an efficient semicoherent method and promising candidates are then refined through a fully coherent analysis. We also find that for the first stage of a blind search incoherent harmonic summing of powers is not worthwhile at fixed computing cost for typical gamma-ray pulsars. Further enhancing sensitivity, we present efficiency-improved interpolation techniques for the semicoherent search stage. Via realistic simulations we demonstrate that overall these optimizations can significantly lower the minimum detectable pulsed fraction by almost 50% at the same computational expense.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

COSMIC EMULATION: FAST PREDICTIONS FOR THE GALAXY POWER SPECTRUM

The halo occupation distribution (HOD) approach has proven to be an effective method for modeling galaxy clustering and bias. In this approach, galaxies of a given type are probabilistically assigned to individual halos in N-body simulations. In this paper, we present a fast emulator for predicting the fully nonlinear galaxy–galaxy auto and galaxy–dark matter cross power spectrum and correlation function over a range of freely specifiable HOD modeling parameters. The emulator is constructed using results from 100 HOD models run on a large ΛCDM N-body simulation, with Gaussian Process interpolation applied to a PCA-based representation of the galaxy power spectrum. The total error is currently ∼1% in the auto correlations and ∼2% in the cross correlations from z = 1 to z = 0, over the considered parameter range. We use the emulator to investigate the accuracy of various analytic prescriptions for the galaxy power spectrum, parametric dependencies in the HOD model, and the behavior of galaxy bias as a function of HOD parameters. Additionally, we obtain fully nonlinear predictions for tangential shear correlations induced by galaxy–galaxy lensing from our galaxy–dark matter cross power spectrum emulator. All emulation products are publicly available at http://www.hep.anl.gov/cosmology/CosmicEmu/emu.html.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

CONSTRUCTING A FLEXIBLE LIKELIHOOD FUNCTION FOR SPECTROSCOPIC INFERENCE

We present a modular, extensible likelihood framework for spectroscopic inference based on synthetic model spectra. The subtraction of an imperfect model from a continuously sampled spectrum introduces covariance between adjacent datapoints (pixels) into the residual spectrum. For the high signal-to-noise data with large spectral range that is commonly employed in stellar astrophysics, that covariant structure can lead to dramatically underestimated parameter uncertainties (and, in some cases, biases). We construct a likelihood function that accounts for the structure of the covariance matrix, utilizing the machinery of Gaussian process kernels. This framework specifically addresses the common problem of mismatches in model spectral line strengths (with respect to data) due to intrinsic model imperfections (e.g., in the atomic/molecular databases or opacity prescriptions) by developing a novel local covariance kernel formalism that identifies and self-consistently downweights pathological spectral line “outliers.” By fitting many spectra in a hierarchical manner, these local kernels provide a mechanism to learn about and build data-driven corrections to synthetic spectral libraries. An open-source software implementation of this approach is available at http://iancze.github.io/Starfish, including a sophisticated probabilistic scheme for spectral interpolation when using model libraries that are sparsely sampled in the stellar parameters. We demonstrate some salient features of the framework by fitting the high-resolution V-band spectrum of WASP-14, an F5 dwarf with a transiting exoplanet, and the moderate-resolution K-band spectrum of Gliese 51, an M5 field dwarf.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

AN ALGORITHM FOR RADIATION MAGNETOHYDRODYNAMICS BASED ON SOLVING THE TIME-DEPENDENT TRANSFER EQUATION

We describe a new algorithm for solving the coupled frequency-integrated transfer equation and the equations of magnetohydrodynamics in the regime that light-crossing time is only marginally shorter than dynamical timescales. The transfer equation is solved in the mixed frame, including velocity-dependent source terms accurate to O(v/c). An operator split approach is used to compute the specific intensity along discrete rays, with upwind monotonic interpolation used along each ray to update the transport terms, and implicit methods used to compute the scattering and absorption source terms. Conservative differencing is used for the transport terms, which ensures the specific intensity (as well as energy and momentum) are conserved along each ray to round-off error. The use of implicit methods for the source terms ensures the method is stable even if the source terms are very stiff. To couple the solution of the transfer equation to the MHD algorithms in the ATHENA code, we perform direct quadrature of the specific intensity over angles to compute the energy and momentum source terms. We present the results of a variety of tests of the method, such as calculating the structure of a non-LTE atmosphere, an advective diffusion test, linear wave convergence tests, and the well-known shadow test. We use new semi-analytic solutions for radiation modified shocks to demonstrate the ability of our algorithm to capture the effects of an anisotropic radiation field accurately. Since the method uses explicit differencing of the spatial operators, it shows excellent weak scaling on parallel computers.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

SEPARATION OF THE RIBBON FROM GLOBALLY DISTRIBUTED ENERGETIC NEUTRAL ATOM FLUX USING THE FIRST FIVE YEARS OF IBEX OBSERVATIONS

The Interstellar Boundary Explorer (IBEX) observes the IBEX ribbon, which stretches across much of the sky observed in energetic neutral atoms (ENAs). The ribbon covers a narrow (∼20°-50°) region that is believed to be roughly perpendicular to the interstellar magnetic field. Superimposed on the IBEX ribbon is the globally distributed flux that is controlled by the processes and properties of the heliosheath. This is a second study that utilizes a previously developed technique to separate ENA emissions in the ribbon from the globally distributed flux. A transparency mask is applied over the ribbon and regions of high emissions. We then solve for the globally distributed flux using an interpolation scheme. Previously, ribbon separation techniques were applied to the first year of IBEX-Hi data at and above 0.71 keV. Here we extend the separation analysis down to 0.2 keV and to five years of IBEX data enabling first maps of the ribbon and the globally distributed flux across the full sky of ENA emissions. Our analysis shows the broadening of the ribbon peak at energies below 0.71 keV and demonstrates the apparent deformation of the ribbon in the nose and heliotail. We show global asymmetries of the heliosheath, including both deflection of the heliotail and differing widths of the lobes, in context of the direction, draping, and compression of the heliospheric magnetic field. We discuss implications of the ribbon maps for the wide array of concepts that attempt to explain the ribbon's origin. Thus, we present the five-year separation of the IBEX ribbon from the globally distributed flux in preparation for a formal IBEX data release of ribbon and globally distributed flux maps to the heliophysics community.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

RICH: OPEN-SOURCE HYDRODYNAMIC SIMULATION ON A MOVING VORONOI MESH

We present here RICH, a state-of-the-art two-dimensional hydrodynamic code based on Godunov’s method, on an unstructured moving mesh (the acronym stands for Racah Institute Computational Hydrodynamics). This code is largely based on the code AREPO. It differs from AREPO in the interpolation and time-advancement schemeS as well as a novel parallelization scheme based on Voronoi tessellation. Using our code, we study the pros and cons of a moving mesh (in comparison to a static mesh). We also compare its accuracy to other codes. Specifically, we show that our implementation of external sources and time-advancement scheme is more accurate and robust than is AREPO when the mesh is allowed to move. We performed a parameter study of the cell rounding mechanism (Lloyd iterations) and its effects. We find that in most cases a moving mesh gives better results than a static mesh, but it is not universally true. In the case where matter moves in one way and a sound wave is traveling in the other way (such that relative to the grid the wave is not moving) a static mesh gives better results than a moving mesh. We perform an analytic analysis for finite difference schemes that reveals that a Lagrangian simulation is better than a Eulerian simulation in the case of a highly supersonic flow. Moreover, we show that Voronoi-based moving mesh schemes suffer from an error, which is resolution independent, due to inconsistencies between the flux calculation and the change in the area of a cell. Our code is publicly available as open source and designed in an object-oriented, user-friendly way that facilitates incorporation of new algorithms and physical processes.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

Evolution of entanglement entropy in the early universe

We investigate the entropy evolution in the early universe by computing the change of the entanglement entropy in Freedmann-Robertson-Walker quantum cosmology in the presence of particle horizon. The matter is modeled by a Chaplygin gas so as to provide a smooth interpolation between inflationary and radiation epochs, rendering the evolution of entropy from early time to late time trackable. We found that soon after the onset of the inflation, the total entanglement entropy rapidly decreases to a minimum. It then rises monotonically in the remainder of the inflation epoch as well as the radiation epoch. Our result is in qualitative agreement with the area law of Ryu and Takayanagi including the logarithmic correction. We comment on the possible implication of our finding to the cosmological entropy problem.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

A general reconstruction of the recent expansion history of the universe

Distance measurements are currently the most powerful tool to study the expansion history of the universe without specifying its matter content nor any theory of gravitation. Assuming only an isotropic, homogeneous and flat universe, in this work we introduce a model-independent method to reconstruct directly the deceleration function via a piecewise function. Including a penalty factor, we are able to vary continuously the complexity of the deceleration function from a linear case to an arbitrary (n+1)-knots spline interpolation. We carry out a Monte Carlo (MC) analysis to determine the best penalty factor, evaluating the bias-variance trade-off, given the uncertainties of the SDSS-II and SNLS supernova combined sample (JLA), compilations of baryon acoustic oscillation (BAO) and H(z) data. The bias-variance analysis is done for three fiducial models with different features in the deceleration curve. We perform the MC analysis generating mock catalogs and computing their best-fit. For each fiducial model, we test different reconstructions using, in each case, more than 10{sup 4} catalogs in a total of about 5× 10{sup 5}. This investigation proved to be essential in determining the best reconstruction to study these data. We show that, evaluating a single fiducial model, the conclusions about the bias-variance ratio are misleading. We determine the reconstruction method in which the bias represents at most 10% of the total uncertainty. In all statistical analyses, we fit the coefficients of the deceleration function along with four nuisance parameters of the supernova astrophysical model. For the full sample, we also fit H{sub 0} and the sound horizon r{sub s}(z{sub d}) at the drag redshift. The bias-variance trade-off analysis shows that, apart from the deceleration function, all other estimators are unbiased. Finally, we apply the Ensemble Sampler Markov Chain Monte Carlo (ESMCMC) method to explore the posterior of the deceleration function up to redshift 1.3 (using only JLA) and 2.3 (JLA+BAO+H(z)). We obtain that the standard cosmological model agrees within 3σ level with the reconstructed results in the whole studied redshift intervals. Since our method is calibrated to minimize the bias, the error bars of the reconstructed functions are a good approximation for the total uncertainty.

79 ASTRONOMY AND ASTROPHYSICS↗