Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “experimental evolution”

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 523 records · Page 29

Checking nonflow assumptions and results via PHENIX published correlations in p+p, p+Au, d+Au, and 3 He+Au at √ s NN=200 GeV

Recently the PHENIX Collaboration has made available two-particle correlation Fourier coefficients for multiple detector combinations in minimum bias p+p and 0–5% central p+Au, d+Au, and 3 He+Au collisions at √ s NN = 200 GeV [Phys. Rev. C 105, 024901 (2022)]. Using these coefficients for three sets of two-particle correlations, azimuthal anisotropy coefficients v 2 and v 3 are extracted for midrapidity charged hadrons as a function of transverse momentum. In this paper, we use the available coefficients to explore various nonflow hypotheses as well as to compare the results with theoretical model calculations. The nonflow methods fail basic closure tests with ampt and pythia/angantyr, particularly when including correlations with particles in the low multiplicity light-projectile going direction. In data, the nonflow adjusted v 2 results are modestly lower in p+Au and the adjusted v 3 results are more significantly higher in p+Au and d+Au. However, the resulting higher values for the ratio v 3 /v 2 in p+Au at RHIC compared to p+Pb at the LHC is additional evidence for a significant overcorrection. Incorporating these additional checks, the conclusion that these flow coefficients are dominated by initial geometry coupled with final-state interactions (e.g., hydrodynamic expansion of quark-gluon plasma) remains true, and explanations based solely on initial-state glasma are ruled out. The detailed balance between intrinsic and fluctuation-driven geometry and the exact role of weakly versus strongly coupled prehydrodynamic evolution remains an open question for triangular flow, requiring further theoretical and experimental investigation.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

From Reproducible Edge–Cloud Experimentation to Real-World Practice: The E2Clab Experience

Reproducibility is already difficult in distributed systems; on the computing continuum, it becomes substantially harder. Applications that span sensing devices, edge and fog resources, and cloud platforms must be evaluated across heterogeneous hardware, variable network conditions, cross-layer orchestration decisions, and long-running workflow lifecycles. We use E2Clab as a case study to examine these challenges and their implications for experimental methodology. We explain why reproducible experimentation is harder on the continuum, then revisit E2Clab as an initial response based on explicit modeling of infrastructure, workflow lifecycle, and artifacts. Lastly, we discuss how its evolution toward more realistic application settings can be understood through the lens of Translational Computer Science. We argue that reproducible continuum experimentation requires methods that are rigorous enough for research while remaining adaptable to real-world practice.

42 ENGINEERING↗

Deep Koopman Neural Network for Analyzing High-Energy-Density Simulations of Electrical Wire Explosions

Megaampere-scale electrical wire experiments (EWEs) provide a platform for studying magnetohydrodynamic (MHD) instability growth in magneto-inertial fusion (MIF) devices. Even when nonlinear simulations of these experiments can digitally reproduce much of the experimentally observed instability growth, interpreting the results and understanding mode growth and evolution can be non-trivial. As a first step toward providing better interpretation of these simulation features, this work investigates the use of a deep neural network that uses Koopman operator theory to analyze the dynamics of pulsed-power-driven explosions of EWEs. This deep neural network is trained on 1-D resistive MHD simulations of EWEs. This neural network learns to transform the nonlinear data into a lower-dimensional representation where the time dynamics are linear. Layers of this neural network are shown to learn features of the simulations, including the locations of shock waves and different physical regimes of the simulation. Using the learned features, the network can compress a time state of the simulation consisting of 5120 data point into a 36-parameter lower-dimensional latent space embedding. Furthermore, these embeddings are shown to be clustered in the latent space by initial radius and time state.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Multi-modelling predictions show high uncertainty of required carbon input changes to reach a 4‰ target

Soils store vast amounts of carbon (C) on land, and increasing soil organic carbon (SOC) stocks in already managed soils such as croplands may be one way to remove C from the atmosphere, thereby limiting subsequent warming. The main objective of this study was to estimate the amount of additional C input needed to annually increase SOC stocks by 4‰ at 16 long-term agricultural experiments in Europe, including exogenous organic matter (EOM) additions. We used an ensemble of six SOC models and ran them under two configurations: (1) with default parametrization and (2) with parameters calibrated site-by-site to fit the evolution of SOC stocks in the control treatments (without EOM). We compared model simulations and analysed the factors generating variability across models. The calibrated ensemble was able to reproduce the SOC stock evolution in the unfertilised control treatments. We found that, on average, the experimental sites needed an additional 1.5 ± 1.2 Mg C ha -1 year -1 to increase SOC stocks by 4‰ per year over 30 years, compared to the C input in the control treatments (multi-model median ± median standard deviation across sites). That is, a 119% increase compared to the control. While mean annual temperature, initial SOC stocks and initial C input had a significant effect on the variability of the predicted C input in the default configuration (i.e., the relative standard deviation of the predicted C input from the mean), only water-related variables (i.e., mean annual precipitation and potential evapotranspiration) explained the divergence between models when calibrated. Here, our work highlights the challenge of increasing SOC stocks in agriculture and accentuates the need to increasingly lean on multi-model ensembles when predicting SOC stock trends and related processes. To increase the reliability of SOC models under future climate change, we suggest model developers to better constrain the effect of water-related variables on SOC decomposition.

