Search NASA⌕ Search

SEARCH · Search NASA

Results for “Corrected equation”

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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Nonresonant two-photon x-ray absorption in Cu

We present a real-space Green's function theory and calculations of two-photon x-ray absorption (TPA). Our focus is on nonresonant 𝐾-shell TPA in metallic Cu, which has been observed experimentally at intense x-ray free electron laser (XFEL) sources. The theory is based on an independent particle Green's function treatment of the Kramers-Heisenberg equation and an approximation for the sum over nonresonant intermediate states in terms of a static quadrupole transition operator. XFEL effects are modeled by a partially depleted 𝑑 band. This approach is shown to give results for 𝐾-shell TPA in quantitative agreement with XFEL experiment and with a Bethe-Salpeter equation approach. Furthermore, we also briefly discuss many-body corrections and TPA sum rules.

Approximation methods for many-body systems↗
DOE OSTI2021

Analysis of Lobe Power Calculator and Indication System with Physics and Cycle Based Models

The Advanced Test Reactor (ATR) at INL measures reactor power through various methods, two of them being thermal and Nitrogen-16 (N-16) activity. Water power calculator (WPC) is a thermal power system that measures flow and temperature to determine thermal quadrant powers. The N-16 system utilizes a beta detector that outputs Nitrogen activation levels to calculate lobe power through an algorithm called Lobe Power Calculation and Indication system (LPCIS). The LPCIS utilizes the N-16 system and the WPC system to determine reactor core power levels. The WPC provides accurate calculations of quadrant and total reactor thermal power. With the use of WPC measurements, thermal-to-N-16 (T2N) power ratios are produced to determine if the two indication systems agree on core power. Relative magnitude equations are used to utilize N-16 coefficients and multipliers to improve the indications of the LPCIS. These correct indications are crucial for maintaining safety limits because operators rely on this information for decision making. Currently, the LPCIS system uses linear equations and matrices to calculate lobe power through multipliers and coefficients. Advancements in technology and system upgrades have increased the accuracy of power readings by making the system more dynamic. The new proposed coefficient and multiplier method implements a cycle specific and physics-based model intended for changing coefficients during operation. This calibration experiment focused on power splitting, outer shim and neck shim, as well as fuel burning into the reactor digital acquisition system (RDAS) weighting factors. Results demonstrated that the physics learning method yields a smaller error margin inside of the desired power range for the data set 166-A. Continuing to improve the physics-based model will help improve the power accuracy of the LPCIS system.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2024

Diffusion-Limited Kinetics in Reactive Systems

A proper representation of chemical kinetics is vital to understanding, modeling, and optimizing many important chemical processes. In liquid and surface phases, where diffusion is slow, the rate at which the reactants diffuse together limits the overall rate of many elementary reactions. Commonly, the textbook Smoluchowski theory is utilized to estimate effective rate coefficients in the liquid phase. On surfaces, modelers commonly resort to much more complex and expensive Kinetic Monte Carlo (KMC) simulations. Here, in this study, we extend the Smoluchowski model to allow the diffusing species to undergo chemical reactions and derive analytical formulas for the diffusion-limited rate coefficients for 3D, 2D, and 2D/3D interface cases. With these equations, we are able to demonstrate that when species react faster than they diffuse they can react orders of magnitude faster than predicted by Smoluchowski theory, through what we term “the reactive transport effect”. We validate the derived steady-state equations against particle Monte Carlo (PMC) simulations, KMC simulations, and non-steady-state solutions. Furthermore, using PMC and KMC simulations, we propose corrections that agree with all limits and the computed data for the 2D and 2D/3D interface steady-state equations, accounting for unique limitations in the associated derived equations. Additionally, we derive equations to handle couplings between diffusion-limited rate coefficients in reaction networks. We believe these equations should make it possible to run much more accurate mean-field simulations of liquids, surfaces, and liquid–surface interfaces accounting for diffusion limitations and the reactive transport effect.

Johnson, Matthew Sean↗
DOE OSTI2020

Improvements to CTF Code Verification and Unit Testing (FY2020)

