Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “function approximation”

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,693 records · Page 94
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

THE MASS DISTRIBUTION FUNCTION OF PLANETS

The distribution of orbital period ratios of adjacent planets in extrasolar planetary systems discovered by the Kepler space telescope exhibits a peak near ∼1.5–2, a long tail of larger period ratios, and a steep drop-off in the number of systems with period ratios below ∼1.5. We find from these data that the dimensionless orbital separations have an approximately log-normal distribution. Using Hill’s criterion for the dynamical stability of two planets, we find an upper bound on planet masses such that the most common planet mass does not exceed 10{sup −3.2}m{sub ∗}, or about two-thirds of Jupiter’s mass for solar-mass stars. Assuming that the mass ratio and the dynamical separation (orbital spacings in units of mutual Hill radius) of adjacent planets are independent random variates, and adopting empirical distributions for these, we use Hill’s criterion in a statistical way to estimate the planet mass distribution function from the observed distribution of orbital separations. We find that the planet mass function is peaked in logarithm of mass, with a peak value and standard deviation of logm/M{sub ⊕} of ∼(0.6−1.0) and ∼(1.1−1.2), respectively.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Paraxial, Thin-Lens Analysis of Fixed-Tune, Non-Scaling FFAs with Two Magnets per Cell

Fixed-field accelerators (FFAs) have closed orbits that change as a function of beam momentum. It is sometimes useful to avoid various resonances by keeping the tunes constant as these orbits change. A well-known example is the scaling FFA where the entire beam orbit and optics are geometrically scaled as a function of momentum. However, this is stricter than necessary: cells with three or more lenses can have fixed tunes even when the focussing strengths in the lenses change. Recently, Dejan Trbojevic has found a pair of nonlinear magnets that produce fixed tunes and fixed beta functions while not being a scaling FFA. This improves on a more approximate solution in. Notably a scaling FFA would require one magnet to be entirely reverse-bending, whereas the Trbojevic solution does not do this and instead resembles an intermediate point between the traditional scaling field profiles and the nonscaling profiles centred around a momentum with equal and positive fields. This suggests there are at least three levels of stringency that one can apply to a fixed-tune FFA design: 1. Fixed cell tunes (in both planes) as a function of momentum; 2. Fixed optics (beta functions) as a function of momentum; 3. Similarity of all orbits via a scaling symmetry law ↔ traditional scaling FFA. This note studies the interesting case #2 above (#3 being fully characterised by the orbit at a single energy) in the simplest possible example: a cell of two thin lenses in the small angle (paraxial) approximation.

43 PARTICLE ACCELERATORS↗
DOE OSTI2021

Sequential Learning of Active Subspaces

This study shows that in recent years, active subspace methods (ASMs) have become a popular means of performing subspace sensitivity analysis on black-box functions. Naively applied, however, ASMs require gradient evaluations of the target function. In the event of noisy, expensive, or stochastic simulators, evaluating gradients via finite differencing may be infeasible. In such cases, often a surrogate model is employed, on which finite differencing is performed. When the surrogate model is a Gaussian process (GP), we show that the ASM estimator is available in closed form, rendering the finite-difference approximation unnecessary. We use our closed-form solution to develop acquisition functions focused on sequential learning tailored to sensitivity analysis on top of ASMs. We also show that the traditional ASM estimator may be viewed as a method of moments estimator for a certain class of GPs. We demonstrate how uncertainty on GP hyperparameters may be propagated to uncertainty on the sensitivity analysis, allowing model-based confidence intervals on the active subspace. Our methodological developments are illustrated on several examples.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2015

Non-local bias in the halo bispectrum with primordial non-Gaussianity

