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 847 records · Page 47
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Runaway electron seed formation at reactor-relevant temperature

Systematic variation of the pre-disruption core electron temperature (T e ) from 1 to 12 keV using an internal transport barrier scenario reveals a dramatic increase in the production of ‘seed’ runaway electrons (REs), ultimately accessing near-complete conversion of the pre-disruption current into sub-MeV RE current. Injected Ar pellets are observed to ablate more intensely and promptly as T e rises. At high T e , the observed ablation exceeds predictions from published thermal ablation models. Simultaneously, the thermal quench (TQ) is observed to significantly shorten with increasing T e —a surprising result. While the reason for the shorter TQ is not yet understood, candidate mechanisms include: insufficiently accurate thermal ablation models, enhanced ablation driven by the seed RE population, or significant parallel heat transport along stochastic fields. Kinetic modeling that self-consistently treats the plasma cooling via radiation, the induced electric field, and the formation of the seed RE is performed. Including the combined effect of the inherent dependence of hot-tail RE seeding on T e together with the shortened TQ, modeling recovers the progression towards near-complete conversion of the pre-disruption current to RE current as T e rises. Measurement of the HXR spectrum during the early current quench (CQ) reveals a trend of decreasing energy with pre-disruption T e . At the very highest T e (≈ 12 keV), ≈ 100% conversion of the thermal current to runaway current is found. The energy of this peculiar RE beam is inferred to be sub-MeV as it emits vanishingly few MeV hard x-rays (HXRs). These measurements demonstrate novel TQ dynamics as Te is varied and illustrate the limitations of treating the RE seed formation problem without considering the inter-related dependencies of the pellet ablation, radiative energy loss, and resultant variations of the TQ duration. Finally, if the observed shortening of the TQ with increasing T e extends to fusion-grade plasmas, than their propensity to form large quantities of RE seeds at high T e may be far worse than previously thought. Positively, the high T e scenario in DIII-D produces REs so prodigiously that it can serve as a meaningful new platform for demonstrating RE avoidance techniques.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Feasibility of measuring the speed of sound of the quark-gluon plasma from the multiplicity and mean 𝑝 𝑇 of ultracentral heavy-ion collisions

The mean transverse momentum ⟨𝑝 𝑇 ⟩ of hadrons has been observed experimentally and in numerical simulations to have a power-law dependence on the hadronic multiplicity 𝑁 in ultracentral relativistic heavy-ion collisions: ⟨𝑝 𝑇 ⟩∝𝑁 𝑏 UC . It has been put forward that this exponent 𝑏 UC is the speed of sound of quark-gluon plasma measured at a temperature determined from ⟨𝑝 𝑇 ⟩. We study step by step the connection between (i) the energy and entropy of hydrodynamic simulations and (ii) experimentally measurable observables. We show that an argument based on energy and entropy should yield an exponent equal to the pressure over energy density 𝑃/ɛ, rather than the speed of sound 𝑐$_s^2$; however, we also observe that ⟨𝑝 𝑇 ⟩ and 𝑁 are not sufficiently accurate proxies for the energy and entropy to make this possible in practice. From simulations, we find that the exponent 𝑏 UC is significantly different whether the “effective volume” is strictly constant or not, a condition that cannot be enforced experimentally. Additional tests using a modified equation of state find that the exponent 𝑏 UC exhibits a variable degree of correlations with the speed of sound and with 𝑃/ɛ, but is not an accurate measurement of either quantity in general.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Non-Gaussian distribution of displacements for Lennard-Jones particles in equilibrium

Most mesoscale simulation methods assume Gaussian distributions of velocity-like quantities. These quantities are not true velocities, however, but rather time-averaged velocities or displacements of particles. Here, we show that there is a large range of coarse-graining scales where the assumption of a Gaussian distribution of these displacements fails, and a more complex distribution is required to adequately express these distribution functions of displacements.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Non-isometry, state dependence and holography

We establish an equivalence between non-isometry of quantum codes and state dependence of operator reconstruction, and discuss implications of this equivalence for holographic duality. Specifically, we define quantitative measures of non-isometry and state dependence and describe bounds relating these quantities. In the context of holography we show that, assuming known gravitational path integral results for overlaps between semiclassical states, non-isometric bulk-to-boundary maps with a trivial kernel are approximately isometric and bulk reconstruction approximately state-independent. In contrast, non-isometric maps with a non-empty kernel always lead to state-dependent reconstruction. We also show that if a global bulk-to-boundary map is non-isometric, then there exists a region in the bulk which is causally disconnected from the boundary. Finally, we conjecture that, under certain physical assumptions for the definition of the Hilbert space of effective field theory in AdS space, the presence of a global horizon implies a non-isometric global bulk-to-boundary map.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

