Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Statistical Physics”

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,459 records · Page 81
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

Environment-assisted non-Markovian speedup dynamics control

For the damped Jaynes–Cummings model in the weak system–environment coupling regime, a scheme to manipulate the non-Markovian dynamics has been previously proposed in the controllable environment. We here mainly focus on the two-Lorentzian model and the band-gap model, which are composed by adding a new Lorentzian spectrum to the original single-Lorentzian spectrum. By manipulating the parameters of the new added Lorentzian spectrum, i.e., changing the weights of the two Lorentzians, and tuning the spectral width of the new second Lorentzian, the dynamical crossover from Markovian to non-Markovian is achieved in the weak system–environment coupling regime. In addition, according to the controllable non-Markovian dynamics behavior, the roles of the parameters of the environment on the speedup evolution of the system are discussed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

A method for determining optimal electric range by considering electric vehicle lightweighting on perceived ownership cost

The limited driving range due to high costs and low energy densities of batteries constraints the battery electric vehicle (BEV) market growth. Lightweighting in theory can reduce energy consumption rate and extend the driving range. The knowledge gap is to quantitatively link the cost-effectiveness of light-weight technologies with range extension and consumer acceptance of BEVs. In this study, a physics-based energy consumption model of BEVs is constructed and associated with a statistics-based model on the basis of travel surveys. A perceived cost of ownership (PCO) is then developed by adding intangible costs to traditional total cost of ownership models. We estimate, at the disaggregate vehicle model and driver level and the aggregate market level, 1) the extended range due to lightweighting for a given battery size; and 2) the optimal electric range based on lightweighting decisions. The cost-effectiveness of lightweighting for BEV range extension is found to vary with income-dependent daily range limitation value, driving patterns and lightweighting technology costs. In general, adopting lightweighting in BEVs is more cost-effective for consumers with higher daily limitation value, as well as for those with higher driving intensity or suitable daily driving patterns. When the lightweighting involves a higher vehicle production cost, less lightweighting could reduce the overall PCO for BEV owners. 4 of the selected top ten BEV models are found to benefit from additional 2.09%–4.45% lightweighting. Finally, the method built in this study can guide automakers in planning R&D investments in battery and lightweighting technologies.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2015

Adaptive error covariances estimation methods for ensemble Kalman filters

This paper presents a computationally fast algorithm for estimating, both, the system and observation noise covariances of nonlinear dynamics, that can be used in an ensemble Kalman filtering framework. The new method is a modification of Belanger's recursive method, to avoid an expensive computational cost in inverting error covariance matrices of product of innovation processes of different lags when the number of observations becomes large. When we use only product of innovation processes up to one-lag, the computational cost is indeed comparable to a recently proposed method by Berry–Sauer's. However, our method is more flexible since it allows for using information from product of innovation processes of more than one-lag. Extensive numerical comparisons between the proposed method and both the original Belanger's and Berry–Sauer's schemes are shown in various examples, ranging from low-dimensional linear and nonlinear systems of SDEs and 40-dimensional stochastically forced Lorenz-96 model. Our numerical results suggest that the proposed scheme is as accurate as the original Belanger's scheme on low-dimensional problems and has a wider range of more accurate estimates compared to Berry–Sauer's method on L-96 example.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

Current-Induced Forces for Nonadiabatic Molecular Dynamics

We present general first-principles derivation of the expression for current-induced forces. The expression is applicable in nonequilibrium molecular systems with arbitrary intramolecular interactions and for any electron-nuclei coupling. It provides a controlled consistent way to account for quantum effects of nuclear motion around a classical trajectory, accounts for electronic non-Markov character of the friction tensor, and suggests a method for treatment beyond strictly adiabatic approximation. Lastly, we show the connection of the expression with previous studies and discuss effective ways to evaluate the friction tensor.

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

Improving the Quasi‐Biennial Oscillation via a Surrogate‐Accelerated Multi‐Objective Optimization

