Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Physics Regularization”

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 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Sub-system self-consistency in coupled cluster theory

Here, in this article, we provide numerical evidence indicating that the single-reference coupled-cluster (CC) energies can be calculated alternatively to their copybook definition. We demonstrate that the CC energy can be reconstructed by diagonalizing the effective Hamiltonians describing correlated sub-systems of the many-body system. In the extreme case, we provide numerical evidence that the CC energy can be reproduced through the diagonalization of the effective Hamiltonian describing sub-system composed of a single electron. These properties of the CC formalism can be exploited to design protocols to define effective interactions in sub-systems used as probes to calculate the energy of the entire system and introduce a new type of self-consistency for approximate CC approaches.

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

KATRIN: status and prospects for the neutrino mass and beyond

Abstract The Karlsruhe Tritium Neutrino (KATRIN) experiment is designed to measure a high-precision integral spectrum of the endpoint region of T 2 β decay, with the primary goal of probing the absolute mass scale of the neutrino. After a first tritium commissioning campaign in 2018, the experiment has been regularly running since 2019, and in its first two measurement campaigns has already achieved a sub-eV sensitivity. After 1000 days of data-taking, KATRIN’s design sensitivity is 0.2 eV at the 90% confidence level. In this white paper we describe the current status of KATRIN; explore prospects for measuring the neutrino mass and other physics observables, including sterile neutrinos and other beyond-Standard-Model hypotheses; and discuss research-and-development projects that may further improve the KATRIN sensitivity.

beyond standard model↗
DOE OSTI2023

Completeness of the NASA/IPAC Extragalactic Database (NED) Local Volume Sample

