Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “asymptotic analysis”

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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Systematic Local Simulations of Fast Neutrino Flavor Conversions with Scattering Effects

We investigate the dynamics of fast neutrino flavor conversions (FFCs) in the one-dimensional (1D) and zero-dimensional (0D) models, in which spatial advection is considered and ignored, respectively. In this study, we employ snapshots obtained by our self-consistent, realistic Boltzmann-neutrino-radiation-hydrodynamics simulations. We show that the FFC growth rate is considerably larger in the 1D model than in the 0D model, as expected from the previous linear analysis results. We find that the momentum space dimension does not significantly influence the neutrino transition probability in 1D models. On the other hand, in the 0D model without collisions, the FFC depends on the momentum space, and the azimuthal angle dependence breaks the periodicity of the FFC. Our study demonstrates that collisional instability can lead to further flavor conversions on a long timescale in 1D models after the asymptotic state of FFC has been reached. Such an effect should be taken into consideration when the fast and collisional flavor instabilities coexist.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Small instantons in weakly-gauged holographic models

Small instantons can play an important role in Yang-Mills theories whose gauge couplings are sizeable at small distances. An interesting class of theories where this could occur is in weakly-gauged holographic models (dual to Yang-Mills theories interacting with strongly-coupled CFTs), since gauge couplings are indeed enhanced towards the UV boundary of the 5D AdS space. However, contrary to expectations, we show that small instantons in these non-asymptotically-free models are highly suppressed and ineffective. This is due to the conservation of topological charge that forbids instantons to be localized near the UV boundary. Despite this fact we find non-trivial UV localized instanton-anti-instanton solutions of the Yang-Mills equations where the topological charges annihilate in the AdS bulk. These analytic solutions arise from a 5D conformal transformation of the uplifted 4D instanton. Our analysis therefore reveals unexpected nonperturbative configurations of Yang-Mills theories when they interact with strongly-coupled CFTs.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Scalar modular bootstrap and zeros of the Riemann zeta function

Using the technology of harmonic analysis, we derive a crossing equation that acts only on the scalar primary operators of any two-dimensional conformal field theory with U(1) c symmetry. From this crossing equation, we derive bounds on the scalar gap of all such theories. Rather remarkably, our crossing equation contains information about all nontrivial zeros of the Riemann zeta function. As a result, we rephrase the Riemann hypothesis purely as a statement about the asymptotic density of scalar operators in certain two-dimensional conformal field theories. We discuss generalizations to theories with only Virasoro symmetry.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Implications of inertial subrange scaling for stably stratified mixing

In this work, we investigate the effects of the turbulent dynamic range on active scalar mixing in stably stratified turbulence by adapting the theoretical passive scalar modelling and demonstrating their usefulness through consideration of the results of direct numerical simulations of statistically stationary homogeneous stratified and sheared turbulence. By analysis of inertial and inertial–convective subrange scalings, we show that the relationship between the active scalar and turbulence time scales is predicted by the ratio of the Kolmogorov and Obukhov–Corrsin constants, provided mean flow parameters permit the two subrange scalings to be appropriate approximations. We use the resulting relationship between time scales to parameterise an appropriate turbulent mixing coefficient $\varGamma \equiv \chi /\epsilon$, defined here as the ratio of available potential energy ($E_p$) and turbulent kinetic energy ($E_k$) dissipation rates. With the analysis presented here, we show that $\varGamma$ can be estimated by $E_p,E_k$ and a universal constant provided an appropriate Reynolds number is sufficiently high. This large Reynolds number regime appears here to occur at $ {{Re_b}} \equiv \epsilon / \nu N^{2} \gtrapprox 300$ where $\nu$ is the kinematic viscosity and $N$ is the characteristic buoyancy frequency. We propose a model framework for irreversible diapycnal mixing with robust theoretical parametrisation and asymptotic behaviour in this high-$ {{Re_b}}$ limit.

42 ENGINEERING↗
DOE OSTI2023

A data-driven method for modelling dissipation rates in stratified turbulence

