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,477 records · Page 82
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

SDSS-IV MaStar: Stellar parameter determination with continuum-supplemented full-spectrum fitting

We present a stellar parameter catalog built to accompany the MaStar Stellar Library, which is a comprehensive collection of empirical, medium-resolution stellar spectra. We constructed this parameter catalog by using a multicomponent χ 2 fitting approach to match MaStar spectra to models generated by interpolating the ATLAS9-based BOSZ model spectra. The total χ 2 for a given model is defined as the sum of components constructed to characterize narrow-band features of observed spectra (e.g., absorption lines) and the broadband continuum shape separately. Extinction and systematics due to flux calibration were taken into account in the fitting. The χ 2 distribution for a given region of model space was sampled using a Markov chain Monte Carlo (MCMC) algorithm, the data from which were then used to extract atmospheric parameter estimates (T eff , log g, [Fe/H], and [α/Fe]), their corresponding uncertainties, and direct extinction measurements. Two methods were used to extract parameters and uncertainties: one that accepts the MCMC’s prescribed minimum-χ 2 result, and one that uses Bayesian inference to compute a likelihood-weighted mean from the χ 2 distribution sampled by the MCMC. Results were evaluated for internal consistency using repeat observations where available and by comparing them with external data sets (e.g., APOGEE-2 and Gaia DR2). Our spectral-fitting exercise reveals possible deficiencies in current theoretical model spectra, illustrating the potential power of MaStar spectra for helping to improve the models. This paper represents an update to the parameters that were originally presented with SDSS-IV DR17. The MaStar parameter catalog containing our BestFit results is available on the SDSS-IV DR17 website as part of version 2 of the MaStar stellar parameter value-added catalog.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Improved temperature dependence of rate coefficients for rotational state-to-state transitions in H 2 O + H 2 O collisions

Aims. We present an improved database of temperature-dependent rate coefficients for rotational state-to-state transitions in H 2 O + H 2 O collisions. The database includes 231 transitions between the lower para-states of H 2 O and 210 transitions between its lower ortho-states (up toj= 7) and can be employed in cometary and planetary applications up to the temperature of 1000 K. Methods. We developed and applied a new general method that allows the generation of rate coefficients for excitation and quenching processes that automatically satisfy the principle of microscopic reversibility and also helps to cover the range of low collision energies by interpolation of cross sections between the process threshold and the computed data points. Results. We find that in the range of intermediate temperatures, 150

Astronomy & Astrophysics↗
DOE OSTI2024

Performance of high-order Godunov-type methods in simulations of astrophysical low Mach number flows

High-order Godunov methods for gas dynamics have become a standard tool for simulating different classes of astrophysical flows. Their accuracy is mostly determined by the spatial interpolant used to reconstruct the pair of Riemann states at cell interfaces and by the Riemann solver that computes the interface fluxes. In most Godunov-type methods, these two steps can be treated independently, so that many different schemes can in principle be built from the same numerical framework. Because astrophysical simulations often test out the limits of what is feasible with the computational resources available, it is essential to find the scheme that produces the numerical solution with the desired accuracy at the lowest computational cost. However, establishing the best combination of numerical options in a Godunov-type method to be used for simulating a complex hydrodynamic problem is a nontrivial task. In fact, formally more accurate schemes do not always outperform simpler and more diffusive methods, especially if sharp gradients are present in the flow. For this work, we used our fully compressible Seven-League Hydro (SLH) code to test the accuracy of six reconstruction methods and three approximate Riemann solvers on two- and three-dimensional (2D and 3D) problems involving subsonic flows only. We considered Mach numbers in the range from 10 −3 to 10 −1 , which are characteristic of many stellar and geophysical flows. In particular, we considered a well-posed, 2D, Kelvin–Helmholtz instability problem and a 3D turbulent convection zone that excites internal gravity waves in an overlying stable layer. Although the different combinations of numerical methods converge to the same solution with increasing grid resolution for most of the quantities analyzed here, we find that (i) there is a spread of almost four orders of magnitude in computational cost per fixed accuracy between the methods tested in this study, with the most performant method being a combination of a low-dissipation Riemann solver and a sextic reconstruction scheme; (ii) the low-dissipation solver always outperforms conventional Riemann solvers on a fixed grid when the reconstruction scheme is kept the same; (iii) in simulations of turbulent flows, increasing the order of spatial reconstruction reduces the characteristic dissipation length scale achieved on a given grid even if the overall scheme is only second order accurate; (iv) reconstruction methods based on slope-limiting techniques tend to generate artificial, high-frequency acoustic waves during the evolution of the flow; and (v) unlimited reconstruction methods introduce oscillations in the thermal stratification near the convective boundary, where the entropy gradient is steep.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