4 per 1000 initiative↗

Formation kinetics and thermodynamic stability of the Aurivillius compounds in Bi 4 Ti 3 O 12 –BiFeO 3 system

Abstract The Aurivillius compounds in the Bi 2 O 3 –Fe 2 O 3 –TiO 2 system, combining ferroelectric, semiconducting, and ferromagnetic properties, have attracted particular interest. Formation kinetics and thermodynamic stability are the fundamental knowledge needed for modeling and predicting the temporal microstructure and property evolution during materials processing but have not yet been addressed by quantitative experimental measurement. This article focuses on the Bi n +1 Fe n –3 Ti 3 O 3 n +3 Aurivillius compounds on the Bi 4 Ti 3 O 12 –BiFeO 3 tie‐line to elucidate the mechanisms and thermodynamic controls responsible for phase formation of compounds with various perovskite‐like layers. Five high‐purity Aurivillius compounds Bi 4 Ti 3 O 12 , Bi 5 FeTi 3 O 15 , Bi 6 Fe 2 Ti 3 O 18 , Bi 7 Fe 3 Ti 3 O 21 , and Bi 8 Fe 4 Ti 3 O 24 with integer n = 3–7 values were synthesized and their phase transformation properties and enthalpies of formation were studied by X‐ray diffraction in situ, high temperature differential scanning calorimetry, and high temperature oxide melt solution calorimetry. Thermodynamic stability of the compounds decreases with increasing n , and formation kinetics gradually slow down, demonstrating the inherent difficulty to synthesize pure Aurivillius compounds with n larger than 8. This difficulty was confirmed by an impurity phase coexisting with Bi 9 Fe 5 Ti 3 O 27 .

Gong, Weiping↗

Real-time (nanoseconds) determination of liquid phase growth during shock-induced melting

Melting of solids is a fundamental natural phenomenon whose pressure dependence has been of interest for nearly a century. However, the temporal evolution of the molten phase under pressure has eluded measurements because of experimental challenges. By using the shock front as a fiducial, we investigated the time-dependent growth of the molten phase in shock-compressed germanium. In situ x-ray diffraction measurements at different times (1 to 6 nanoseconds) behind the shock front quantified the real-time growth of the liquid phase at several peak stresses. These results show that the characteristic time for melting in shock-compressed germanium decreases from ~7.2 nanoseconds at 35 gigapascals to less than 1 nanosecond at 42 gigapascals. Our melting kinetics results suggest the need to consider heterogeneous nucleation as a mechanism for shock-induced melting and provide an approach to measuring melting kinetics in shock-compressed solids.

36 MATERIALS SCIENCE↗

Light transport with weak angular dependence in fog

Random scattering and absorption of light by tiny particles in aerosols, like fog, reduce situational awareness and cause unacceptable down-time for critical systems or operations. Computationally efficient light transport models are desired for computational imaging to improve remote sensing capabilities in degraded optical environments. To this end, we have developed a model based on a weak angular dependence approximation to the Boltzmann or radiative transfer equation that appears to be applicable in both the moderate and highly scattering regimes, thereby covering the applicability domain of both the small angle and diffusion approximations. An analytic solution was derived and validated using experimental data acquired at the Sandia National Laboratory Fog Chamber facility. The evolution of the fog particle density and size distribution were measured and used to determine macroscopic absorption and scattering properties using Mie theory. A three-band (0.532, 1.55, and 9.68 µ m) transmissometer with lock-in amplifiers enabled changes in fog density of over an order of magnitude to be measured due to the increased transmission at higher wavelengths, covering both the moderate and highly scattering regimes. The meteorological optical range parameter is shown to be about 0.6 times the transport mean free path length, suggesting an improved physical interpretation of this parameter.

74 ATOMIC AND MOLECULAR PHYSICS↗

