Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Physical quantities”

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 919 records · Page 51

Chiral properties of the nucleon interpolating current and θ-dependent observables

We revisit the chiral properties of nucleon interpolating currents, and show that of the two leading order currents j 1 and j 2 , only two linear combinations j 1 ± j 2 transform covariantly under the anomalous U⁢(1) A symmetry. As a result, calculations of quantities which vanish by symmetry in the chiral limit may produce unphysical results if carried out with different linear combinations of the currents. This includes observables such as electric dipole moments, induced by the quantum chromodynamics (QCD) parameter θ, and the θ-dependence of the nucleon mass. For completeness, we also exhibit the leading order results for nucleon electric dipole moments (d n,p ) induced by θ, and the nucleon magnetic moments (μ n,p ), when calculated using QCD sum rules for both the covariant choices of the nucleon interpolating current. The results in each channel, conveniently expressed as the ratios, d n,p /μ n,p , are numerically consistent, and reflect the required physical dependence on θ.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

A system identification approach for non-intrusive reduced order modeling of radiation-induced photocurrents

In this study, development of compact photocurrent models is currently dominated by analytical techniques that rely on physical assumptions to render the governing equations solvable in a closed form. Violation of these assumptions can reduce the accuracy of the models and/or limit their scope. In this paper we show that system identification of nonlinear state-space systems can serve as an alternative numerical basis for non-intrusive reduced order modeling of photocurrent effects. To that end we develop a compact gray box photocurrent model (GBPM) by using a state-space representation with a low-dimensional latent state equation that mimics a mathematical model for the response of an idealized class of devices to ionizing radiation. In so doing we obtain a model that learns the dynamics of a quantity of interest directly from its measurements without requiring snapshots of the internal device state or its discretized model, and can be inferred from very small data sets. To demonstrate the approach we train the GBPM using a small experimental data set for a Z5236 Zener diode and a small synthetic data set obtained by simulating a synthetic pn-junction device. We then compare the GBPMs with black box models trained on the same data and show that performance of the latter is limited by the size of the data set, while the former are able to achieve excellent performance in both the reproductive and the predictive regimes.

97 MATHEMATICS AND COMPUTING↗

InteraChem: Virtual Reality Visualizer for Reactive Interactive Molecular Dynamics

Interactive molecular dynamics in virtual reality (IMD-VR) simulations provide a digital molecular playground for students as an alternative or complement to traditional molecular modeling kits or 2D illustrations. Previous IMD-VR studies have used molecular mechanics to enable simulations of macromolecules such as proteins and nanostructures for the classroom setting with considerable success. Here, we present the InteraChem molecular visualizer, intended for reactive IMD-VR simulation using semiempirical and ab initio methods. InteraChem visualizes not only the molecular geometry but also (1) isosurfaces such as molecular orbitals and electrostatic potentials and (2) two-dimensional graphs of time-varying simulation quantities such as kinetic/potential energy, internal coordinates, and user-applied force. Additionally, InteraChem employs speech recognition to facilitate user interaction and introduces a novel “atom happiness” visualization using emojis to indicate the energetic feasibility of a particular bonding arrangement. Furthermore, we include a set of accompanying exercises that we have used to teach chemical reactivity in small molecular systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Milestone in predicting core plasma turbulence: successful multi-channel validation of the gyrokinetic code GENE

On the basis of several recent breakthroughs in fusion research, many activities have been launched around the world to develop fusion power plants on the fastest possible time scale. In this context, high-fidelity simulations of the plasma behavior on large supercomputers provide one of the main pathways to accelerating progress by guiding crucial design decisions. When it comes to determining the energy confinement time of a magnetic confinement fusion device, which is a key quantity of interest, gyrokinetic turbulence simulations are considered the approach of choice – but the question, whether they are really able to reliably predict the plasma behavior is still open. The present study addresses this important issue by means of careful comparisons between state-of-the-art gyrokinetic turbulence simulations with the GENE code and experimental observations in the ASDEX Upgrade tokamak for an unprecedented number of simultaneous plasma observables.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A description of structured waves in shock compressed particulate composites