Primordial non-Gaussianity can lead to a scale-dependent bias in the density of collapsed halos relative to the underlying matter density. The galaxy power spectrum already provides constraints on local-type primordial non-Gaussianity complementary those from the cosmic microwave background (CMB), while the bispectrum contains additional shape information and has the potential to outperform CMB constraints in future. We develop the bias model for the halo density contrast in the presence of local-type primordial non-Gaussianity, deriving a bivariate expansion up to second order in terms of the local linear matter density contrast and the local gravitational potential in Lagrangian coordinates. Nonlinear evolution of the matter density introduces a non-local tidal term in the halo model. Furthermore, the presence of local-type non-Gaussianity in the Lagrangian frame leads to a novel non-local convective term in the Eulerian frame, that is proportional to the displacement field when going beyond the spherical collapse approximation. We use an extended Press-Schechter approach to evaluate the halo mass function and thus the halo bispectrum. We show that including these non-local terms in the halo bispectra can lead to corrections of up to 25% for some configurations, on large scales or at high redshift.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Investigating the minimum post-burst energy required to function an exploding bridgewire detonator

The minimum energy required after bridge-burst to function a commercially available RP-1 exploding bridgewire detonator with a gold bridgewire was investigated using two capacitor discharge units (CDUs) making use of state-of-the-art on/off high-voltage high-current switches. Both switches are capable of shutting off the current flow in ≈100ns and the conduction time was varied systematically to limit the duration of electrical energy delivered from the CDU until a detonation go/no-go threshold was reached. It was found that the minimal extra electrical energy required for a soft-fire (a statistically marginal go) was approximately half that required for a hard-fire (a definite go). Furthermore, the measured energy required to function the RP-1 detonator from the 11.5 kV CDU was approximately 20 mJ greater than for the 600 V CDU, and it is postulated that this was due to arcing across the detonator terminals and an associated current shunt that was not efficient at depositing energy conducive to detonator functioning. Measurements of the bridgewire temperatures bursting in air showed that the maximum temperature reached was consistent (≈18kK) when sufficient energy was delivered that the detonator functioned (either soft- or hard-fire), but that the temperature for a sub-threshold bridgewire burst was significantly cooler at approximately 8 kK. Furthermore, it was demonstrated that the intensity of the light emitted at burst scaled with deposited energy. A related finding is that the light emission from the bridgewire peaks some tens of nanoseconds after both peak voltage and peak electrical power.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
DOE OSTI2015

Beyond Atomic Sizes and Hume-Rothery Rules: Understanding and Predicting High-Entropy Alloys

High-entropy alloys constitute a new class of materials that provide an excellent combination of strength, ductility, thermal stability, and oxidation resistance. Although they have attracted extensive attention due to their potential applications, little is known about why these compounds are stable or how to predict which combination of elements will form a single phase. Here, we present a review of the latest research done on these alloys focusing on the theoretical models devised during the last decade. We discuss semiempirical methods based on the Hume-Rothery rules and stability criteria based on enthalpies of mixing and size mismatch. To provide insights into the electronic and magnetic properties of high-entropy alloys, we show the results of first-principles calculations of the electronic structure of the disordered solid-solution phase based on both Korringa Kohn Rostoker coherent potential approximation and large supercell models of example face-centered cubic and body-centered cubic systems. Furthermore, we discuss in detail a model based on enthalpy considerations that can predict which elemental combinations are most likely to form a single-phase high-entropy alloy. The enthalpies are evaluated via first-principles high-throughput density functional theory calculations of the energies of formation of binary compounds, and therefore it requires no experimental or empirically derived input. Finally, the model correctly accounts for the specific combinations of metallic elements that are known to form single-phase alloys while rejecting similar combinations that have been tried and shown not to be single phase.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Two-point sum-rules in three-dimensional Yang-Mills theory

We compute the stress-tensor two-point function in three-dimensional Yang-Mills theory to three-loops in perturbation theory. Using its calculable shape at high momenta, we test the notion that its Borel transform is saturated at low energies by the lowest glueball state(s). This assumption provides relatively stable estimates for the mass of the lightest glueball that we compare with lattice simulations. We also provide estimates for the coupling of the lightest glueball to the stress tensor. Along the way, we comment on the extent that such estimates are non-rigorous. Lastly, we discuss the possibility of applying the sum-rule analysis to two-point functions of higher-spin operators and obtain a crude approximation for the glueball couplings to these operators.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2016

Reflections on conformal spectra

