Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “perturbation 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 811 records · Page 45
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Strong QCD Insights from Excited Nucleon Structure Studies with CLAS and CLAS12

Studies of the spectrum of hadrons and their structure in experiments with electromagnetic probes offer unique insight into many facets of the strong interaction in the regime of large quark-gluon running coupling, i.e. the regime of strong QCD. The experimental program within Hall B at Jefferson Laboratory based on data acquired with the CLAS spectrometer using electron and photon beams with energies up to 6 GeV has already considerably extended the scope of research in hadron physics in joint efforts between experiment and phenomenological data analysis. Impressive progress in relating the hadron structure observables inferred from the data to the strong QCD mechanisms underlying hadron mass generation has been achieved in the past decade. Furthermore, these results will be considerably extended with data from the experimental program with the new CLAS12 spectrometer that has begun data taking using electron beams with energies up to 11 GeV. With this extended kinematic reach the structure of nucleon resonances will be probed at the highest photon virtualities ever achieved in the studies of exclusive electroproduction, which will allow for the exploration of the distance scale where>98% of light hadron mass emerges from QCD in the transition of the strong interaction from the regime of quark-gluon confinement to perturbative QCD.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Non-Gaussianity constraints with anisotropic μ distortion measurements from Planck

ABSTRACT Primordial non-Gaussianity can source μ-distortion anisotropies that are correlated with the large-scale temperature and polarization signals of the cosmic microwave background (CMB). A measurement of μT and μE correlations can therefore be used to constrain it on wavelengths of perturbations not directly probed by the standard CMB anisotropies. We carry out a first rigorous search for μ-distortion anisotropies with Planck data, applying the well-tested constrained ILC component-separation method combined with the needlet framework. We correlate the reconstructed μ map with the CMB anisotropies to derive constraints on the amplitude fNL of the local form bispectrum, specifically on the squeezed configurations with effective wavenumbers $k_s \simeq {740}\, \mathrm{Mpc^{-1}}$ and $k_L \simeq {0.05}\, \mathrm{Mpc^{-1}}$, improving previously estimated constraints by more than an order of magnitude. This enhancement is owing to the fact that we are able to use the full multipole information by carefully controlling biases and systematic effects in the analysis. We also for the first time incorporate constraints from measurements of μE correlations, which further tighten the limits. A combination of the derived Planck μT and μE power spectra yields |fNL| ≲ 6800 (95 per cent c.l.) on this highly squeezed bispectrum. This is only ≃3 times weaker than the anticipated constraint from Litebird. Furthermore we show that a combination of Litebird with Planck can improve the expected future constraint by ${\simeq}20{{\%}}$. These limits can be used to constrain multifield inflation models and primordial black hole formation scenarios, thus providing a promising novel avenue forward in CMB cosmology.

Rotti, Aditya (ORCID:0000000341245226)↗
DOE OSTI2023

Studying the hadron structure with PANDA and CLAS using machine learning techniques

The hadron spectroscopy and structure are currently very active fields of research to study the non-perturbative regime of quantum chronodynamics. The first one studies the complex structure of excited hadrons by looking at their decay products, while the latter uses lepton scattering on nucleons. Both methods require reconstruction algorithms with great efficiency and good particle identification and background rejection rates. This work aims to provide these by either improving the existing methods or developing new ones. The first part of this document presents a feasibility study of a predicted hybrid charmonium state for the PANDA experiment. Lattice QCD calculations predict the ground state hybrid charmonium to be a spin exotic with quantum numbers of JP C = 1?+ at a mass of around 4.3 GeV with a width to be around 20 MeV. A machine learning based data analysis scheme is proposed to further improve the signal efficiency and the background reduction, alongside with improvements of the analysis software (PandaRoot), that are vital for this study. These improvements include a reworked clustering algorithm for the electromagnetic calorimeter (EMC) and an optimized monte carlo matching for neutral particles. The second part of this document is about studying the proton structure. A multidimensional study of the structure function ratio Fsin(?)LU /FUU has been performed for K±, based on the measurement of beam-spin asymmetries. It uses the high statistics data recorded with the CLAS12 spectrometer at Jefferson Laboratory. Fsin(?)LU is a twist-3 quantity that provides information about the quark gluon correlations in the proton. This document will present for the first time a simultaneous analysis of two kaon channels over a large kinematic range of z, xB , PT and Q2 with virtualities Q2 ranging from 1 GeV2 up to 8 GeV2 using machine learning techniques for improved particle identification.

Kripko, Aron↗
DOE OSTI2025

The structure of quark mass corrections in the gg → HH amplitude at high-energy