Dynamic compression of composite materials is of scientific interest because the mechanical mismatch between internal phases challenges continuum theories. Typical assumptions about steady wave propagation and quasi-instantaneous state changes require reexamination along with the need for time-dependent models. To that end, data and models are presented here for the shock compression of an idealized particulate composite. To serve as a generic representative of this material class, a polymer matrix was filled with tungsten particles, ranging from 1 to 50 vol. %. This creates a simple microstructure containing randomly scattered particles with an extreme impedance mismatch to the binding matrix. These materials were parallel plate impact loaded by Al flyers traveling at 1.8–5.0 km/s. Velocimetry provided records of the equilibrium state and the compression wave structure for each case with trends quantified by an empirical fit. The same quantities were also studied as a function of the wave's propagation distance. A homogenized viscoelastic model then made it possible to progress from cataloging trends to making predictions. Starting from a Mie–Gruneisen equation of state, additional time varying terms were added to capture the transient response. After calibration, accurate predictions of the steady wave structure were possible.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Neutrino oscillation phenomenology in the standard model and beyond

Over the last 20 years measurements of neutrino oscillation parameters have become very precise. In the standard neutrino oscillation picture most of the parameters are measured at the percent level. In this thesis we study neutrino oscillations in the standard picture and beyond. We analyze data from all types of neutrino oscillation experiments to obtain a global picture of neutrino oscillations. The remaining unknowns in the standard picture are the value of the CP-violating phase $\delta$, the octant of the atmospheric angle $\theta_{23}$ and the neutrino mass ordering. We discuss the current status of these unknowns and also comment on how well future experiments will do in measuring these quantities, where we discuss new facilities which will observe accelerator, atmossperic and reactor neutrinos. Using oscillation data it is possible to bound a possible violation of the CPT symmetry in the neutrino sector. This is done by performing a fit to neutrino and antineutrino oscillation data separately. We also address the capability of the future experiment DUNE to this scenario. We show that DUNE could improve some of the current bounds considerably. We discuss the sensitivity of DUNE together with the next generation reactor neutrino experiment JUNO to measure quasi-Dirac neutrino oscillations. We find that the future experiments will be able to bound the important quantities, mostly unbounded by current data. Finally, we discuss atmospheric neutrinos at the next generation neutrino telescope ORCA. We show that ORCA will put the best bounds on the invisible neutrino decay of $\nu_3$, coming from neutrino oscillation experiments. We also find that, a scenario like this would not affect the standard measurements of ORCA, namely the determination of the atmospheric neutrino oscillation parameters and the measurement of the neutrino mass ordering.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Deep Neural Network Informed Markov Chain Monte Carlo Methods

In subsurface flow modeling, quantifying the uncertainty of model parameters and the corresponding uncertainly on output quantities is a crucial task for groundwater management. Markov chain Monte Carlo (MCMC) methods can take advantage of observed data to estimate parameters in a Bayesian setting. However, MCMC can be slow to converge and produce highly correlated samples when the dimensions of the parameters is high. Using gradients for the posterior distribution can help samplers explore the parameter space more efficiently, but obtaining gradients can be computationally challenging.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Inverse energy transfer in decaying, three-dimensional, non-helical magnetic turbulence due to magnetic reconnection

ABSTRACT It has been recently shown numerically that there exists an inverse transfer of magnetic energy in decaying, non-helical, magnetically dominated, magnetohydrodynamic turbulence in three dimensions (3D). We suggest that magnetic reconnection is the underlying physical mechanism responsible for this inverse transfer. In the two-dimensional (2D) case, the inverse transfer is easily inferred to be due to smaller magnetic islands merging to form larger ones via reconnection. We find that the scaling behaviour is similar between the 2D and 3D cases, i.e. the magnetic energy evolves as t−1, and the magnetic power spectrum follows a slope of k−2. We show that on normalizing time by the magnetic reconnection time-scale, the evolution curves of the magnetic field in systems with different Lundquist numbers collapse on to one another. Furthermore, transfer function plots show signatures of magnetic reconnection driving the inverse transfer. We also discuss the conserved quantities in the system and show that the behaviour of these quantities is similar between the 2D and 3D simulations, thus making the case that the dynamics in 3D could be approximately explained by what we understand in 2D. Lastly, we also conduct simulations where the magnetic field is subdominant to the flow. Here, too, we find an inverse transfer of magnetic energy in 3D. In these simulations, the magnetic energy evolves as t−1.4 and, interestingly, a dynamo effect is observed.

79 ASTRONOMY AND ASTROPHYSICS↗