In 2010, the U.S. Department of Energy created its first Energy Innovation Hub, which focuses on improving Light Water Reactors (LWRs) through Modeling and Simulation. This hub, named the Consortium for the Advanced Simulation of LWRs (CASL), attempts to characterize and understand LWR behavior under normal operating conditions and use any gained insights to improve their efficiency. In collaboration with North Carolina State University (NCSU), CASL has worked extensively on the thermal-hydraulic subchannel code Coolant Boiling in Rod Arrays—Three Field (COBRA-TF). The NCSU/CASL version of COBRA-TF has been rebranded as CTF. This document focuses on code verification test problems that ensure CTF converges to the correct answer for the intended application. The suite of code verification tests are mapped to the underlying conservation equations of CTF, and significant gaps are addressed. Convergence behavior and numerical errors are quantified for each of the tests. Tests that converge at the correct rate to the corresponding analytic solution are incorporated into the CTF automated regression suite. A new verification utility is created for this purpose, which enables code verification by generalizing the process. For problems that do not behave correctly, the results are reported but the problem is not included in the regression suite. In addition to verification studies, this document also quantifies the existing tests of constitutive models. A few existing gaps are addressed by adding new unit tests.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2020

FY20 Improvements to CTF Code Verification and Unit Testing

In 2010, the U.S. Department of Energy created its first Energy Innovation Hub, which focuses on improving Light Water Reactors (LWRs) through Modeling and Simulation. This hub, named the Consortium for the Advanced Simulation of LWRs (CASL), attempts to characterize and understand LWR behavior under normal operating conditions and use any gained insights to improve their efficiency. In collaboration with North Carolina State University (NCSU), CASL has worked extensively on the thermal-hydraulic subchannel code Coolant Boiling in Rod Arrays–Three Field (COBRA-TF). The NCSU/CASL version of COBRA-TF has been rebranded as CTF. This document focuses on code verification test problems that ensure CTF converges to the correct answer for the intended application. The suite of code verification tests are mapped to the underlying conservation equations of CTF, and significant gaps are addressed. Convergence behavior and numerical errors are quantified for each of the tests. Tests that converge at the correct rate to the corresponding analytic solution are incorporated into the CTF automated regression suite. A new verification utility is created for this purpose, which enables code verification by generalizing the process. For problems that do not behave correctly, the results are reported but the problem is not included in the regression suite. In addition to verification studies, this document also quantifies the existing tests of constitutive models. A few existing gaps are addressed by adding new unit tests.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2023

Consistent truncations in higher derivative supergravity

We consider the torus reduction of heterotic supergravity in the presence of four-derivative corrections. In particular, the reduction on T n generically leads to a half-maximal supergravity coupled to n vector multiplets, and we show that it is consistent to truncate out said vector multiplets. This is done by the analysis of both the bosonic equations of motion and the Killing spinor equations. As an application of the consistent truncation, we examine the four-derivative corrected BPS black string that reduces to a black hole in minimal nine-dimensional supergravity.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Full versus quasiparticle self-consistency in vertex-corrected GW approaches

Here, using seven semiconductors/insulators with band gaps covering the range from 1 eV to 10 eV we systematically explore the performance of two different variants of self-consistency associated with famous Hedin's system of equations: the full self-consistency and the so called quasi-particle approximation to it. The pros and cons of these two variants of self-consistency are sufficiently well documented in literature for the simplest GW approximation to the Hedin's equations. Our study, therefore, aims primarily at the level of theory beyond GW approximation, i.e. at the level of theory which includes vertex corrections. Whereas quasi-particle self-consistency has certain advantages at GW level (well known fact), the situation becomes quite different when vertex corrections are included. In the variant with full self-consistency, vertex corrections (both for polarizability and for self energy) systematically reduce the calculated band gaps making them closer to the experimental values. In the variant with quasi-particle self-consistency, however, an inclusion of the same diagrams has considerably larger effect and calculated band gaps become severely underestimated. Different effect of vertex corrections in two variants of self-consistency can be related to the Z-factor cancellation which plays positive role in quasi-particle self-consistency at GW level of theory but appears to be destructive for the quasi-particle approximation when higher order diagrams are included. Second result of our study is that we were able to reproduce the results obtained with the Questaal code using our FlapwMBPT code when the same variant of self-consistency (quasi-particle) and the same level of vertex corrections (for polarizability only, static approximation for screened interaction, and Tamm-Dancoff approximation for the Bethe-Salpeter equation) are used.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Enhanced contribution of the pairing gap to the QCD equation of state at large isospin chemical potential

In this paper QCD at large isospin density is studied, which is known to be in the superfluid state with Cooper pairs carrying the same quantum number as pions. The gap equation derived from the perturbation theory up to the next-to-leading-order corrections is solved. The pairing gap at large isospin chemical potential is found to be enhanced compared to the color-superconducting gap at large baryon chemical potential due to the 2 difference in the exponent arising from the stronger attraction in one-gluon exchange in the singlet channel. Then, using the gap function, the contribution of the condensation energy of the superfluid state to the QCD equation of state is evaluated. At isospin chemical potential of a few GeV, where the lattice QCD and the perturbative QCD can be both applied, the effect of the condensation energy becomes dominant even compared to the next-to-leading order corrections to the pressure in the perturbation theory. It resolves the discrepancy between the recent lattice QCD results and the perturbative QCD result. Published by the American Physical Society 2024

