Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “density distribution”

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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Assessment of WRF (v 4.2.1) dynamically downscaled precipitation on subdaily and daily timescales over CONUS

Abstract. This study analyzes the quality of simulated historical precipitation across the contiguous United States (CONUS) in a 12 km Weather Research and Forecasting model version 4.2.1 (WRF v 4.2.1)-based dynamical downscaling of the fifth-generation ECMWF atmospheric reanalysis (ERA5). This work addresses the following questions. First, how well are the 3 and 24 h precipitation characteristics (diurnal and annual cycles, precipitation frequency, annual and seasonal mean and maximum precipitation, and distribution of seasonal maximum precipitation) represented in the downscaled simulation, compared to ERA5? And second, how does the performance of the simulated WRF precipitation vary across seasons, regions, and timescales? Performance is measured against the National Centers for Environmental Prediction/Environmental Modeling Center (NCEP/EMC) 4 km Stage IV and Oregon State University Parameter-Elevation Regressions on Independent Slopes Model (PRISM) data on 3 and 24 h timescales, respectively. Our analysis suggests that the 12 km WRF exhibits biases typically found in other WRF simulations, including those at convection-permitting scales. In particular, WRF simulates both the timing and magnitude of the summer diurnal precipitation peak as well as ERA5 over most of the CONUS, except for a delayed diurnal peak over the Great Plains. As compared to ERA5, both the month and the magnitude of the precipitation peak annual cycle are remarkably improved in the downscaled WRF simulation. WRF slightly overestimates 3 and 24 h precipitation maximum over the CONUS, in contrast to ERA5, which generally underestimates these quantities mainly over the eastern half of the CONUS. Notably, WRF better captures the probability density distribution (PDF) of 3 and 24 h annual and seasonal maximum precipitation. WRF exhibits seasonally dependent precipitation biases across the CONUS, while ERA5's biases are relatively consistent year round over most of the CONUS. These results suggest that dynamical downscaling to a higher resolution improves upon some precipitation metrics but is susceptible to common regional climate model biases. Consequently, if used as input data for domain-specific models, we suggest moderate bias correction be applied to the dynamically downscaled product.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Molecular Structure and Conformational Design of Donor‐Acceptor Conjugated Polymers to Enable Predictable Optoelectronic Property

Abstract Tuning the optoelectronic properties of donor‐acceptor conjugated polymers (D‐A CPs) is of great importance in designing various organic optoelectronic devices. However, there remains a critical challenge in precise control of bandgap through synthetic approach, since the chain conformation also affects molecular orbital energy levels. Here, D‐A CPs with different acceptor units are explored that show an opposite trend in energy band gaps with the increasing length of oligothiophene donor units. By investigating their chain conformation and molecular orbital energy, it is found that the molecular orbital energy alignment between donor and acceptor units plays a crucial role in dictating the final optical bandgap of D‐A CPs. For polymers with staggered orbital energy alignment, the higher HOMO with increasing oligothiophene length leads to a narrowing of the optical bandgap despite decreased chain rigidity. On the other hand, for polymers with sandwiched orbital energy alignment, the increased band gap with increasing oligothiophene length originates from the reduction of bandwidth due to more localized charge density distribution. Thus, this work provides a molecular understanding of the role of backbone building blocks on the chain conformation and bandgaps of D‐A CPs for organic optoelectronic devices through the conformation design and segment orbital energy alignment.

Cao, Zhiqiang↗
DOE OSTI2022

Lewis Acid Site Assisted Bifunctional Activity of Tin Doped Gallium Oxide and Its Application in Rechargeable Zn‐Air Batteries

The enhanced safety, superior energy, and power density of rechargeable metal‐air batteries make them ideal energy storage systems for application in energy grids and electric vehicles. However, the absence of a cost‐effective and stable bifunctional catalyst that can replace expensive platinum (Pt)‐based catalyst to promote oxygen reduction reaction (ORR) and oxygen evolution reaction (OER) at the air cathode hinders their broader adaptation. Here, in this work, it is demonstrated that Tin (Sn) doped β‐gallium oxide (β‐Ga 2 O 3 ) in the bulk form can efficiently catalyze ORR and OER and, hence, be applied as the cathode in Zn‐air batteries. The Sn‐doped β‐Ga 2 O 3 sample with 15% Sn (Sn x =0.15 ‐Ga 2 O 3 ) displayed exceptional catalytic activity for a bulk, non‐noble metal‐based catalyst. When used as a cathode, the excellent electrocatalytic bifunctional activity of Sn x =0.15 ‐Ga 2 O 3 leads to a prototype Zn‐air battery with a high‐power density of 138 mW cm −2 and improved cycling stability compared to devices with benchmark Pt‐based cathode. The combined experimental and theoretical exploration revealed that the Lewis acid sites in β‐Ga 2 O 3 aid in regulating the electron density distribution on the Sn‐doped sites, optimize the adsorption energies of reaction intermediates, and facilitate the formation of critical reaction intermediate (O*), leading to enhanced electrocatalytic activity.