W$^+$W$^-$ boson pair production in proton-proton collisions at $\sqrt{s} =$ 13 TeV

A measurement of the W+W- boson pair production cross section in proton-proton collisions at s=13 TeV is presented. The data used in this study are collected with the CMS detector at the CERN LHC and correspond to an integrated luminosity of 35.9 fb-1. The W+W- candidate events are selected by requiring two oppositely charged leptons (electrons or muons). Two methods for reducing background contributions are employed. In the first one, a sequence of requirements on kinematic quantities is applied allowing a measurement of the total production cross section, 117.6±6.8 pb, which agrees well with the theoretical prediction. Fiducial cross sections are also reported for events with zero or one jet, and the change in the zero-jet fiducial cross section with the jet transverse momentum threshold is measured. Normalized differential cross sections are reported within the fiducial region. A second method for suppressing background contributions employs two random forest classifiers. The analysis based on this method includes a measurement of the total production cross section and also a measurement of the normalized jet multiplicity distribution in W+W- events. Finally, a dilepton invariant mass distribution is used to probe for physics beyond the standard model in the context of an effective field theory, and constraints on the presence of dimension-6 operators are derived.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Analytic continuation of noisy data using Adams Bashforth residual neural network

We propose a data-driven learning framework for the analytic continuation problem in numerical quantum many-body physics. Designing an accurate and efficient framework for the analytic continuation of imaginary time using computational data is a grand challenge that has hindered meaningful links with experimental data. The standard Maximum Entropy (MaxEnt)-based method is limited by the quality of the computational data and the availability of prior information. Also, the MaxEnt is not able to solve the inversion problem under high level of noise in the data. Here we introduce a novel learning model for the analytic continuation problem using a Adams-Bashforth residual neural network (AB-ResNet). Additionally, the advantage of this deep learning network is that it is model independent and, therefore, does not require prior information concerning the quantity of interest given by the spectral function. More importantly, the ResNet-based model achieves higher accuracy than MaxEnt for data with higher level of noise. Finally, numerical examples show that the developed AB-ResNet is able to recover the spectral function with accuracy comparable to MaxEnt where the noise level is relatively small.

97 MATHEMATICS AND COMPUTING↗

Toroidal magnetic field ripple in the presence of misaligned toroidal field coils on the COMPASS-U tokamak

Homogeneity of the toroidal magnetic field is a critical aspect of the tokamak fusion device. It is affected by the design of toroidal magnetic field coils and their manufacturing and assembly misalignments. The magnetic field of the COMPASS-U project is analysed using toroidal field ripple quantity. Here ANSYS modelling of bitter-plate-designed toroidal field coils is performed to prepare an accurate simplistic representation of the coils using a virtual current line. Toroidal field ripple is analysed using a Gaussian distribution of possible misalignments. The worst-case toroidal field ripple corresponding to the technical drawings only slightly exceeds the chosen limit of 0.5%.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Molecular Iodine Interactions with Fe, Ni, Cr, and Stainless Steel Alloys

The adsorption behavior of molecular iodine is important for understanding the spread of radioiodine in a nuclear accident. Prior experiments indicate that, in addition to the interaction with Fe, molecular iodine [i.e., I 2(g) ] also interacts with the next most abundant components of austenitic stainless steel (i.e., Ni, and Cr) at room temperature. Here, we investigate iodine adsorption on Fe, Ni, and Cr while focusing on understanding the variables affecting adsorption as well as the iodine compounds that are formed during adsorption. Scanning electron microscopy and energy-dispersive X-ray spectroscopy were used to characterize the surfaces of exposed metal particles and aid in the understanding of the morphology and chemistry of iodine interactions with the substrates. Inductively coupled plasma optical emission spectroscopy was used to detect low levels of metal iodides and X-ray photoelectron spectroscopy was used to confirm the formation of the metal iodides. The role of environmental factors (e.g., humidity and oxygen content) for iodine adsorption on metal substrates is addressed. The individual metals demonstrated formation of metal iodides for Fe and Ni particles from interaction with I 2(g) . The formation of metal iodides may indicate the affinity of iodine for the respective metal. In this study, the iodine affinities ranked Fe > Ni > Cr as determined by the quantity of chemisorbed iodine. This trend is also supported by the distributions and proportions of metals in the corrosion product of the stainless steels. The exposures without oxygen and humidity indicate the potential of a multistep iodine adsorption process where iodine first attacks the oxide layer and then chemisorbs to the exposed metal.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Identification of key oxidative intermediates and the function of chromium dopants in PKU-8: catalytic dehydrogenation of sec-alcohols with tert-butylhydroperoxide