Astronomy & Astrophysics↗
DOE OSTI2022

A hybrid nodal-staggered pseudo-spectral electromagnetic particle-in-cell method with finite-order centering

Electromagnetic particle-in-cell (PIC) codes are widely used to perform computer simulations of a variety of physical systems, including fusion plasmas, astrophysical plasmas, plasma wakefield particle accelerators, and secondary photon sources driven by ultra-intense lasers. In a PIC code, Maxwell's equations are solved on a grid with a numerical method of choice. This article focuses on pseudo-spectral analytical time-domain (PSATD) algorithms and presents a novel hybrid PSATD PIC scheme that combines the respective advantages of standard nodal and staggered methods. The novelty of the hybrid scheme consists in using finite-order centering of grid quantities between nodal and staggered grids, in order to combine the solution of Maxwell's equations on a staggered grid with the deposition of charges and currents and the gathering of electromagnetic forces on a nodal grid. The correctness and performance of the novel hybrid scheme are assessed by means of numerical tests that employ different classes of PSATD equations in a variety of physical scenarios, ranging from the modeling of electron-positron pair creation in vacuum to the simulation of laser-driven and particle beam-driven plasma wakefield acceleration. It is shown that the novel hybrid scheme offers significant numerical and computational advantages, compared to purely nodal or staggered methods, for all the test cases presented.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

A phase-field method for boiling heat transfer

Here we present a phase field method for heat transfer in two-phase flow with boiling. The vapor/liquid interface evolution is modeled by the Cahn-Hilliard equation. The phase change rate is determined by accounting for the heat conduction balance on either vapor or liquid side of the interfacial area, depending on which side the temperature is assumed to be maintained at the saturation temperature during boiling. The velocity correction scheme proposed by Dong & Shen [27] is extended to solve the Navier-Stokes equations for a non-solenoidal velocity field, and the entropy viscosity method is employed for stabilization. The phase change model is verified by two-dimensional simulations of a vapor bubble growing in super-heated liquid and in film boiling. In both cases, mesh independence of the results is systematically performed. Subsequently, the method is applied to predict the growth of three-dimensional vapor bubble in a rectangular microchannel with boiling flow, achieving good agreement with experimental measurements and available simulation results using the level-set method. The numerical experiments demonstrate that the required mesh resolution for the phase field method is comparable with that of volume of fluid (VoF) and level-set methods.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Numerical algorithms for water waves with background flow over obstacles and topography

Abstract We present two accurate and efficient algorithms for solving the incompressible, irrotational Euler equations with a free surface in two dimensions with background flow over a periodic, multiply connected fluid domain that includes stationary obstacles and variable bottom topography. One approach is formulated in terms of the surface velocity potential while the other evolves the vortex sheet strength. Both methods employ layer potentials in the form of periodized Cauchy integrals to compute the normal velocity of the free surface, are compatible with arbitrary parameterizations of the free surface and boundaries, and allow for circulation around each obstacle, which leads to multiple-valued velocity potentials but single-valued stream functions. We prove that the resulting second-kind Fredholm integral equations are invertible, possibly after a physically motivated finite-rank correction. In an angle-arclength setting, we show how to avoid curve reconstruction errors that are incompatible with spatial periodicity. We use the proposed methods to study gravity-capillary waves generated by flow around several elliptical obstacles above a flat or variable bottom boundary. In each case, the free surface eventually self-intersects in a splash singularity or collides with a boundary. We also show how to evaluate the velocity and pressure with spectral accuracy throughout the fluid, including near the free surface and solid boundaries. To assess the accuracy of the time evolution, we monitor energy conservation and the decay of Fourier modes and compare the numerical results of the two methods to each other. We implement several solvers for the discretized linear systems and compare their performance. The fastest approach employs a graphics processing unit (GPU) to construct the matrices and carry out iterations of the generalized minimal residual method (GMRES).

Ambrose, David M.↗
DOE OSTI2026

Recent progress on coarse graining simulations