We present a deep probabilistic convolutional neural network (PCNN) model for predicting local values of small-scale mixing properties in stratified turbulent flows, namely the dissipation rates of turbulent kinetic energy and density variance, $\varepsilon$ and $\chi$ . Inputs to the PCNN are vertical columns of velocity and density gradients, motivated by data typically available from microstructure profilers in the ocean. The architecture is designed to enable the model to capture several characteristic features of stratified turbulence, in particular the dependence of small-scale isotropy on the buoyancy Reynolds number $Re_b:=\varepsilon /(\nu N^2)$ , where $\nu$ is the kinematic viscosity and $N$ is the background buoyancy frequency, the correlation between suitably locally averaged density gradients and turbulence intensity and the importance of capturing the tails of the probability distribution functions of values of dissipation. Empirically modified versions of commonly used isotropic models for $\varepsilon$ and $\chi$ that depend only on vertical derivatives of density and velocity are proposed based on the asymptotic regimes $Re_b\ll 1$ and $Re_b\gg 1$ , and serve as an instructive benchmark for comparison with the data-driven approach. When trained and tested on a simulation of stratified decaying turbulence which accesses a range of turbulent regimes (associated with differing values of $Re_b$ ), the PCNN outperforms assumptions of isotropy significantly as $Re_b$ decreases, and additionally demonstrates improvements over the fitted empirical models. A differential sensitivity analysis of the PCNN facilitates a comparison with the theoretical models and provides a physical interpretation of the features enabling it to make improved predictions.

42 ENGINEERING↗
DOE OSTI2020

Population Dynamics and Temporal Trends of Bull Trout in the East Fork Salmon River, Idaho

Abstract Because of their long-term listing under the Endangered Species Act, much interest has been placed on estimating population vital rates for Bull Trout Salvelinus confluentus, but the biotic and abiotic factors that influence the interannual variability in those vital rates have rarely been evaluated. We used mark–recapture data to estimate fish growth, survival, and trends in abundance for fluvial adult Bull Trout in the East Fork Salmon River, Idaho. Over an 8-year period, a total of 1,205 individual Bull Trout were collected at a weir on the East Fork Salmon River (29 km upstream of its confluence with the Salmon River) during June–September, of which 420 were recaptures from prior years. Bull Trout varied in length from 215 to 756 mm and achieved a slightly larger asymptotic length and a slightly lower rate of growth relative to other fluvial and adfluvial Bull Trout populations. Apparent survival averaged 0.42 across all years, which was similar to previous studies estimating apparent survival for Bull Trout. The number of emigrating anadromous salmonid smolts in the upper Salmon River basin positively influenced East Fork Salmon River Bull Trout growth and survival, and survival was higher in years with lower annual discharge. Assessment of population growth via linear regression analysis indicated that the Bull Trout population was increasing during the study period (λ = 1.08; 95% confidence interval = 1.03–1.14). Our findings highlight the ecological link between the abundance of wild and hatchery Chinook Salmon Oncorhynchus tshawytscha, Sockeye Salmon O. nerka, and steelhead O. mykiss smolts and growth and survival of adult Bull Trout in systems where these species occur in sympatry.

Roth, Curtis J.↗
DOE OSTI2023

Machine Learning-Based Model Predictive Control of Two-Time-Scale Systems

In this study, we present a general form of nonlinear two-time-scale systems, where singular perturbation analysis is used to separate the dynamics of the slow and fast subsystems. Machine learning techniques are utilized to approximate the dynamics of both subsystems. Specifically, a recurrent neural network (RNN) and a feedforward neural network (FNN) are used to predict the slow and fast state vectors, respectively. Moreover, we investigate the generalization error bounds for these machine learning models approximating the dynamics of two-time-scale systems. Next, under the assumption that the fast states are asymptotically stable, our focus shifts toward designing a Lyapunov-based model predictive control (LMPC) scheme that exclusively employs the RNN to predict the dynamics of the slow states. Additionally, we derive sufficient conditions to guarantee the closed-loop stability of the system under the sample-and-hold implementation of the controller. A nonlinear chemical process example is used to demonstrate the theory. In particular, two RNN models are constructed: one to model the full two-time-scale system and the other to predict solely the slow state vector. Both models are integrated within the LMPC scheme, and we compare their closed-loop performance while assessing the computational time required to execute the LMPC optimization problem.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

