Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Weak”

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 1,261 records · Page 70
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

ALESQP: An Augmented Lagrangian Equality-Constrained SQP Method for Optimization with General Constraints

Here we present a new algorithm for infinite-dimensional optimization with general constraints, called ALESQP. In short, ALESQP is an augmented Lagrangian method that penalizes inequality constraints and solves equality-constrained nonlinear optimization subproblems at every iteration. The subproblems are solved using a matrix-free trust-region sequential quadratic programming (SQP) method that takes advantage of iterative, i.e., inexact linear solvers, and is suitable for large-scale applications. A key feature of ALESQP is a constraint decomposition strategy that allows it to exploit problem-specific variable scalings and inner products. We analyze convergence of ALESQP under different assumptions. We show that strong accumulation points are stationary. Consequently, in finite dimensions ALESQP converges to a stationary point. In infinite dimensions we establish that weak accumulation points are feasible in many practical situations. Under additional assumptions we show that weak accumulation points are stationary. We present several infinite-dimensional examples where ALESQP shows remarkable discretization-independent performance in all of its iterative components, requiring a modest number of iterations to meet constraint tolerances at the level of machine precision. Also, we demonstrate a fully matrix-free solution of an infinite-dimensional problem with nonlinear inequality constraints.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

FROSch Preconditioners for Land Ice Simulations of Greenland and Antarctica

Numerical simulations of Greenland and Antarctic ice sheets involve the solution of large-scale highly nonlinear systems of equations on complex shallow geometries. This work is concerned with the construction of Schwarz preconditioners for the solution of the associated tangent problems, which are challenging for solvers mainly because of the strong anisotropy of the meshes and wildly changing boundary conditions that can lead to poorly constrained problems on large portions of the domain. In this study, two-level generalized Dryja--Smith--Widlund (GDSW)--type Schwarz preconditioners are applied to different land ice problems, i.e., a velocity problem, a temperature problem, as well as the coupling of the former two problems. We employ the message passing interface (MPI)--parallel implementation of multilevel Schwarz preconditioners provided by the package FROSch (fast and robust Schwarz) from the Trilinos library. The strength of the proposed preconditioner is that it yields out-of-the-box scalable and robust preconditioners for the single physics problems. To the best of our knowledge, this is the first time two-level Schwarz preconditioners have been applied to the ice sheet problem and a scalable preconditioner has been used for the coupled problem. The preconditioner for the coupled problem differs from previous monolithic GDSW preconditioners in the sense that decoupled extension operators are used to compute the values in the interior of the subdomains. Additionally, several approaches for improving the performance, such as reuse strategies and shared memory OpenMP parallelization, are explored. In our numerical study we target both uniform meshes of varying resolution for the Antarctic ice sheet as well as nonuniform meshes for the Greenland ice sheet. We present several weak and strong scaling studies confirming the robustness of the approach and the parallel scalability of the FROSch implementation. Among the highlights of the numerical results are a weak scaling study for up to 32,K processor cores (8,K MPI ranks and 4 OpenMP threads) and 566,M degrees of freedom for the velocity problem as well as a strong scaling study for up to 4,K processor cores (and MPI ranks) and 68,M degrees of freedom for the coupled problem.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

An Asymptotically Compatible Coupling Formulation for Nonlocal Interface Problems with Jumps

Here, we introduce a mathematically rigorous formulation for a nonlocal interface problem with jumps and propose an asymptotically compatible finite element discretization for the weak form of the interface problem. After proving the well-posedness of the weak form, we demonstrate that solutions to the nonlocal interface problem converge to the corresponding local counterpart when the nonlocal data are appropriately prescribed. Several numerical tests in one and two dimensions show the applicability of our technique, its numerical convergence to exact nonlocal solutions, its convergence to the local limit when the horizons vanish, and its robustness with respect to the patch test.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Thermodynamics conditions of matter in the neutrino decoupling region during neutron star mergers