Here, we use modular invariance and crossing symmetry of conformal field theory to reveal approximate reflection symmetries in the spectral decompositions of the partition function in two dimensions in the limit of large central charge and of the four-point function in any dimension in the limit of large scaling dimensions Δ 0 of external operators. We use these symmetries to motivate universal upper bounds on the spectrum and the operator product expansion coefficients, which we then derive by independent techniques. Some of the bounds for four-point functions are valid for finite Δ 0 as well as for large Δ 0 . We discuss a similar symmetry in a large spacetime dimension limit. Finally, we comment on the analogue of the Cardy formula and sparse light spectrum condition for the four-point function.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2019

Theoretical Model Study of Interplay of Coulomb Interaction and Electron-Phonon Interaction in the Thermal Properties of Monolayer Graphene

We propose here a tight-binding (TB) model Hamiltonian for monolayer graphene-on-substrate describing the nearest-neighbor-hopping, on-site Coulomb interaction on the sub-lattices and the electron-phonon interaction under the high-frequency limit of phonon vibration. Applying Lang-Firsov canonical transformation, the electron and phonon systems are decoupled in the atomic Hamiltonian, such that the effective Coulomb interaction and effective nearest-neighbor-hopping integral respectively appear as U-tilde =U−2t{sub 1}λ and t-tilde {sub 1}=t{sub 1}e{sup ((−t{sub 1}λ)/(ω{sub 0}))}, where U, t{sub 1}, λ and ω{sub 0} are respectively Coulomb energy, nearest-neighbor-hopping integral, electron-phonon (e-ph) coupling and phonon frequency. The effective Coulomb interaction in the Hamiltonian is considered within mean-field approximation. The Hamiltonian is solved by Zubarev’s Green’s function technique. The temperature-dependent electronic entropy and specific heat are calculated from the free energy of graphene system and are computed numerically. The temperature-dependent electronic specific heat exhibits a charge gap peak at room temperature arising due to the effect of Coulomb interaction and electron-phonon interaction. The evolution of these peaks in specific heat is investigated by varying the model parameters of the system.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2019

Nonlinear dynamics from linear quantum evolutions

Linear dynamics restricted to invariant submanifolds generally gives rise to nonlinear dynamics. Submanifolds in the quantum framework may emerge for several reasons: one could be interested in specific properties possessed by a given family of states, either as a consequence of experimental constraints or inside an approximation scheme. In this work we investigate such issues in connection with a one parameter group ϕ{sub t} of transformations on a Hilbert space, H, defining the unitary evolutions of a chosen quantum system. Two procedures will be presented: the first one consists in the restriction of the vector field associated with the Schrödinger equation to a submanifold invariant under the flow ϕ{sub t}. The second one makes use of the Lagrangian formalism and can be extended also to non-invariant submanifolds, even if in such a case the resulting dynamics is only an approximation of the flow ϕ{sub t}. Such a result, therefore, should be conceived as a generalization of the variational method already employed for stationary problems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Using hybrid implicit Monte Carlo diffusion to simulate gray radiation hydrodynamics

This work describes how to couple a hybrid Implicit Monte Carlo Diffusion (HIMCD) method with a Lagrangian hydrodynamics code to evaluate the coupled radiation hydrodynamics equations. This HIMCD method dynamically applies Implicit Monte Carlo Diffusion (IMD) [1] to regions of a problem that are opaque and diffusive while applying standard Implicit Monte Carlo (IMC) [2] to regions where the diffusion approximation is invalid. We show that this method significantly improves the computational efficiency as compared to a standard IMC/Hydrodynamics solver, when optically thick diffusive material is present, while maintaining accuracy. Two test cases are used to demonstrate the accuracy and performance of HIMCD as compared to IMC and IMD. The first is the Lowrie semi-analytic diffusive shock [3]. The second is a simple test case where the source radiation streams through optically thin material and heats a thick diffusive region of material causing it to rapidly expand. We found that HIMCD proves to be accurate, robust, and computationally efficient for these test problems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