The leading and next-to-leading order QCD predictions for Higgs boson pair production at hadron colliders suffer from a significant mass renormalisation scheme uncertainty related to the choice of the top-quark mass. The functional dependence of the result on the value of the intermediate quark mass can be understood in the high-energy limit using the Method of Regions and the tools of Soft-Collinear Effective Theory. In this work, we study the origin of the sizeable logarithmic mass corrections in the gg → HH amplitudes at leading and next-to-leading power in the limit s, |t|, |u| ≫ ${m}_t^2$ ≫ ${m}_H^2$. We argue that the mass corrections follow a predictable factorised pattern that can be exploited to simplify their computation. We present results for the leading power leading logarithmic corrections, our analysis leads to a significant reduction in the theoretical uncertainty of the double Higgs production amplitudes at scattering energies ≳ 1 TeV due to the top-quark mass scheme.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

On a variational formulation of the weakly nonlinear magnetic Rayleigh–Taylor instability

The magnetic-Rayleigh–Taylor (MRT) instability is a ubiquitous phenomenon that occurs in magnetically-driven Z-pinch implosions. It is important to understand this instability since it can decrease the performance of such implosions. In this work, I present a theoretical model for the weakly nonlinear MRT instability. I obtain such a model by asymptotically expanding an action principle, whose Lagrangian leads to the fully nonlinear MRT equations. After introducing a suitable choice of coordinates, I show that the theory can be cast as a Hamiltonian system, whose Hamiltonian is calculated up to the sixth order in a perturbation parameter. The resulting theory captures the harmonic generation of MRT modes. It is shown that the amplitude at which the linear magnetic-Rayleigh–Taylor instability exponential growth saturates depends on the stabilization effect of the magnetic-field tension. Overall, the theory provides an intuitive interpretation of the weakly nonlinear MRT instability and provides a systematic approach for studying this instability in more complex settings.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Molecular structure determination: Equilibrium structure of pyrimidine (m-C 4 H 4 N 2 ) from rotational spectroscopy (reSE) and high-level ab initio calculation (re) agree within the uncertainty of experimental measurement

The pure rotational spectrum of pyrimidine (m-C 4 H 4 N 2 ), the meta-substituted dinitrogen analog of benzene, has been studied in the millimeter-wave region from 235 GHz to 360 GHz. The rotational spectrum of the ground vibrational state has been assigned and fit to yield accurate rotational and distortion constants. Over 1700 distinct transitions were identified for the normal isotopologue in its ground vibrational state and least-squares fit to a partial sextic S-reduced Hamiltonian. Transitions for all four singly substituted 13 C and 15 N isotopologues were observed at natural abundance and were likewise fit. Deuterium-enriched samples of pyrimidine were synthesized, giving access to all eleven possible deuterium-substituted isotopologues, ten of which were previously unreported. Experimental values of rotational constants and computed values of vibration–rotation interaction constants and electron-mass corrections were used to determine semi-experimental equilibrium structures (r e SE ) of pyrimidine. The r e SE structure obtained using coupled-cluster with single, double, and perturbative triple excitations [CCSD(T)] corrections shows exceptional agreement with the re structure computed at the CCSD(T)/cc-pCV5Z level (≤0.0002 Å in bond distance and ≤0.03° in bond angle). Of the various computational methods examined, CCSD(T)/cc-pCV5Z is the only method for which the computed value of each geometric parameter lies within the statistical experimental uncertainty (2σ) of the corresponding semi-experimental coordinate. The exceptionally high accuracy and precision of the structure determination is a consequence of the large number of isotopologues measured, the precision and extent of the experimental frequency measurements, and the sophisticated theoretical treatment of the effects of vibration–rotation coupling and electron mass. Taken together, these demanding experimental and computational studies establish the capabilities of modern structural analysis for a prototypical monocyclic aromatic compound.

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

Growing neoclassical tearing modes seeded via transient-induced-multimode interactions

Nonlinear extended magnetohydrodynamics simulations demonstrating seeding of neoclassical tearing modes (NTMs) via magnetohydrodynamics-transient-induced multimode interactions are presented. Simulations of NTMs are enabled by two recent NIMROD code developments: the implementation of heuristic neoclassical stresses and the application of transient magnetic perturbations (MPs) at the boundary. NTMs are driven unstable by the inherently pressure driven kinetic bootstrap current, which arises due to collisional viscosity between passing and trapped electrons. These simulations use heuristic closures that model the neoclassical electron and ion stresses. NTM growth requires a seed island, which is generated by a transiently applied MP in simulations. The capability is demonstrated using kinetic-based reconstructions with flow of a DIII-D ITER Baseline Scenario discharge (La Haye et al., in Proceedings IAEA FEC, 2020). The applied MP seeds a 2/1 NTM that grows in two phases: a slow growth phase followed by a faster robust growth phase. Additionally, an evolving sequence of higher order core modes are excited at first. Power transfer analysis shows that nonlinear interactions between the core modes and the 2/1 helps drive the initial slow growth. Once the induced 2/1 magnetic island reaches a critical width, the NTM transitions to faster robust growth, which is well described by the nonlinear modified Rutherford equation. Furthermore, this work highlights the role of nonlinear mode coupling in seeding NTMs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Thermal modeling of three lakes within the continuous permafrost zone in Alaska using the LAKE 2.0 model