Pore-Scale Oil Connectivity and Displacement by Controlled-Ionic-Composition Waterflooding Using Synchrotron X-Ray Microtomography

Controlled-ionic-composition waterflooding is an economic and effective method to improve oil recovery in carbonate oil reservoirs. Recent studies show controlling the salinity and ionic composition of injection water can alter the wettability of carbonate mineral surfaces. The pore-scale oil connectivity and displacement by controlled-ionic-composition waterflooding in heterogeneous carbonate reservoirs, especially at the early stage, is still unclear. The goal of this study is to examine the role of ion concentrations and types in the oil displacement efficiency and investigate the impact of the waterflooding on the pore-scale oil displacement using the national synchrotron facility. A carbonate rock sample was flooded with synthetic high-salinity water and other water solutions with different sulfate concentrations. The waterflooding processes were visualized with synchrotron X-ray microtomography to follow the evolution of pore-scale oil/brine interactions at typical field flow rates. Experimental results show that the water with lower sulfate concentration and higher salinity did not change the wettability of the pore surfaces. Higher sulfate ion concentrations in the water, in contrast, altered the wettability of carbonate pore surfaces from oil-wet to neutral-wet within the first few minutes of waterflooding. Finally, novel insight was gained on the ability of water with high-sulfate concentration to displace oil in the small pores and through abundant oil channels, which could consequently lead to higher oil recovery from the carbonate rock.

04 OIL SHALES AND TAR SANDS↗

High Performance SOFCs with a Superior Stability for Reliable and Durable Power Systems

Next generation of fuel cells, electrolyzers, and batteries requires higher power, faster kinetics, and larger energy density, which necessitate the use of compositionally complex oxides to achieve multifunctionalities and activity. These compositionally complex oxides may change their phases and structures during an electrochemical process – a so-called “electrochemically driven phase transformation”. The origin for such a phase change has remained obscure. More importantly, there is a need to develop high performance solid oxide fuel cells with an enhanced stability. In this work, the La 0.6 Sr 0.4 Co 0.2 Fe 0.8 O 3 (LSCF) cathode surface is modified by infiltration of Pr 6 O 11 and the power density at 0.8V and 750 °C is improved by 21%. In addition, by replacing the traditional barrier layer Gd 0.2 Ce 0.8 O 1.9 with mixed conducting Pr 0.1 Gd 0.1 Ce 0.8 O 1.9 , the power density increases by 38%. The different mechanism of promotions was investigated by electrochemical impedance spectroscopy. The ohmic resistance is dramatically reduced by applying the PGCO interlayer, and the distribution of relaxation time was used to analyze the mechanism for which the polarization resistance was decreased attributing to the mixed conduction nature in PrO x . An increase of power density at 0.8 V of 0.358 W/cm 2 (71%) is achieved with the implementation of both surface modification and buffer layer engineering. An experimental study and a theoretical analysis were then carried out on phase evolution in praseodymium nickelates. Nickelate-based electrodes show up to 60× greater phase transformation during operation when compared to thermally annealed ones. Theoretical analysis suggests that the presence of a reduced oxygen partial pressure at the interface between the oxygen electrode and the electrolyte is the origin for the phase change in an oxygen electrode. Guided by the theory, an addition of the electronic conduction in the interface layer leads to the significant suppression of phase change, while improving cell performance and performance stability. When an oxygen electrode is under polarization, the oxygen partial pressure at the interface between the oxygen electrode and the electrolyte is lower than that of incoming oxidant. Under a high polarization, the environment at the aforementioned interface may lead to phase transformation of the oxygen electrode. The local oxygen partial pressure is determined by the transport properties at the interfaces. An addition of the electronic conduction in the interface layer, for instance using (Pr,Gd)-doped ceria to replace Gd-doped ceria, results in improved cell performance and performance stability, while the phase transformation is significantly suppressed. This work provides a fundamental understanding of the origin for phase transformation in oxygen electrodes during operation and use this knowledge to develop a high-performance electrode that exhibits improved performance stability.

30 DIRECT ENERGY CONVERSION↗

Predicting Failure Using Deep Learning SAND Report