A model and variance reduction method for computing statistical outputs of stochastic elliptic partial differential equations

We present a model and variance reduction method for the fast and reliable computation of statistical outputs of stochastic elliptic partial differential equations. Our method consists of three main ingredients: (1) the hybridizable discontinuous Galerkin (HDG) discretization of elliptic partial differential equations (PDEs), which allows us to obtain high-order accurate solutions of the governing PDE; (2) the reduced basis method for a new HDG discretization of the underlying PDE to enable real-time solution of the parameterized PDE in the presence of stochastic parameters; and (3) a multilevel variance reduction method that exploits the statistical correlation among the different reduced basis approximations and the high-fidelity HDG discretization to accelerate the convergence of the Monte Carlo simulations. The multilevel variance reduction method provides efficient computation of the statistical outputs by shifting most of the computational burden from the high-fidelity HDG approximation to the reduced basis approximations. Furthermore, we develop a posteriori error estimates for our approximations of the statistical outputs. Based on these error estimates, we propose an algorithm for optimally choosing both the dimensions of the reduced basis approximations and the sizes of Monte Carlo samples to achieve a given error tolerance. We provide numerical examples to demonstrate the performance of the proposed method.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Parameterization of the Stoner-Wohlfarth model of magnetic hysteresis

The Stoner-Wohlfarth (SW) is the most used model of magnetic hysteresis, but its computation is time-consuming. We approximate piecewise this model by easy-to-compute analytic functions. Overall, our parametrization is suitable for fast quantitative evaluations and fitting experimental data, which we exemplify. Applicability of this parametrization is as broad as that of the SW model, which can be extended to include more parameters.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Toward Quantum Chemical Free Energy Simulations of Platinum Nanoparticles on Titania Support

Platinum nanoparticles (Pt-NPs) supported on titania surfaces are costly but indispensable heterogeneous catalysts because of their highly effective and selective catalytic properties. Therefore, it is vital to understand their physicochemical processes during catalysis to optimize their use and to further develop better catalysts. However, simulating these dynamic processes is challenging due to the need for a reliable quantum chemical method to describe chemical bond breaking and bond formation during the processes but, at the same time, fast enough to sample a large number of configurations required to compute the corresponding free energy surfaces. Density functional theory (DFT) is often used to explore Pt-NPs; nonetheless, it is usually limited to some minimum-energy reaction pathways on static potential energy surfaces because of its high computational cost. In this work, we report a combination of the density functional tight binding (DFTB) method as a fast but reliable approximation to DFT, the steered molecular dynamics (SMD) technique, and the Jarzynski equality to construct free energy surfaces of the temperature-dependent diffusion and growth of platinum particles on a titania surface. In particular, we present the parametrization for Pt-X (X = Pt, Ti, or O) interactions in the framework of the second-order DFTB method, using a previous parametrization for titania as a basis. The optimized parameter set was used to simulate the surface diffusion of a single platinum atom (Pt 1 ) and the growth of Pt 6 from Pt 5 and Pt 1 on the rutile (110) surface at three different temperatures (T = 400, 600, 800 K). The free energy profile was constructed by using over a hundred SMD trajectories for each process. We found that increasing the temperature has a minimal effect on the formation free energy; nevertheless, it significantly reduces the free energy barrier of Pt atom migration on the TiO 2 surface and the transition state (TS) of its deposition. In a concluding remark, the methodology opens the pathway to quantum chemical free energy simulations of Pt-NPs’ temperature-dependent growth and other transformation processes on the titania support.

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

Uncloaking the thermodynamics of the studtite to metastudtite shear-induced transformation

