Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Moment equations”

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.

224 records · Page 13
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Towards grounding nuclear physics in QCD

Exascale computing could soon enable a predictive theory of nuclear structure and reactions rooted in the Standard Model, with quantifiable and systematically improvable uncertainties. Such a predictive theory will help exploit experiments that use nucleons and nuclei as laboratories for testing the Standard Model and its limitations. Examples include direct dark matter detection, neutrinoless double beta decay, and searches for permanent electric dipole moments of the neutron and atoms. It will also help connect QCD to the properties of cold neutron stars and hot supernova cores. We discuss how a quantitative bridge between QCD and the properties of nuclei and nuclear matter will require a synthesis of lattice QCD (especially as applied to two- and three- nucleon interactions), effective field theory, and ab initio methods for solving the nuclear many-body problem. While there are significant challenges that must be addressed in developing this triad of theoretical tools, the rapid advance of computing is accelerating progress. In particular, we focus this review on the anticipated advances from lattice QCD and how these advances will impact few-body effective theories of nuclear physics by providing critical input, such as constraints on unknown low-energy constants of the effective (field) theories. We also review particular challenges that must be overcome for the successful application of lattice QCD for low-energy nuclear physics. We describe progress in developing few-body effective (field) theories of nuclear physics, with an emphasis on HOBET, a non-relativistic effective theory of nuclear physics, which is less common in the literature. We use the examples of neutrinoless double beta decay and the nuclear-matter equation of state to illustrate how the coupling of lattice QCD to effective theory might impact our understanding of symmetries and exotic astrophysical environments.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Improvements to a class of hybrid methods for radiation transport: Nyström reconstruction and defect correction methods

In this study, two modifications are introduced for improving the accuracy, versatility, and robustness of a class of hybrid methods for radiation transport. In general, such methods are constructed by splitting the radiative flux into collided and uncollided components to which low- and high-resolution angular approximations are applied, respectively. In this work we focus on discrete ordinates discretizations of high and low order. The first modification we introduce changes the way in which the collided component is mapped into the uncollided component at the end of each time step in a simulation. The new mapping is a Nyström-type reconstruction that is applicable to arbitrary discrete ordinates quadratures, is guaranteed to preserve positivity of the solution provided that all ordinate weights are positive, is significantly more accurate than previous methods, and can be readily extended to other discretizations such as moment methods, finite element methods, and diffusion approximations. The second modification leverages integral deferred correction (IDC) to iteratively correct for the splitting error introduced by the inconsistency in angular discretization between the collided and uncollided components, in addition to improving the accuracy of the low-order temporal error that is treated by traditional IDC methods. Numerical tests in one- and two-dimensional geometries are used to demonstrate the increased accuracy and efficiency of the proposed modifications. It is found that the two techniques combined yield methods with solution accuracy and memory requirements comparable to that of monolithic discrete ordinates methods while reducing runtime by as much as a factor of between two and ten, depending on the problem.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Exact and locally implicit source term solvers for multifluid-Maxwell systems

Recently, a family of models that couple multifluid systems to the full Maxwell equations have been used in laboratory, space, and astrophysical plasma modeling. These models are more complete descriptions of the plasma than reduced models like magnetohydrodynamic (MHD) since they are derived more closely from the full kinetic Vlasov-Maxwell system, without assumptions like quasi-neutrality, negligible electron mass, etc. Thus these models naturally retain non-ideal MHD effects like electron inertia, Hall term, pressure anisotropy/nongyrotropy, displacement current, among others. One obstacle to broader application of these model is that an explicit treatment of their source terms leads to the need to resolve rapid processes like plasma oscillation and electron cyclotron motion, even when these are not important. In this paper, we suggest two ways to address this issue. First, we derive the analytic solutions to the source update equations, which can be implemented as a practical, but less generic solver. We then develop a time-centered, locally implicit algorithm to update the source terms, allowing stepping over the fast kinetic time-scales. For a plasma with S species, the locally implicit algorithm involves inverting a local (3 S + 3) × (3 S + 3) matrix only, thus is very efficient. The performance can be further increased by using the direct update formulas to skip null calculations. In this paper, we present benchmarks illustrating the exact energy-conservation of the locally implicit solver, as well as its efficiency and robustness for both small-scale, idealized problems and largescale, complex systems. The locally implicit algorithm can be also easily extended to include other local sources, like collisions and ionization, which are difficult to solve analytically.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