An optimization method for chaotic turbulent flow

Evidence indicates that quantities-of-interest in some turbulent flows can be controlled despite the overall chaotic dynamics. It is typically thought that this is via relatively deterministic, larger-scale components of the turbulence. However, finding such controls, if they exist, is challenging because chaos causes sensitivity gradients to explode and the search space to become intractably non-convex. This challenge is analyzed, and a penalty method is introduced to cope with it. In the new approach, the time domain is broken into segments approximately matching the chaos time scales, so that the solution within each segment is both physical and relatively deterministic. The initial condition of each segment is included in an adjoint-based gradient optimization, which temporarily introduces artificial Δq discontinuities in the overall solution. The optimization then proceeds in stages with increasing penalization of Δq. The method is developed and illustrated for a logistic map, the Lorenz Equation, and an advection augmented Kuramoto–Sivashinsky Equation. These examples show how the Δq temporarily increases the search scale prior to the strong Δq → 0 penalization that recovers a physical solution. It is then applied to turbulent Kolmogorov flow, for which it also far outperforms a standard adjoint-based gradient search. Finally, the utility of such an optimized chaotic solution is discussed.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Development of a deep potential model for F and CF 2 etching of Si and SiO 2

An understanding of plasma-surface interactions at increasingly smaller scales is invaluable for the development of novel technologies and processing techniques. Molecular dynamics (MD) simulations can provide insights into atomic-scale interactions, though they are restricted by the availability of interatomic potentials. Machine learning methods, such as Deep Potential Molecular Dynamics (DeepMD), provide a systematic framework for the development of accurate and flexible ab initio-based models. In this work, we develop DeepMD models for the ion-enhanced etching of Si and SiO 2 by F and CF 2 radicals. We employ an active learning process to expand the data set on which the model is trained and demonstrate its effect on the model accuracy. The DeepMD results are compared to data from classical MD simulations and experiments. Physical sputtering yields of SiO 2 by Ar + ions show good agreement with previous simulation results using conventional interatomic potentials, though the predicted depth profiles are different. Etching yields are calculated as a function of ion energy and neutral to ion flux ratio for the Ar + ion-enhanced etching of SiO 2 and Si by F atoms, as well as for etching of SiO 2 by CF 2 radicals, showing reasonable agreement with experimental data. Finally, an ion-enhanced surface kinetic model is fitted to the DeepMD etch yields, and the fitted parameters are compared to quantities computed directly from DeepMD simulations. This study illustrates how molecular dynamics simulations using machine learning potentials can provide an accurate model of etching processes relevant to device manufacturing.