The interplay between thermodynamics and mechanical properties in the transformation of studtite, (UO 2 )(O 2 )(H 2 O) 2 ·2H 2 O, into metastudtite, (UO 2 )(O 2 )(H 2 O) 2 , two important corrosion phases observed on the surface of uranium dioxide exposed to water, is revealed using density functional perturbation theory. Phonon calculations within the quasi-harmonic approximation predict that the standard entropy change for the (UO 2 )(O 2 )(H 2 O) 2 ·2H 2 O → (UO 2 )(O 2 )(H 2 O) 2 + 2H 2 O reaction is ΔS 0 = +80 J·mol –1 ·K –1 for the production of water in the liquid state and +389 J·mol–1·K–1 for water vapor. Similar to bulk H 2 O(l), the bulk modulus of (UO 2 )(O 2 )(H 2 O) 2 ·2H 2 O increases with temperature, contrasting with (UO 2 )(O 2 )(H 2 O) 2 which features the typical Anderson–Gruneisen temperature dependence of oxide solids. Upon removal of interstitial H 2 O in studtite, the most important changes in the shear modulus, the parameter limiting the mechanical stability, arise in the planes normal to chain propagation directions. Lastly, the present findings have important implications for the dehydration of other hygroscopic materials.

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

Stabilizing lattice oxygen redox in layered sodium transition metal oxide through spin singlet state

Abstract Reversible lattice oxygen redox reactions offer the potential to enhance energy density and lower battery cathode costs. However, their widespread adoption faces obstacles like substantial voltage hysteresis and poor stability. The current research addresses these challenges by achieving a non-hysteresis, long-term stable oxygen redox reaction in the P3-type Na 2/3 Cu 1/3 Mn 2/3 O 2 . Here we show this is accomplished by forming spin singlet states during charge and discharge. Detailed analysis, including in-situ X-ray diffraction, shows highly reversible structural changes during cycling. In addition, local CuO 6 Jahn-Teller distortions persist throughout, with dynamic Cu-O bond length variations. In-situ hard X-ray absorption and ex-situ soft X-ray absorption study, along with density function theory calculations, reveal two distinct charge compensation mechanisms at approximately 3.66 V and 3.99 V plateaus. Notably, we observe a Zhang-Rice-like singlet state during 3.99 V charging, offering an alternative charge compensation mechanism to stabilize the active oxygen redox reaction.

25 ENERGY STORAGE↗
DOE OSTI2024

Structural diversity and clustering of bacterial flagellar outer domains

Supercoiled flagellar filaments function as mechanical propellers within the bacterial flagellum complex, playing a crucial role in motility. Flagellin, the building block of the filament, features a conserved inner D0/D1 core domain across different bacterial species. In contrast, approximately half of the flagellins possess additional, highly divergent outer domain(s), suggesting varied functional potential. In this study, we report atomic structures of flagellar filaments from three distinct bacterial species: Cupriavidus gilardii , Stenotrophomonas maltophilia , and Geovibrio thiophilus . Our findings reveal that the flagella from the facultative anaerobic G. thiophilus possesses a significantly more negatively charged surface, potentially enabling adhesion to positively charged minerals. Furthermore, we analyze all AlphaFold predicted structures for annotated bacterial flagellins, categorizing the flagellin outer domains into 682 structural clusters. This classification provides insights into the prevalence and experimental verification of these outer domains. Remarkably, two of the flagellar structures reported herein belong to a distinct cluster, indicating additional opportunities on the study of the functional diversity of flagellar outer domains. Our findings underscore the complexity of bacterial flagellins and open up possibilities for future studies into their varied roles beyond motility.

Science & Technology - Other Topics↗
DOE OSTI2015

QTAIM studies of [Li(DMSO){sub n}]{sup +} and [Al(DMSO){sub n}]{sup 3+}

Geometry optimization of [Li(DMSO){sub n}]{sup +} and [Al(DMSO){sub n}]{sup 3+} complexes with n = 1 - 6 in DMSO (dimethyl sulfoxide, (CH{sub 3}){sub 2}SO) solutions is performed using DFT treatment with B3LYP hybrid functional and cc-pVDZ basis sets. Solvent effects are approximated within Integral Equation Formalism Polarisable Continuum Model. Electron structure of individual model systems is investigated in terms of Quantum Theory of Atoms-in-Molecule topological analysis of electron density. Metal-DMSO bonding through oxygen atom is preferred. The most probable coordination numbers n = 4 for Li{sup +} and n = 6 for Al{sup 3+} complexes with DMSO are concluded.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