N -representability violations in truncated equation-of-motion coupled-cluster methods

One-electron reduced density matrices (1RDMs) from equation-of-motion (EOM) coupled-cluster with single and double excitations (CCSD) calculations are analyzed to assess their N-representability (i.e., whether they are derivable from a physical N-electron state). We identify EOM-CCSD stationary states whose 1RDMs violate either ensemble-state N-representability conditions or pure-state conditions known as generalized Pauli constraints. As such, these 1RDMs do not correspond to any physical N-electron state. Unphysical states are also encountered in the course of time-dependent EOM-CC simulations; when an external field drives transitions between a pair of stationary states with pure-state N-representable 1RDMs, the 1RDM of the time-dependent state can violate ensemble-state conditions. Furthermore, these observations point to potential challenges in interpreting the results of time-dependent EOM-CCSD simulations.

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

Quantitative radiography for determining density fluctuations in HED experiments

We have developed a method to extract density fluctuation measurements from x-ray radiographs of high-energy density (HED) instability growth and turbulence experiments. We use this information to calculate density fluctuation statistics for constraining the performance of turbulent mix models in HED systems. The density calculation combines image filtering, removal of systemic effects such as backlighter variation, calculation of transmission across multiple materials, and use of tracer materials to generate an approximate single-material density field. From the density map, we calculate both average density and a variance-like moment b (density-specific-volume covariance), which we compare to our models. We infer both quantities from a single image, which is significantly more information than the historic single scalar mix width measurements. We also develop a method of analyzing simulation outputs that incorporate both the density fluctuation metric from a turbulence model and the bulk material maps from the hydrodynamic code. This analysis helps address the question of how to initialize the simulations for best comparison to data from systems with large separations of scale in the mixing perturbation initial condition. We find that our data analysis method yields 1D average density and b curves with similar morphology and amplitudes as those from preliminary simulation comparisons.

47 OTHER INSTRUMENTATION↗
DOE OSTI2026

Geometric decoherence time in Lindbladian dynamics

The onset of decoherence in open many-body systems lacks a dynamical timescale grounded in the loss of bipartite entanglement. Here, we introduce the geometric decoherence time, defined as the earliest moment the monotone relation between logarithmic negativity and Rényi-$\frac{1}{2}$ entropy—exactly equal across any bipartition for pure states—breaks down under open-system evolution, signaling entropy growth without accompanying entanglement growth. We establish this criterion in both single-particle Gaussian dynamics and many-body Lindbladian evolution. We show that quantum mutual information provides a complementary long-time diagnostic: Its asymptotic vanishing is equivalent to factorization of the steady state across the bipartition, a condition strictly stronger than separability, and whenever a product steady state is approached exponentially in trace norm, negativity and mutual information share the same decay rate. In the presence of a strong symmetry, this tracking can fail—residual classical correlations can survive after entanglement has vanished. In the Kitaev chain with balanced gain and loss, we derive a closed-form solution and show that the topological phase sustains longer coherence times than the trivial phase at identical dissipation, with a local minimum at the chiral-symmetric point. In the interacting XXZ chain, exact many-body evolution shows that local 𝑍 dephasing preserves residual classical correlations, whereas gain and loss restore the mutual-information tracking of negativity. Furthermore, our results establish the geometric decoherence time as a dynamical scale tracking the onset of decoherence.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2025

A New Measurement of the Neutron Electric Form Factor with the Super Bigbite Spectrometer Apparatus

