Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Classical force fields”

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 235 records · Page 13

Strong-field QED limitations on TeV-class plasma wakefield accelerators

We demonstrate that quantum and classical radiation effects can become non-negligible for TeV-class beams propagating through plasma channels typical of staged plasma accelerators. Although the quantum nonlinearity parameter χ e remains small under currently envisioned experimental conditions, the cumulative influence of radiation over long acceleration distances can lead to significant modifications to the beam’s energy spread, emittance, and polarization. Our analytic models, validated by particle-in-cell simulations, highlight that for standard Gaussian beams, the orbit-induced energy spread dominates over quantum stochastic effects but can be mitigated by tailoring the beam profile, for example, through ring-shaped transverse distributions. In regimes where the radiation reaction approaches the accelerating force, the emittance may be cooled, forming distinctive ring-shaped phase-space structures. Finally, we analyze the influence of the radiation effect on spin transport inside the wakefield. These findings underscore the importance of considering both classical and quantum radiation dynamics in the design and optimization of future high-intensity plasma accelerators.

accelerator↗

A numerical model of nonmigrating diurnal tides between the surface and 65 km

Observations of surface pressure and middle atmosphere temperatures and winds indicate that a substantial nonmigrating component is present in the diurnal tide. The nonmigrating tides, which propagate with a zonal phase speed that is different from the earth's rotation, are attributed to the diurnal heating of geographically fixed sources. In this study we utilize a classical tidal model to examine the propagation characteristics of diurnal tides. The global fields of tropospheric sensible, radiative, and latent heating used to drive the model are supplied from summer and winter diurnal climatologies of the National Center for Atmospheric Research (NCAR) Community Climate Model (CCM2). A novel aspect of this study is the focus on the relative importance of the nonmigrating components. The classical model successfully reproduces many observed features of the low-latitude diurnal surface pressure tides. In the middle atmosphere, the simulated migrating (or sun-synchronous) tide shows qualitative agreement with November-March Limb Infrared Monitor of the Stratosphere (LIMS) observations. Tropospheric solar heating is clearly the dominant driving force for the migrating tide, with secondary contributions from boundary-layer sensible heating and tropospheric latent heat release. The leading modes of the zonal mean tide are also driven chiefly by tropospheric solar heating. The higher-order modes of the zonal mean and eastward propagating tides may be attributed to the joint effects of tropospheric solar heating, sensible heating, and latent heat release. The LIMS and other data reveal features that cannot be explained or examined within the context of the classical model used in the present study. These include upward phase propagation, vertical attenuation, and temporal variations in the migrating diurnal tide.

Lieberman, Ruth S.↗

Accurate ab initio quartic force fields for the ions HCO(+) and HOC(+)

The quartic force fields of HCO(+) and HOC(+) have been computed using augmented coupled cluster methods and basis sets of spdf and spdfg quality. Calculations on HCN, CO, and N2 have been performed to assist in calibrating the computed results. Going from an spdf to an spdfg basis shortens triple bonds by about 0.004 A, and increases the corresponding harmonic frequency by 10-20/cm, leaving bond distances about 0.003 A too long and triple bond stretching frequencies about 5/cm too low. Accurate estimates for the bond distances, fundamental frequencies, and thermochemical quantities are given. HOC(+) lies 37.8 +/- 0.5 kcal/mol (0 K) above HCO(+); the classical barrier height for proton exchange is 76.7 +/- 1.0 kcal/mol.

Martin, J. M. L.↗

A chain stretch-based gradient-enhanced model for damage and fracture in elastomers