We introduce the NASA/IPAC Extragalactic Database (NED) Local Volume Sample (NED-LVS), a subset of ~1.9 million objects with distances out to 1000 Mpc. We use UV and IR fluxes available in NED from all-sky surveys to derive physical properties, and estimate the completeness relative to the expected local luminosity density. The completeness relative to near-IR luminosities (which traces a galaxy's stellar mass) is roughly 100% at D < 30 Mpc and remains moderate (70%) out to 300 Mpc. For brighter galaxies (≳L * ), NED-LVS is ~100% complete out to ~400 Mpc. When compared to other local Universe samples (GLADE and HECATE), all three are ~100% complete below 30 Mpc. At distances beyond ~80 Mpc, NED-LVS is more complete than both GLADE and HECATE by ~10%–20%. NED-LVS is the underlying sample for the NED gravitational-wave follow-up service (NED-GWF), which provides prioritized lists of host candidates for GW events within minutes of alerts issued by the LIGO–Virgo–KAGRA collaboration. We test the prioritization of galaxies in the volume of GW170817 by three physical properties, where we find that both stellar mass and inverse specific star formation rate place the correct host galaxy in the top 10. In addition, NED-LVS can be used for a wide variety of other astrophysical studies: galaxy evolution, star formation, large-scale structure, galaxy environments, and more. The data in NED are updated regularly, and NED-LVS will be updated concurrently. Consequently, NED-LVS will continue to provide an increasingly complete sample of galaxies for a multitude of astrophysical research areas for years to come.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

The Palomar Transient Factory Core-collapse Supernova Host-galaxy Sample. I. Host-galaxy Distribution Functions and Environment Dependence of Core-collapse Supernovae

Several thousand core-collapse supernovae (CCSNe) of different flavors have been discovered so far. However, identifying their progenitors has remained an outstanding open question in astrophysics. Studies of SN host galaxies have proven to be powerful in providing constraints on the progenitor populations. In this paper, we present all CCSNe detected between 2009 and 2017 by the Palomar Transient Factory. This sample includes 888 SNe of 12 distinct classes out to redshift z ≈ 1. We present the photometric properties of their host galaxies from the far-ultraviolet to the mid-infrared and model the host-galaxy spectral energy distributions to derive physical properties. The galaxy mass function of Type Ic, Ib, IIb, II, and IIn SNe ranges from 105 to 1011.5 M o˙, probing the entire mass range of star-forming galaxies down to the least-massive star-forming galaxies known. Moreover, the galaxy mass distributions are consistent with models of star-formation-weighted mass functions. Regular CCSNe are hence direct tracers of star formation. Small but notable differences exist between some of the SN classes. Type Ib/c SNe prefer galaxies with slightly higher masses (i.e., higher metallicities) and star formation rates than Type IIb and II SNe. These differences are less pronounced than previously thought. H-poor superluminous supernovae (SLSNe) and SNe Ic-BL are scarce in galaxies above 1010 M o˙. Their progenitors require environments with metallicities of < 0.4 and < 1 solar, respectively. In addition, the hosts of H-poor SLSNe are dominated by a younger stellar population than all other classes of CCSNe. Our findings corroborate the notion that low metallicity and young age play an important role in the formation of SLSN progenitors.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Domain-Dependent Electronic Properties of 2D Copper Boride Synthesized from Reversible Subsurface Diffusion on Cu(111)

Emerging boron-based 2D materials hold properties with potential applications ranging from electronic devices to catalysis and energy storage. Using physical vapor deposition, 2D copper borides are formed on Cu(111) at elevated temperatures. However, the as-prepared surfaces often contain boron clusters that compromise structural homogeneity, and the dynamics of their formation is not yet fully understood. Through in situ low-energy electron microscopy (LEEM), Auger electron spectroscopy (AES), and scanning tunneling microscopy/spectroscopy (STM/STS), we show that cluster-free 2D CuB x is synthesized via reversible subsurface diffusion of boron atoms. The resolved regular atomic structure aligns with the previously calculated Cu 8 B 14 model. In a rare case, possibly due to surface strain and a change of orientation, a second domain of CuB x was observed with a pronounced change in electronic properties, yet electronic states remain delocalized in both domains. In conclusion, these findings demonstrate the sensitivity of electronic properties to domain structure and highlight a potential opportunity to tune 2D CuB x for diverse applications.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2026

Fast solvers for tokamak fluid models with PETSc

Multigrid (MG) is widely recognized as a highly effective solver for the model problem, the Laplacian, but textbook MG fails on most problems of interest. MG methods have been applied to complex, real-world applications with careful consideration of the physical model and discretization. In this work we develop the first step in applying MG methods to science and engineering relevant magnetohydrodynamics (MHD) tokamak models in the M3D-C1 (https://m3dc1.pppl.gov) fusion energy science code. The semi-implicit time integrator in M3D-C1 is composed of many linear solves. The implicit advance of the momentum equation is the most challenging and is the focus of this work. The current production solver in M3D-C1 is a block Jacobi (BJ) preconditioner within a Krylov solver, where blocks group degrees of freedom on planes of constant toroidal coordinate. BJ convergence degrades as the number of planes increases due to the spectral properties of the matrix preconditioned with BJ. The partially magnetic field-aligned, regular toroidal grid structure in M3D-C1 is amenable to semi-coarsening geometric MG in the toroidal direction. This paper develops such a solver and demonstrates competitive performance on a runaway electron model of a SPARC (https://cfs.energy/technology/sparc) disruption, and superior robustness on a stellarator model on which the BJ solver fails to converge.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

A generalized Selberg zeta function for flat space cosmologies

Flat space cosmologies (FSCs) are time dependent solutions of three-dimensional (3D) gravity with a vanishing cosmological constant. They can be constructed from a discrete quotient of empty 3D flat spacetime and are also called shifted-boost orbifolds. Using this quotient structure, we build a new and generalized Selberg zeta function for FSCs, and show that it is directly related to the scalar 1-loop partition function. We then propose an extension of this formalism applicable to more general quotient manifolds $\mathcal{M}$/ℤ, based on representation theory of fields propagating on this background. Our prescription constitutes a novel and expedient method for calculating regularized 1-loop determinants, without resorting to the heat kernel. We compute quasinormal modes in the FSC using the zeroes of a Selberg zeta function, and match them to known results.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Holographic thermal correlators: a tale of Fuchsian ODEs and integration contours

We analyze real-time thermal correlation functions of conserved currents in holographic field theories using the grSK geometry, which provides a contour prescription for their evaluation. We demonstrate its efficacy, arguing that there are situations involving components of conserved currents, or derivative interactions, where such a prescription is, in fact, essential. To this end, we first undertake a careful analysis of the linearized wave equations in AdS black hole backgrounds and identify the branch points of the solutions as a function of (complexified) frequency and momentum. All the equations we study are Fuchsian with only regular singular points that for the most part are associated with the geometric features of the background. Special features, e.g., the appearance of apparent singular points at the horizon, whence outgoing solutions end up being analytic, arise at higher codimension loci in parameter space. Using the grSK geometry, we demonstrate that these apparent singularities do not correspond to any interesting physical features in higher-point functions. We also argue that the Schwinger-Keldysh collapse and KMS conditions, implemented by the grSK geometry, continue to hold even in the presence of such singularities. For charged black holes above a critical charge, we furthermore demonstrate that the energy density operator does not possess an exponentially growing mode, associated with ‘pole-skipping’, from one such apparent singularity. Our analysis suggests that the connection between the scrambling physics of black holes and energy transport has, at best, a limited domain of validity.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

t 1 -noise elimination by continuous chemical shift anisotropy refocusing

We report that due to their high gyromagnetic ratio, there is considerable interest in measuring distances and correlations involving protons, but such measurements are compounded by the simultaneous recoupling of chemical shift anisotropy (CSA). This secondary recoupling adds additional modulations to the signal intensities that ultimately lead to t 1 -noise and signal decay. Recently, Venkatesh et al. demonstrated that the addition of CSA refocusing periods during 1 H-X dipolar recoupling led to sequences with far higher stability and performance. Herein, we describe a related effort and develop a symmetry-based recoupling sequence that continually refocuses the 1 H CSA. This sequence shows superior performance to the regular and t1-noise eliminated D-HMQC sequences in the case of spin-1/2 nuclei and comparable performance to the later for half-integer quadrupoles.

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

Modification of Lie's transform perturbation theory for charged particle motion in a magnetic field

It is pointed out that the conventional Lie transform perturbation theory for the guiding center motion of charged particles in a magnetic field needs to be modified for ordering inconsistency. There are two reasons. First, the ordering difference between the temporal variation of gyrophase and that of the other phase space coordinates needs to be taken into account. Second, it is also important to note that the parametric limit of the derivative of a function is not equivalent to the derivative of the limit function. When these facts are taken into account, the near identity transformation rule for one form related to the Lagrangian is modified. With the modified near identity transformation rule, the drift motion of charged particles can be described in the first order, instead of the second order and beyond through a tedious expansion process as in the conventional formulation. Here, this resolves the discrepancy between the direct and Lie transform treatments in the Lagrangian perturbation theory for charged particle motion in a magnetic field.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Differential equations for cosmological correlators

Cosmological fluctuations retain a memory of the physics that generated them in their spatial correlations. The strength of correlations varies smoothly as a function of external kinematics, which is encoded in differential equations satisfied by cosmological correlation functions. In this work, we provide a broader perspective on the origin and structure of these differential equations. As a concrete example, we study conformally coupled scalar fields in a power-law cosmology. The wavefunction coefficients in this model have integral representations, with the integrands being the product of the corresponding flat-space results and “twist factors” that depend on the cosmological evolution. Similar twisted integrals arise for loop amplitudes in dimensional regularization, and their recent study has led to the discovery of rich mathematical structures and powerful new tools for computing multi-loop Feynman integrals in quantum field theory. The integrals of interest in cosmology are also part of a finite-dimensional basis of master integrals, which satisfy a system of first-order differential equations. We develop a formalism to derive these differential equations for arbitrary tree graphs. The results can be represented in graphical form by associating the singularities of the differential equations with a set of graph tubings. Upon differentiation, these tubings grow in a local and predictive fashion. In fact, a few remarkably simple rules allow us to predict — by hand — the equations for all tree graphs. While the rules of this “kinematic flow” are defined purely in terms of data on the boundary of the spacetime, they reflect the physics of bulk time evolution. We also study the analogous structures in tr ϕ 3 theory, and see some glimpses of hidden structure in the sum over planar graphs. This suggests that there is an autonomous combinatorial or geometric construction from which cosmological correlations, and the associated spacetime, emerge.

Cosmological models↗
DOE OSTI2021

Highly efficient Pt catalyst on newly designed CeO 2 -ZrO 2 -Al 2 O 3 support for catalytic removal of pollutants from vehicle exhaust

Pt-CeO 2 catalysts have been widely studied for the vehicle emission control. Designing novel CeO 2 based supports with improved physical-chemical properties has become a research hotspot to further promote the catalytic performance and stability of Pt-CeO 2 catalysts. In this work, through utilizing a unique, two-step incipient wetness impregnation (T-IWI) method for ceria-zirconia-alumina (CZA-T) support preparation, a Pt single site catalyst (Pt/CZA-T) with excellent thermal stability was synthesized. Higher oxidation activity and Oxygen storage capacity (OSC) were achieved on activated Pt/CZA-T, comparing to Pt catalysts on regular CeO 2 /Al 2 O 3 (Pt/CA) and one-step prepared CeZrO x /Al 2 O 3 (Pt/CZA). Via the modification of hydrophilic/hydrophobic properties of γ-Al 2 O 3 by this unique T-IWI method, finer Ce 0.9 Zr 0.1 O 2 particles with higher density of surface defects were formed on CZA-T, on which a higher Pt dispersion and stronger Pt-O-Ce interaction were achieved. Upon activation, smaller, well-dispersed Pt clusters on CZA-T were generated. It was concluded that the CO oxidation performance and OSC were highly related to the size of Pt clusters on different supports that we have developed. Finally, more Pt sites located at Pt cluster-CeZrO x interfaces, which were the real active sites, were responsible for the highest OSC function and CO oxidation activity of activated Pt/CZA-T catalyst.

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

Numerical Investigation of Fluid Flow and Space Charge in Liquid Argon Time Projection Chamber (LArTPC) Detectors

Overview This project focused on developing a high-fidelity numerical framework to simulate the multiphysics environment within Liquid Argon Time Projection Chamber (LArTPC) detectors. The primary objective was to characterize the complex interplay between ion transport, background fluid dynamics, and electric field distortions—a critical factor for the calibration and sensitivity of next-generation High Energy Physics experiments, such as DUNE. Technical Achievements The research successfully yielded a hybrid numerical space-charge solver utilizing a Cell-Centered Finite Volume Method (FVM) for ion transport coupled with a Finite Element Method (FEM) for electric potential. Key accomplishments include: • Verification & Validation: The 3-D solver was rigorously verified against 1-D analytical solutions, demonstrating high numerical accuracy in predicting space-charge-induced field deviations. • Field Distortion Analysis: 3D simulations revealed that space charge effects introduce significant non-uniformities in the electric field. Critically, the research identified that background LAr flow velocities, when comparable to ion drift velocities, markedly exacerbate these distortions. • Technology Transfer: The resulting source code and comprehensive user manuals were successfully transferred to collaborators at Fermilab, providing a portable computational tool for the broader scientific community. Challenges and Future Directions While the space-charge solver achieved all performance metrics, the integrated fluid dynamics modeling encountered convergence challenges stemming from the extreme 200-fold disparity in length scales between the detector's 37 mm inlet pipes and the 8-meter global domain. To address this, the project has identified a clear technical pivot toward Hierarchical Geometric Adaptive Mesh Refinement (HG-AMR). By implementing an h-type refinement strategy with hanging nodes, future iterations of this solver will be capable of resolving localized high-gradient inlet flows without the prohibitive computational costs of regular grids. This advancement, combined with data-driven uncertainty quantification based on MicroBooNE-style calibration, will enable the precise modeling of detector responses in large-scale cryogenic environments where direct measurement remains difficult. Impact The computational tools developed under this award provide a foundation for enhancing the energy resolution and spatial reconstruction of noble liquid detectors. By bridging the gap between theoretical fluid dynamics and experimental field calibration, this work supports the DOE’s mission to advance the frontiers of neutrino physics and dark matter detection.

42 ENGINEERING↗
DOE OSTI2026

Overview of the MAST Upgrade physics programme: testing novel concepts at low aspect ratio to inform future devices

The research programme performed on the Mega Amp Spherical Tokamak (MAST) Upgrade experiment has made significant advances in developing the physics understanding of low aspect ratio tokamaks in support of the operation of ITER and design of fusion powerplants. High performance plasma scenarios have been developed to facilitate a broad programme of experiments, in which confinement is constrained by the presence of m/n = 2/1 modes that cause substantial losses of fast ions. The onset of these modes coincides with the q = 2 surface residing in a local minimum in the toroidal current density profile. The maximum electron temperature at the pedestal top, T e,ped is limited with gas fuelling to ∼350 eV to maintain regular ELMs; higher T e,ped results in a transition to a non-stationary ELM-free regime. The operational space of spherical tokamaks has been expanded into small and ELM-free regimes. Strong shaping of the last closed flux surface can induce a transition from large to small ELMs, and ELM suppression with resonant magnetic perturbations has been observed for the first time in a low aspect ratio tokamak. Negative triangularity shaping has induced a transition from ELMy H-mode to a high-performance L-mode regime for the first time in a low aspect ratio tokamak. In studies of fast ion confinement, losses of fast particles due to Global Alfvén Eigenmodes have been identified. Interactions between fast ions generated by off-axis neutral beam injection and thermal neutrals can result in significant losses of fast ions. Experiments with on- and off-axis neutral beam injection exhibit a flux pumping mechanism, where the central safety factor is held to ∼1 in the absence of sawteeth. In studies of pedestal physics, it has been found that elevated main chamber neutral pressures result in an increase in the electron density and reduction in the temperature at the pedestal top. Advances in understanding plasma exhaust include the integration of a high-performance plasma core with detached outer divertors in the X-point target configuration. A newly commissioned lower divertor cryopump reduces the lower divertor neutral pressure by up to 50%, with minimal effect on the main chamber or upper divertor. New measurements and SOLPS-ITER simulations emphasise the importance of plasma–neutral interactions on divertor detachment in the conditions accessible in experiments. Real-time control of the ionisation front location in both divertor chambers independently has been demonstrated in double null experiments, enabled by the tightly baffled divertor chambers.

MAST Upgrade↗
DOE OSTI2021

Limited-angle computed tomography with deep image and physics priors

Computed tomography is a well-established x-ray imaging technique to reconstruct the three-dimensional structure of objects. It has been used extensively in a variety of fields, from diagnostic imaging to materials and biological sciences. One major challenge in some applications, such as in electron or x-ray tomography systems, is that the projections cannot be gathered over all the angles due to the sample holder setup or shape of the sample. This results in an ill-posed problem called the limited angle reconstruction problem. Typical image reconstruction in this setup leads to distortion and artifacts, thereby hindering a quantitative evaluation of the results. To address this challenge, we use a generative model to effectively constrain the solution of a physics-based approach. Our approach is self-training that can iteratively learn the nonlinear mapping from partial projections to the scanned object. Because our approach combines the data likelihood and image prior terms into a single deep network, it is computationally tractable and improves performance through an end-to-end training. We also complement our approach with total-variation regularization to handle high-frequency noise in reconstructions and implement a solver based on alternating direction method of multipliers. We present numerical results for various degrees of missing angle range and noise levels, which demonstrate the effectiveness of the proposed approach.

47 OTHER INSTRUMENTATION↗
DOE OSTI2024

Approximate entropy analysis for nonlinear beam dynamics

In this paper, we applied an approximate entropy (ApEn) analysis to nonlinear beam dynamics. Due to the presence of strong nonlinear magnets, the degree of chaos in beam motion gradually increases proportionally with its amplitude. Such chaos can be quantitatively characterized with ApEn of beam turn-by-turn readings. ApEn analysis is a technique used to quantify the amount of regularity and the unpredictability of fluctuations with respect to time-series data. The ApEn, when applied as a chaos indicator, can then be used for nonlinear lattice optimization and analysis. The National Synchrotron Light Source II (NSLS-II) electron storage ring lattice was used as an example of a real-world application of our technique. Published by the American Physical Society 2024

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Solving Coupled Cluster Equations by the Newton Krylov Method

We describe using the Newton Krylov method to solve the coupled cluster equation. The method uses a Krylov iterative method to compute the Newton correction to the approximate coupled cluster amplitude. The multiplication of the Jacobian with a vector, which is required in each step of a Krylov iterative method such as the Generalized Minimum Residual (GMRES) method, is carried out through a finite difference approximation, and requires an additional residual evaluation. The overall cost of the method is determined by the sum of the inner Krylov and outer Newton iterations. We discuss the termination criterion used for the inner iteration and show how to apply pre-conditioners to accelerate convergence. We will also examine the use of regularization technique to improve the stability of convergence and compare the method with the widely used direct inversion of iterative subspace (DIIS) methods through numerical examples.

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

Optically induced umklapp shift currents in striped cuprates

Motivated by recent experiments that observed low-frequency second-order optical responses in doped striped superconductors, here we investigate the nonlinear electrodynamics of systems exhibiting a charge density wave (CDW) order parameter. Due to the Bragg scattering off the CDW order, an incoming spatially homogeneous electric field in addition to zero momentum current generates umklapp currents that are modulated in space at momenta of the reciprocal CDW lattice. In particular, here we predict and microscopically evaluate the umklapp shift current, a finite momentum analog of the regular shift current which represents the second-order optical process that downconverts homogeneous AC electric field into low-frequency, zero momentum current. Specifically, we evaluate real-time response functions within mean-field theory via the Keldysh technique and use the Peierls substitution to compute observables at finite momenta in lattice models. We find that systems with certain lattice symmetries (such as inversion symmetry), where the regular shift current is disallowed, may give rise to the umklapp one. We apply our framework to investigate lattice symmetries in layered materials with helical-like stripes and show that both types of shift currents provide insight into the nature of intertwined phases of matter. Finally, we discuss the relation of our findings to recent experiments in striped superconductors. Published by the American Physical Society 2024

Materials Science↗