Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “FORMULATION”

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 757 records · Page 42
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

Formulation of the Rice-Walsh equation of state based on Shock Hugoniot data for porous metals

The dimensionless material parameter R introduced by Wu and Jing into the Rice-Walsh equation of state (EOS) has been deduced from the LASL shock Hugoniot data for porous Al and Cu. It was found that the parameter R/p decays smoothly with shock pressure p and displays small experimental scatter in the high pressure region. This finding led to the conclusion that the parameter has only a weak temperature dependence and is well approximated by a function of pressure alone, and the Grüneisen parameter should be temperature dependent under compression. The thermodynamic formulation of the Rice-Walsh EOS for Al and Cu was realized using the empirically determined function R(p) for each material and their known shock Hugoniot. It was then possible to reproduce porous shock Hugoniot for these metals. For most degrees of porosity, agreement between the porous data and the calculated Hugoniots using the empirical function described was very good. However, slight discrepancies were seen for Hugoniots with very high porosity. Two new thermal variables were introduced after further analysis, which enabled the calculation of the cold compression curve for these metals. The Grüneisen parameters along full-density and porous Hugoniot curve were calculated using a thermodynamic identity connecting R and the Grüneisen parameter. It was shown that the Grüneisen parameter is strongly temperature dependent. The present analysis suggested that the Rice-Walsh type EOS is a preferable choice for the analysis with its simple form, pressure-dependent empirical Wu-Jing parameter, and its compatibility with porous shock data.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

Action-angle formulation of generalized, orbit-based, fast-ion diagnostic weight functions

Due to the usually complicated and anisotropic nature of the fast-ion distribution function, diagnostic velocity-space weight functions, which indicate the sensitivity of a diagnostic to different fast-ion velocities, are used to facilitate the analysis of experimental data. Additionally, when velocity-space weight functions are discretized, a linear equation relating the fast-ion density and the expected diagnostic signal is formed. In a technique known as velocity-space tomography, many measurements can be combined to create an ill-conditioned system of linear equations that can be solved using various computational methods. However, when velocity-space weight functions (which by definition ignore spatial dependencies) are used, velocity-space tomography is restricted, both by the accuracy of its forward model and also by the availability of spatially overlapping diagnostic measurements. In this work, we extend velocity-space weight functions to a full 6D generalized coordinate system and then show how to reduce them to a 3D orbit-space without loss of generality using an action-angle formulation. Moreover, we show how diagnostic orbit-weight functions can be used to infer the full fast-ion distribution function, i.e., orbit tomography. In depth derivations of orbit weight functions for the neutron, neutral particle analyzer, and fast-ion D-α diagnostics are also shown

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2018

Accurate, rapid identification of dislocation lines in coherent diffractive imaging via a min-max optimization formulation

Defects such as dislocations impact materials properties and their response during external stimuli. Imaging these defects in their native operating conditions to establish the structure-function relationship and, ultimately, to improve performance via defect engineering has remained a considerable challenge for both electron-based and x-ray-based imaging techniques. While Bragg coherent x-ray diffractive imaging (BCDI) is successful in many cases, nuances in identifying the dislocations has left manual identification as the preferred method. Derivative-based methods are also used, but they can be inaccurate and are computationally inefficient. Here we demonstrate a derivative-free method that is both more accurate and more computationally efficient than either derivative-or human-based methods for identifying 3D dislocation lines in nanocrystal images produced by BCDI. We formulate the problem as a min-max optimization problem and show exceptional accuracy for experimental images. We demonstrate a 227x speedup for a typical experimental dataset with higher accuracy over current methods. We discuss the possibility of using this algorithm as part of a sparsity-based phase retrieval process. Here, we also provide MATLAB code for use by other researchers.

47 OTHER INSTRUMENTATION↗
DOE OSTI2018

Perturbative variational formulation of the Vlasov-Maxwell equations