In this work we investigate the thermodynamics conditions at which neutrinos decouple from matter in neutron star merger remnants by post-processing results of merger simulations. We find that the matter density and the neutrino energies are the most relevant quantities in determining the decoupling surface location. For mean energy neutrinos (~9, 15 and 24 MeV for ν e , $\bar{ν}$ e and ν μ,τ , respectively) the transition between diffusion and free-streaming conditions occurs around 10 11 g cm -3 for all neutrino species. Weak and thermal equilibrium freeze-out occurs deeper (several 10 12 g cm -3 ) for heavy-flavor neutrinos than for $\bar{ν}$ e and ν e (≳10 11 g cm -3 ). Decoupling temperatures are broadly in agreement with the average neutrino energies, with softer equations of state characterized by ~1 MeV larger decoupling temperatures. Neutrinos streaming at infinity with different energies come from different remnant parts. While low-energy neutrinos (~3 MeV) decouple at ρ~10 13 g cm -3 , T~10 MeV and Y e ≲0.1 close to weak equilibrium, high-energy ones (~50 MeV ) decouple from the disk at ρ~10 9 g cm -3 , T~2 MeV and Y e ≳0.25. The presence of a massive NS or a BH influences the neutrino thermalization. While in the former case decoupling surfaces are present for all relevant energies, the lower maximum density ( ≲10 12 g cm -3 ) in BH-torus systems does not allow softer neutrinos to thermalize and diffuse.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Wigner SU(4) symmetry, clustering, and the spectrum of $$^{12}$$C

Abstract We present lattice calculations of the low-lying spectrum of $$^{12}$$ 12 C using a simple nucleon–nucleon interaction that is independent of spin and isospin and therefore invariant under Wigner’s SU(4) symmetry. We find strong signals for all excited states up to $$\sim 15$$ ∼ 15 MeV above the ground state, and explore the structure of each state using a large variety of $$\alpha $$ α cluster and harmonic oscillator trial states, projected onto given irreducible representations of the cubic group. We are able to verify earlier findings for the $$\alpha $$ α clustering in the Hoyle state and the second $$2^+$$ 2 + state of $$^{12}$$ 12 C. The success of these calculations to describe the full low-lying energy spectrum using spin-independent interactions suggest that either the spin-orbit interactions are somewhat weak in the $$^{12}$$ 12 C system, or the effects of $$\alpha $$ α clustering are diminishing their influence. This is in agreement with previous findings from ab initio shell model calculations.

Shen, Shihang (ORCID:0000000280516466)↗
DOE OSTI2022

On the stability of the open-string QED neutron and dark matter

We study the stability of a hypothetical QED neutron, which consists of a color-singlet system of two d quarks and a u quark interacting with the QED interaction. As a quark cannot be isolated, the intrinsic motion of the three quarks in the lowest-energy state may lie predominantly in 1 + 1 dimensions, as in a d-u-d open string. The attractive d-u and u-d QED interactions may overcome the weaker repulsive d-d QED interaction to bind the three quarks together. We examine the QED neutron in a phenomenological three-body problem in 1 + 1 dimensions with an effective interaction extracted from Schwinger’s exact QED solution in 1 + 1 dimensions. The phenomenological model in a variational calculation yields a stable QED neutron at 44.5 MeV. The analogous QED proton with two u quarks and a d quark has been found to be too repulsive to be stable and does not have a bound or continuum state, onto which the QED neutron can decay via the weak interaction. Consequently, the QED neutron is stable against the weak decay, has a long lifetime, and is in fact a QED dark neutron. It may be produced following the deconfinement-to-confinement phase transition of the quark gluon plasma in high-energy heavy-ion collisions. Because of the long lifetime of the QED dark neutron, self-gravitating assemblies of QED dark neutrons or dark antineutrons may be good candidates for a part of the primordial dark matter produced during the phase transition of the quark gluon plasma in the evolution of the early Universe.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

The s process in massive stars, a benchmark for neutron capture reaction rates

A clear definition of the contribution from the slow neutron-capture process (s process) to the solar abundances between Fe and the Sr-Zr region is a crucial challenge for nuclear astrophysics. Robust s-process predictions are necessary to disentangle the contribution from other stellar processes producing elements in the same mass region. Nuclear uncertainties are affecting s-process calculations, but most of the needed nuclear input are accessible to present nuclear experiments or they will be in the near future. Neutron-capture rates have a great impact on the s process in massive stars, which is a fundamental source for the solar abundances of the lighter s-process elements heavier than Fe (weak s-process component). In this work we present a new nuclear sensitivity study to explore the impact on the s process in massive stars of 86 neutron-capture rates, including all the reactions between C and Si and between Fe and Zr. We derive the impact of the rates at the end of the He-burning core and at the end of the C-burning shell, where the 22 Ne(α,n) 25 Mg reaction is is the main neutron source. We confirm the relevance of the light isotopes capturing neutrons in competition with the Fe seeds as a crucial feature of the s process in massive stars. For heavy isotopes we study the propagation of the neutron-capture uncertainties, finding a clear difference of the impact of Fe and Co isotope rates with respect to the rates of heavier stable isotopes. The local uncertainty propagation due to the neutron-capture rates at the s-process branching points is also considered, discussing the example of 85 Kr. The complete results of our study for all the 86 neutron-capture rates are available online. Finally, we present the impact on the weak s process of the neutron-capture rates included in the new ASTRAL library (v0.2).

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Lifetimes of low-lying levels in 158 Gd