Lewis acid sites↗
DOE OSTI2021

A Time-Dependent Directional Damage Theory for Brittle Rocks Considering the Kinetics of Microcrack Growth

In this work, a novel microcrack damage theory for describing the time-dependent behavior of brittle rocks is proposed. Instead of using scalar or tensorial variables to approximate the effect of microcracks, the concept of directional damage is introduced by upscaling the directional distribution density of microcracks through the noninteracting homogenization scheme. The contribution of an individual crack on global strain is first derived from the crack-opening displacements including sliding, closure, and dilation. The model is then cast in the free energy form, and the range of parameters ensuring compliance with the second law of thermodynamics is derived rigorously. Considering that the subcritical propagation of microcracks is the dominant mechanism behind the time-dependent deformation and failure of brittle rocks, the driving force for directional damage is defined as the corresponding energy release rate, and the effect of time is introduced through a damage evolution law inspired from the kinetics of subcritical crack growth. The model is applied to simulate the behavior of basalt rock under various loading conditions. In addition to predicting the stress–strain–time response, the model offers enhanced resolution in representing the anisotropic characteristics of microcrack-induced damage in brittle rocks and their time-dependent evolutions.

58 GEOSCIENCES↗
DOE OSTI2020

Finite Element Analysis and Techno-economic Modeling of Solar Silicon Molten Salt Electrolysis

A new process is presented for low-cost one-step production of pure solid silicon from natural quartzite by molten salt electrolysis. At a process temperature of 1100°C, a techno-economic model including detailed mass and energy balances estimates energy consumption below 15 kWh/kg, with operating cost of $\$1.74$/kg and capital cost around $\$10,500$ per t/a (tonne annually) of production capacity for a 160,000 t/a plant. Use of an inert solid oxide membrane anode results in a pure oxygen by-product and no direct emissions. As a result, finite element analysis estimates the current density distribution and total current to inform the design of slab-shaped solid silicon cathodes.

14 SOLAR ENERGY↗
DOE OSTI2024

Physics-coupled data-driven design of high-temperature alloys

We present a materials design loop, which streamlines physics-coupled machine learning (ML) surrogate models to discover new alloy chemistries with improved properties. The efficacy is demonstrated by discovering a high-temperature alumina-forming austenitic (AFA) stainless steel with enhanced creep, followed by experimental validation. The ML models have been trained using a well-curated, highly consistent experimental dataset augmented with synthetic microstructural features from a computational thermodynamic approach. We have populated a large number of hypothetical AFA alloys to explore the high-dimensional composition space and have predicted their creep properties by providing the same synthetic input features obtained from the trained ML models. Uncertainties from the ML training were taken as thresholds for truncating predicted results to identify alloys with improved or deteriorated creep. Individual elemental compositions have been determined via probability density distribution analysis from the group of alloys at the top and bottom of the predicted creep values for further virtual and experimental validations. In conclusion, we anticipate that this workflow can be applied to screen desired conditions, such as chemistry and processing parameters, in high-dimensional space through physics-guided data analytics.

Alloy design↗
DOE OSTI2021

Theoretical predictive screening of noble-metal-containing M 3 AuC 2 (M = Ti, V, and Cr) MAX phases