Accurate simulation of the quasi-biennial oscillation (QBO) is challenging due to uncertainties in representing convectively generated gravity waves. We develop an end-to-end uncertainty quantification workflow that calibrates these gravity wave processes in E3SM for a realistic QBO. Central to our approach is a domain knowledge-informed, compressed representation of high-dimensional spatio-temporal wind fields. By employing a parsimonious statistical model that learns the fundamental frequency from complex observations, we extract interpretable and physically meaningful quantities capturing key attributes. Building on this, we train a probabilistic surrogate model that approximates the fundamental characteristics of the QBO as functions of critical physics parameters governing gravity wave generation. Leveraging the Karhunen–Loève decomposition, our surrogate efficiently represents these characteristics as a set of orthogonal features, capturing cross-correlations among multiple physics quantities evaluated at different pressure levels and enabling rapid surrogate-based inference at a fraction of the computational cost of full-scale simulations. Finally, we analyze the inverse problem using a multi-objective approach. Our study reveals a tension between amplitude and period that constrains the QBO representation, precluding a single optimal solution. To navigate this, we quantify the bi-criteria trade-off and generate a set of Pareto optimal parameter values that balance the conflicting objectives. This integrated workflow improves the fidelity of QBO simulations and offers a versatile template for uncertainty quantification in complex geophysical models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Uptake of N 2 O 5 by aqueous aerosol unveiled using chemically accurate many-body potentials

The reactive uptake of N 2 O 5 to aqueous aerosol is a major loss channel for nitrogen oxides in the troposphere. Despite its importance, a quantitative picture of the uptake mechanism is missing. Here we use molecular dynamics simulations with a data-driven many-body model of coupled-cluster accuracy to quantify thermodynamics and kinetics of solvation and adsorption of N 2 O 5 in water. The free energy profile highlights that N 2 O 5 is selectively adsorbed to the liquid–vapor interface and weakly solvated. Accommodation into bulk water occurs slowly, competing with evaporation upon adsorption from gas phase. Leveraging the quantitative accuracy of the model, we parameterize and solve a reaction–diffusion equation to determine hydrolysis rates consistent with experimental observations. We find a short reaction–diffusion length, indicating that the uptake is dominated by interfacial features. The parameters deduced here, including solubility, accommodation coefficient, and hydrolysis rate, afford a foundation for which to consider the reactive loss of N 2 O 5 in more complex solutions.

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

Coupled Langmuir oscillations in 2-dimensional quantum plasmas

In this work, we present a hydrodynamic model to study the coupled quantum electron plasma oscillations (QEPO) for two dimensional (2D) degenerate plasmas, which incorporates all the essential quantum ingredients such as the statistical degeneracy pressure, electron-exchange, and electron quantum diffraction effect. Effects of diverse physical aspects like the electronic band-dispersion effect, the electron exchange-correlations and the quantum Bohm-potential as well as other important plasma parameters such as the coupling parameter (plasma separation) and the plasma electron number-densities on the linear response of the coupled system are investigated. By studying three different 2D plasma coupling types, namely, graphene-graphene, graphene-metalfilm, and metalfilm-metalfilm coupling configurations, it is remarked that the collective quantum effects can influence the coupled modes quite differently, depending on the type of the plasma configuration. It is also found that the slow and fast QEPO frequency modes respond very differently to the change in plasma parameters. Current findings can help in understanding of the coupled density oscillations in multilayer graphene, graphene-based heterojunctions, or nanofabricated integrated circuits.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2014

Relativistic self-focusing of ultra-high intensity X-ray laser beams in warm quantum plasma with upward density profile