Kounis-Melas, Andreas [Princeton Univ., NJ (United↗
DOE OSTI2020

Variability in the Massive Open Cluster NGC 1817 from K2: A Rich Population of Asteroseismic Red Clump, Eclipsing Binary, and Main-sequence Pulsating Stars

We present a survey of variable stars detected in K2 Campaign 13 within the massive intermediate-age (∼1 Gyr) open cluster NGC 1817. We identify a complete sample of 44 red clump stars in the cluster, and have measured asteroseismic quantities (ν {sub max} and/or Δν) for 29 of them. Five stars showed suppressed dipole modes, and the occurrence rates indicate that mode suppression is unaffected by evolution through core helium burning. A subset of the giants in NGC 1817 (and in the similarly aged cluster NGC 6811) have ν {sub max} and Δν values at or near the maximum observed for core helium–burning stars, indicating they have core masses near the minimum for fully nondegenerate helium ignition. Further asteroseismic study of these stars can constrain the minimum helium core mass in red clump stars and the physics that determines this limit. Two giant stars show photometric variations on timescales similar to previously measured spectroscopic orbits. Thirteen systems in the field show eclipses, but only five are probable cluster members. We identify 32 δ Sct pulsators, 27 γ Dor candidates, and 7 hybrids that are probable cluster members, with most being new detections. We used the ensemble properties of the δ Sct stars to identify stars with possible radial pulsation modes. Among the oddities we have uncovered are: an eccentric orbit for a short-period binary containing a δ Sct pulsating star; a rare subgiant within the Hertzsprung gap showing δ Sct pulsations; and two hot γ Dor pulsating star candidates.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Topical Report on Water Balance Model Development and Initial Application

This topical report summarizes results from developing a generic dynamic water balance tool which can be used by coal-fired power plant operators to evaluate FGD water management strategies. EPRI is using a statistical software model created on the GoldSim platform and data from coal-fired energy plants to determine the fraction of power plant water needs that can be satisfied by electrodialysis reversal treatment and, correspondingly, the discharge volume and/or byproducts that still need to be managed. The generic model is designed to accept a number of easily obtainable process variables and generate probabilities of resultant water characteristics. The model is also constructed as a universal process flow; one that can be adjusted to accommodate most major components of an FGD blowdown handling system. The generic model is dynamic and utilizes user-inputted operational variation in water quality and quantity to best provide a probabilistic representation of potential performance outcomes. In order to account for complexities of water composition and treatment, a separate chemical modeling software is used to evaluate water chemistry and unit operation performance. The thermodynamic modeling software, OLI Studio, is employed to convert ion data from laboratory analyses into balanced molecular concentrations that are classified as dissolved and suspended solids. The dynamic water model is designed to import a total of 11 species commonly found in FGD blowdown water as molecular concentrations.Water management unit operations are used in the model and categorized by either physical unit operation blocks or chemical unit operations blocks. As the name implies, the physical water treatment blocks are unit operations that are exclusively physical changes to the system (e.g., separations and volume changes). The chemical water treatment blocks are unit operations where chemical changes occur and require the use of OLI Studio to characterize subsequent streams. Water and material balances from three unique coal-fired power plants were obtained and preliminary modeling on the dynamic water balance tool was performed using the extracted data from these power plants. These three power plants were chosen due to their unique configurations and water profiles that represented a variety of operational changes. Results from using the dynamic tool with the three plant configurations are presented and discussed in the context of the existing tool and planned future capabilities.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2023

The hazard components of representative key risks. The physical climate perspective

The framework of Representative Key Risks (RKRs) has been adopted by the Intergovernmental Panel on Climate Change Working Group II in order to categorize, assess and communicate a wide range of regional and sectoral key risks, expected to become severe as a consequence of the potentially detrimental convergence of changing climate conditions, exposure and vulnerability of human and natural systems, and responses through (mal)adaptation measures. Other papers in this special issue treat each of eight RKRs by holistically assessing their current status and future evolution as a result of this convergence. However, their treatment cannot always organize such assessment according to a systematic gradation of climatic changes, often times having to extrapolate the big-picture evolution of risk from either qualitative effects of “low” medium” and “high” warming, or limited/focused analysis of the consequences of particular mitigation choices (e.g., benefits of limiting warming to 1.5 or 2C). In this study we offer as systematic a representation as possible – given current literature and assessments – of the future evolution of the hazard components of RKRs. We draw on the assessment of changes in climatic impact-drivers relevant to RKRs described in the 6th Assessment Report by Working Group 1 (WGI) supplemented when needed by more recent literature. We identify the relevant hazards for each RKR, based upon the WGII authors’ assessment, and we report on their current state and expected future changes in magnitude, intensity and/or frequency, linking these changes to Global Warming Levels to the extent possible. For some of these quantities - like regional trends in oceanic and atmospheric temperature and precipitation, some heat and precipitation extremes, permafrost thaw and Northern Hemisphere snow cover - a strong and quantitative relationship with global average temperature change has been assessed. For others - like frequency and intensity of tropical cyclones and extra-tropical storms, and fire weather - the link to increasing global temperatures can only be described qualitatively. For some processes - like the behavior of ice-sheets, or changes in circulation dynamics - large uncertainties about the effects of different levels of warming remain, and for a few others - like ocean pH and air pollution - the composition of the scenario of anthropogenic emissions is most relevant, rather than the warming reached. In almost all cases, however, the basic message remains that every small increment in CO2 concentration in the atmosphere and associated warming will bring changes in climate phenomena that will contribute to increasing risk of impacts on human and natural systems when these changes happen in conjunction with these systems’ conditions of exposure and vulnerability, and in lack of effective adaptation. Our picture of the evolution of RKR-relevant climatic impact-drivers complements and enriches the treatment of RKRs in the other papers in at least two ways: by filling in an often only cursory or limited representation of the physical climate aspects driving impacts, and by providing a fuller representation of their future possible evolution which represents an important component – if never the only one – of the potential future evolution of risk severity.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Uncertainty Quantification Framework for Predicting Material Response with Large Number of Parameters: Application to Creep Prediction in Ferritic-Martensitic Steels Using Combined Crystal Plasticity and Grain Boundary Models

This paper presents an uncertainty quantification (UQ) framework for the physics-based model prediction of material response with a large number of parameters. The application problem presented in this work is that of predicting creep in Grade 91 steel at 600°C. The material response is defined with a physically based microstructural model with constitutive equations emulating several observed phenomena in Grade 91 and embodied into an explicit geometry mesoscale finite element model for prior austenite grains and grain boundaries. Creep within the grains and in grain boundaries are represented by crystal plasticity for dislocation motion and a physics-based model for cavity growth and nucleation, respectively. The creep behavior of this material is influenced by several parameters, some of which have a wide range of variation based on experimental data. UQ combined with microstructural modeling can discover the core microstructural causes of experimental variability, leading to improved materials with lower variability in critical long-term material properties. In this study, we investigate the model's uncertainty to identify material properties that may be modified during production to increase creep life and analyze different components of the crystal plasticity model for improvements. For this purpose, a quantity of interest is defined as time to minimum creep rate, which correlates well to the creep failure of the material. A deep neural network model was trained and validated to be used as a surrogate for the finite element model. Then, a variance-based sensitivity analysis is performed on the surrogate model to find the Sobol indices of the input parameters in respect to the output quantity of interest. The Sobol indices are used to reduce the dimensionality of the model. Generalized polynomial chaos expansion is used on the reduced basis models to propagate the uncertainty from the input parameters to the quantity of interest using the deep neural network surrogate model. These results are benchmarked against uncertainty propagation using Monte Carlo simulations. In conclusion, the UQ performed through the reduced basis model captures almost all the uncertainty in the model with significantly fewer simulations, making it possible to perform the UQ directly via simulations with the finite element model rather than surrogate machine-learned models.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Energy dependent calculations of fission product, prompt, and delayed neutron yields for neutron induced fission on 235 U, 238 U, and 239 Pu

In this article, we perform energy-dependent calculations of independent and cumulative fission product yields for 235 U, 238 U, and 239 Pu in the first chance fission region. Starting with the primary fission fragment distributions taken from available experimental data and analytical functions based on assumptions for the excitation energy and spin-parity distributions, the Hauser-Feshbach statistical decay treatment for fission fragment de-excitation is applied to more than 1,000 fission fragments for the incident neutron energies up to 5 MeV. The calculated independent yields are then used as an input of β-decay calculations to produce the cumulative yield, and summation calculations are performed. Model parameters in these procedures are adjusted by applying the Bayesian technique at the thermal energy for 235 U and 239 Pu and in the fast energy range for 238 U. The calculated fission observable quantities, such as the energy-dependent cumulative yields, and prompt and delayed neutron yields, are compared with available experimental data. We also study the possible impact of the second chance fission opening on the energy dependence of the delayed neutron yield by on the opening of second chance fission on the energy dependence calculation.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Fluctuations and Conservation Laws

We report on recent progress concerning effects of global conservation laws on cumulants of conserved quantities. Specifically, we will relate for an arbitrary equation of state cumulants of a conserved charge measured in a subvolume of a thermal system with the corresponding grandcanonical susceptibilities, taking into account exact global conservation of that charge. Applications to actual measurement at the RHIC and LHC as well as extensions to multiple conserved charges will be discussed.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Linear embedding of nonlinear dynamical systems and prospects for efficient quantum algorithms

The simulation of large nonlinear dynamical systems, including systems generated by discretization of hyperbolic partial differential equations, can be computationally demanding. Such systems are important in both fluid and kinetic computational plasma physics. This motivates exploring whether a future error-corrected quantum computer could perform these simulations more efficiently than any classical computer. In this work, we describe a method for mapping any finite nonlinear dynamical system to an infinite linear dynamical system (embedding) and detail three specific cases of this method that correspond to previously studied mappings. Then we explore an approach for approximating the resulting infinite linear system with finite linear systems (truncation). Using a number of qubits only logarithmic in the number of variables of the nonlinear system, a quantum computer could simulate truncated systems to approximate output quantities if the nonlinearity is sufficiently weak. Other aspects of the computational efficiency of the three detailed embedding strategies are also discussed.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Multi-fidelity thermal modeling of laser powder bed additive manufacturing

Laser powder bed fusion (LPBF) Additive manufacturing (AM) has attracted interest as an agile method of building production metal parts to reduce design-build-test cycle times for systems. However, predicting part performance is difficult due to inherent process variabilities. This makes qualification challenging. Computational process models have attempted to address some of these challenges, including mesoscale, full physics models and reduced fidelity conduction models. The goal of this work is credible multi-fidelity modeling of the LPBF process by investigating methods for estimating the error between models of two different fidelities. Two methods of error estimation are investigated, adjoint-based error estimation and Bayesian calibration. Adjoint-based error estimation is found to effectively bounding the error between the two models, but with very conservative bounds, making predictions highly uncertain. Bayesian parameter calibration applied to conduction model heat source parameters is found to effectively bound the observed error between the models for melt pool morphology quantities of interest. However, the calibrations do not effectively bound the error in heat distribution.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Thermodynamics–Structure–Dynamics Correlations and Nonuniversal Effects in the Elastically Collective Activated Hopping Theory of Glass-Forming Liquids

We employ the microscopic Elastically Collective Nonlinear Langevin Equation (ECNLE) theory of activated dynamics in combination with crystal-avoiding simulations to study four inter-related questions for metastable monodisperse hard sphere fluids. The first is how significantly improved integral equation theory structural input (Modified-Verlet (MV) closure) changes the dynamical predictions of ECNLE theory. The main consequence is a modest enhancement of the importance of the collective elastic barrier relative to its local cage contribution, which increases the alpha relaxation time and fragility relative to prior results based on the Percus–Yevick closure. Second, ECNLE-MV theory predictions for the alpha time and self-diffusion constant in the metastable regime are quantitatively compared to our new simulations. The small adjustment of a numerical prefactor that enters the collective elastic barrier leads to quantitative agreement over three decades. Third, using the more accurate MV structural input, ECNLE theory is shown to predict thermodynamics–structure–dynamics “correlations” based on various long and short wavelength scalar properties all related to static two-point collective density fluctuations. The logarithm of the alpha relaxation time scales as a power law with these scalar metrics with an exponent that is significantly lower in the less dense noncooperative activated regime compared to the very dense highly cooperative regime. However, the discovered correlation of activated relaxation with a thermodynamic property (dimensionless compressibility) is not causal in ECNLE theory, but rather reflects a strong connection between the local structural quantities that quantify kinetic constraints in the theory with the amplitude of long wavelength density fluctuations. Fourth, the consequences of chemically specific nonuniversalities associated with the onset condition and relative importance of collective elasticity are studied. In conclusion, the predicted thermodynamics–structure–dynamics correlations are found to be robust, albeit with nontrivial shifts of the onset condition.

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

Validation of a Global Geospace Model With a Systems Science Approach Based on Canonical Correlation Analysis

A systems science approach based on canonical correlation analysis (CCA) is applied as a new, behavioral way to validate global geospace models. The biggest novelty of the technique is that it validates models at a system level, whereby a side‐by‐side comparison is performed of CCA applied to a 30‐day observational and the corresponding simulation data sets comprising quiet, moderate and active times. The simulation used the Multiscale Atmosphere‐Geospace Environment (MAGE) model. It is shown that (a) CCA must be combined with sensitivity analysis to be effective, (b) the MAGE model generally reproduces the observed behavior (more so for quieter time intervals), quantified by the intercorrelations between different variables and (c) the technique identifies the SuperMAG SML index as a quantity for which refinements of the model are needed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Equation-of-motion coupled-cluster method with double electron-attaching operators: Theory, implementation, and benchmarks

We report a production-level implementation of the equation-of-motion (EOM) coupled-cluster (CC) method with double electron-attaching (DEA) EOM operators of 2 p and 3 p 1 h types, EOM-DEA-CCSD. This ansatz, suitable for treating electronic structure patterns that can be described as two-electrons-in-many orbitals, represents a useful addition to the EOM-CC family of methods. We analyze the performance of EOM-DEA-CCSD for energy differences and molecular properties. By considering reduced quantities, such as state and transition one-particle density matrices, we compare EOM-DEA-CCSD wave functions with wave functions computed by other EOM-CCSD methods. The benchmarks illustrate that EOM-DEA-CCSD is capable of treating diradicals, bond-breaking, and some types of conical intersections.

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

Diffusive hydrodynamics from integrability breaking

We describe here the crossover from generalized to conventional hydrodynamics in nearly integrable systems. Integrable systems have infinitely many conserved quantities, which spread ballistically, in general. When integrability is broken, only a few of these conserved quantities survive. The remaining conserved quantities are generically transported diffusively; we further derive a compact and general diffusion equation for these. The diffusion constant depends on the matrix elements of the integrability-breaking perturbation; for a certain class of integrability-breaking perturbations, including long-range interactions, the diffusion constant can be expressed entirely in terms of generalized hydrodynamic data.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