Similar to quasi-brittle materials, it has been recently shown that elastomers can exhibit a macroscopically diffuse damage zone that accompanies the fracture process. In this study, we introduce a stretch-based gradient-enhanced damage (GED) model that allows the fracture to localize and also captures the development of a physically diffuse damage zone. This capability contrasts with the paradigm of the phase field method for fracture, where a sharp crack is numerically approximated in a diffuse manner. Capturing fracture localization and diffuse damage in our approach is achieved by considering nonlocal effects that encompass network topology, heterogeneity, and imperfections. These considerations motivate the use of a statistical damage function dependent upon the nonlocal deformation state. From this model, fracture toughness is realized as an output. While GED models have been classically utilized for damage modeling of structural engineering materials (e.g., concrete), they face challenges when trying to capture the cascade from damage to fracture, often leading to damage zone broadening (de Borst and Verhoosel, 2016). This deficiency contributed to the popularity of the phase-field method over the GED model for elastomers and other quasi-brittle materials. Other groups have proceeded with damage-based GED formulations that prove identical to the phase-field method (Lorentz et al., 2012), but these inherit the aforementioned limitations. To address this issue in a thermodynamically consistent framework, we implement two modeling features (a nonlocal driving force bound and a simple relaxation function) specifically designed to capture the evolution of a physically meaningful damage field and the simultaneous localization of fracture, thereby overcoming a longstanding obstacle in the development of these nonlocal strain- or stretch-based approaches. Here, we discuss several numerical examples to understand the features of the approach at the limit of incompressibility, and compare them to the phase-field method as a benchmark for the macroscopic response and fracture energy predictions.

Elastomers↗

Application of Quantum Computing to Biochemical Systems: A Look to the Future

Chemistry is considered as one of the more promising applications to science of near-term quantum computing. Recent work in transitioning classical algorithms to a quantum computer has led to great strides in improving quantum algorithms and illustrating their quantum advantage. Because of the limitations of near-term quantum computers, the most effective strategies split the work over classical and quantum computers. There is a proven set of methods in computational chemistry and materials physics that has used this same idea of splitting a complex physical system into parts that are treated at different levels of theory to obtain solutions for the complete physical system for which a brute force solution with a single method is not feasible. These methods are variously known as embedding, multi-scale, and fragment techniques and methods. We review these methods and then propose the embedding approach as a method for describing complex biochemical systems, with the parts not only treated with different levels of theory, but computed with hybrid classical and quantum algorithms. Such strategies are critical if one wants to expand the focus to biochemical molecules that contain active regions that cannot be properly explained with traditional algorithms on classical computers. While we do not solve this problem here, we provide an overview of where the field is going to enable such problems to be tackled in the future.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Hierarchical Reinforcement Learning of a Short-Range Bond-Order Potential for Silica: Analytic Embedding of Coordination with Classical Efficiency

Reinforcement learning (RL) has recently emerged as a data-efficient strategy to parametrize short-range interatomic potentials. Building on our past RL optimization of pairwise silica models, we extend the framework to a bond-order (Tersoff-type) potential that provides an analytic embedding of local coordination through a three-body term. A hierarchical RL workflow combining continuous-action Monte Carlo Tree Search and property-based rewards efficiently explores the 26-dimensional parameter space, sequentially optimizing lattice parameters, densities, angles, and cohesive energies of 21 silica polymorphs. The resulting models, Q-Tersoff and ML-Tersoff, reproduce the energetic ordering of low-energy phases and capture the angular correlations and amorphous structure factors of silica with improved fidelity over pairwise force fields, while remaining orders of magnitude faster than high-dimensional machine-learned potentials. Both models underperform for elastic constants and high-energy frameworks, delineating the limits of the current analytic form. The approach establishes a general and interpretable route to angle-aware, short-range potentials that bridge physics-based and machine-learned descriptions of silicate materials.

36 MATERIALS SCIENCE↗

Magnetic charge propagation upon a 3D artificial spin-ice