On the hot surface ignition of a wall-stagnating spray flame

The ignition of fuel sprays due to interaction with hot surfaces is an important phenomenon in the safety analysis of many engineering systems. We perform a parametric study of the hot surface ignition (HSI) of a fuel spray approaching a heated surface caused by the accidental leakage of a fuel line. Here, to this end, we employ a one-dimensional Eulerian-Eulerian formulation with a non-equilibrium evaporation model and a realistic chemical mechanism to describe $n$-dodecane fuel chemistry. We first describe and analyze the phenomenology of the unsteady processes leading to ignition using non-dimensionalized quantities. Through consideration of the temporal development of the most reactive mixture, we demonstrate that ignition occurs at a fuel-lean composition in a premixed region near the hot surface. Using non-dimensional parameters identified from the governing equations, we perform a parametric study of the time, location and local mixture composition at ignition and determine the ignition limits. We then identify the most important parametric sensitivities for physical analysis using a data-driven classification method. Our analysis demonstrates a contraction of the ignition limits with increased Stokes number and a regime of parametric insensitivity of igniting mixture composition. We also show that at high Damköhler numbers, the ignition location conforms to the parametric behavior of the thermal boundary layer, whereas at low Damköhler numbers approaching the ignition limit it reaches a near-unity value of the quenching Peclet number. We then compare the demonstrated parametric dependencies to the results of the quasi-steady asymptotic ignition literature, showing that our results are consistent with those obtained analytically within the limitations imposed by the simplified formulation of the latter.

42 ENGINEERING↗
DOE OSTI2022

Running Primal-Dual Gradient Method for Time-Varying Nonconvex Problems

This paper focuses on a time-varying constrained nonconvex optimization problem, and considers the synthesis and analysis of online regularized primal-dual gradient methods to track a Karush-Kuhn-Tucker (KKT) trajectory. The proposed regularized primal-dual gradient method is implemented in a running fashion, in the sense that the underlying optimization problem changes during the execution of the algorithms. In order to study its performance, we first derive its continuous-time limit as a system of differential inclusions. We then study sufficient conditions for tracking a KKT trajectory, and also derive asymptotic bounds for the tracking error (as a function of the time-variability of a KKT trajectory). Further, we provide a set of sufficient conditions for the KKT trajectories not to bifurcate or merge, and also investigate the optimal choice of the parameters of the algorithm. Illustrative numerical results for a time-varying nonconvex problem are provided.

differential inclusion↗
DOE OSTI2022

Asymptotically Free AdS 3 /CFT 2

We propose a new AdS 3 /CFT 2 duality, in which the bulk string theory has a target spacetime AdS 3 times a squashed three-sphere $S$$^{3}_{b}$, and the dual CFT 2 is a symmetric product of sigma models on R Φ × $S$$^{3}_{b}$, deformed by a Φ-dependent Z 2 twist operator. The duality maps the asymptotic region of AdS 3 to the region Φ → ∞, where the twist interaction in the CFT 2 turns off. The AdS 3 backgrounds in question have R AdS < ℓ s , and so lie on the string side of the string/black hole correspondence transition. As a consequence, the high energy density of states consists of a string gas in AdS 3 rather than an ensemble of BTZ black holes. Furthermore, this property allows us to derive the dual CFT 2 by a systematic analysis of the worldsheet string theory on AdS 3 .

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Generalized parton distributions and gravitational form factors at large momentum transfer

Within the soft collinear effective theory (SCET), we derive a factorization theorem which resums Sudakov logarithms (α s ln 2 ( –t)) n to all orders in the quark-in-quark generalized parton distribution (GPD) at large momentum transfer t, and perform a consistency check to one-loop. We show that the same Sudakov factor appears in the ‘Feynman’ contribution to the GPDs of the nucleon. Our result enables the resummation of all the large logarithms ln Q 2 and ln 2 t in exclusive processes with two hard scales Λ$^{2}_{QCD}$ ≪ |t| ≪ Q 2 . We also present a SCET power counting analysis of the Feynman contributions to the GPDs and show that the x-dependence of GPDs factorizes at large-t with controlled corrections. This in particular implies that any ratio of GPD moments such as the electromagnetic and gravitational form factors (GFF) is perturbatively calculable in this approximation. Furthermore, we identify a novel order α s power-law t-dependence in the GPD and the D-type GFF that will dominate over the standard order α$^{2}_{s}$ ‘leading twist’ asymptotic contribution in the phenomenologically relevant region of t.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