The perturbative variational formulation of the Vlasov-Maxwell equations is presented up to the third order in the perturbation analysis. From the second and third-order Lagrangian densities, the first-order and second-order Vlasov-Maxwell equations are expressed in gauge-invariant and gauge-independent forms, respectively. Upon deriving the reduced second-order Vlasov-Maxwell Lagrangian for the linear nonadiabatic gyrokinetic Vlasov-Maxwell equations, the reduced Lagrangian densities for the linear drift-wave equation and the linear hybrid kinetic-magnetohydrodynamic (MHD) equations are derived, with their associated wave-action conservation laws obtained by the Noether method. The exact wave-action conservation law for the linear hybrid kinetic-MHD equations is written explicitly. To conclude, a new form of the third-order Vlasov-Maxwell Lagrangian is derived in which ponderomotive effects play a crucial role.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Formulation of 8-moment plasma transport with application to the Nernst effect

8-moment plasma models using two different distribution functions are used to study the Nernst effect and heat transport in dense plasma. These models are presented in hyperbolic form in contrast to traditional parabolic systems derived from perturbing the distribution function, as in Braginskii [Rev. Plasma Phys. 1, 205 (1965)]. The hyperbolic moment formulation can be solved implicitly in time with straightforward and fast local solvers. The numerical implementation of 8-moment models with the relaxation method in the PERSEUS code is also presented. To test 8-moment PERSEUS compared to Braginskii’s transport equations, a verification test for the Nernst thermomagnetic wave by Velikovich et al. [Phys. Plasmas 26, 112702 (2019)] is performed that confirms the presence of the same physics, but with slight differences in the transport coefficients, which are tabulated in the limits of high and low magnetization.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Analytic evaluation of energy first derivatives for spin–orbit coupled-cluster singles and doubles augmented with noniterative triples method: General formulation and an implementation for first-order properties

The paper reports a formulation of analytic energy first derivatives for the coupled-cluster singles and doubles augmented with noniterative triples [CCSD(T)] method with spin–orbit coupling included at the orbital level and an implementation for evaluation of first-order properties. Here, the new development extends the applicability of relativistic coupled cluster methods for heavy-element chemistry and spectroscopy. Benchmark calculations of first-order properties including dipole moments and electric-field gradients using the corresponding exact two-component property integrals are presented for heavy-element containing molecules to demonstrate the usefulness of the present analytic scheme.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Spin-free formulation of the multireference driven similarity renormalization group: A benchmark study of first-row diatomic molecules and spin-crossover energetics

We report a spin-free formulation of the multireference (MR) driven similarity renormalization group (DSRG) based on the ensemble normal ordering of Mukherjee and Kutzelnigg (1997). This ensemble averages over all microstates of a given total spin quantum number, and therefore, it is invariant with respect to SU(2) transformations. As such, all equations may be reformulated in terms of spin-free quantities and they closely resemble those of spin-adapted closed-shell coupled cluster (CC) theory. The current implementation is used to assess the accuracy of various truncated MR-DSRG methods (perturbation theory up to third order and iterative methods with single and double excitations) in computing the constants of 33 first-row diatomic molecules. The accuracy trends for these first-row diatomics are consistent with our previous benchmark on a small subset of closed-shell diatomic molecules. We then present the first MR-DSRG application on transition-metal complexes by computing the spin splittings of the [Fe(H 2 O) 6 ] 2+ and [Fe(NH 3 ) 6 ] 2+ molecules. A focal point analysis (FPA) shows that third-order perturbative corrections are essential to achieve reasonably converged energetics. The FPA based on the linearized MR-DSRG theory with one- and two-body operators and up to a quintuple- ζ basis set predicts the spin splittings of [Fe(H 2 O) 6 ] 2+ and [Fe(NH 3 ) 6 ] 2+ to be –35.7 and –17.1 kcal mol –1 , respectively, showing good agreement with the results of local CC theory with singles, doubles, and perturbative triples.

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

A modified Villain formulation of fractons and other exotic theories

We reformulate known exotic theories (including theories of fractons) on a Euclidean spacetime lattice. We write them using the Villain approach and then we modify them to a convenient range of parameters. The new lattice models are closer to the continuum limit than the original lattice versions. In particular, they exhibit many of the recently found properties of the continuum theories including emergent global symmetries and surprising dualities. Also, these new models provide a clear and rigorous formulation to the continuum models and their singularities. In appendices, we use this approach to review well-studied lattice models and their continuum limits. These include the XY-model, the $\mathbb{Z}_N$ clock-model, and various gauge theories in diverse dimensions. Here, this presentation clarifies the relation between the condensed-matter and the high-energy views of these systems. It emphasizes the role of symmetries associated with the topology of field space, duality, and various anomalies.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Efficient formulation of multitime generalized quantum master equations: Taming the cost of simulating 2D spectra