Magnetic charge propagation in spin-ice materials has yielded a paradigm-shift in science, allowing the symmetry between electricity and magnetism to be studied. Recent work is now suggesting the spin-ice surface may be important in mediating the ordering and associated phase space in such materials. Here, we detail a 3D artificial spin-ice, which captures the exact geometry of bulk systems, allowing magnetic charge dynamics to be directly visualized upon the surface. Using magnetic force microscopy, we observe vastly different magnetic charge dynamics along two principal directions. For a field applied along the surface termination, local energetics force magnetic charges to nucleate over a larger characteristic distance, reducing their magnetic Coulomb interaction and producing uncorrelated monopoles. In contrast, applying a field transverse to the surface termination yields highly correlated monopole-antimonopole pairs. Detailed simulations suggest it is the difference in effective chemical potential as well as the energy landscape experienced during dynamics that yields the striking differences in monopole transport.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Yang-Lee edge singularity triggered entanglement transition

Here we show that a class of $\mathscr{PT}$ symmetric non-Hermitian Hamiltonians realizing the Yang-Lee edge singularity exhibits an entanglement transition in the long-time steady state evolved under the Hamiltonian. Such a transition is induced by a level crossing triggered by the critical point associated with the Yang-Lee singularity and hence is first order in nature. At the transition, the entanglement entropy of the steady state jumps discontinuously from a volume-law to an area-law scaling. We exemplify this mechanism using a one-dimensional transverse field Ising model with additional imaginary fields, as well as the spin-1 Blume-Capel model and the three-state Potts model. We further make a connection to the forced-measurement induced entanglement transition in a Floquet nonunitary circuit subject to continuous measurements followed by post-selections. Our results demonstrate a new mechanism for entanglement transitions in non-Hermitian systems harboring a critical point.

36 MATERIALS SCIENCE↗

Semi-Classical Models for Virtual Antiparticle Pairs

Virtual particle-antiparticle pairs of massive elementary particle& are predicted in Quantum Field Theory (QFT) to appear from the vacuum and annihilate each other again within their Heisenberg lifetimes h/4mc(exp 2). In this work, semiclassical models of this process - for the cases of massive leptons, quarks, and the massive weak bosons W and Z - are constructed. It is shown that the dynamical lifetime of the particle- antiparticle system in each case equals the Heisenberg lifetime to good approximation, and obeys appropriate quantization conditions on the field fluctuation action. In other words, the dynamical lifetime of the semiclassical model agrees with QED and QCD to good approximation. But the formula for the dynamical lifetime in each model includes the force strength coupling constant (e in the lepton case, alpha(sup s) (q(exp 2)) in the quark cases), while the Heisenberg lifetime formula does not. Observing the agreement of the Heisenberg and dynamical lifetimes, we may derive the QED and QCD coupling constants in terms of h, c, and numerical factors only.

Batchelor, David↗

Small-scale properties from exascale computations of turbulence on a $\mathbf{32\,768^3}$ periodic cube

To study the physics of small-scale properties of homogeneous isotropic turbulence at increasingly high Reynolds numbers, direct numerical simulation results have been obtained for forced isotropic turbulence at Taylor-scale Reynolds number R λ = 2500 on a 32 768 3 three-dimensional periodic domain using a GPU pseudo-spectral code on a 1.1 exaflop GPU supercomputer (Frontier). These simulations employ the multi-resolution independent simulation (MRIS) technique (Yeung & Ravikumar 2020, Phys. Rev. Fluids, vol. 5, 110517) where ensemble averaging is performed over multiple short segments initiated from velocity fields at modest resolution, and subsequently taken to higher resolution in both space and time. Reynolds numbers are increased by reducing the viscosity with the large-scale forcing parameters unchanged. Although MRIS segments at the highest resolution for each Reynolds number last for only a few Kolmogorov time scales, small-scale physics in the dissipation range is well captured – for instance, in the probability density functions and higher moments of the dissipation rate and enstrophy density, which appear to show monotonic trends persisting well beyond the Reynolds number range in prior works in the literature. Attainment of range of length and time scales consistent with classical scaling also reinforces the potential utility of the present high-resolution data for studies of short-time-scale turbulence physics at high Reynolds numbers where full-length simulations spanning many large-eddy time scales are still not accessible. A single snapshot of the 32 768 3 data is publicly available for further analyses via the Johns Hopkins Turbulence Database.