Using first-principles density functional theory (DFT), the ground state physical properties of the newly synthesized noble-metal-containing Ti 3 AuC 2 MAX phase have been investigated. The effect of transition-element Ti replacement (V and Cr) on physical properties, including structural, elastic, electronic, thermal, and optical are presented. The optimized lattice parameters of Ti 3 AuC 2 MAX phase are in good agreement with the experimental values and decrease when we replace Ti with V and then Cr. The magnetic properties of M 3 AuC 2 are predicted by using generalized gradient approximation (GGA) incorporating onsite electron correlation parameter (Hubbard U parameter) within DFT. Out of three studied MAX phases, V 3 AuC 2 and Cr 3 AuC 2 phases exhibit magnetism. The thermodynamical stability of M 3 AuC 2 (M = Ti, V, and Cr) phases is discussed by the formation enthalpy with respect to their most competing phases. The electronic structure reveals that these phases have metallic nature and are electrically anisotropic. The M-element replacement has an effect on the bonding properties of M 3 AuC 2 . The bonding between M-C is in the order of Cr 3 AuC 2 > V 3 AuC 2 > Ti 3 AuC 2 according to their peak positions and heights of density of states in the occupied site, which is also confirmed from the charge density distribution. The elastic properties indicate that M 3 AuC 2 phases are ductile, machinable, less stiff, and better resistant to thermal shock and are elastically anisotropic. To understand the properties of M 3 AuC 2 for the extreme environment, we employed a quasi-harmonic Debye model at the pressure and temperature range of 0–50 GPa and 0–1600 K, respectively. Additionally, we evaluated thermal conductivity and melting temperature to further understand the potential of these MAX phases in coating applications at elevated temperature. In relation to optical properties, these MAX phases have a reasonable absorption coefficient in the visible as well as in the ultra-violet regions. The reflectivity of M 3 AuC 2 phases is polarization dependent and of 45 % against the visible region, which reveals its potential as coating material to minimize solar heating.

312 MAX phases↗
DOE OSTI2021

Effects of fiber dosage, loading orientation and stress on frequency response of enhanced Carbon Nano-Fiber Aggregates

Nanomaterial additives have been broadly used in cement and concrete-based sensors to measure the stress and strain in the structure. Most of the studies focused on measuring the resistivity of sensors with Direct Current (DC) or relatively low-frequency. Here in this work, the impedance responses of enhanced Carbon Nano-Fiber Aggregates (CNFAs) at different frequencies are rigorously studied and utilized as a tool for real-time structural health monitoring (SHM). The CNFA impedance measurements at different measurement frequencies provide an extra dimension (frequency dimension) to allow the characterization of different phenomena and multimodal measurements simultaneously. First, a comprehensive physical model is established to understand the frequency responses of the CNFA impedance. Second, COMSOL Multiphysics is used to simulate the frequency responses of the CNFA impedance, which provides more insights into the electric field and current density distribution. The impedance of CNFAs decreased with the increased dosage of carbon nanofibers (CNFs). Third, experimental studies are reported in detail. The electrical impedance variation (EZV) of the CNFA in parallel orientation is 151% higher than that of perpendicular orientation at the frequency of 2612 Hz and uniaxial compression stress of 4.65 MPa. The strain-EZV curve obtained from the response spectrum is linear for frequencies ranging from 98 Hz to 463.9 kHz. The strain-EZV curve for 5625 Hz has a linear fit with the gauge factor of 147.78. At 1000 Hz, the CNFA exhibited a repeatable behavior up to 9.35 MPa and detection limit up to the stress of 18.62 MPa.

36 MATERIALS SCIENCE↗
DOE OSTI2020

CAKE: Consistent Automatic Kinetic Equilibrium reconstruction

Magnetohydrodynamic equilibrium reconstructions with kinetic constraints are an essential input of many tokamak stability and transport analysis work- flows. The reconstructions themselves also contain useful information about pressure and current density distributions within the plasma which can be used to drive other calculations. While standardized magnetic-only, and partially kinetic equilibrium reconstruction is widely available and used, re- constructions for well diagnosed machines such as DIII-D, that incorporates a large and heterogeneous set of internal diagnostic measurements yet robustly reaching high numerical accuracy have been a time intensive, manual process. The Consistent Automatic Kinetic Equilibrium reconstruction tool (CAKE) has been developed and implemented at DIII-D for producing low error, kinetically constrained reconstructions without human intervention, and to identify strategies for doing so that would be broadly applicable. Here, the techniques employed to allow CAKE to handle data of limited quality and availability, as well as the full automation of the workflow may assist future efforts to form kinetic constraint profiles in real time plasma control applications.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Fission gas retention of densely packed uranium carbonitride tristructural-isotropic fuel particles in a 3D printed SiC matrix