The results of a numerical study of high-intensity X-ray laser beam interaction with warm quantum plasma (WQP) are presented. By means of an upward ramp density profile combined with quantum factors specially the Fermi velocity, we have demonstrated significant relativistic self-focusing (RSF) of a Gaussian electromagnetic beam in the WQP where the Fermi temperature term in the dielectric function is important. For this purpose, we have considered the quantum hydrodynamics model that modifies refractive index of inhomogeneous WQPs with the inclusion of quantum correction through the quantum statistical and diffraction effects in the relativistic regime. Also, to better illustration of the physical difference between warm and cold quantum plasmas and their effect on the RSF, we have derived the envelope equation governing the spot size of X-ray laser beam in Q-plasmas. In addition to the upward ramp density profile, we have found that the quantum effects would be caused much higher oscillation and better focusing of X-ray laser beam in the WQP compared to that of cold quantum case. Our computational results reveal the importance of the use of electrons density profile and Fermi speed in enhancing self-focusing of laser beam.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

Effect of secondary emission on the argon plasma afterglow with large dust density

A zero-dimensional, space-averaged model for argon plasma afterglow with large dust density is developed. In the model, three groups of electrons in the plasma afterglow are assumed: (i) thermal electrons with Maxwellian distribution, (ii) energetic electrons generated by metastable-metastable collisions (metastable pooling), and (iii) secondary electrons generated at collisions of ions with the electrodes, which have sufficiently large negative voltages in the afterglow. The model calculates the time-dependencies for electron densities in plasma afterglow based on experimental decay times for metastable density and electrode bias. The effect of secondary emission on electron density in the afterglow is estimated by varying secondary emission yields. It is found that this effect is less important than metastable pooling. The case of dust-free plasma afterglow is considered also, and it is found that in the afterglow the effect of secondary emission may be more important than metastable pooling. The secondary emission may increase thermal electron density n{sub e} in dust-free and dusty plasma afterglows on a few ten percentages. The calculated time dependencies for n{sub e} in dust-free and dusty plasma afterglows describe well the experimental results.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

Formation of collisional sheath in electronegative plasma with two species of positive ions

Sheath formation is investigated for electronegative plasma in presence of two species of positive ions in collisional environment. The gas under consideration is a mixture of oxygen and argon. Argon is the considered as having fixed volume and impact of collision is studied with increasing pressure of oxygen. Fluid equations are solved for three species namely, the two positive ions and a negative ion. Electrons are considered to follow Boltzmann distribution. Collision is modeled by constant mean free path model and has been used as a parameter. It has been found that collision enhances the sheath formation. The negative ion core is nearly unaffected by the presence of collision and is governed by the electric potential. The negative flux field is, however, affected by the presence of collision and shows a steady behavior in front of the wall. The two positive ions are heavily affected by the presence of collision and the modeling is such that their equilibrium densities can be estimated by solving simultaneous cubic equations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

An analytic formula for the relativistic incoherent Thomson backscattering spectrum for a drifting bi-Maxwellian plasma

An analytic formula has been derived for the relativistic incoherent Thomson backscattering spectrum for a drifting anisotropic plasma when the scattering vector is parallel to the drifting direction. The shape of the scattering spectrum is insensitive to the electron temperature perpendicular to the scattering vector, but its amplitude may be modulated. As a result, while the measured temperature correctly represents the electron distribution parallel to the scattering vector, the electron density may be underestimated when the perpendicular temperature is higher than the parallel temperature. Since the scattering spectrum in shorter wavelengths is greatly enhanced by the existence of drift, the diagnostics might be used to measure local electron current density in fusion plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

Ionic liquids—Conventional solvent mixtures, structurally different but dynamically similar