intermittency↗

Fluid Physics in a Fluctuating Acceleration Environment

Our program of research aims at developing a stochastic description of the residual acceleration field onboard spacecraft (g-jitter) to describe in quantitative detail its effect on fluid motion. Our main premise is that such a statistical description is necessary in those cases in which the characteristic time scales of the process under investigation are long compared with the correlation time of g-jitter. Although a clear separation between time scales makes this approach feasible, there remain several difficulties of practical nature: (i), g-jitter time series are not statistically stationary but rather show definite dependences on factors such as active or rest crew periods; (ii), it is very difficult to extract reliably the low frequency range of the power spectrum of the acceleration field. This range controls the magnitude of diffusive processes; and (iii), models used to date are Gaussian, but there is evidence that large amplitude disturbances occur much more frequently than a Gaussian distribution would predict. The lack of stationarity does not constitute a severe limitation in practice, since the intensity of the stochastic components changes very slowly during space missions (perhaps over times of the order of hours). A separate analysis of large amplitude disturbances has not been undertaken yet, but it does not seem difficult a priori to devise models that may describe this range better than a Gaussian distribution. The effect of low frequency components, on the other hand, is more difficult to ascertain, partly due to the difficulty associated with measuring them, and partly because they may be indistinguishable from slowly changing averages. This latter effect is further complicated by the lack of statistical stationarity of the time series. Recent work has focused on the effect of stochastic modulation on the onset of oscillatory instabilities as an example of resonant interaction between the driving acceleration and normal modes of the system, and on cavity flow as an example of how an oscillatory response under periodic driving becomes diffusive if the forcing is random instead. This paper describes three different topics that illustrate behavior that is peculiar to a stochastic acceleration field. In the first case, we show that g-jitter can induce effective attractive or repulsive forces between a pair of spherical particles that are suspended in an incompressible fluid of different density provided that the momentum diffusion length is larger than the interparticle separation (as in the case in most colloidal suspensions). Second, a stochastic modulation of the control parameter in the vicinity of a pitchfork or supercritical bifurcation is known not to affect the location of the threshold. We show, however, that resonance between the modulation and linearly stable modes close to onset can lead to a shift in threshold. Finally, we discuss the classical problem of vorticity diffusion away from a plane boundary that is being vibrated along its own plane. Periodic motion with zero average vorticity production results in an exponential decay of the vorticity away from the boundary. Random vibration, on the other hand, results in power law decay away from the boundary even if vorticity production averages to zero.

Drolet, Francois↗

Turbomachine Sealing and Secondary Flows: Review of Rotordynamics Issues in Inherently Unsteady Flow Systems With Small Clearances - Part 2

Today's computational methods enable the determination of forces in complex systems, but without field validation data, or feedback, there is a high risk of failure when the design envelope is challenged. The data of Childs and Bently and field data reported in NASA Conference Proceedings serve as sources of design information for the development of these computational codes. Over time all turbomachines degrade and instabilities often develop, requiring responsible, accurate, turbomachine diagnostics with proper decisions to prevent failures. Tam et al. (numerical) and Bently and Muszynska (analytical) models corroborate and implicate that destabilizing factors are related through increases in the fluid-force average circumferential velocity. The stability threshold can be controlled by external swirl and swirl brakes and increases in radial fluid film stiffness (e.g., hydrostatic and ambient pressures) to enhance rotor stability. Also cited are drum rotor self-excited oscillations, where the classic fix is to add a split or severed damper ring or cylindrical damper drum, and the Benkert-Wachter work that engendered swirl brake concepts. For a smooth-operating, reliable, long-lived machine, designers must pay very close attention to sealing dynamics and diagnostic methods. Correcting the seals enabled the space shuttle main engine high-pressure fuel turbopump (SSME HPFTP) to operate successfully.

Hendricks, R. C.↗