The Transformational Challenge Reactor (TCR) fuel form was designed to contain large, densely packed uranium carbonitride (UCN) tristructural-isotropic (TRISO) fuel particles within a 3D printed SiC matrix, increasing the uranium density compared to conventional TRISO fuel forms and offering full geometric freedom for core design. Here, this work summarizes initial low-burnup, high-power irradiation testing of TCR fuel materials, including loose UCN TRISO particles and integral fuel compacts with ~55% TRISO particles by volume, to evaluate fission gas retention. Fission gasses were fully retained in all loose particle tests and in integral compacts irradiated at low (<250 °C) surface temperatures. Initial testing at higher (~700–750 °C) fuel surface temperatures showed fission gas release (FGR) and complete fracture of three compacts, but no FGR was observed in later high temperature tests (~300–750 °C) of both fueled compacts and loose TRISO particles. Calculated thermal stresses in the failed compacts were far less than the measured strength of the SiC matrix and the stresses in some failed compacts were less than those in compacts that did not show FGR. Thermal stress-induced matrix cracks also would not cause complete fracture because the tensile stresses transition to compression in the higher temperature regions. Therefore, fuel failure was likely not caused by thermal stresses and may have been related to leakage currents from the electrical heaters and erratic fuel surface temperatures that were only observed in the test for which failure was observed. In any case, the matrix cracks propagated through the coatings of TRISO particles located in the high-density matrix regions on the peripheries of the compacts, resulting in measurable fission gas release. The discussion focuses on the importance of understanding matrix density distributions and the particle-matrix interface properties to prevent matrix cracks from causing TRISO particle failures.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

Micro-beam bending of FCC bicrystals: A comparison between defect dynamics simulations and experiments

To understand the role of the grain boundary (GB) in plasticity at small scale, a concurrently coupled mesoscale plasticity model was developed to simulate micro-bending of bicrystalline micron-sized beams. By coupling dislocation dynamics (DD) with a finite element model (FEM), a novel defect dynamics model provides the means to investigate intricate interactions between dislocations and GBs under various loading conditions. Our simulations of micro-bending agree well with corresponding micro-bending experiments, and they show that mechanical response of bicrystals could have not only hardening but also softening depending on the characters of the GB. In addition, changing the location of the GB in the microbeams results in different mechanical responses; GBs located at the neutral plane show softening compared to single crystals, while inclined GBs located halfway along the length of the beam show little effect. In conclusion, simulation results could provide a clear picture on detailed dislocation-GB interactions, and quantitative resolved shear stress analysis supplemented by dislocation density distribution is used to analyze the mechanical response of bicrystalline samples.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Effects of dissipative baryon current in heavy-ion collisions at RHIC-BES energies

The CLVisc (3+1)D viscous hydrodynamic model is extended to include the equation of net baryon conservation and the Israel-Stewart-like equations for dissipative baryon current. Using the NEOSB equation of state, we simulate the dynamical evolution and collectivity of the quark-gluon plasma with finite chemical potential, assuming smooth energy density and net baryon density distributions at the initial proper time. Numerical results are shown for the impact of net-baryon dissipation on particle yields and p T spectra in heavy-ion collisions at beam energy scan energies.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

A machine learning study on spinodal clumping in heavy ion collisions

Possible observables of baryon number clustering due to the instabilities occurring at a first order QCD phase transition are discussed. The dynamical formation of baryon clusters at a QCD phase transition can be described by numerical fluid dynamics, augmented with a gradient term and an equation of state with a mechanically unstable region. It is shown that the dynamical description of this phase transition, in nuclear collisions, will lead to the formation of dense baryon clusters at the phase boundary. State-of-the-art machine learning methods find that the coordinate space clumping leaves characteristic imprints on the spatial net density distribution in almost every event. On the other hand the momentum distributions do not show any clear event-by-event features. Lastly, it is shown that the 'third order' cumulant, the skewness, shows a peak at the beam energy where the system, created in the heavy ion collision, reaches the deconfinement phase transition.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Leptodermic corrections to the TOV equations and nuclear astrophysics within the effective surface approximation

The macroscopic model for a neutron star (NS) as a liquid drop at the equilibrium is used to extend the Tolman-Oppenheimer-Volkoff (TOV) equations taking into account the gradient terms responsible for the system surface. The parameters of the Schwarzschild metric in the spherical case are found with these surface corrections to the known leading (zero) order of the leptodermic approximation a/R << 1, where a is the NS effective-surface (ES) thickness, and R is the effective NS radius. The energy density $\mathscr{E}$ is considered in a general form including the functions of the particle number density and of its gradient terms. The macroscopic gravitational component $Φ$(ρ) of the energy density is taken into account in the simplest form as expansion in powers of $ρ$ – $\overline{ρ}$, where $\overline{ρ}$ is the saturation density, up to second order, in terms of its contributions to the separation particle energy and incompressibility. Density distributions ρ across the NS ES in the normal direction to the ES, which are derived in the simple analytical form at the same leading approximation, was used for the derivation of the modified TOV (MTOV) equations by accounting for their NS surface corrections. As a result, the MTOV equations are analytically solved at first order and the results are compared with the standard TOV approach of the zero order.