Modern 4-wave mixing spectroscopies are expensive to obtain experimentally and computationally. In certain cases, the unfavorable scaling of quantum dynamics problems can be improved using a generalized quantum master equation (GQME) approach. However, the inclusion of multiple (light–matter) interactions complicates the equation of motion and leads to seemingly unavoidable cubic scaling in time. In this paper, we present a formulation that greatly simplifies and reduces the computational cost of previous work that extended the GQME framework to treat arbitrary numbers of quantum measurements. Specifically, we remove the time derivatives of quantum correlation functions from the modified Mori–Nakajima–Zwanzig framework by switching to a discrete-convolution implementation inspired by the transfer tensor approach. We then demonstrate the method’s capabilities by simulating 2D electronic spectra for the excitation-energy-transfer dimer model. In our method, the resolution of data can be arbitrarily coarsened, especially along the t2 axis, which mirrors how the data are obtained experimentally. Even in a modest case, this demands O(103) fewer data points. We are further able to decompose the spectra into one-, two-, and three-time correlations, showing how and when the system enters a Markovian regime where further measurements are unnecessary to predict future spectra and the scaling becomes quadratic. This offers the ability to generate long-time spectra using only short-time data, enabling access to timescales previously beyond the reach of standard methodologies.

Chemistry↗
DOE OSTI2024

Lagrangian formulation of nuclear–electronic orbital Ehrenfest dynamics with real-time TDDFT for extended periodic systems

Here, we present a Lagrangian-based implementation of Ehrenfest dynamics with nuclear–electronic orbital (NEO) theory and real-time time-dependent density functional theory for extended periodic systems. In addition to a quantum dynamical treatment of electrons and selected protons, this approach allows for the classical movement of all other nuclei to be taken into account in simulations of condensed matter systems. Furthermore, we introduce a Lagrangian formulation for the traveling proton basis approach and propose new schemes to enhance its application for extended periodic systems. Validation and proof-of-principle applications are performed on electronically excited proton transfer in the o-hydroxybenzaldehyde molecule with explicit solvating water molecules. These simulations demonstrate the importance of solvation dynamics and a quantum treatment of transferring protons. This work broadens the applicability of the NEO Ehrenfest dynamics approach for studying complex heterogeneous systems in the condensed phase.

Calculus of variations↗
DOE OSTI2026

An algebraic convolution formulation for multiple-scattering correction in small-angle neutron scattering

Multiple scattering in small-angle neutron scattering (SANS) redistributes spectral weight and distorts structural interpretation, particularly for thick or strongly scattering samples. We develop a finite-dimensional spectral desmearing framework that corrects multiple scattering without resorting to integral transforms or model-dependent extrapolation. The primary intensity is expanded in an orthonormal basis adapted to the isotropic transverse-momentum measure, under which convolution reduces to a recursive tensor contraction, allowing the Poisson-weighted multiple-scattering series to be evaluated directly in a finite-dimensional basis representation. This formulation yields a stable forward–inverse mapping between apparent and primary spectra. Numerical tests demonstrate convergence under repeated convolution and accurate recovery of the single-scattering intensity. Application to SANS measurements collected at multiple neutron facilities, including the Spallation Neutron Source, the High Flux Isotope Reactor, and the Institut Laue-Langevin, shows the quantitative reconstruction of the underlying primary spectrum across a wide range of transmission conditions, including strongly attenuating samples. Here, the method provides a stable, model-agnostic framework for multiple-scattering correction in SANS and enables consistent structural interpretation across instruments and scattering regimes.