An asymptotic-preserving 2D-2P relativistic Drift-Kinetic-Equation solver for runaway electron simulations in axisymmetric tokamaks

We propose an asymptotic-preserving (AP), uniformly convergent numerical scheme for the relativistic collisional Drift-Kinetic Equation (rDKE) to simulate runaway electrons in axisymmetric toroidal magnetic field geometries typical of tokamak devices. The approach is derived from an exact Green's function solution with numerical approximations of quantifiable impact, and results in a simple, two-step operator-split algorithm, consisting of a collisional Eulerian step, and a Lagrangian orbit-integration step with analytically prescribed kernels. The AP character of the approach is demonstrated by analysis of the dominant numerical errors, as well as by numerical experiments. We demonstrate the ability of the algorithm to provide accurate answers regardless of plasma collisionality on a circular axisymmetric tokamak geometry.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

The Average Shape of Sea Ice Ridge Keels

Through analysis of over 64,000 ridge profiles identified from moored upward-looking sonars, we identify a well-defined average shape of pressure ridge keels that is concave or cusped, not triangular as widely assumed in other literature. On the basis of this average shape, we put forward a new, dimensional-definition, of a pressure ridge cross-section that follows a negative exponential form and allows an average pressure ridge cross section to be constructed with knowledge or choice of a single parameter. The horizonal asymptote of the profile represents the draft of the “background” ice in which the ridge is embedded. The draft of this background ice scales with keel depth and is typically greater than can be accounted through thermodynamic growth, indicating that ridges tend to be embedded in fields of mechanically thickened rubble. Using a variational ridge model we simulated keel shapes for a range of ridge building conditions. Here, the model results agree well with the observations and indicate the cusped shape of an average ridge profile arises from the varying angle of horizontal shear in the ice cover when ridges form. The modeling results also explain the elongated tail of the area draft distribution of the pack.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Unified wavefront singularity characterization of three-dimensional elastodynamic time-domain half-space Green's function under impulsive boundary and internal loads

Founded on a novel analytical formulation that led to a rigorous yet compact path-integral representation of the time-domain elastodynamic half-space Green's function, a unified analysis of the possible occurrence of different singular wavefront behaviour in the response under arbitrary impulsive internal or surface point loads at arbitrary source-receiver locations is presented. With the decomposition of the general solution into distinct initiating and reflected wave group integrals that share a common factored format and simple contour definitions, the mathematical framework is shown to allow a straightforward identification of the specific conditions and the particular wave groups that are responsible for the singular wavefront phenomena without resorting to advanced analytic function theories or asymptotic methods. Analytic characterizations of the nature, strength and direction of all intrinsic singular wavefront behaviours of the three-dimensional Green's function in three canonical cases of source-receiver configurations are given in a dual integral-closed form format to facilitate their theoretical understanding as well as computational applications. Graphical illustrations of their variation with the source-receiver configuration and the medium's Poisson's ratio together with relevant comparison and clarifications of some classical treatments are included.

Green’s functions↗
DOE OSTI2023

Effective field theory analysis of the Coulomb breakup of the one-neutron halo nucleus $$^{19}$$C

Abstract We analyse the Coulomb breakup of $$^{19}$$ 19 C measured at 67 A MeV at RIKEN. We use the Coulomb-Corrected Eikonal (CCE) approximation to model the reaction and describe the one-neutron halo nucleus $$^{19}$$ 19 C within Halo Effective Field Theory (Halo EFT). At leading order we obtain a fair reproduction of the measured cross section as a function of energy and angle. The description is insensitive to the choice of optical potential, as long as it accurately represents the size of $${}^{18}$$ 18 C. It is also insensitive to the interior of the $${}^{19}$$ 19 C wave function. Comparison between theory and experiment thus enables us to infer asymptotic properties of the ground state of $$^{19}$$ 19 C: these data put constraints on the one-neutron separation energy of this nucleus and, for a given binding energy, can be used to extract an asymptotic normalisation coefficient (ANC). These results are confirmed by CCE calculations employing next-to-leading order Halo EFT descriptions of $${}^{19}$$ 19 C: at this order the results for the Coulomb breakup cross section are completely insensitive to the choice of the regulator. Accordingly, this reaction can be used to constrain the one-neutron separation energy and ANC of $${}^{19}$$ 19 C.