Sound Transmission through Cylindrical Shell Structures Excited by Boundary Layer Pressure Fluctuations

This paper examines sound transmission into two concentric cylindrical sandwich shells subject to turbulent flow on the exterior surface of the outer shell. The interior of the shells is filled with fluid medium and there is an airgap between the shells in the annular space. The description of the pressure field is based on the cross-spectral density formulation of Corcos, Maestrello, and Efimtsov models of the turbulent boundary layer. The classical thin shell theory and the first-order shear deformation theory are applied for the inner and outer shells, respectively. Modal expansion and the Galerkin approach are used to obtain closed-form solutions for the shell displacements and the radiation and transmission pressures in the cavities including both the annular space and the interior. The average spectral density of the structural responses and the transmitted interior pressures are expressed explicitly in terms of the summation of the cross-spectral density of generalized force induced by the boundary layer turbulence. The effects of acoustic and hydrodynamic coincidences on the spectral density are observed. Numerical examples are presented to illustrate the method for both subsonic and supersonic flows.

Tang, Yvette Y.↗

Logarithmic profile of temperature in sheared and unstably stratified atmospheric boundary layers

The impact of buoyancy on the mean velocity, temperature, and scalar concentration profiles in the lower atmosphere is typically investigated within the framework of Monin-Obukhov similarity theory (MOST). MOST is the theoretical foundation for parametrizing surface-atmosphere exchanges in nearly all weather, climate, and hydrological models.According to MOST, the classic logarithmic profiles of mean velocity and temperature break down as the buoyancy effects become important. However, recent studies on turbulent Rayleigh-Bénard convection and natural convection along vertical walls suggest that the mean temperature in the near-surface region still follows a logarithmic profile.Motivated by these new results, we study the mean potential temperature profile in sheared and unstably stratified atmospheric boundary layers using direct numerical simulations and field observations. We find that the mean potential temperature profile remains logarithmic across a wide range of stability parameters, which characterizes the relative importance of buoyancy versus shear effects. Compared to MOST, our results suggest that the buoyancy force does not modify the logarithmic nature of the mean potential temperature profile, but instead modulates its slope, which is no longer universal and differs from 1/κ,whereκisthe von Kármán constant. Furthermore, this study provides another perspective on scalar turbulence in the atmospheric boundary layer

58 GEOSCIENCES↗

On the origin of ground-state vacuum-field catalysis: Equilibrium consideration

Recent experiments suggest that vibrational strong coupling (VSC) may significantly modify ground-state chemical reactions and their rates even without external pumping. The intrinsic mechanism of this “vacuum-field catalysis” remains largely unclear. Generally, modifications of thermal reactions in the ground electronic states can be caused by equilibrium or non-equilibrium effects. The former are associated with modifications of the reactant equilibrium distribution as expressed by the transition state theory of chemical reaction rates, while the latter stem from the dynamics of reaching and leaving transition state configurations. In this work, we examine how VSC can affect chemical reactions rates in a cavity environment according to transition state theory. Our approach is to examine the effect of coupling to cavity mode(s) on the potential of mean force (PMF) associated with the reaction coordinate. Within the context of classical nuclei and classical photons and also assuming no charge overlap between molecules, we find that while the PMF can be affected by the cavity environment, this effect is negligible for the usual micron-length cavities used to examine VSC situations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Deep Solar Activity Minimum 2007-2009: Solar Wind Properties and Major Effects on the Terrestrial Magnetosphere

