Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Formal Methods”

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 667 records · Page 37

A hybrid calorimetry-simulation model of mixing enthalpy for molten salt

Calorimetric determination of enthalpies of mixing (ΔH mix ) in multicomponent molten salts is often interpreted using empirical models that lack physically meaningful parameters. However, for improving pyrochemical separation of spent nuclear fuel, where lanthanides are major fission products and critical elements, a deeper thermodynamic understanding of the link between excess thermodynamic properties and solvation structure is critically needed. In this work, we implement a hybrid and physics-informed framework, MIVM+Calorimetry+AIMD, which integrates experimentally measured ΔH mix (via high temperature drop calorimetry) with solvation structures from ab initio molecular dynamics (AIMD). This approach is demonstrated using LaCl 3 mixed with eutectic LiCl-KCl (58 mol% – 42 mol%) at 873 K and 1133 K. MIVM-derived parameters enable extrapolation of excess Gibbs energy and La 3+ activity across compositions. In contrast, direct ΔH mix predictions from AIMD and polarizable ion model simulations deviate significantly. By incorporating experimentally benchmarked solvation structures into an interpretable thermodynamic model, the MIVM+Calorimetry+AIMD formalism achieves higher accuracy and generalizable method for studying molten salts, offering a robust path for understanding and optimizing molten salt chemistry relevant to nuclear fuel cycles and separation science.