Accurate prediction of ductile failure is critical to Sandia’s NW mission, but the models are computationally heavy. The costs of including high-fidelity physics and mechanics that are germane to the failure mechanisms are often too burdensome for analysts either because of the person-hours it requires to input them or because of the additional computational time, or both. In an effort to deliver analysts a tool for representing these phenomena with minimal impact to their existing workflow, our project sought to develop modern data-driven methods that would add microstructural information to business-as-usual calculations and expedite failure predictions. The goal is a tool that receives as input a structural model with stress and strain fields, as well as a machine-learned model, and output predictions of structural response in time, including failure. As such, our project spent substantial time performing high-fidelity, three-dimensional experiments to elucidate materials mechanisms of void nucleation and evolution. We developed crystal-plasticity finite-element models from the experimental observations to enrich the findings with fields not readily measured. We developed engineering length-scale simulations of replicated test specimens to understand how the engineering fields evolve in the presence of fine-scale defects. Finally, we developed deep learning convolutional neural networks, and graph-based neural networks to encode the findings of the experiments and simulations and make forward predictions in time for structural performance. This project demonstrated the power of data-driven methods for model development, which have the potential to vastly increase both the accuracy and speed of failure predictions. These benefits and the methods necessary to develop them are highlighted in this report. However, many challenges remain to implementing these in real applications, and these are discussed along with potential methods for overcoming them.

97 MATHEMATICS AND COMPUTING↗

Evolution of eddy viscosity in the wake of a wind turbine

The eddy viscosity hypothesis is a popular method in wind turbine wake modeling for estimating turbulent Reynolds stresses. We document the downstream evolution of eddy viscosity in the wake of a wind turbine from experimental and large-eddy-simulation data. Wake eddy viscosity is isolated from its surroundings by subtracting the inflow profile, and the driving forces are identified in each wake region. Eddy viscosity varies in response to changes in turbine geometry and nacelle misalignment with larger turbines generating stronger velocity gradients and shear stresses. We propose a model for eddy viscosity based on a Rayleigh distribution. Model parameters are obtained from scaling the eddy viscosity hypothesis and demonstrate satisfactory agreement with the reference data. The model is implemented in the curled wake formulation in the FLOw Redirection and Induction in Steady State (FLORIS) framework and assessed through comparisons with the previous formulation. Our approach produced more accurate flow field estimates with lower total error for the majority of cases.

17 WIND ENERGY↗

Research and applications modules /RAM/ - A family of payload carriers.

RAM is a family of payload carriers, pressurized, unpressurized and pressurizable, that operate in conjunction with the Space Shuttle and associated space systems to provide modular facilities to conduct manned or man-tended research in low earth orbit. RAM family elements are discussed as the next step in evolution from the Sortie Laboratory required to provide a significant increase in experimental capability, flexibility, and economy necessary to meet the needs of an expanding program of research and applications activities. Examples of RAM missions and operations as sortie mission laboratories and observatories, and free-flying man-tended observatories are illustrated.

Withee, W. W.↗

Comparison of theory and experiment for homogeneous turbulence with shear

Solutions for uniformly-sheared turbulence, in which the interaction of the turbulence with the mean shear dominates the turbulent self-interaction, are compared with experiment. An anisotropic spectral tensor, which appears general enough to represent the initial experimental turbulence, is used for the initial condition in the calculations. The evolution of one-point turbulence components and microscales, as well as two-point velocity correlations, are considered. In most cases the agreement with experiment is good. The theory correctly predicts the presence of a negative region for two-point longitudinal-velocity correlations only for point separations in the direction normal to the flow and the mean gradient.

Deissler, R. G.↗

Modern developments in shear flow control with swirl

Passive and active control of swirling turbulent jets is experimentally investigated. Initial swirl distribution is shown to dominate the free jet evolution in the passive mode. Vortex breakdown, a manifestation of high intensity swirl, was achieved at below critical swirl number (S = 0.48) by reducing the vortex core diameter. The response of a swirling turbulent jet to single frequency, plane wave acoustic excitation was shown to depend strongly on the swirl number, excitation Strouhal number, amplitude of the excitation wave, and core turbulence in a low speed cold jet. A 10 percent reduction of the mean centerline velocity at x/D = 9.0 (and a corresponding increase in the shear layer momentum thickness) was achieved by large amplitude internal plane wave acoustic excitation. Helical instability waves of negative azimuthal wave numbers exhibit larger amplification rates than the plane waves in swirling free jets, according to hydrodynamic stability theory. Consequently, an active swirling shear layer control is proposed to include the generation of helical instability waves of arbitrary helicity and the promotion of modal interaction, through multifrequency forcing.

Farokhi, Saeed↗

A viscoplastic model with application to LiF-22 percent CaF2 hypereutectic salt