We discuss the temporal variations and frequency distributions of solar wind and interplanetary magnetic field parameters during the solar minimum of 2007 - 2009 from measurements returned by the IMPACT and PLASTIC instruments on STEREO-A.We find that the density and total field strength were significantly weaker than in the previous minimum. The Alfven Mach number was higher than typical. This reflects the weakness of magnetohydrodynamic (MHD) forces, and has a direct effect on the solar wind-magnetosphere interactions.We then discuss two major aspects that this weak solar activity had on the magnetosphere, using data from Wind and ground-based observations: i) the dayside contribution to the cross-polar cap potential (CPCP), and ii) the shapes of the magnetopause and bow shock. For i) we find a low interplanetary electric field of 1.3+/-0.9 mV/m and a CPCP of 37.3+/-20.2 kV. The auroral activity is closely correlated to the prevalent stream-stream interactions. We suggest that the Alfven wave trains in the fast streams and Kelvin-Helmholtz instability were the predominant agents mediating the transfer of solar wind momentum and energy to the magnetosphere during this three-year period. For ii) we determine 328 magnetopause and 271 bow shock crossings made by Geotail, Cluster 1, and the THEMIS B and C spacecraft during a three-month interval when the daily averages of the magnetic and kinetic energy densities attained their lowest value during the three years under survey.We use the same numerical approach as in Fairfield's empirical model and compare our findings with three magnetopause models. The stand-off distance of the subsolar magnetopause and bow shock were 11.8 R(sub E) and 14.35 R(sub E), respectively. When comparing with Fairfield's classic result, we find that the subsolar magnetosheath is thinner by approx. 1 R(sub E). This is mainly due to the low dynamic pressure which results in a sunward shift of the magnetopause. The magnetopause is more flared than in Fairfield's model. By contrast the bow shock is less flared, and the latter is the result of weaker MHD forces.

solar minimum↗

Discrete stochastic model of point defect-dislocation interaction for simulating dislocation climb

Dislocation climb is an important high temperature process in metals plasticity, responsible for the phenomena such as creep, swelling, or hardening. Climb is defined by the ability of dislocations to leave their original glide plane by interacting with point defects. As such, dislocation climb is controlled by point defect diffusion/absorption/emission, all of which involve thermal activation. The existing thermodynamically consistent models for climb are generally formulated in a continuum framework, through the definition of effective defect fluxes and climb propensities in response to thermodynamic driving forces. However, the point-wise discrete nature of vacancies (and/or self-interstitials) confers a highly discrete nature to the climb dynamics, which is also strongly affected by elastic forces. The combination of discreteness, thermal activation, and elasticity makes this process too challenging for direct atomistic methods such as molecular dynamics. Here we develop a kinetic Monte Carlo model that captures vacancy generation and transport kinetics acting in conjuction with the evolving elastic fields provided by discrete dislocation dynamics simulations. The two models are coupled via the applied stresses and stress gradients generated by dislocation structures at vacancy locations. Our simulations reveal two surprising results. First that climb is dominated by vacancy emission even when the background vacancy concentration is much higher than the equilibrium one. And, second, that climb velocities might be much faster than otherwise believed when one uses the classical theories of climb. These effects are due to the locality of vacancy-dislocation processes, which are not captured in classical treatments that assume smooth vacancy fluxes and homogeneous concentrations. We apply the method to study elementary climb processes in body-centered cubic iron and furnish climb mobility functions to be used in parametric dislocation dynamics and/or crystal plasticity simulations. Here, we apply the technique to study non-conservative plastic bypass of spherical precipitates by edge dislocations and point out the differences between our discrete approach and existing continuum formulations.

36 MATERIALS SCIENCE↗

Direct numerical simulation of a combusting droplet with convection

The evaporation and combustion of a single droplet under forced and natural convection was studied numerically from first principles using a numerical scheme that solves the time-dependent multiphase and multispecies Navier-Stokes equations and tracks the sharp gas-liquid interface cutting across an arbitrary Eulerian grid. The flow fields both inside and outside of the droplet are resolved in a unified fashion. Additional governing equations model the interphase mass, energy, and momentum exchange. Test cases involving iso-octane, n-hexane, and n-propanol droplets show reasonable comparison rate, and flame stand-off distance. The partially validated code is, thus, readied to be applied to more demanding droplet combustion situations where substantial drop deformation render classical models inadequate.

Liang, Pak-Yan↗