Abstract. Lakes in the Arctic are important reservoirs of heat with much lower albedo in summer and greater absorption of solar radiation than surrounding tundra vegetation. In the winter, lakes that do not freeze to their bed have a mean annual bed temperature >0 ∘C in an otherwise frozen landscape. Under climate warming scenarios, we expect Arctic lakes to accelerate thawing of underlying permafrost due to warming water temperatures in the summer and winter. Previous studies of Arctic lakes have focused on ice cover and thickness, the ice decay process, catchment hydrology, lake water balance, and eddy covariance measurements, but little work has been done in the Arctic to model lake heat balance. We applied the LAKE 2.0 model to simulate water temperatures in three Arctic lakes in northern Alaska over several years and tested the sensitivity of the model to several perturbations of input meteorological variables (precipitation, shortwave radiation, and air temperature) and several model parameters (water vertical resolution, sediment vertical resolution, depth of soil column, and temporal resolution). The LAKE 2.0 model is a one-dimensional model that explicitly solves vertical profiles of water state variables on a grid. We used a combination of meteorological data from local and remote weather stations, as well as data derived from remote sensing, to drive the model. We validated modeled water temperatures with data of observed lake water temperatures at several depths over several years for each lake. Our validation of the LAKE 2.0 model is a necessary step toward modeling changes in Arctic lake ice regimes, lake heat balance, and thermal interactions with permafrost. The sensitivity analysis shows us that lake water temperature is not highly sensitive to small changes in air temperature or precipitation, while changes in shortwave radiation and large changes in precipitation produced larger effects. Snow depth and lake ice strongly affect water temperatures during the frozen season, which dominates the annual thermal regime of Arctic lakes. These findings suggest that reductions in lake ice thickness and duration could lead to more heat storage by lakes and enhanced permafrost degradation.

Clark, Jason A. (ORCID:0000000248105261)↗
DOE OSTI2020

High speed PLIF study of the Richtmyer–Meshkov instability upon re-shock

The Richtmyer–Meshkov instability (RMI) of a twice-shocked gas interface is researched using high-speed planar-laser induced fluorescence (PLIF) in the Wisconsin Shock Tube Laboratory’s vertical shock tube. The initial condition (IC) is a shear layer with broadband diffuse perturbations at the interface between a helium-acetone mixture and argon. This IC is accelerated by a shock of nominal strength M = 1.8, and then accelerated again by the transmitted shock that reflects off the end wall of the tube. An estimate of the light gas mole fraction is extracted from high-speed imaging using an iterative process that accounts for the nonlinear temperature dependence of the acetone’s fluorescence quantum yield (FQY) and absorption cross-section. A vorticity deposition model for the initial growth rate after reshock is compared with the Mikaelian model for re-shock. Previously used in literature, the number of generations is shown to naturally arise from a normalisation of the scalar transport equation. A self-similar analysis is then further performed using the mole fraction data to explore the evolution of the RMI after reshock and the higher order moments of the light gas mole fraction are compared with a proposed model.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Investigation of resonant layer response in electron viscosity regime

We present a supplementary study of previous work in Waybright and Park [Phys. Plasmas 31, 022502 (2024)] which demonstrates a substantial effect of electron viscosity on the resonant layer response to non-axisymmetric magnetic perturbations. A main refinement is to include a curl element of electron viscosity in the generalized Ohm's law. The refinement reveals a resonant layer response in the electron viscosity regime corresponding to slowly rotating and highly viscous plasmas.

Asymptotic analysis↗
DOE OSTI2026

An overview of 3D field optimization for control of transport and edge instabilities on KSTAR