Catalytic oxidation using green oxidants has played an important role in modern chemical engineering, however the in-situ generated active groups and the related catalytic mechanism need to be further understood. Combined with several measurements, such as isotopic tracers, molecular probes and cycle voltammetry, chromium ions were observed to undergo a Cr 3+ –Cr 2+ –Cr 3+ redox cycle. Simultaneously, t-BuOO* groups are generated in a large quantity and seem to function as the most critical active species to enhance the dehydrogenation of sec-alcohols when using TBHP as oxidant. DFT calculations give an energetically favorable pathway for TBHP activation, agreeing well with the obtained catalytic results. Here, a reaction scheme was finally proposed to interpret the probable mechanistic process.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Symmetry of Cosmological Observables, a Mirror World Dark Sector, and the Hubble Constant

We find that a uniform scaling of the gravitational free-fall rates and photon-electron scattering rate leaves most dimensionless cosmological observables nearly invariant. This result opens up a new approach to reconciling cosmic microwave background and large-scale structure observations with high values of the Hubble constant H0: Find a cosmological model in which the scaling transformation can be realized without violating any measurements of quantities not protected by the symmetry. A “mirror world” dark sector allows for effective scaling of the gravitational free-fall rates while respecting the measured mean photon density today. Further model building might bring consistency with the two constraints not yet satisfied: the inferred primordial abundances of deuterium and helium.

79 ASTRONOMY AND ASTROPHYSICS↗

Learning Forecasts of Rare Stratospheric Transitions from Short Simulations

Abstract Rare events arising in nonlinear atmospheric dynamics remain hard to predict and attribute. We address the problem of forecasting rare events in a prototypical example, sudden stratospheric warmings (SSWs). Approximately once every other winter, the boreal stratospheric polar vortex rapidly breaks down, shifting midlatitude surface weather patterns for months. We focus on two key quantities of interest: the probability of an SSW occurring, and the expected lead time if it does occur, as functions of initial condition. These optimal forecasts concretely measure the event’s progress. Direct numerical simulation can estimate them in principle but is prohibitively expensive in practice: each rare event requires a long integration to observe, and the cost of each integration grows with model complexity. We describe an alternative approach using integrations that are short compared to the time scale of the warming event. We compute the probability and lead time efficiently by solving equations involving the transition operator, which encodes all information about the dynamics. We relate these optimal forecasts to a small number of interpretable physical variables, suggesting optimal measurements for forecasting. We illustrate the methodology on a prototype SSW model developed by Holton and Mass and modified by stochastic forcing. While highly idealized, this model captures the essential nonlinear dynamics of SSWs and exhibits the key forecasting challenge: the dramatic separation in time scales between a single event and the return time between successive events. Our methodology is designed to fully exploit high-dimensional data from models and observations, and has the potential to identify detailed predictors of many complex rare events in meteorology.

Finkel, Justin↗

Helium Bubble Growth and Retention Characteristics in Aging Palladium Tritide

The helium-related characteristics of aged palladium tritide, including tritium loss, the formation of He bubbles from tritium decay-induced He atoms, and He thermal desorption behavior are computed using a continuum model of bubble evolution that incorporates the log-normal bubble spacing distribution deduced from 3 He nuclear magnetic resonance (NMR) measurements. This spacing distribution produces significant differences between mean (expectation) values of the bubble size and pressure and the average (median) values of these aging characteristics obtained by a previous calculation. The new calculations find much of the He is retained in a large quantity of smaller bubbles with higher bubble pressures, improving the overall stability of the bubbles and extending the age of He retention. By contrast, the integrated tritide material characteristics of swelling and tritium pressure are found relatively insensitive to details of the bubble spacing distribution. Inter-bubble fracture is predicted to begin with small, closely-spaced bubbles and progress to include larger bubbles near the onset of rapid He release. The critical concentration for this onset increases almost linearly with material tensile strength and decreases with shear modulus. Release of He from the aging solid is modeled as a growing network of fractured inter-bubble ligaments. For small particle or small grain material, the physical size of the fracture cluster, relative to the material dimensions, becomes important. In contrast with normal aging, inter-bubble fracture occurring during thermal desorption appears to begin with large, widely-spaced bubbles, indicating elevated temperature techniques may be of limited use for evaluating development of the bubble fracture network.