Uncertainty analysis in geospatial merit matrix–based hydropower resource assessment

Hydraulic head and mean annual streamflow, two main input parameters in hydropower resource assessment, are not measured at every point along the stream. Translation and interpolation are used to derive these parameters, resulting in uncertainties. This study estimates the uncertainties and their effects on model output parameters: the total potential power and the number of potential locations (stream-reach). These parameters are quantified through Monte Carlo Simulation (MCS) linking with a geospatial merit matrix based hydropower resource assessment (GMM-HRA) Model. The methodology is applied to flat, mild, and steep terrains. Results show that the uncertainty associated with the hydraulic head is within 20% for mild and steep terrains, and the uncertainty associated with streamflow is around 16% for all three terrains. Output uncertainty increases as input uncertainty increases. However, output uncertainty is around 10% to 20% of the input uncertainty, demonstrating the robustness of the GMM-HRA model. Hydraulic head is more sensitive to output parameters in steep terrain than in flat and mild terrains. Furthermore, mean annual streamflow is more sensitive to output parameters in flat terrain.

13 HYDRO ENERGY↗
DOE OSTI2014

Virtual edge illumination and one dimensional beam tracking for absorption, refraction, and scattering retrieval

We propose two different approaches to retrieve x-ray absorption, refraction, and scattering signals using a one dimensional scan and a high resolution detector. The first method can be easily implemented in existing procedures developed for edge illumination to retrieve absorption and refraction signals, giving comparable image quality while reducing exposure time and delivered dose. The second method tracks the variations of the beam intensity profile on the detector through a multi-Gaussian interpolation, allowing the additional retrieval of the scattering signal.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

Investigations of α-helix↔β-sheet transition pathways in a miniprotein using the finite-temperature string method

A parallel implementation of the finite-temperature string method is described, which takes into account the invariance of coordinates with respect to rigid-body motions. The method is applied to the complex α-helix↔β-sheet transition in a β-hairpin miniprotein in implicit solvent, which exhibits much of the complexity of conformational changes in proteins. Two transition paths are considered, one derived from a linear interpolant between the endpoint structures and the other derived from a targeted dynamics simulation. Two methods for computing the conformational free energy (FE) along the string are compared, a restrained method, and a tessellation method introduced by E. Vanden-Eijnden and M. Venturoli [J. Chem. Phys. 130, 194103 (2009)]. It is found that obtaining meaningful free energy profiles using the present atom-based coordinates requires restricting sampling to a vicinity of the converged path, where the hyperplanar approximation to the isocommittor surface is sufficiently accurate. This sampling restriction can be easily achieved using restraints or constraints. The endpoint FE differences computed from the FE profiles are validated by comparison with previous calculations using a path-independent confinement method. The FE profiles are decomposed into the enthalpic and entropic contributions, and it is shown that the entropy difference contribution can be as large as 10 kcal/mol for intermediate regions along the path, compared to 15–20 kcal/mol for the enthalpy contribution. This result demonstrates that enthalpic barriers for transitions are offset by entropic contributions arising from the existence of different paths across a barrier. The possibility of using systematically coarse-grained representations of amino acids, in the spirit of multiple interaction site residue models, is proposed as a means to avoid ad hoc sampling restrictions to narrow transition tubes.

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