An international team from several laboratories and universities has made key advances over the last few years in the control of plasma transport and edge instabilities with applied 3D fields in the KSTAR tokamak to optimize long pulse operation scenarios. This overview begins with the optimization of both core and edge resonant magnetic perturbations (RMP) to improve fast ion confinement to avoid excessive limiter heat loads due to fast ion losses and successful modeling of the experimental results. Integrated and advanced plasma control techniques with machine learning (ML) and adaptive control were then used to optimize the 3D field spectrum in real-time to control edge localized modes (ELMs) while avoiding core locked modes that could disrupt the plasma. Accelerating the offline model of 3D fields with a surrogate ML model can optimize ELM suppression in the edge while limiting the impact of the applied RMP fields deeper in the plasma core in real-time. In addition, the impact of the 3D fields on the divertor heat load has been modeled and compared with experimental measurements. An analysis of a multi-machine database including KSTAR has been performed to better understand the metrics for the observed RMP thresholds for ELM suppression and the resulting plasma performance. Predictive modeling of the operational space for ELM suppression and density pumpout due to RMP has shown the importance of magnetic islands in the plasma edge and their impact on plasma turbulence. This research has culminated in the development of successful long pulse operational scenarios on KSTAR while attempting to overcome challenges of the new tungsten divertor.

3D fields↗
DOE OSTI2020

Microscopic theory of onset of decaging and bond-breaking activated dynamics in ultradense fluids with strong short-range attractions

We theoretically study thermally activated “in cage” elementary dynamical processes that precede full structural relaxation in ultra-dense particle liquids interacting via strong short range attractive forces. The analysis is based on a microscopic theory formulated at the particle trajectory level built on the dynamic free energy concept and an explicit treatment of how attractive forces control the formation and lifetime of physical bonds. Mean time scales for bond breaking, the early stage of cage escape, and non-Fickian displacement by a fixed amount are analyzed in the repulsive glass, bonded repulsive (attractive) glass, fluid, and dense gel regimes. The theory predicts a strong length-scale-dependent growth of these time scales with attractive force strength at fixed packing fraction, a much weaker slowing down with density at fixed attraction strength, and a strong decoupling of the shorter bond breaking time with the other two time scales that are controlled mainly by perturbed steric caging. All results are in good accord with simulations, and additional testable predictions are made. In conclusion, the classic statistical mechanical projection approximation of replacing all bare attractive and repulsive forces with a single effective force determined by pair structure incurs major errors for describing processes associated with thermally activated escape from transiently localized states.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Singling out modified gravity parameters and data sets reveals a dichotomy between Planck and lensing

ABSTRACT An important route to testing general relativity (GR) at cosmological scales is usually done by constraining modified gravity (MG) parameters added to the Einstein perturbed equations. Most studies have analysed so far constraints on pairs of MG parameters, but here, we explore constraints on one parameter at a time while fixing the other at its GR value. This allows us to analyse various models while benefiting from a stronger constraining power from the data. We also explore which specific data sets are in tension with GR. We find that models with (μ = 1, η) and (μ, η = 1) exhibit a 3.9σ and 3.8σ departure from GR when using Planck18 + Supernovae type Ia (SNe) + Baryon acoustic oscillation (BAO), while (μ, η) shows a tension of 3.4σ. We find no tension with GR for models with the MG parameter Σ fixed to its GR value. Using a Bayesian model selection analysis, we find that some one-parameter MG models are moderately favoured over ΛCDM when using all data set combinations except Planck cosmic microwave background lensing and dark energy survey data. Namely, Planck18 shows a moderate tension with GR that only increases when adding any combination of redshift space distortion, SNe, or BAO. However, adding lensing diminishes or removes these tensions, which can be attributed to the ability of lensing in constraining the MG parameter Σ. The two overall groups of data sets are found to have a dichotomy when performing consistency tests with GR, which may be due to systematic effects, lack of constraining power, or modelling. These findings warrant further investigation using more precise data from ongoing and future surveys.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Characterization and interaction of geometric and contact/impact nonlinearities in dynamical systems

Here, we study how a contact/impact nonlinearity interacts with a geometric cubic nonlinearity in an oscillator system. Specific focus is shown to the effects on bifurcation behavior and secondary resonances (i.e., super- and sub-harmonic resonances). The effects of the individual nonlinearities are first explored for comparison, and then the influences of the combined nonlinearities, varying one parameter at a time, are analyzed and discussed. Nonlinear characterization is then performed on an arbitrary system configuration to study super- and sub-harmonic resonances and grazing contacts or bifurcations. Both the cubic and contact nonlinearities cause a drop in amplitude and shift up in frequency for the primary resonance, and they activate high-amplitude subharmonic resonance regions. The nonlinearities seem to never destructively interfere. The contact nonlinearity generally affects the system’s superharmonic resonance behavior more, particularly with regard to the occurrence of grazing contacts and the activation of many bifurcations in the system’s response. The subharmonic resonance behavior is more strongly affected by the cubic nonlinearity and is prone to multistable behavior. Perturbation theory proved useful for determining when the cubic nonlinearity would be dominant compared to the contact nonlinearity. The limiting behaviors of the contact stiffness and freeplay gap size indicate the cubic nonlinearity is dominant overall. It is demonstrated that the presence of contact may result in the activation of several bifurcations. In addition, it is proved that the system’s subharmonic resonance region is prone to multistable dynamical responses having distinct magnitudes.