Physics↗
DOE OSTI2025

Convergence Analysis of the Alternating Anderson–Picard Method for Nonlinear Fixed-Point Problems

Anderson acceleration (AA) has been widely used to solve nonlinear fixed-point problems due to its rapid convergence. This work focuses on a variant of AA in which multiple Picard iterations are performed between each AA step, referred to as the Alternating Anderson–Picard (AAP) method. Furthermore, despite introducing more “slow” Picard iterations, this method has been shown to be efficient and even more robust in both linear and nonlinear cases. However, there is a lack of theoretical analysis for AAP in the nonlinear case. In this paper, we address this gap by establishing the equivalence between AAP and a multisecant-GMRES method that uses GMRES to solve a multisecant linear system at each iteration. From this perspective, we show that AAP “converges” to the Newton-GMRES method. Specifically, as the residual approaches zero, the multisecant matrix, the approximate Jacobian inverse, the search direction, and the optimization gain of AAP converge to their counterparts in the Newton-GMRES method. These connections provide insights for analyzing the asymptotic convergence properties of AAP. Consequently, we show that AAP is locally 𝑞-linear convergent and provide an upper bound for the convergence factor of AAP. To validate the theoretical results, numerical examples are provided.

Anderson acceleration↗
DOE OSTI2022

The RADIOSTAR Project

Radioactive nuclei are the key to understanding the circumstances of the birth of our Sun because meteoritic analysis has proven that many of them were present at that time. Their origin, however, has been so far elusive. The ERC-CoG-2016 RADIOSTAR project is dedicated to investigating the production of radioactive nuclei by nuclear reactions inside stars, their evolution in the Milky Way Galaxy, and their presence in molecular clouds. So far, we have discovered that: (i) radioactive nuclei produced by s l o w ( 107 Pd and 182 Hf) and r a p i d ( 129 I and 247 Cm) neutron captures originated from stellar sources &mdash;asymptotic giant branch (AGB) stars and compact binary mergers, respectively&mdash;within the galactic environment that predated the formation of the molecular cloud where the Sun was born; (ii) the time that elapsed from the birth of the cloud to the birth of the Sun was of the order of 10 7 years, and (iii) the abundances of the very short-lived nuclei 26 Al, 36 Cl, and 41 Ca can be explained by massive star winds in single or binary systems, if these winds directly polluted the early Solar System. Our current and future work, as required to finalise the picture of the origin of radioactive nuclei in the Solar System, involves studying the possible origin of radioactive nuclei in the early Solar System from core-collapse supernovae, investigating the production of 107 Pd in massive star winds, modelling the transport and mixing of radioactive nuclei in the galactic and molecular cloud medium, and calculating the galactic chemical evolution of 53 Mn and 60 Fe and of the p -process isotopes 92 Nb and 146 Sm.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Scalar fields in 3D asymptotically flat higher-spin gravity

Abstract In this work we construct a novel associative algebra and use it to define a theory of higher-spin gravity in (2 + 1)-dimensional asymptotically flat spacetimes. Our construction is based on a quotient of the universal enveloping algebra of i s l ( 2 , R ) with respect to the ideal generated by its Casimir elements, the mass squared M 2 and the three-dimensional analogue of the square of the Pauli–Lubanski vector S and propose to call the resulting associative algebra i h s ( M 2 , S ) . We provide a definition of its generators and even though we are not yet able to provide the complete set of multiplication rules of this algebra our analysis allows us to study many interesting and relevant sub-structures of i h s ( M 2 , S ) . We then show how to consistently couple a scalar field to an i h s ( M 2 , S ) higher-spin gauge theory.

Physics↗