Stochastic many-body perturbation theory for anharmonic molecular vibrations

A new quantum Monte Carlo (QMC) method for anharmonic vibrational zero-point energies and transition frequencies is developed, which combines the diagrammatic vibrational many-body perturbation theory based on the Dyson equation with Monte Carlo integration. The infinite sums of the diagrammatic and thus size-consistent first- and second-order anharmonic corrections to the energy and self-energy are expressed as sums of a few m- or 2m-dimensional integrals of wave functions and a potential energy surface (PES) (m is the vibrational degrees of freedom). Each of these integrals is computed as the integrand (including the value of the PES) divided by the value of a judiciously chosen weight function evaluated on demand at geometries distributed randomly but according to the weight function via the Metropolis algorithm. In this way, the method completely avoids cumbersome evaluation and storage of high-order force constants necessary in the original formulation of the vibrational perturbation theory; it furthermore allows even higher-order force constants essentially up to an infinite order to be taken into account in a scalable, memory-efficient algorithm. The diagrammatic contributions to the frequency-dependent self-energies that are stochastically evaluated at discrete frequencies can be reliably interpolated, allowing the self-consistent solutions to the Dyson equation to be obtained. This method, therefore, can compute directly and stochastically the transition frequencies of fundamentals and overtones as well as their relative intensities as pole strengths, without fixed-node errors that plague some QMC. It is shown that, for an identical PES, the new method reproduces the correct deterministic values of the energies and frequencies within a few cm{sup −1} and pole strengths within a few thousandths. With the values of a PES evaluated on the fly at random geometries, the new method captures a noticeably greater proportion of anharmonic effects.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

Electrical and thermal conductivities in dense plasmas

Expressions for the electrical and thermal conductivities in dense plasmas are derived combining the Chester-Thellung-Kubo-Greenwood approach and the Kramers approximation. The infrared divergence is removed assuming a Drude-like behaviour. An analytical expression is obtained for the Lorenz number that interpolates between the cold solid-state and the hot plasma phases. An expression for the electrical resistivity is proposed using the Ziman-Evans formula, from which the thermal conductivity can be deduced using the analytical expression for the Lorenz number. The present method can be used to estimate electrical and thermal conductivities of mixtures. Comparisons with experiment and quantum molecular dynamics simulations are done.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

Structural and electronic properties of Si{sub 1–x}Ge{sub x} alloy nanowires

We present first-principles density-functional calculations of Si{sub 1–x}Ge{sub x} alloy nanowires. We show that given the composition of the alloy, the structural properties of the nanowires can be predicted with great accuracy by means of Vegard's law, linearly interpolating the values of a pure Si and a pure Ge nanowire of the same diameter. The same holds, to some extent, also for electronic properties such as the band-gap. We also assess to what extend the band-gap varies as a function of disorder, i.e., how it changes for different random realization of a given concentration. These results make possible to tailor the desired properties of SiGe alloy nanowires starting directly from the data relative to the pristine wires.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2014

Data approximation using a blending type spline construction

Generalized expo-rational B-splines (GERBS) is a blending type spline construction where local functions at each knot are blended together by C{sup k}-smooth basis functions. One way of approximating discrete regular data using GERBS is by partitioning the data set into subsets and fit a local function to each subset. Partitioning and fitting strategies can be devised such that important or interesting data points are interpolated in order to preserve certain features. We present a method for fitting discrete data using a tensor product GERBS construction. The method is based on detection of feature points using differential geometry. Derivatives, which are necessary for feature point detection and used to construct local surface patches, are approximated from the discrete data using finite differences.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

A boundary integral formalism for stochastic ray tracing in billiards

Determining the flow of rays or non-interacting particles driven by a force or velocity field is fundamental to modelling many physical processes. These include particle flows arising in fluid mechanics and ray flows arising in the geometrical optics limit of linear wave equations. In many practical applications, the driving field is not known exactly and the dynamics are determined only up to a degree of uncertainty. This paper presents a boundary integral framework for propagating flows including uncertainties, which is shown to systematically interpolate between a deterministic and a completely random description of the trajectory propagation. A simple but efficient discretisation approach is applied to model uncertain billiard dynamics in an integrable rectangular domain.