42 ENGINEERING↗
DOE OSTI2025

Microbial Community Analysis & Functional Evaluation in Soils

The overall objective of this proposal was to develop technologies to alter the composition and function of important members of microbial communities. In particular, the overall objective of the microbial community editing portion of the proposal focuses on developing foundational tools and understanding required to predict, alter and design grass rhizosphere communities impacting DOE missions. Specifically, the project is centered on the Microbial Community Analysis & Functional Evaluation in Soils (m-CAFES) to manipulate microbial consortia associated with plants of interest for the bioenergy sector, under the presumption that bacterial communities can be manipulated to enhance plant health. For tasks of specific interest to us, we are focusing on developing novel Clustered Regularly Interspaced Short Palindromic Repeats (CRISPR) based technologies (primarily focusing on Aim 1) and their delivery modalities (notably subaim 1.2) to edit specific bacterial genomes of interest to enhance their functionalities, and programmably ablate specific undesirable members of bacterial communities for plant health. We are focusing on engineering bacteriophages (bacterial viruses, for subaim 1.2) to carry programmable CRISPR-Cas systems (subaim 1.1) to target (ablate) or alter (edit) genomes of interest. This will enable us to carry out microbial perturbations that will impact community composition and function and ultimately plant growth and health, to enable the next phase of the project by deploying them in situ (subaims 1.3 and 1.4).

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Propagating information content: an example with advection

The mathematical algorithm to derive geophysical information from remote sensing observations is called a retrieval. The mathematics of many retrieval problems are ill-posed, and thus a priori information is used to help constrain the derived geophysical variable to realistic values. One quantity of interest, therefore, is the information content of the observation. Perfect information content in the observation would be achieved if the retrieval were able to capture any perturbation in the desired geophysical variable with the proper magnitude. Many new data products can be derived by combining geophysical variables retrieved from multiple different remote sensors. This paper explores, for the first time, how to derive the information content of these derived products. The approach uses traditional error propagation techniques to derive the uncertainty of the derived field twice, both when the observations are used in the retrieval and also when only the a priori information from each remote sensor is propagated. These two uncertainties are then used to provide an estimate of the information content of the derived geophysical variable. This study demonstrates how to propagate the uncertainties from six different instruments to provide the information content for water vapor and temperature advection. A multi-month analysis demonstrates that, in a mean sense, the information content for temperature advection is nearly unity for all heights below 700 m while, the information content for water vapor advection is somewhat more variable but still larger than 0.6 in the convective boundary layer.

Turner, David D. [National Oceanic and Atmospheric↗
DOE OSTI2021

Flow boiling transient critical heat flux tests with stainless steel and FeCrAl: Transient correlation implementation, model calibration, and sensitivity analysis

In this study, steady-state and transient internal-flow critical heat flux (CHF) experiments were carried out under two flow conditions, at atmospheric pressure, on 316 stainless steel (316-SS) and iron-chromium-aluminum (FeCrAl) tubes. Slow power transients at a low mass flow (300 kg/m 2 - s) with very low subcooling (X e =-0.0054) generated premature CHFs, which were prevented by faster power transients and high mass flow (1000 kg/m 2 - s). The measured transient CHFs increased linearly with increasing power transients, compared with the steady-state CHF. Yet, the wall superheat at the CHFs decreased with faster power transients. Transient CHF correlations highlighting heterogeneous spontaneous nucleation were calibrated to the measured CHFs and compared with other existing correlations. Transient CHF multipliers were acquired from pool and flow boiling empirical CHF correlations that were generated. The multipliers were applied to the RELAP5-3D nuclear system code to analyze the discrepancies between the measured data and the predicted CHF and post-CHF behavior, which improved peak cladding temperature predictions by 24.6%. A variance-based global sensitivity study perturbing the experimental uncertainties and cladding material thermal properties showed the diminishing influence of flow boiling heat transfer with increasing power transients, highlighting the significance of the volumetric heat capacity for the cladding integrity during transients. Transient CHF correlations were applied to the most limiting design basis accident: a hot-zero-power reactivity-initiated accident on a generic pressurized water reactor RELAP5-3D model.

36 MATERIALS SCIENCE↗