Protons and neutrons, collectively known as nucleons, make up the nuclei at the core of atoms which form our world. The nucleon has been under intensive study for over 100 years, and yet we still do not fully understand the internal dynamics which govern properties like its spin or its mass-which contributes to almost all of the visible mass in the universe. These dynamics are governed by quantum chromodynamics(QCD), the predictions of which are experimentally tested at high energy accelerator facilities such as Jefferson Lab. TheGEN-II experiment (E12-09-016) is one such experiment. GEN-II is part of the Super Bigbite Spectrometer (SBS) experimental form factor pro gramme taking place in Hall A at Jefferson Lab, which aims to make precision measurements of the nucleon electromagnetic form factors(EMFFs)at record high values of squared four momentum transfer ¿2. EMFFs describe the electric and magnetic moment distributions within the nucleon. They can be measured through elastic electron scattering off the nucleon, and describe the recoil response of the target nucleon at a given energy scale. GEN-II is a double polarized semi-exclusive beam target asymmetry (BTA) experiment, seeking to measure the electric form factor of the neutron,¿¿ ¿,at three new values of squared four-momentum transfer¿2 =2.92,6.74and9.82GeV². The latter two points being at record high ¿2. The form factor is determined through measuring the BTA of quasielliptical scattering of a neutron from a polarized nuclear target. The experiment utilized the CEBAF accelerator to produce longitudinally polarized electrons up to ~85% polarization, which were scattered off neutrons within a novel polarized helium-3 (³He) target. This new polarized ³He target was employed by building on the technology of its precursors which existed in similar preceding experiments. This target was designed to operate at the high luminosities typical of Hall A, and reached a record breaking combination of polarization and beam intensity known as figure of merit, three times larger than those predecessors. The SBS collaboration designed and constructed two brand new high acceptance spectrometers for these experiments, an electron arm named Bigbite (BB) and a hadron arm named Super Bigbite. Both spectrometers featured a large acceptance EM dipole magnet, and complementary detector systems. The electron arm contained gaseous electronmultipliers(GEMs) which were used for high precision tracking of the scattered electrons, a heavy gas cherenkov (GRINCH) which was used for PID between electrons and pions, a plastic i ii scintillator timing hodoscope to provide high resolution timing of the start of events, and a pair of EM calorimeters (BBCal) which provided energy measurements of detected particles, and provided the experimental trigger. The hadron arm also contained a system of GEMs which will be utilized for future SBS experiments, and a hadron calorimeter designed to provide position, timing and energy measurements of the recoiling nucleon. The calibration of all detector subsystems, beam and target data is discussed, with a focus on novel timing calibrations to the hodoscope and hadron calorimeter. An analysis of selecting quasielliptical events and suppressing background contributions from a number of sources which contaminate the final event sample is given. The largest irreducible backgrounds are found to be from misidentified protons, timing accidentals and inelastic events. The physical asymmetry is measured and used to extract a value for the form factor ratio ¿¿ ¿/¿¿ ¿. High precision ¿2 data for ¿¿ ¿ is used to then extract ¿¿ ¿. This work finds at ¿2 = 2.92 GeV2 that ¿¿ ¿ = 0.0129+0.0019 -0.0020. This result is in statistical agreement with existing fits to world data, and predictions from the constituent quark model and Dyson–Schwinger equations, in this region of ¿2.

Penman, Gary [Univ. of Glasgow, Scotland (United K↗
DOE OSTI2021

Development of five-moment two-fluid modeling for Z-pinch physics

The Z-pinch m = 0 instability as well as its stabilization by radially sheared axial flow is studied using the nonlinear ideal five-moment twofluid (5M2F) model with an extension of that model to include Braginskii heat and momentum transport. Using the ideal 5M2F model, linear growth rate results are compared with prior work using MHD and Hall MHD. At small normalized wavenumber, 1 < k a < 4 , where a is the effective pinch radius, 5M2F results agree with Hall MHD within 20% in scenarios without radially sheared axial flow. With the sheared flow and focusing on ka = 10/3, agreement with Hall MHD is excellent. In the limit of small ion inertial length, results also match with MHD. A comparison with PIC modeling of shear-free m = 0 stability focuses on a plasma scenario based on recent experimental results. In a scan of mode wavenumber, ideal 5M2F results are qualitatively similar to PIC: the growth rate rises to a peak at a moderate wavenumber and declines at a large wavenumber in contrast to MHD results, which show the saturation of the growth rate with the increasing wavenumber rather than a decline. The peak normalized 5M2F growth rate is &#x3B3; &#x3C4; A = 1.5 , where sA is the Alfven transit time across the pinch. The peak occurs at normalized wavenumber ka = 10. For comparison, PIC results have a peak growth of &#x3B3; &#x3C4; A = 0.77 at ka = 5. Including Braginskiibased closure of the 5M2F model does not qualitatively change the ideal results in this particular case. Nonlinear saturation is studied using the 5M2F model with the dissipative Braginskii-based closure in cases with pinch-edge sheared-flow speed equal to half the Alfven speed. Nonlinear mixing due to the sheared flow yields a quasi-steady state after modest losses of pinch ion inventory and pinch thermal energy, approximately 30% and 10%, respectively. 5M2F modeling captures the essential physics of m = 0 instability and offers a computationally tractable route to high-fidelity modeling of 3D Z-pinch behavior, including m = 1 instability.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