42 ENGINEERING↗
DOE OSTI2014

Experimental and theoretical study on the excited-state dynamics of ortho-, meta-, and para-methoxy methylcinnamate

The S{sub 1} state dynamics of methoxy methylcinnamate (MMC) has been investigated under supersonic jet-cooled conditions. The vibrationally resolved S{sub 1}-S{sub 0} absorption spectrum was recorded by laser induced fluorescence and mass-resolved resonant two-photon ionization spectroscopy and separated into conformers by UV-UV hole-burning (UV-UV HB) spectroscopy. The S{sub 1} lifetime measurements revealed different dynamics of para-methoxy methylcinnamate from ortho-methoxy methylcinnamate and meta-methoxy methylcinnamate (hereafter, abbreviated as p-, o-, and m-MMCs, respectively). The lifetimes of o-MMC and m-MMC are on the nanosecond time scale and exhibit little tendency of excess energy dependence. On the other hand, p-MMC decays much faster and its lifetime is conformer and excess energy dependent. In addition, the p-MMC-H{sub 2}O complex was studied to explore the effect of hydration on the S{sub 1} state dynamics of p-MMC, and it was found that the hydration significantly accelerates the nonradiative decay. Quantum chemical calculation was employed to search the major decay route from S{sub 1}(ππ{sup ∗}) for three MMCs and p-MMC-H{sub 2}O in terms of (i) trans → cis isomerization and (ii) internal conversion to the {sup 1}nπ{sup ∗} state. In o-MMC and m-MMC, the large energy barrier is created for the nonradiative decay along (i) the double-bond twisting coordinate (∼1000 cm{sup −1}) in S{sub 1} as well as (ii) the linear interpolating internal coordinate (∼1000 cm{sup −1}) from S{sub 1} to {sup 1}nπ{sup ∗} states. The calculation on p-MMC decay dynamics suggests that both (i) and (ii) are available due to small energy barrier, i.e., 160 cm{sup −1} by the double-bond twisting and 390 cm{sup −1} by the potential energy crossing. The hydration of p-MMC raises the energy barrier of the IC route to the S{sub 1}/{sup 1}nπ{sup ∗} conical intersection, convincing that the direct isomerization is more likely to occur.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2014

Communication: Does force spectroscopy of biomolecules probe their intrinsic dynamic properties?

In single-molecule pulling experiments, the molecule of interest is attached to a much larger object such as an atomic force microscope tip or a micrometer sized bead. The measured dynamics of molecular transitions is therefore affected by the hydrodynamic drag on the pulling instrument itself. By considering the transitions within the combined system (the molecule and the instrument), it is shown here that two distinct physical regimes exist: when the intrinsic stiffness of the molecule is greater than that of the linker connecting the molecule to the pulling setup then the pulling experiment probes the intrinsic dynamics of the molecule with only relatively small (and quantifiable) corrections resulting from the pulling setup. In contrast, when the stiffness of the linker exceeds that of the molecule, the molecular transition in question involves concerted motion of the molecule and the pulling setup and the hydrodynamic drag on the pulling instrument becomes the dominant source of friction along the molecular reaction coordinate. An analytical formula interpolating between these two cases is further derived. These results explain recent conflicting observations where some single-molecule pulling measurements report anomalously low diffusion coefficients along molecular reaction coordinates while others do not.

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

Anyons in the operational formalism