Goncharov, Vitaliy G. [Washington State Univ., Pul↗

Formation mechanism of guided resonances and bound states in the continuum in photonic crystal slabs

We develop a formalism, based on the mode expansion method, to describe the guided resonances and bound states in the continuum (BICs) in photonic crystal slabs with one-dimensional periodicity. This approach provides analytic insights to the formation mechanisms of these states: the guided resonances arise from the transverse Fabry–Pérot condition, and the divergence of the resonance lifetimes at the BICs is explained by a destructive interference of radiation from different propagating components inside the slab. As a result, we show BICs at the center and on the edge of the Brillouin zone protected by symmetry, BICs at generic wave vectors not protected by symmetry, and the annihilation of BICs at low-symmetry wave vectors.

36 MATERIALS SCIENCE↗

Control of simulation of a particle accelerator as a spatially distributed system

The work deals with the construction of strategy for finding optimal control of systems for spatially distributed systems on an example of problems in physics beams. The suggested approach is based on use of exact and approximate invariants as a controlling mechanism for control of modelling process of beam lines. The classical problem of optimal control of similar systems as an optimal control problem is formulated in a functional space. This approach is based on algebraic methods Lie and the corresponding matrix formalism that allows you to build not only the formal procedures for finding optimal solutions, but also provide the necessary constructive computational procedures.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Fano effect dominance over Coulomb blockade in transport properties of parallel coupled quantum dot system

Theoretical study of the Coulomb blockade effect on transport properties (Transmission Probability and I-V characteristics) for varied configuration of coupled quantum dot system has been studied by using Non Equilibrium Green Function(NEGF) formalism and Equation of Motion(EOM) method in the presence of magnetic flux. The self consistent approach and intra-dot Coulomb interaction is being taken into account. As the key parameters of the coupled quantum dot system such as dot-lead coupling, inter-dot tunneling and magnetic flux threading through the system can be tuned, the effect of asymmetry parameter and magnetic flux on this tuning is being explored in Coulomb blockade regime. The presence of the Coulomb blockade due to on-dot Coulomb interaction decreases the width of transmission peak at energy level ε + U and by adjusting the magnetic flux the swapping effect in the Fano peaks in asymmetric and symmetric parallel configuration sustains despite strong Coulomb blockade effect.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Hierarchy of forward-backward stochastic Schrödinger equation

Driven by the impetus to simulate quantum dynamics in photosynthetic complexes or even larger molecular aggregates, we have established a hierarchy of forward-backward stochastic Schrödinger equation in the light of stochastic unravelling of the symmetric part of the influence functional in the path-integral formalism of reduced density operator. The method is numerically exact and is suited for Debye-Drude spectral density, Ohmic spectral density with an algebraic or exponential cutoff, as well as discrete vibrational modes. The power of this method is verified by performing the calculations of time-dependent population differences in the valuable spin-boson model from zero to high temperatures. By simulating excitation energy transfer dynamics of the realistic full FMO trimer, some important features are revealed.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Microscopic Theory of Nuclear Fission: A Review

This paper reviews how nuclear fission is described within nuclear density functional theory. A distinction should be made between spontaneous fission, where half-lives are the main observables and quantum tunnelling the essential concept, and induced fission, where the focus is on fragment properties and explicitly time-dependent approaches are often invoked. Overall, the cornerstone of the density functional theory approach to fission is the energy density functional formalism. The basic tenets of this method, including some well-known tools such as the Hartree–Fock–Bogoliubov (HFB) theory, effective two-body nuclear potentials such as the Skyrme and Gogny force, finite-temperature extensions and beyond mean-field corrections, are presented succinctly. The energy density functional approach is often combined with the hypothesis that the time-scale of the large amplitude collective motion driving the system to fission is slow compared to typical time-scales of nucleons inside the nucleus. In practice, this hypothesis of adiabaticity is implemented by introducing (a few) collective variables and mapping out the many-body Schrödinger equation into a collective Schrödinger-like equation for the nuclear wave-packet. The region of the collective space where the system transitions from one nucleus to two (or more) fragments defines what are called the scission configurations. The inertia tensor that enters the kinetic energy term of the collective Schrödinger-like equation is one of the most essential ingredients of the theory, since it includes the response of the system to small changes in the collective variables. For this reason, the two main approximations used to compute this inertia tensor, the adiabatic time-dependent HFB and the generator coordinate method, are presented in detail, both in their general formulation and in their most common approximations. The collective inertia tensor enters also the Wentzel–Kramers–Brillouin (WKB) formula used to extract spontaneous fission half-lives from multi-dimensional quantum tunnelling probabilities (For the sake of completeness, other approaches to tunnelling based on functional integrals are also briefly discussed, although there are very few applications.) It is also an important component of some of the time-dependent methods that have been used in fission studies. Concerning the latter, both the semi-classical approaches to time-dependent nuclear dynamics and more microscopic theories involving explicit quantum-many-body methods are presented. One of the hallmarks of the microscopic theory of fission is the tremendous amount of computing needed for practical applications. In particular, the successful implementation of the theories presented in this article requires a very precise numerical resolution of the HFB equations for large values of the collective variables. This aspect is often overlooked, and several sections are devoted to discussing the resolution of the HFB equations, especially in the context of very deformed nuclear shapes. In particular, the numerical precision and iterative methods employed to obtain the HFB solution are documented in detail. Finally, a selection of the most recent and representative results obtained for both spontaneous and induced fission is presented, with the goal of emphasizing the coherence of the microscopic approaches employed. In conclusion, although impressive progress has been achieved over the last two decades to understand fission microscopically, much work remains to be done. Several possible lines of research are outlined in the conclusion.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

The Gutzwiller conjugate gradient minimization method for correlated electron systems

In this report we review our recent work on the Gutzwiller conjugate gradient minimization method, an ab initio approach developed for correlated electron systems. The complete formalism has been outlined that allows for a systematic understanding of the method, followed by a discussion of benchmark studies of dimers, one- and two-dimensional single-band Hubbard models. In the end, we present some preliminary results of multi-band Hubbard models and large-basis calculations of F 2 to illustrate our efforts to further reduce the computational complexity.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Phenotypic redshifts with self-organizing maps: A novel method to characterize redshift distributions of source galaxies for weak lensing

Wide-field imaging surveys such as the Dark Energy Survey (DES) rely on coarse measurements of spectral energy distributions in a few filters to estimate the redshift distribution of source galaxies. In this regime, sample variance, shot noise, and selection effects limit the attainable accuracy of redshift calibration and thus of cosmological constraints. We present a new method to combine wide-field, few-filter measurements with catalogues from deep fields with additional filters and sufficiently low photometric noise to break degeneracies in photometric redshifts. The multiband deep field is used as an intermediary between wide-field observations and accurate redshifts, greatly reducing sample variance, shot noise, and selection effects. Our implementation of the method uses self-organizing maps to group galaxies into phenotypes based on their observed fluxes, and is tested using a mock DES catalogue created from N-body simulations. Here, it yields a typical uncertainty on the mean redshift in each of five tomographic bins for an idealized simulation of the DES Year 3 weak-lensing tomographic analysis of σ Δz = 0.007, which is a 60 per cent improvement compared to the Year 1 analysis. Although the implementation of the method is tailored to DES, its formalism can be applied to other large photometric surveys with a similar observing strategy.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Entanglement and seniority

In quantum mechanics entanglement is the most striking phenomenon which has no counterpart in classical systems. Although, different approaches have already been developed to study correlations in the case of indistinguishable particles, the exploration of the so-called mode entanglement is still in its initial stage in nuclear physics. Study of mode-entanglement in the seniority model, derivation of analytic formulas for the one-body reduced density matrix of states with seniority $ν$ = 0, 1, 2, and $ν$ = 3, and also determination of the particle number dependence of the one-body reduced density matrix for arbitrary seniority. In addition, comparison of the predictions of the seniority model with the results of the shell model to gain insight into the structure and correlations of the ground and lowest yrast states. In the seniority model the analytic results are calculated using the quasi-spin formalism. The numerical shell model calculations are carried out in the standard shell model framework using an inert core. The applied realistic effective interactions are derived earlier using the $G$-matrix formalism. The density matrix renormalization group method is also applied in order to directly calculate the mode entropies. In the seniority model simple analytical expressions are given for the mode entropies. The peculiar behavior of the half-filled shells and the seniority-zero states are revealed. Numerical results are presented for the lightest stable calcium isotopes and for 94 Ru nucleus. For 94 Ru, the seniority model accounts for the 0$g$ 9/2 mode entropies, but seniority mixing is important for certain yrast states. Furthermore, interaction induced quantum fluctuations decrease the occupation of the 0$f$ 5/2 , 1$p$ 3/2 and 1$p$ 1/2 shells, and amount in non-negligible mode entropies on these shells, too, clearly outside the scope of the simple (0$g$ 9/2 ) 4 seniority model. The 0$f$ 7/2 shell-based seniority model is more accurate for the Ca isotopes, but seniority mixing is substantial for some 44 Ca yrast states, too. Mode and one-body entanglement entropies are useful tools to investigate the structure of quantum correlations in nuclei.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Gravitational form factors of charmonia

We investigate the gravitational form factors of charmonium. Our method is based on a Hamiltonian formalism on the light front known as basis light-front quantization. The charmonium mass spectrum and light-front wave functions were obtained from diagonalizing an effective Hamiltonian that incorporates confinement from holographic QCD and one-gluon exchange interaction from light-front QCD. We proposed a quantum many-body approach to construct the hadronic matrix elements of the energy momentum tensor T + + and T + − , which are used to extract the gravitational form factors A ( Q 2 ) and D ( Q 2 ) . The obtained form factors satisfy the known constraints, e.g., the von Laue condition. From these quantities, we also extract the energy, pressure and light-front energy distributions of the system. We find that hadrons are multilayer systems. Published by the American Physical Society 2024

Astronomy & Astrophysics↗

Parameterization of the nuclear Hulthén potentials

Within the formalism of supersymmetry-inspired factorization method, a two-term nuclear Hulthén potential has been developed and parameterized to reproduce the nucleon–nucleon scattering phase shifts for P and D partial wave states.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Vidyut3d: A Gpu Accelerated Fluid Solver for Non-Equilibrium Plasmas on Adaptive Grids

We present the numerical methods, programming methodology, verification, and performance assessment of a non-equilibrium plasma fluid solver that can effectively utilize current and upcoming central processing and graphics processing unit (CPU+GPU) architectures, in this work. Our plasma fluid model solves the coupled conservation equations for species transport, electrostatic Poisson and electron temperature on adaptive Cartesian grids. Our solver is written using performance portable adaptive-grid/particle management library, AMReX, and is portable over widely available vendor specific GPU architectures. We present verification of our solver using method of manufactured solutions that indicate formal second order accuracy with central diffusion and fifth-order weighted-essentially-non-oscillatory (WENO) advection scheme. We also verify our solver with published literature on capacitive discharges and atmospheric pressure streamer propagation. We demonstrate the use of our solver on two 3D simulation cases: an atmospheric streamer propagation in Ar-H2 mixtures and a low pressure twin electrode radio frequency reactor. Our performance studies on three different CPU+GPU architectures indicate approximately 150-400X speed-up using AMD and NVIDIA GPUs per time step compared to a single CPU core for a 4 million cell simulation with 15 species.

Sitaraman, Hariswaran↗

Achieving algorithmic resilience for temporal integration through spectral deferred corrections

Spectral deferred corrections (SDC) is an iterative approach for constructing higher-order-accurate numerical approximations of ordinary differential equations. SDC starts with an initial approximation of the solution defined at a set of Gaussian or spectral collocation nodes over a time interval and uses an iterative application of lower-order time discretizations applied to a correction equation to improve the solution at these nodes. Each deferred correction sweep increases the formal order of accuracy of the method up to the limit inherent in the accuracy defined by the collocation points. In this paper, we demonstrate that SDC is well suited to recovering from soft (transient) hardware faults in the data. A strategy where extra correction iterations are used to recover from soft errors and provide algorithmic resilience is proposed. Specifically, in this approach the iteration is continued until the residual (a measure of the error in the approximation) is small relative to the residual of the first correction iteration and changes slowly between successive iterations. Here, we demonstrate the effectiveness of this strategy for both canonical test problems and a comprehensive situation involving a mature scientific application code that solves the reacting Navier-Stokes equations for combustion research.

97 MATHEMATICS AND COMPUTING↗

Estimating Eigenenergies from Quantum Dynamics: A Unified Noise-Resilient Measurement-Driven Approach

Ground state energy estimation in physical, chemical, and materials sciences is one of the most promising applications of quantum computing. In this work, we introduce a new hybrid approach that finds the eigenenergies by collecting real-time measurements and post-processing them using the machinery of dynamic mode decomposition (DMD). From the perspective of quantum dynamics, we establish that our approach can be formally understood as a stable variational method on the function space of observables available from a quantum many-body system. We also provide strong theoretical and numerical evidence that our method converges rapidly even in the presence of a large degree of perturbative noise, and show that the method bears an isomorphism to robust matrix factorization methods developed independently across various scientific communities. Our numerical benchmarks on spin and molecular systems demonstrate an accelerated convergence and a favorable resource reduction over state-of-the-art algorithms. The DMD-centric strategy can systematically mitigate noise and stands out as a leading hybrid quantum-classical eigensolver.

Shen, Yizhi↗

Theory of molecular rate processes in the presence of intense laser radiation

The present paper deals with the influence of intense laser radiation on gas-phase molecular rate processes. Representations of the radiation field, the particle system, and the interaction involving these two entities are discussed from a general rather than abstract point of view. The theoretical methods applied are outlined, and the formalism employed is illustrated by application to a variety of specific processes. Quantum mechanical and semiclassical treatments of representative atom-atom and atom-diatom collision processes in the presence of a field are examined, and examples of bound-continuum processes and heterogeneous catalysis are discussed within the framework of both quantum-mechanical and semiclassical theories.

George, T. F.↗

The fate of the downgoing slab - A study of the moment tensors from body waves of complex deep-focus earthquakes

Two large, complex deep-focus earthquakes are analyzed, using body waves, to determine the relative locations and equivalent forces of the component events. Details of the source mechanism determination are given, with special emphasis on the means of comparing results for double-couple and more general models. The formalism includes a description of a method for calculating the confidence ellipses around the principal axes of the moment tensor. An interpretation of the results suggests that the mechanism of deep-focus earthquakes is controlled by the buoyancy of the slab and the change in viscosity at the 650 km discontinuity.

Strelitz, R. A.↗

Improved perfect-fluid energy-momentum tensor with spin in Einstein-Cartan space-time

The description of the spin given here is classical in that it is intrinsic but not quantized. The approach in this matter is similar to, for example, the work of Bailey and Israel (1973, 1975, 1979), where the fluid particles, which have intrinsic spin, may be galaxies or clusters of galaxies. The elementary particles of these objects and the 'ferromagnetic alignment' of their quantum spins are not resorted to in order to describe a fluid with spin. Physically this means that the equation of motion for the spin tensor is a modified Fermi-Walker transport equation (Misner et al., 1973), arising as a direct result of the inclusion of spin as an intrinsic variable in the thermodynamic description of the internal energy. The variables in this description are classical variables throughout and are not microscopic fields. An improved perfect-fluid energy-momentum tensor that includes spin and torsion is presented. Use is made of a Lagrangian variational principle based on the tetrad formalism of Halbwach (1960) and the method od constraints of Ray (1972).

Ray, J. R.↗

Partial redistribution in high-density, highly ionized plasmas

The conditions for which partial redistribution functions must be used for radiation transport in high-Z, high-density, laser-produced plasmas are examined. A previously developed two-photon formalism based on the model microfield method is used to calculate redistribution functions including electron and ion Stark broadening with ion dynamic effects. The competition between the relaxation rates and spontaneous emission is shown to determine the conditions for partial redistribution. The discussion makes use of microfield fluctuation rates and broadening coefficients which can be determined from simulation calculations. The redistribution function for the Ly-alpha transition of Ar XVIII is presented for typical plasma conditions.

Talin, B.↗