A viscoplastic model for class M (metal-like behavior) materials is presented. One novel feature is its use of internal variables to change the stress exponent of creep (where n is approximately = 5) to that of natural creep (where n = 3), in accordance with experimental observations. Another feature is the introduction of a coupling in the evolution equations of the kinematic and isotropic internal variables, making thermal recovery of the kinematic variable implicit. These features enable the viscoplastic model to reduce to that of steady-state creep in closed form. In addition, the hardening parameters associated with the two internal state variables (one scalar-valued, the other tensor-valued) are considered to be functions of state, instead of being taken as constant-valued. This feature enables each internal variable to represent a much wider spectrum of internal states for the material. The model is applied to a LiF-22 percent CaF2 hypereutectic salt, which is being considered as a thermal energy storage material for space-based solar dynamic power systems.

Freed, A. D.↗

Modern developments in shear flow control with swirl

Passive and active control of swirling turbulent jets is experimentally investigated. Initial swirl distribution is shown to dominate the free jet evolution in the passive mode. Vortex breakdown, a manifestation of high intensity swirl, was achieved at below critical swirl number (S = 0.48) by reducing the vortex core diameter. The response of a swirling turbulent jet to single frequency, plane wave acoustic excitation was shown to depend strongly on the swirl number, excitation Strouhal number, amplitude of the excitation wave, and core turbulence in a low speed cold jet. A 10 percent reduction of the mean centerline velocity at x/D = 9.0 (and a corresponding increase in the shear layer momentum thickness) was achieved by large amplitude internal plane wave acoustic excitation. Helical instability waves of negative azimuthal wave numbers exhibit larger amplification rates than the plane waves in swirling free jets, according to hydrodynamic stability theory. Consequently, an active swirling shear layer control is proposed to include the generation of helical instability waves of arbitrary helicity and the promotion of modal interaction, through multifrequency forcing.

Farokhi, S.↗

Analysis of the progressive failure of brittle matrix composites

This report investigates two of the most common modes of localized failures, namely, periodic fiber-bridged matrix cracks and transverse matrix cracks. A modification of Daniels' bundle theory is combined with Weibull's weakest link theory to model the statistical distribution of the periodic matrix cracking strength for an individual layer. Results of the model predictions are compared with experimental data from the open literature. Extensions to the model are made to account for possible imperfections within the layer (i.e., nonuniform fiber lengths, irregular crack spacing, and degraded in-situ fiber properties), and the results of these studies are presented. A generalized shear-lag analysis is derived which is capable of modeling the development of transverse matrix cracks in material systems having a general multilayer configuration and under states of full in-plane load. A method for computing the effective elastic properties for the damaged layer at the global level is detailed based upon the solution for the effects of the damage at the local level. This methodology is general in nature and is therefore also applicable to (0(sub m)/90(sub n))(sub s) systems. The characteristic stress-strain response for more general cases is shown to be qualitatively correct (experimental data is not available for a quantitative evaluation), and the damage evolution is recorded in terms of the matrix crack density as a function of the applied strain. Probabilistic effects are introduced to account for the statistical nature of the material strengths, thus allowing cumulative distribution curves for the probability of failure to be generated for each of the example laminates. Additionally, Oh and Finney's classic work on fracture location in brittle materials is extended and combined with the shear-lag analysis. The result is an analytical form for predicting the probability density function for the location of the next transverse crack occurrence within a crack bounded region. The results of this study verified qualitatively the validity of assuming a uniform crack spacing (as was done in the shear-lag model).

Thomas, David J.↗

Fundamental mechanisms in premixed flame propagation via vortex-flame interactions: Numerical simulations

The goal of the present study is to assess numerically the ability of single-step and two-step chemical models to describe the main features encountered during the interaction between a two-dimensional vortex pair and a premixed laminar flame. In the two-step mechanism, the reaction kinetics are represented by a first chain branching reaction A + X yields 2X and a second chain termination reaction X + X yields P. This paper presents the fundamental mechanisms occurring during vortex-flame interactions and the relative impact of the major parameters encountered in turbulent premixed flames and suspected of playing a role in quenching mechanism: (1) Influence of stretch is investigated by analyzing the contribution of curvature and tangential strain on the local structure of the flame. The effect of Lewis number on the flame response to a strained field is analyzed. (2) Radiative heat losses which are suspected to be partially or totally responsible for quenching are also investigated. (3) The effect of the diffusion of the radicals is studied using a two-step mechanism in which an intermediate species is present. The parameters of the two-step mechanism are entirely determined from physical arguments. (4) Precise quantitative comparisons between the DNS and the experimental results of Samaniego et al are performed. These comparisons concern the evolution of the minimum heat release rate found along the flame front during the interaction and the distribution of the heat release rate along the flame front.

Mantel, Thierry↗