We focus on coarse graining simulations based on the primary conservation equations, effectively codesigned physics and algorithms, and low-Mach-number corrected (LMC) hydrodynamics. Simulation methods involve LANL’s x-Radiation-Adaptive-Grid-Eulerian Large-Eddy Simulation, Besnard-Harlow-Rauenzahn (BHR) Reynolds-Averaged Navier-Stokes (RANS) approach, and Dynamic BHR – a paradigm bridging RANS and LES. A relevant question addressed relates to whether 3D RANS and RANS/LES hybrids – the industry standards for aerospace and automotive research, are presently relevant for practical variable-density applications involving shocked and accelerated interface instabilities. Furthermore, recent simulations of the GaTECH inclined mixing-layer shock-tube and NIF ICF-capsule experiments are used to demonstrate issues, challenges, and potential for 3D coarse grained LMC simulation strategies for robustly simulating complex transitional and coupled hydrodynamics-multiphysics with coarser resolution. Present LES readiness to provide accurate predictions at scale is demonstrated – whereas 3D RANS and RANS/LES bridging do not appear impactful in this context.

42 ENGINEERING↗
DOE OSTI2022

Symplectic effective field theory for nuclear structure studies

Here, a Symplectic Effective Field Theory that unveils the observed emergence of symplectic symmetry in atomic nuclei is advanced. Specifically, starting from a simple extension of the harmonic-oscillator Lagrangian, an effective field theory applied against symplectic basis states is shown to yield a Hamiltonian system with one fitted parameter. The scale of the system can be determined self consistently as the ratio of the average volume of a nucleus assumed to be spherical to its volume as determined by the average number of oscillator quanta, which is stretched by the fact that the plane-wave solution satisfies the equations of motion at every order without the need for perturbative corrections. As an application of the theory, results for 20 Ne, 22 Ne and 22 Mg are presented that yield energy spectra, B(E2) values, and matter radii in good agreement with experimentally measured results.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Hyperbolic string tadpole

Hyperbolic geometry on the one-bordered torus is numerically uniformized using Liouville theory. This geometry is relevant for the hyperbolic string tadpole vertex describing the one-loop quantum corrections of closed string field theory. We argue that the Lamé equation, upon fixing its accessory parameter via Polyakov conjecture, provides the input for the characterization. The explicit expressions for the Weil-Petersson metric as well as the local coordinates and the associated vertex region for the tadpole vertex are given in terms of classical torus conformal blocks. The relevance of this vertex for vacuum shift computations in string theory is highlighted.

Firat, Atakan Hilmi↗
DOE OSTI2023

Single inclusive hadron production in DIS at small x: next to leading order corrections

We calculate the one-loop corrections to single inclusive hadron production in Deep Inelastic Scattering (DIS) at small x in the forward rapidity region using the Color Glass Condensate formalism. We show that the divergent parts of the next to leading order (NLO) corrections either cancel among each other or lead to x (rapidity) evolution of the leading order (LO) dipole cross section according to the JIMWLK evolution equation and DGLAP evolution of the parton-hadron fragmentation function. The remaining finite parts constitute the NLO (α s ) corrections to the LO single inclusive hadron production cross section in DIS at small x.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Machine learning of superconducting critical temperature from Eliashberg theory

Abstract The Eliashberg theory of superconductivity accounts for the fundamental physics of conventional superconductors, including the retardation of the interaction and the Coulomb pseudopotential, to predict the critical temperature T c . McMillan, Allen, and Dynes derived approximate closed-form expressions for the critical temperature within this theory, which depends on the electron–phonon spectral function α 2 F ( ω ). Here we show that modern machine-learning techniques can substantially improve these formulae, accounting for more general shapes of the α 2 F function. Using symbolic regression and the SISSO framework, together with a database of artificially generated α 2 F functions and numerical solutions of the Eliashberg equations, we derive a formula for T c that performs as well as Allen–Dynes for low- T c superconductors and substantially better for higher- T c ones. This corrects the systematic underestimation of T c while reproducing the physical constraints originally outlined by Allen and Dynes. This equation should replace the Allen–Dynes formula for the prediction of higher-temperature superconductors.

36 MATERIALS SCIENCE↗
DOE OSTI2024

An asymptotic Grad–Shafranov equation for quasisymmetric stellarators

A first-order model is derived for quasisymmetric stellarators where the vacuum field due to coils is dominant, but plasma-current-induced terms are not negligible and can contribute to magnetic differential equations, with $\beta$ of the order of the ratio induced to vacuum fields. Under these assumptions, it is proven that the aspect ratio must be large and a simple expression can be obtained for the lowest-order vacuum field. The first-order correction, which involves both vacuum and current-driven fields, is governed by a Grad–Shafranov equation and the requirement that flux surfaces exist. These two equations are not always consistent, and so this model is generally overconstrained, but special solutions exist that satisfy both equations simultaneously. One family of such solutions is the set of first-order near-axis solutions. Thus, the first-order near-axis model is a subset of the model presented here. Several other solutions outside the scope of the near-axis model are also found. A case study comparing one such solution to a VMEC-generated solution shows good agreement.

Nikulsin, Nikita (ORCID:0000000318611777)↗