74 ATOMIC AND MOLECULAR PHYSICS↗

Blast Effects on Buildings (Final Report)

Lawrence Livermore National Laboratory (LLNL) has conservatively reduced the explosive safety standards suggested by the Small Quantities in Research Laboratories (SQRL) program testing by 60%. For this reason, further research into the detailed effects of small amounts of explosives in typically constructed rooms is needed in order to improve their factor of safety. The team was tasked with designing an experiment to investigate the effect of different variables on drywall under explosive blasts. In order to meet this objective, the team conducted a comprehensive literature review to gain an understanding of industry-standard construction practices and review previous tests conducted by the Army corp of engineers. Since the team used PBXN-5 rather than C4, as in SQRL, an initial shot was conducted to compare damage levels. From those results and the physical constraints of the testing chamber, the team redesigned multiple single-panel drywall frames to capture the entirety of the incurred damage. Proposed designs were narrowed down using a decision matrix. From the study of previous tests and literature review, variables were chosen that the team hypothesized to have an impact on drywall strength. The variables that were tested were chosen from the results of that work and specific variables the sponsor was interested in, and they were paint, humidity/moisture content, and explosive positioning relative to the studs. Detailed plans were made for each variable according to what conditions the team wanted to investigate. For humidity, this involved testing low, ambient, and high conditions by treating the panels in a chamber. Preliminary shots were performed to test the structural integrity of the frame and streamline the test diagnostics which involved a high-speed camera placed behind the drywall, outside of the chamber, and a pressure probe placed behind the drywall, inside the chamber. Once the instrumentation, diagnostics, and frame design were finalized, a quantitative damage criteria matrix was created to categorize the results of the main shot series. In conjunction with evidence from the high-speed video, the achieved damage levels indicate that high moisture content drywall is better able to withstand explosive blasts. Larger stud damage and lower drywall damage occurred when the explosive was located directly in front of a stud. Paint had no noticeable impact on strength. Ultimately, the team conducted a total of 17 tests, leaving the door open for future in-depth research into the impact of humidity.

36 MATERIALS SCIENCE↗

Forecasting Solar-Thermal Systems Performance under Transient Operation Using a Data-Driven Machine Learning Approach Based on the Deep Operator Network Architecture

Modeling and prediction of the dynamic behavior of thermal systems operating under intermittent energy input and variable load requirements represent one of the greatest challenges in the development of efficient and reliable renewable-based power generation technologies. In this work, a data-driven machine learning modeling framework was developed based on a modified version of the Deep Operator Network architecture where the time coordinate in the trunk net is replaced with historical data of the predicting quantity. The modeling framework can be used to accurately predict the performance of renewable-based energy conversion technologies including wind- and solar-based power plants. This novel framework was applied on a solar-thermal system that consists of a solar collection loop using a flat plate collector, a power generation loop comprising an Organic Rankine Cycle, and a thermal energy storage tank connecting both loops. Variable solar irradiance, air temperature, and power load profiles were used by the Deep Operator Network to predict the State-of-Charge and the efficiency of the thermal system for several days. The results were compared with the State-of-Charge and efficiency functions calculated using a physics-based model. For a simple operation scenario, characterized by a clear sky solar irradiance profile and constant load, the standard deviation in the State-of-Charge prediction by Deep Operator Network is below 0.9% during a seven-day prediction time horizon. For the most realistic operation scenario that considers real solar irradiance and a rough load profile, the maximum standard deviation in the predictions for the State-of-Charge and efficiency are below 6.8% and 2.5%, respectively. A comparison between Deep Operator Network and Long Short Term Memory network was also performed. In general, both networks predict very well the State-of-Charge for different data density conditions; however, a higher accuracy, with a standard deviation below 2.0%, is obtained by the Deep Operator Network during three and half days using sparser training data of 20-minute points. The same accuracy for the State-of-Charge prediction with the Long Short Term Memory network is achieved only for 14 h. Average standard deviations for the State-of-Charge prediction of 1.1% with the Deep Operator Network and 1.5% with the Long Short Term Memory network are obtained for a four-day prediction time using a denser training data of 5-minute points.

DeepONet↗