The operational formalism to quantum mechanics seeks to base the theory on a firm foundation of physically well-motivated axioms [1]. It has succeeded in deriving the Feynman rules [2] for general quantum systems. Additional elaborations have applied the same logic to the question of identical particles, confirming the so-called Symmetrization Postulate [3]: that the only two options available are fermions and bosons [4, 5]. However, this seems to run counter to results in two-dimensional systems, which allow for anyons, particles with statistics which interpolate between Fermi-Dirac and Bose-Einstein (see [6] for a review). In this talk we will show that the results in two dimensions can be made compatible with the operational results. That is, we will show that anyonic behavior is a result of the topology of the space in two dimensions [7], and does not depend on the particles being identical; but that nevertheless, if the particles are identical, the resulting system is still anyonic.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Temperature dependent low-field measurements of the magnetocaloric ΔT with sub-mK resolution in small volume and thin film samples

We present temperature dependent ΔT measurements of the magnetocaloric effect in a thin film sample of Gd, employing magnetomodulation and detection of thermal radiation. A bulk sample of the metamagnetic material LaFe{sub 11.05}Co{sub 0.91}Si{sub 1.04} shows a strong broadening of the ΔT peak for increasing field amplitudes between 4 and 45 mT. Bulk Gd in comparison shows only a weak broadening. All investigated samples exhibit a clear quadratic dependence of ΔT on the external field H{sub ext} at the ΔT peak maximum, contrary to earlier predictions. An analytic expression is derived that interpolates between the H{sub ext}{sup 2}-behavior at low and the well-known H{sub ext}{sup 2/3}-behavior at high fields.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Full-dimensional and reduced-dimensional calculations of initial state-selected reaction probabilities studying the H + CH{sub 4} → H{sub 2} + CH{sub 3} reaction on a neural network PES

Initial state-selected reaction probabilities of the H + CH{sub 4} → H{sub 2} + CH{sub 3} reaction are calculated in full and reduced dimensionality on a recent neural network potential [X. Xu, J. Chen, and D. H. Zhang, Chin. J. Chem. Phys. 27, 373 (2014)]. The quantum dynamics calculation employs the quantum transition state concept and the multi-layer multi-configurational time-dependent Hartree approach and rigorously studies the reaction for vanishing total angular momentum (J = 0). The calculations investigate the accuracy of the neutral network potential and study the effect resulting from a reduced-dimensional treatment. Very good agreement is found between the present results obtained on the neural network potential and previous results obtained on a Shepard interpolated potential energy surface. The reduced-dimensional calculations only consider motion in eight degrees of freedom and retain the C{sub 3v} symmetry of the methyl fragment. Considering reaction starting from the vibrational ground state of methane, the reaction probabilities calculated in reduced dimensionality are moderately shifted in energy compared to the full-dimensional ones but otherwise agree rather well. Similar agreement is also found if reaction probabilities averaged over similar types of vibrational excitation of the methane reactant are considered. In contrast, significant differences between reduced and full-dimensional results are found for reaction probabilities starting specifically from symmetric stretching, asymmetric (f{sub 2}-symmetric) stretching, or e-symmetric bending excited states of methane.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

A multilevel stochastic collocation method for SPDEs

We present a multilevel stochastic collocation method that, as do multilevel Monte Carlo methods, uses a hierarchy of spatial approximations to reduce the overall computational complexity when solving partial differential equations with random inputs. For approximation in parameter space, a hierarchy of multi-dimensional interpolants of increasing fidelity are used. Rigorous convergence and computational cost estimates for the new multilevel stochastic collocation method are derived and used to demonstrate its advantages compared to standard single-level stochastic collocation approximations as well as multilevel Monte Carlo methods.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Average M shell fluorescence yields for elements with 70≤Z≤92

The theoretical, experimental and analytical methods for the calculation of average M-shell fluorescence yield (ω{sup ¯}{sub M}) of different elements are very important because of the large number of their applications in various areas of physical chemistry and medical research. In this paper, the bulk of the average M-shell fluorescence yield measurements reported in the literature, covering the period 1955 to 2005 are interpolated by using an analytical function to deduce the empirical average M-shell fluorescence yield in the atomic range of 70≤Z≤92. The results were compared with the theoretical and fitted values reported by other authors. Reasonable agreement was typically obtained between our result and other works.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