Magner, A. G. [Institute for Nuclear Research, Kyi↗
DOE OSTI2021

Estimating the color lifetime of energetic quarks

Using a simple geometric framework with a realistic nuclear density distribution, we fit published HERMES data to determine fundamental properties of hadronization using the nuclear medium as a spatial analyzer. Our approach uses a fit to the transverse momentum broadening observable and the hadronic multiplicity ratio; the simultaneous fit to two different observables strongly constrains the outcome. Using the known sizes of the target nuclei, we extract the color lifetime, finding a z h -dependent range of values from 2 to 8 fm/c for these data. We also extract estimates for the qˆ transport coefficient characterizing the strength of the interaction between the quark and the cold nuclear medium, finding an average value of 0.035±0.011 GeV 2 /fm. With a three-parameter model we obtain satisfactory fits to the data with a goodness-of-fit parameter χ 2 /dof of 1.1 or less. In a secondary fit of the results from that model we independently find a value for the Lund String Model string tension of 1.00±0.05 GeV/fm. We evaluated the sensitivity for extracting quark energy loss and effective in-medium hadronic cross sections using four-parameter variants of the model, finding large uncertainties in both cases. Our results suggest that hadronic interaction of forming hadrons in the nuclear medium is the primary dynamical cause of meson attenuation in the HERMES data, with quark energy loss playing a more minor role.

Electron-ion collider↗
DOE OSTI2022

Probing neutron-skin thickness with free spectator neutrons in ultracentral high-energy isobaric collisions

We show that the yield ratio of free spectator neutrons produced in high-energy 96 Zr+ 96 Zr to 96 Ru+ 96 Ru collisions is a clean probe of the neutron-skin thickness of colliding nuclei and the slope parameter L of the symmetry energy. The idea is demonstrated based on the proton and neutron density distributions via a state-of-the-art Skyrme-Hartree-Fock-Bogolyubov calculation. Among spectator nucleons given by the Glauber model, free spectator neutrons include those from direct production that survive from clusterization as well as those from deexcitation of heavy clusters described by the popular GEMINI model. More free neutrons are produced in collisions of 96 Zr nucleus due to its larger neutron skin, compared to those produced in collisions of 96 Ru nucleus with a smaller neutron skin. The difference of the free spectator neutron yield is further increased with the increasing difference of the neutron-skin thickness between 96 Zr and 96 Ru with a larger L value, and the increase in ultracentral collisions is particularly insensitive to model details and experimental uncertainties. Since the production of free spectator neutrons is not affected by the complicated dynamics in the mid-rapidity region, the ratio of their multiplicities in ultracentral isobaric collisions is a robust observable for constraining the neutron skin and L value.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Porosity evolution under increasing tension in wire-arc additively manufactured aluminum using in-situ micro-computed tomography and convolutional neural network

Internal defects such as porosities are often formed in additively manufactured metal components. The pores nucleate, grow, and coalesce to form cracks under loads, leading to eventual catastrophic failure. Here, in this paper, the full-field porosity evolution, including pore growth and coalescence in a wire-arc additively manufactured (WAAM) aluminum alloy cylinder under tension is observed with in-situ X-ray micro-computed tomography (μCT). The pore size distribution, density, and tensile stress are calculated from the volumetric images analyzed by a convolutional neural network (CNN) algorithm, which provides rapid analysis of 12,950 slice images from μCT volumetric images at the reference state and 13 tensile strains. The results show the quantitative evolution of the growth and coalescence of macropores under tension. A strong correlation is found between the local pore volume fraction and the true tensile stress when the tensile strain is larger than 5%.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Fractal analysis on Ag 2 O thin film using a data-driven approach

The synthesis of fractal Ag oxide (Ag 2 O) on the surface of Ag thin film has been achieved at room temperature by using Synchrotron X-ray irradiation. We have performed an automated quantitative analysis of a batch of 1879 fractal Ag 2 O patterns in a scanning electron microscopy (SEM) image within a radius of 2000 mm outward from the center of the X-ray beam. The morphology is similar to that of the diffusion-limited cluster aggregation (DLCA) model. The fractal dimension (D) of Ag 2 O is between 1.7 and 1.5 from the center to the edge. The area distribution density of fractal Ag 2 O follows a quadratic function with radius R. It is found that the branches’ number of fractal Ag 2 O is a key factor affecting the fractal dimension. The more branches the fractal has, the greater the D is. This is the first time that Ag fractal has been investigated by combining automated data analysis methods with batch experimental data. This data-driven approach provides a new research perspective for rationally regulating materials’ fractal morphology and performance.

36 MATERIALS SCIENCE↗