In more ways than one, pure ionic liquids (ILs) can be seen as mixtures. By definition they are comprised of cationic and anionic components and they also possess dual charge and apolar characteristics. We recently uncovered interesting dynamical behavior [Araque et al., J. Phys. Chem. B 119(23), 7015 (2015)] that can be ascribed to this duality. For small neutral solutes local friction can be high in certain regions and much lower in others. It is only reasonable to ask whether this interesting behavior is unique to ILs or is also common in certain conventional solvent mixtures such as dimethylsufoxide/glycerol for which the viscosity can be tuned to be similar. We make the case that the latter scenario is correct and that whereas viscous conventional solvent mixtures are structurally very different from ILs, dynamically they are not. From the perspective of a solute that is small, both ILs and viscous conventional solvent mixtures display frictionally stiff and soft regions associated with cage and jump diffusive regimes. In the case of ILs these are associated with charge-enhanced and charge-depleted liquid regions, whereas in the case of the conventional solvents by the distinct frictional properties of the two components.

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

Determination of deuterium–tritium critical burn-up parameter by four temperature theory

Conditions for thermonuclear burn-up of an equimolar mixture of deuterium-tritium in non-equilibrium plasma have been investigated by four temperature theory. The photon distribution shape significantly affects the nature of thermonuclear burn. In three temperature model, the photon distribution is Planckian but in four temperature theory the photon distribution has a pure Planck form below a certain cut-off energy and then for photon energy above this cut-off energy makes a transition to Bose-Einstein distribution with a finite chemical potential. The objective was to develop four temperature theory in a plasma to calculate the critical burn up parameter which depends upon initial density, the plasma components initial temperatures, and hot spot size. All the obtained results from four temperature theory model are compared with 3 temperature model. It is shown that the values of critical burn-up parameter calculated by four temperature theory are smaller than those of three temperature model.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

Dynamics of photoexcited Ba{sup +} cations in {sup 4}He nanodroplets

We present a joint experimental and theoretical study on the desolvation of Ba{sup +} cations in {sup 4}He nanodroplets excited via the 6p ← 6s transition. The experiments reveal an efficient desolvation process yielding mainly bare Ba{sup +} cations and Ba{sup +}He{sub n} exciplexes with n = 1 and 2. The speed distributions of the ions are well described by Maxwell-Boltzmann distributions with temperatures ranging from 60 to 178 K depending on the excitation frequency and Ba{sup +} He{sub n} exciplex size. These results have been analyzed by calculations based on a time-dependent density functional description for the helium droplet combined with classical dynamics for the Ba{sup +}. In agreement with experiment, the calculations reveal the dynamical formation of exciplexes following excitation of the Ba{sup +} cation. In contrast to experimental observation, the calculations do not reveal desolvation of excited Ba{sup +} cations or exciplexes, even when relaxation pathways to lower lying states are included.

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

Structures of glass-forming liquids by x-ray scattering: Glycerol, xylitol, and D-sorbitol

In this work, synchrotron X-ray scattering has been used to investigate three liquid polyalcohols of different sizes (glycerol, xylitol and D-sorbitol) from above the glass transition temperatures $T_g$ to below. We focus on two structural orders: the association of the polar OH groups by hydrogen bonds (HBs) and the packing of the non-polar hydrocarbon groups. We find that the two structural orders evolve very differently reflecting the different natures of bonding. Upon cooling from 400 K, the O$\cdots$O correlation at 2.8 Å increases significantly in all three systems, indicating more HBs, until kinetic arrest at $T_g$; the increase is well described by an equilibrium between bonded and non-bonded OH with ΔH = 9.1 kJ/mol and ΔS = 13.4 J/mol/K. When heated above $T_g$, glycerol loses the fewest HBs per OH, followed by xylitol and by D-sorbitol, in the same order the number of OH groups per molecule increases (3, 5, and 6). The pair correlation functions of all three liquids show exponentially damped density modulations of wavelength 4.5 Å, which are associated with the main scattering peak and with the intermolecular C$\cdots$C correlation. In this respect, glycerol is the most ordered with the most persistent density ripples, followed by D-sorbitol and by xylitol. Heating above $T_g$ causes faster damping of the density ripples with the rate of change being the slowest in xylitol, followed by glycerol and by D-sorbitol. Given the different dynamic fragility of the three liquids (glycerol being the strongest and D-sorbitol the most fragile), we relate our results to the current theories of the structural origin for the difference. We find that the fragility difference is better understood on the basis of the thermal stability of HB clusters than that of the structure associated with the main scattering peak.

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