Tung, Chi-Huan [Oak Ridge National Laboratory (ORN↗
DOE OSTI2016

An introduction to optimal power flow: Theory, formulation, and examples

The set of optimization problems in electric power systems engineering known collectively as Optimal Power Flow (OPF) is one of the most practically important and well-researched subfields of constrained nonlinear optimization. OPF has enjoyed a rich history of research, innovation, and publication since its debut five decades ago. Nevertheless, entry into OPF research is a daunting task for the uninitiated--both due to the sheer volume of literature and because OPF's ubiquity within the electric power systems community has led authors to assume a great deal of prior knowledge that readers unfamiliar with electric power systems may not possess. This article provides an introduction to OPF from an operations research perspective; it describes a complete and concise basis of knowledge for beginning OPF research. The discussion is tailored for the operations researcher who has experience with nonlinear optimization but little knowledge of electrical engineering. Topics covered include power systems modeling, the power flow equations, typical OPF formulations, and common OPF extensions.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2011

POLARIZED LINE FORMATION IN MULTI-DIMENSIONAL MEDIA. IV. A FOURIER DECOMPOSITION TECHNIQUE TO FORMULATE THE TRANSFER EQUATION WITH ANGLE-DEPENDENT PARTIAL FREQUENCY REDISTRIBUTION

To explain the linear polarization observed in spatially resolved structures in the solar atmosphere, the solution of polarized radiative transfer (RT) equation in multi-dimensional (multi-D) geometries is essential. For strong resonance lines, partial frequency redistribution (PRD) effects also become important. In a series of papers, we have been investigating the nature of Stokes profiles formed in multi-D media including PRD in line scattering. For numerical simplicity, so far we have restricted our attention to the particular case of PRD functions which are averaged over all the incident and scattered directions. In this paper, we formulate the polarized RT equation in multi-D media that takes into account the Hanle effect with angle-dependent PRD functions. We generalize here to the multi-D case the method for Fourier series expansion of angle-dependent PRD functions originally developed for RT in one-dimensional geometry. We show that the Stokes source vector S = (S{sub I} , S{sub Q} , S{sub U} ){sup T} and the Stokes vector I = (I, Q, U){sup T} can be expanded in terms of infinite sets of components S-tilde{sup (k)}, I-tilde{sup (k)}, respectively, k in [0, +{infinity}). We show that the components S-tilde{sup (k)} become independent of the azimuthal angle ({psi}) of the scattered ray, whereas the components I-tilde{sup (k)} remain dependent on {psi} due to the nature of RT in multi-D geometry. We also establish that S-tilde{sup (k)} and I-tilde{sup (k)} satisfy a simple transfer equation, which can be solved by any iterative method such as an approximate Lambda iteration or a Bi-Conjugate Gradient-type projection method provided we truncate the Fourier series to have a finite number of terms.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

EFFECTS OF DUST GROWTH AND SETTLING ON THE IONIZATION BY RADIONUCLIDES. I. FORMULATION AND RESULTS IN A QUIESCENT STATE OF PROTOPLANETARY DISKS

We investigate the evolution of the ionization rates by the decay of radionuclides in protoplanetary disks at the early stage of planet formation where size growth and settling of dust particles proceed extensively. Because most of the nuclides to ionize gas, such as short-lived nuclide {sup 26}Al and long-lived one {sup 40}K, are refractory elements, they are contained in the solid material of dust particles. Thus, the ionization by these nuclides is affected by the following three processes: (1) the change of the relative abundance of dust particles due to the settling toward the midplane of the disk, (2) the energy loss of emitted energetic particles inside the solid material of dust particles, and (3) the absorption of energetic particles by the other dust particles located nearby. In this series of papers we comprehensively investigate the basic physical processes, calculate the settling and size growth of dust particles numerically, and clarify the evolution of the ionization rates relative to their initial values in various disk models at this stage. In this paper we investigate the energy-loss processes concerning dust particles, formulate the coalescence equation for settling particles, and apply them to quiescent disk models that are similar to the solar nebula. For simplicity, dust particles are assumed to be compact spheres that remain perfect sticking for mutual collisions. Because the settling of dust particles is not appreciable in the first 10{sup 3} yr, the ionization rate varies little except in the outermost part near the disk surface. As the settling proceeds, the rate around the midplane increases considerably. The maximum ionization rates by {sup 26}Al in the minimum mass solar nebula are about 100, 51, and 14 times larger than their initial values for the orbits R = 0.5, 1, and 5 AU, respectively, which are close to or exceed the ionization rate by cosmic ray in the interstellar medium. The rates by {sup 40}K also increase by factors of about 36, 19, and 5 at the same orbits. In the inner orbital regions, these rates exceed the rates by the attenuated cosmic rays by an order of magnitude. The rates in the residual parts decrease extensively as time goes by, because amounts of the floating dust particles decrease continuously.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Magnetic potential based formulation for linear and non-linear 3D RF sheath simulation

Abstract This paper reports a new numerical scheme to simulate the radio-frequency (RF) induced RF sheath, which is suitable for a large 3D simulation. In the RF sheath boundary model, the tangential component of the electric field ( E t ) is given by the gradient of a scalar electric field potential. We introduce two additional scalar potentials for the tangential components of the magnetic field, which effectively impose the normal electric displacement ( D n ) on the plasma sheath boundary condition via in-homogeneous Neumann boundary condition and constrain the tangential electric field on the surface as curl-free ( ∇ × E t = 0 ). In our approach, the non-linear sheath impedance is formulated as a natural extension of the large thickness (or asymptotic) sheath limit ( D n = 0 ), allowing for handling both asymptotic and non-linear regimes seamlessly. The new scheme is implemented using the Petra-M finite element method analysis framework and is verified with simulations in the literature. The significance of non-linearity is discussed in various plasma conditions. An application of this scheme to asymptotic RF sheath simulation on the WEST ICRF antenna side limiters is also discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Nonequilibrium formulation of varying-temperature bit erasure

Landauer's principle states that erasing a bit of information at fixed temperature T costs at least $k$ B $T$ ln 2 units of work. Here we investigate erasure at varying temperature, to which Landauer's result does not apply. Here we formulate bit erasure as a stochastic nonequilibrium process involving a compression of configuration space, with physical and logical states associated in a symmetric way. Erasure starts and ends at temperature T, but temperature can otherwise vary with time in an arbitrary way. Defined in this way, erasure is governed by a set of nonequilibrium fluctuation relations that show that varying-temperature erasure can done with less work than $k$ B $T$ ln 2. As a result, erasure and the complementary process of bit randomization can be combined to form a work-producing engine cycle.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Energy and momentum conservation in the Euler–Poincaré formulation of local Vlasov–Maxwell-type systems

The action principle by Low for the classic Vlasov–Maxwell system contains a mix of Eulerian and Lagrangian variables. This renders the Noether analysis of reparametrization symmetries inconvenient, especially since the well-known energy- and momentum-conservation laws for the system are expressed in terms of Eulerian variables only. While an Euler–Poincaré formulation of Vlasov–Maxwell-type systems, effectively starting with Low's action and using constrained variations for the Eulerian description of particle motion, has been known for a while Cendra et al, it is hard to come by a documented derivation of the related energy- and momentum-conservation laws in the spirit of the Euler–Poincaré machinery. To our knowledge only one such derivation exists in the literature so far, dealing with the so-called guiding-center Vlasov–Darwin system Sugama et al. The present exposition discusses a generic class of local Vlasov–Maxwell-type systems, with a conscious choice of adopting the language of differential geometry to exploit the Euler–Poincaré framework to its full extent. After reviewing the transition from a Lagrangian picture to an Eulerian one, we demonstrate how symmetries generated by isometries in space lead to conservation laws for linear- and angular-momentum density and how symmetry by time translation produces a conservation law for energy density. Here, we also discuss what happens if no symmetries exist. Finally, two explicit examples will be given—the classic Vlasov–Maxwell and the drift-kinetic Vlasov–Maxwell—and the results expressed in the language of regular vector calculus for familiarity.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2012

MAGNETIC BRAKING FORMULATION FOR SUN-LIKE STARS: DEPENDENCE ON DIPOLE FIELD STRENGTH AND ROTATION RATE

We use two-dimensional axisymmetric magnetohydrodynamic simulations to compute steady-state solutions for solar-like stellar winds from rotating stars with dipolar magnetic fields. Our parameter study includes 50 simulations covering a wide range of relative magnetic field strengths and rotation rates, extending from the slow- and approaching the fast-magnetic-rotator regimes. Using the simulations to compute the angular momentum loss, we derive a semi-analytic formulation for the external torque on the star that fits all of the simulations to a precision of a few percent. This formula provides a simple method for computing the magnetic braking of Sun-like stars due to magnetized stellar winds, which properly includes the dependence on the strength of the magnetic field, mass loss rate, stellar radius, surface gravity, and spin rate, and which is valid for both slow and fast rotators.

79 ASTRONOMY AND ASTROPHYSICS↗