The low-lying structure of the well-deformed nucleus 158 Gd has been revisited to elucidate the nature of the low-lying states in 158 Gd. Earlier (p, t) studies identified numerous 0 + states below 4.3 MeV, prompting questions about whether these states correspond to collective vibrations or shape coexistence. New and previously reported (n, n' γ) measurements are combined, including γ-γ coincidences, excitation functions, and angular distributions, to extract lifetimes and transition probabilities for 44 excited states up to 2.7 MeV, including 32 previously unmeasured levels. Our results confirm or revise γ-ray placements and provide detailed transition strengths, revealing both weakly collective and strongly enhanced B(E2) and B(E1) transition probabilities. In particular, a tentative 0 + state at 2437.8 keV exhibits a strong interband B(E2) transition, which may be a candidate for a possible two-phonon (ββ) excitation. Systematic comparisons with neighboring Gd isotopes, Hartree–Fock–Bogoliubov, and interacting-boson model predictions suggest that the first excited 0 + state in 158 Gd is predicted to be a β-vibration, although it is weakly collective.

158Gd↗
DOE OSTI2019

Is the magnitude of the Peccei–Quinn scale set by the landscape?

Rather general considerations of the string theory landscape imply a mild statistical draw towards large soft SUSY breaking terms tempered by the requirement of proper electroweak symmetry breaking where SUSY contributions to the weak scale are not too far from m weak ~ 100 GeV. Such a picture leads to the prediction that m h ~ 125 GeV while most sparticles are beyond current LHC reach. Here we explore the possibility that the magnitude of the Peccei-Quinn (PQ) scale f a is also set by string landscape considerations within the framework of a compelling SUSY axion model. First, we examine the case where the PQ symmetry arises as an accidental approximate global symmetry from a more fundamental gravity-safe Z$R\atop{24}$$ symmetry and where the SUSY mu parameter arises from a Kim-Nilles operator. The pull towards large soft terms then also pulls the PQ scale as large as possible. Unless this is tempered by rather severe (unknown) cosmological or anthropic bounds on the density of dark matter, then we would expect a far greater abundance of dark matter than is observed. This conclusion cannot be negated by adopting a tiny axion misalignment angle θ i because WIMPs are also overproduced at large f a . Hence, we conclude that setting the PQ scale via anthropics is highly unlikely. Instead, requiring soft SUSY breaking terms of order the gravity-mediation scale m 3/2 ~ 10 -100 TeV places the mixed axion-neutralino dark matter abundance into the intermediate scale sweet zone where f a ~ 10 11 -10 12 GeV. We compare our analysis to the more general case of a generic SUSY DFSZ axion model with uniform selection on theta_i but leading to the measured dark matter abundance: this approach leads to a preference for f a ~ 10 12 GeV.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Differential cross-section measurements for the electroweak production of dijets in association with a Z boson in proton–proton collisions at ATLAS

Differential cross-section measurements are presented for the electroweak production of two jets in association with a Z boson. These measurements are sensitive to the vector-boson fusion production mechanism and provide a fundamental test of the gauge structure of the Standard Model. The analysis is performed using proton–proton collision data collected by ATLAS at √s=13 TeV and with an integrated luminosity of 139 fb -1 . The differential cross-sections are measured in the Z→ℓ + ℓ - decay channel (ℓ=e,μ) as a function of four observables: the dijet invariant mass, the rapidity interval spanned by the two jets, the signed azimuthal angle between the two jets, and the transverse momentum of the dilepton pair. The data are corrected for the effects of detector inefficiency and resolution and are sufficiently precise to distinguish between different state-of-the-art theoretical predictions calculated using Powheg+Pythia8, Herwig7+Vbfnlo and Sherpa 2.2. The differential cross-sections are used to search for anomalous weak-boson self-interactions using a dimension-six effective field theory. The measurement of the signed azimuthal angle between the two jets is found to be particularly sensitive to the interference between the Standard Model and dimension-six scattering amplitudes and provides a direct test of charge-conjugation and parity invariance in the weak-boson self-interactions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Neutrinos in Stellar Astrophysics (in The Encyclopedia of Cosmology)

The physics of the mysterious and stealthy neutrino is at the heart of many phenomena in the cosmos. These particles interact with matter and with each other through the aptly named weak interaction. At typical astrophysical energies the weak interaction is some twenty orders of magnitude weaker than the electromagnetic interaction. However, in the early universe and in collapsing stars neutrinos can more than make up for their feeble interaction strength with huge numbers. Neutrinos can dominate the dynamics in these sites and set the conditions that govern the synthesis of the elements. Here we journey through the history of the discovery of these particles and describe their role in stellar evolution and collapse, the big bang, and multi-messenger astrophysics. Neutrino physics is at the frontier of elementary particle physics, nuclear physics, astrophysics and cosmology. All of these fields overlap in the neutrino story.

Haxton, Wick↗
DOE OSTI2023

Emergent Halperin–Saslow mode, gauge glass and quenched disorders in quantum Ising magnet TmMgGaO 4

Quenched disorders could bring novel quantum excitations and models to certain quantum magnets. Motivated by recent experiments on the quantum Ising magnet TmMgGaO 4 , we explore the effects of the quenched disorder and the interlayer coupling in this triangular lattice Ising antiferromagnet. It is pointed out that the weak quenched (nonmagnetic) disorder would convert the emergent 2D Berezinskii–Kosterlitz–Thouless (BKT) phase and the critical region into a U(1) gauge glass. There will be an emergent Halperin–Saslow mode associated with this gauge glass. Using the Imry-Ma’s renormalization group result, we explain the fate of the finite-field [Formula: see text] symmetry breaking transition at the low temperatures. The ferromagnetic interlayer coupling would suppress the BKT phase and generate a tiny ferromagnetism. With quenched disorders, this interlayer coupling changes the 2D gauge glass into a 3D gauge glass, and the Halperin–Saslow mode persists. This work merely focuses on addressing a phase regime in terms of emergent U(1) gauge glass behaviors and hope to inspire future works and thoughts in weakly disordered frustrated magnets in general.

Physics↗
DOE OSTI2021

On a nonlocal Cahn–Hilliard model permitting sharp interfaces

A nonlocal Cahn–Hilliard model with a non-smooth potential of double-well obstacle type that promotes sharp interfaces in the solution is presented. To capture long-range interactions between particles, a nonlocal Ginzburg–Landau energy functional is defined which recovers the classical (local) model as the extent of nonlocal interactions vanish. In contrast to the local Cahn–Hilliard problem that always leads to diffuse interfaces, the proposed nonlocal model can lead to a strict separation into pure phases of the substance. In this work, the lack of smoothness of the potential is essential to guarantee the aforementioned sharp-interface property. Mathematically, this introduces additional inequality constraints that, in a weak formulation, lead to a coupled system of variational inequalities which at each time instance can be restated as a constrained optimization problem. We prove the well-posedness and regularity of the semi-discrete and continuous in time weak solutions, and derive the conditions under which pure phases are admitted. Moreover, we develop discretizations of the problem based on finite element methods and implicit–explicit time-stepping methods that can be realized efficiently. Finally, we illustrate our theoretical findings through several numerical experiments in one and two spatial dimensions that highlight the differences in features of local and nonlocal solutions and also the sharp interface properties of the nonlocal model.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Sphaleron transition rates and the chiral magnetic effect

The chiral magnetic effect is a novel quantum phenomenon proposed for high-energy nuclear collisions but which has yet to be observed. We quantify the axial charge relaxation time, due to sphalerons, which enters in simulations of this effect. An extrapolation of weak coupling calculations of the sphaleron rate yields rather different relaxation times than strong coupling AdS/CFT calculations. The AdS/CFT relaxation time is the larger one of the two by an order of magnitude, but the weak coupling relaxation time may not be reliable because it is only marginally bigger than the microscopic thermalization time. The role of quark masses has yet to be accurately assessed.

Physics↗
DOE OSTI2020

Performance Evaluation of a Two-Dimensional Flood Model on Heterogeneous High-Performance Computing Architectures

This paper describes the implementation of a two-dimensional hydrodynamic flood model with two different numerical schemes on heterogeneous high-performance computing architectures. Both schemes were able to solve the nonlinear hyperbolic shallow water equations using an explicit upwind first-order approach on finite differences and finite volumes, respectively, and were conducted using MPI and CUDA. Four different test cases were simulated on the Summit supercomputer at Oak Ridge National Laboratory. Both numerical schemes scaled up to 128 nodes (768 GPUs) with a maximum 98.2x speedup of over 1 GPU. The lowest run time for the 10 day Hurricane Harvey event simulation at 5 meter resolution (272 million grid cells) was 50 minutes. GPUDirect communication proved to be more convenient than the standard communication strategy. Both strong and weak scaling are shown.

Sharif, Md Bulbul↗
DOE OSTI2024

Dark Matter Candidates of a Very Low Mass

We review dark matter (DM) candidates of a very low mass appearing in the window below the traditional weakly interacting massive particle (m χ ≲ 10 GeV) and extending down to m χ ≳ 1 meV, somewhat below the mass limit at which DM becomes wavelike. Such candidates are motivated by hidden sectors such as hidden valleys, which feature hidden forces and rich dynamics, but have evaded traditional accelerator searches for New Physics because of their relatively weak coupling to the Standard Model (SM). Such sectors can still be detected through dedicated low-energy colliders, which, through their intense beams, can have sensitivity to smaller couplings, or through astrophysical observations of the evolution of DM halos and stellar structures, which, through the Universe's epochs, can be sensitive to small DM interactions. We also consider mechanisms whereby the DM abundance is fixed through the interaction with the SM, which directly motivates the search for light DM in terrestrial experiments. The bulk of this review is dedicated to the new ideas that have been proposed for direct detection of such DM candidates of a low mass through nuclear recoils, electronic excitations, or collective modes such as phonons and magnons. The rich tapestry of materials and modes in the condensed matter landscape is reviewed along with specific prospects for detection.

Physics↗
DOE OSTI2024

Spectroscopy in Nanoscopic Cavities: Models and Recent Experiments

The ability of nanophotonic cavities to confine and store light to nanoscale dimensions has important implications for enhancing molecular, excitonic, phononic, and plasmonic optical responses. Spectroscopic signatures of processes that are ordinarily exceedingly weak such as pure absorption and Raman scattering have been brought to the single-particle limit of detection, while new emergent polaritonic states of optical matter have been realized through coupling material and photonic cavity degrees of freedom across a wide range of experimentally accessible interaction strengths. In this review, we discuss both optical and electron beam spectroscopies of cavity-coupled material systems in weak, strong, and ultrastrong coupling regimes, providing a theoretical basis for understanding the physics inherent to each while highlighting recent experimental advances and exciting future directions.

Chemistry↗
DOE OSTI2019

Reduced Model for Direct Induction Startup Scenario Development on MAST-U and NSTX-U

A reduced semi-empirical model using time-dependent axisymmetric vacuum field calculations is used to develop the prefill and feed-forward coil current targets required for reliable direct induction (DI) startup on the new MA-class spherical tokamaks, MAST-U and NSTX-U. The calculations are constrained by operational limits unique to each device, such as the geometry of the conductive elements and active coils, power supply specifications and coil heating and stress limits. The calculations are also constrained by semi-empirical models for sufficient breakdown, current drive, equilibrium and stability of the plasma developed from a shared database. A large database of DI startup on NSTX and NSTX-U is leveraged to quantify the requirements for achieving a reliable breakdown (Ip ~ 20 kA). It is observed that without pre-ionization, STs access the large E/P regime at modest loop voltage (Vloop) where the electrons in the weakly ionized plasma are continually accelerating along the open field lines. This ensures a rapid (order millisecond) breakdown of the neutral gas, even without pre-ionization or high-quality field nulls. The timescale of the initial increase in Ip on NSTX is reproduced in the reduced model provided a mechanism for impeding the applied electric field is included. Most discharges that fail in the startup phase are due to an inconsistency in the evolution of the plasma current (Ip) and equilibrium field or loss of vertical stability during the burn-through phase. The requirements for the self-consistent evolution of the fields in the weakly and full-ionized plasma states are derived from demonstrated DI startup on NSTX, NSTX-U and MAST. The predictive calculations completed for MAST-U and NSTX-U illustrate that the maximum Ip ramp rate (dIp/dt) in the early startup phase is limited by the voltage limits on the poloidal field coils on MAST-U and passive vertical stability on NSTX-U.

Direct induction↗