Modified Fermi’s golden rule rate expressions

Fermi’s golden rule (FGR) serves as the basis for many expressions of spectroscopic observables and quantum transition rates. The utility of FGR has been demonstrated through decades of experimental confirmation. However, there still remain important cases where the evaluation of a FGR rate is ambiguous or ill-defined. Examples are cases where the rate has divergent terms due to the sparsity in the density of final states or time dependent fluctuations of system Hamiltonians. Strictly speaking, assumptions of FGR are no longer valid for such cases. However, it is still possible to define modified FGR rate expressions that are useful as effective rates. The resulting modified FGR rate expressions resolve a long standing ambiguity often encountered in using FGR and offer more reliable ways to model general rate processes. Finally, simple model calculations illustrate the utility and implications of new rate expressions.

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

CSI 2264: simultaneous optical and infrared light curves of young disk-bearing stars in NGC 2264 with CoRoT and Spitzer—evidence for multiple origins of variability

We present the Coordinated Synoptic Investigation of NGC 2264, a continuous 30 day multi-wavelength photometric monitoring campaign on more than 1000 young cluster members using 16 telescopes. The unprecedented combination of multi-wavelength, high-precision, high-cadence, and long-duration data opens a new window into the time domain behavior of young stellar objects. Here we provide an overview of the observations, focusing on results from Spitzer and CoRoT. The highlight of this work is detailed analysis of 162 classical T Tauri stars for which we can probe optical and mid-infrared flux variations to 1% amplitudes and sub-hour timescales. We present a morphological variability census and then use metrics of periodicity, stochasticity, and symmetry to statistically separate the light curves into seven distinct classes, which we suggest represent different physical processes and geometric effects. We provide distributions of the characteristic timescales and amplitudes and assess the fractional representation within each class. The largest category (>20%) are optical 'dippers' with discrete fading events lasting ∼1-5 days. The degree of correlation between the optical and infrared light curves is positive but weak; notably, the independently assigned optical and infrared morphology classes tend to be different for the same object. Assessment of flux variation behavior with respect to (circum)stellar properties reveals correlations of variability parameters with Hα emission and with effective temperature. Overall, our results point to multiple origins of young star variability, including circumstellar obscuration events, hot spots on the star and/or disk, accretion bursts, and rapid structural changes in the inner disk.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

SELF-CONSISTENT ION CYCLOTRON ANISOTROPY-BETA RELATION FOR SOLAR WIND PROTONS

We derive a set of self-consistent marginally stable states for a system of ion-cyclotron waves propagating parallel to the large-scale magnetic field through a homogeneous proton-electron plasma. The proton distributions and the wave dispersions are related through the condition that no further ion-cyclotron resonant particle scattering or wave growth/damping may take place. The thermal anisotropy of the protons in these states therefore defines the threshold value for triggering the proton-cyclotron anisotropy instability. A number of recent papers have noted that the anisotropy of solar wind protons at 1 AU does not seem to be limited by the proton-cyclotron anisotropy threshold, even at low plasma beta. However, this puzzle seems to be due solely to the estimation of this anisotropy threshold under the assumption that the protons have a bi-Maxwellian distribution. We note that bi-Maxwellian distributions are never marginally stable to the resonant cyclotron interaction, so these estimates do not represent physically valid thresholds. The threshold anisotropies obtained from our marginally stable states are much larger, as a function of proton parallel beta, than the bi-Maxwellian estimates, and we show that the measured data remains below these more rigorous thresholds. Thus, the results of this paper resolve the apparent contradiction presented by the solar wind anisotropy observations at 1 AU: the bi-Maxwellian anisotropies are not rigorous thresholds, and so do not limit the proton distributions in the solar wind.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