Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Conditional 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 451 records · Page 25
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Phase field simulation of anode microstructure evolution of solid oxide fuel cell through Ni(OH) 2 diffusion

Microstructure evolution in a solid oxide fuel cell anode is strongly affected by operating condition. The detailed mesoscale evolution mechanism under operating condition is still under debate. Here we develop a phase-field model to simulate microstructure evolution through formation and diffusion of gaseous Ni(OH)2. We studied the coarsening kinetics and redistribution of Ni under different steam distributions and compare it to that under pure hydrogen condition. Overall, the results suggest that although the presence of a steady gradient of steam leads to redistribution of Ni, Ni(OH) 2 formation and diffusion do not significantly change the Ni coarsening rate under strictly humid conditions, contrary to commonly reported hypotheses. It is concluded that competing mechanisms other than Ni(OH) 2 diffusion must be responsible for the experimentally observed Ni redistribution and enhanced coarsening under humid conditions.

25 ENERGY STORAGE↗
DOE OSTI2026

Multipoint Aerostructural Optimization of Wind Turbine Rotors Using a Coupled Blade‐Resolved Aerostructural Solver

Physics‐based design optimization workflows thread the needle between computational cost limitations and simulation complexity, often compromising between modeling detail and the range of operating design conditions. Multipoint aerostructural optimization of wind turbine rotors has so far been confined to low‐fidelity analyses or to high‐fidelity studies with simplified structural models, leaving the most complex design trade‐offs unexplored. We close this gap by performing the first tightly coupled gradient‐based multipoint aerostructural rotor optimization using 3D aerodynamic and structural solvers with discrete coupled adjoints. The optimizer simultaneously varies blade planform, airfoil shapes, and structural thickness through more than 270 design variables, minimizing a weighted combination of rotor mass and power across multiple wind speeds. Applied to a modified DTU 10‐MW benchmark under conservative structural and aerodynamic constraints, our multipoint optimization reduces rotor mass by up to 36% and increases power by 12%–15% across the main operating conditions; biasing the objective toward power yields power gains up to 18% and a 17% mass reduction. For a nominal wind distribution, 3‐point rotor designs accounting for low RPM and high thrust conditions capture dominant trade‐offs and outperform single‐point designs. Adding two off‐design points changes individual‐condition power by less than 3% but leaves the weighted average within 0.5%, and the mass‐power bias has a stronger effect on the final design than the operating‐point weighting itself. Our framework extends naturally to richer load cases and site‐specific wind distributions, providing a basis for high‐fidelity multipoint design earlier in industrial workflows.

17 WIND ENERGY↗
DOE OSTI2024

An implicit solution for Asay foil trajectories generated by separable, sustained-production ejecta source models

We present a simple implicit solution for the time-dependent trajectory of a thin Asay foil ejecta diagnostic for the general case where the impinging ejecta cloud is generated by a source function characterized by an arbitrary (sustained) time dependence and a time-independent (stationary) particle velocity distribution. In the limit that the source function time dependence becomes a delta function, this solution—which is amenable to rapid numerical calculations of arbitrary accuracy—exactly recovers a previously published solution for the special case of instantaneous ejecta production. We also derive simple expressions for the free-surface arrival (catch-up) time as well as the true ejecta areal mass accumulation on the accelerating foil and place bounds on the level of error incurred when applying instant-production mass solutions to a sustained-production trajectory. We demonstrate these solutions with example calculations for hypothetical source functions spanning a wide range of ejecta production durations, velocity distributions, and temporal behaviors. These calculations demonstrate how the foil trajectory is often insensitive to the temporal dependence of the source function, instead being dominated by the velocity distribution. We quantify this insensitivity using a “compatibility score” metric. Under certain conditions, one may capitalize upon this insensitivity to obtain a good approximation of the second integral of the velocity distribution from the observed foil trajectory.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Modeling of H2 Dispersion at ARIES

Hydrogen is a versatile and clean energy carrier that can be produced from various renewable sources such as wind, solar, and hydropower and help decarbonize electricity grids, industry, and transportation. Using the Hydrogen Research Facility under Advanced Research on Integrated Energy Systems (ARIES) at the National Renewable Energy Laboratory's (NREL) Flatirons campus as a test bench, the study examines the feasibility, useability, and value of using computational fluid dynamics (CFD) techniques to model hydrogen dispersion. The ARIES facility was chosen because controlled hydrogen releases can be performed at a rate of 27 kg-H2/hr. Site-specific atmospheric and weather condition data such as wind speed and temperature were used as inputs to the model. The results show statistical distributions and ranges of hydrogen concentrations at locations throughout the domain. Wind conditions are found to significantly impact the release behavior, including the hydrogen cloud's direction and concentrations. At low wind speeds (below 1 mph), hydrogen forms a cloud and at higher wind speeds (> 2-4 mph) hydrogen plume stretches in the direction of wind momentum. From >100 simulations for ARIES site-specific conditions, statistical quantities combined with a clustering algorithm were used to propose sensor location at various elevations from ground.

dispersion↗
DOE OSTI2026

Plutonium Retention by Crystalline Silicotitanate under Hyperalkaline Conditions Relevant to Tank-Side Cesium-Removal at the Hanford Site

Crystalline silicotitanate (CST) is used in Hanford’s Tank-Side Cesium-Removal (TSCR) process to selectively remove Cs-137 from highly caustic, nitrate-rich tank supernatants. Recent testing with actual waste samples suggests that CST can also retain measurable plutonium (Pu), which could affect radiological classification and disposal pathways for spent CST. To quantify this behavior, Pu partitioning to CST was studied under Hanford-relevant conditions using batch-contact experiments in a representative simulant (2 M NaNO3, 0.7 M NaOH). Isotherm data were measured and distribution ratios calculated, with Cs+ uptake used as benchmark. Under low-carbonate conditions, Pu was retained strongly by CST in systems initially contacted with either PuO2 nanoparticles (Pu(IV)) or aqueous Pu(VI), with distribution ratios of ~2,200–3,700 mL/g, generally exceeding those for Cs+ (~400–1,000 mL/g). Increasing carbonate concentration strongly reduced PuO2 nanoparticle retention; at [Na2CO3] = 1 M, distribution ratios decreased by up to one order of magnitude to roughly 100–300 mL/g. Electron microscopy suggests that Pu retention involves a combination of mechanisms such as PuO2 NP aggregation induced by CST leachate components, and association with CST bead surfaces.

Neumann, J.↗
DOE OSTI2025

Solving high-dimensional inverse problems using amortized likelihood-free inference with noisy and incomplete data

Here, we present a likelihood-free probabilistic inversion method based on normalizing flows for high-dimensional inverse problems. The proposed method is composed of two complementary networks: a summary network for data compression and an inference network for parameter estimation. The summary network encodes raw observations into a fixed-size vector of summary features, while the inference network generates samples of the approximate posterior distribution of the model parameters based on these summary features. The posterior samples are produced in a deep generative fashion by sampling from a latent Gaussian distribution and passing these samples through an invertible transformation. We construct this invertible transformation by sequentially alternating conditional invertible neural network and conditional neural spline flow layers. The summary and inference networks are trained simultaneously. We apply the proposed method to an inversion problem in groundwater hydrology to estimate the posterior distribution of the log-conductivity field conditioned on spatially sparse time-series observations of the system’s hydraulic head responses. The conductivity field is represented with 706 degrees of freedom in the considered problem. Comparison with the likelihood-based iterative ensemble smoother PEST-IES method demonstrates that the proposed method accurately estimates the parameter posterior distribution and the observations’ predictive posterior distribution at a fraction of the inference time of PEST-IES.

conditional invertible neural network↗
DOE OSTI2023

Modeling of Stress Distribution in Molten Salt Thermal Energy Storage Tanks for In-Service Central Receiver Power Plants

Failures in molten nitrate salts thermal energy storage tanks (TES) have been occurring in several concentrating solar power (CSP) plants around the world after a few months or years of operation. These failures are mainly related to a combination of high stress, corrosion, large deformation, and thermal cycling. Most of these failures are a consequence of the infancy of this technology; technical standards for design, fabrication, implementation, and operation of molten salt tanks operating at temperatures up to 565 degrees Celsius do not exist and most of current in-service tanks have been designed using the recommendations from the American Petroleum Institute (API) 650 and the American Society of Mechanical Engineers ASME section II standards, which are applicable for TES tanks operating at lower temperatures. We have developed a 3-dimensional model of a representative molten salt hot tank design that uses typical CSP plant operation conditions and allows determining the temperature and stress distribution in the tank. The model also considers the initial floor residual stresses and distortion after its welding fabrication. We have identified critical thermal gradients in the tank floor and shell and locations in the tank where the stresses are high and could become susceptible points for failure.

concentrating solar power↗
DOE OSTI2021

Closures and Simulation for Thermal Radiation Transport in Stochastic Media with Nonlinear Temperature Dependence

Because of the practical challenge of rendering very complex realistic spatial structures for numerical work, it is common practice to resort to characterizing such media as stochastic mixtures of materials, ideally parametrized with low order statistics such as the mean, variance, and correlation functions of the now random material properties. This enables realizations of the medium to be repeatedly generated and radiation transport computations to, in principle, be performed for a large ensemble of these realizations to obtain a statistically well-characterized radiation field. Statistical post-processing yields desired quantities such as conditional and unconditional mean radiation flux and probability distributions of transmitted radiation. However, such computations prove expensive for all but the simplest stochastic geometries and are most suited for benchmarking approximate models. The most common approximations lead to homogenized media so that transport computations are required only on a single medium realization but by construct provide only limited statistical information on the radiation field. Almost all approximate approaches to this problem attempt to develop equations for low order moments of the radiation intensity (mean, second moment, correlation function) but inevitably encounter a closure problem: the equation for any statistical moment will contain terms depending on unknown higher-order moments. Thus, the challenge shifts to one of developing closure relations that relate the unknown moments to the lower order moments. Under very special conditions, an exact closure can be derived but in general closures are heuristically stated constitutive relations. Also, closure approaches depend on whether the mixing statistics are spatially and/or temporally continuous as in fluctuating turbulent fields, or discontinuous as in randomly mixed solid chunks of material. Thus, unconditional averaging is generally applied in the former case but conditional averaging is more appropriate when the mixing is discontinuous. In this work, the emphasis is on binary statistical mixtures of immiscible fluids as as such the quantities of interest are averages (flux, temperature) conditioned on the material type.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Hardness Prediction by Incorporating Heat Transfer and Molten Pool Fluid Flow in a Multi-pass, Multi-layer Weld for Onsite Repair of Grade 91 steel

Introduction: In the current fleet of fossil-fired power plants, creep strength enhanced ferritic steels (CSEF) are used to sustain the harsh service conditions. Enhanced properties of Grade 91 steel result from tempered martensite with a fine distribution of MX and M23C6 carbides. Grade 91 steel is subjected to onsite welding repair to remedy their degradation due to extreme service condition. Knowledge of weld repairability of these steels, such as as-welded hardness distribution, is essential to establishing sound repair procedures. Experimental trial and error tests can consume a lot of time as many welding variables need to be studied. For numerical modelling, most of the multi-pass multi-layer models are based on finite element method, which are limited to solve the heat conduction equation and ignore convective heat transfer due to melt flow. Moreover, the mesh has to be pre-built based on a known or assumed weld cross-section geometry. These finite element based models thus have limited predictive capability as defects are not considered and nugget size are pre-assumed. This research aims at developing a thermal and microstructure evolution model incorporating molten pool dynamics in a multi-pass multi-layer material deposition to predict the as-welded hardness distribution. Technical Approach: All the thermal, physical, and metallurgical properties of Grade 91 as a function of temperature are collected from the literature and inputted into the thermo-fluid model based on Flow-3D, a computational fluid dynamics software. A multi-pass, multi-layer material deposition is simulated where the melting of filler wire into the molten pool is directly considered based on the volume of fluid (VOF) method. The flow behaviour of the molten pool is used to understand the formation of deposition geometry and defects. The temperature profiles during the multi-pass, multi-layer welding are calculated. The results computed using the new model are compared against the experimental data of fusion zone geometry and thermal cycles. Hardness prediction in the heat-affected zone (HAZ) are made using Johnson-Mehl-Avrami (JMA) equation for solid-state phase transformation kinetics. The JMA parameters are extracted from the experimental data available in the literature. For comparison, a standard finite element heat conduction model is also developed to predict the thermal cycles and hardness distribution in the multi-pass, multi-layer weld. Expected Result: Results obtained using the molten pool dynamic simulation versus the finite element heat conduction model are compared. Specifically, the effects of convective heat transfer on the accuracy of the calculated thermal history, bead shape and size, and HAZ hardness distribution are examined.

Modeling, Grade 91 steel, Multi-pass multi-layer d↗
DOE OSTI2024

Space-Time Structured Plasma Waves

Electrostatic waves play a critical role in nearly every branch of plasma physics from fusion to advanced accelerators, to astro, solar, and ionospheric physics. The properties of planar electrostatic waves are fully determined by the plasma conditions, such as density, temperature, ionization state, or details of the distribution functions. Furthermore we demonstrate that electrostatic wave packets structured with space-time correlations can have properties that are independent of the plasma conditions. For instance, an appropriately structured electrostatic wave packet can travel at any group velocity, even backward with respect to its phase fronts, while maintaining a localized energy density. These linear, propagation-invariant wave packets can be constructed with or without orbital angular momentum by superposing natural modes of the plasma and can be ponderomotively excited by space-time structured laser pulses like the flying focus.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Hollow-structured Ni-N-C catalysts for highly selective CO 2 electroreduction

Atomically dispersed single-atom catalysts have emerged as promising non-precious catalyst alternatives to expensive Ag and Au catalysts for electrochemical CO 2 reduction reaction (CO 2 RR). In particular, nickel-nitrogen-carbon (Ni-N-C) catalysts have demonstrated a high faradaic efficiency (FE) toward CO formation at low overpotentials. Nonetheless, the exact nature of Ni active sites under CO 2 RR remains elusive and conventional Ni-N-C catalysts are limited by microporosity and low density of Ni single atoms, hindering performance in CO 2 electrolyzers. Here, we report the synthesis of hollow-structured Ni-N-C ( hs -Ni-N-C) catalysts via a post-synthesis modification (PSM) strategy using partial ligand exchange of 2-methylimidazole with 3-amino-1,2,4-triazole. This approach enables the formation of a hollow structure, resulting in more than a twofold increase in Ni atom density compared to regular Ni-N-C (r-Ni-N-C). In a zero-gap CO 2 electrolyzer, the optimized hs -Ni-N-C allows for achieving an FE CO of 97% at a current density of > 100 mA cm⁻ 2 , while maintaining high CO selectivity with stable performance over 100 h at 2.5 V. hs-Ni-N-C shows a more than sevenfold increase in the CO partial current density relative to r-Ni-N-C resulting from the combined effects of a higher density of Ni single-atom sites, improved kinetics, and lower transport resistance under the operating conditions, as indicated by electrochemical impedance spectra and distribution of relaxation times analysis. Operando high energy-resolution X-ray absorption spectroscopy (XAS) reveals that atop-bonded CO on Ni single sites induces dynamic transformations of the Ni–N coordination environment, leading to a symmetric coordination structure of hs -Ni-N-C. Under CO 2 RR, the catalysts undergo a more pronounced structural change and form a minor fraction of Ni nanoparticles. Density functional theory calculations are consistent with the XAS results and provide molecular insights showing that the interplay between protonation and CO adsorption leads to adsorbate-induced restructuring of the Ni single atom. This work demonstrates the synergistic role of hollow structure and high-density Ni atoms in governing CO 2 RR selectivity and provides mechanistic insights into the structural dynamics of single-atom catalysts under operating conditions.

36 MATERIALS SCIENCE↗
DOE OSTI2020

A numerical study of the meso-structure variability in the compaction process of prepreg platelet molded composites

Compaction deformations in a thermoplastic composite system molded from chopped prepreg platelets were analyzed to obtain statistical distributions of local distortions that define the final consolidated meso-structure. This study provides an insight into and allows quantification of the intrinsic geometric variability of the composite meso-structure that controls the variability of the effective composite mechanical properties. A two-phase viscoelastic finite element model was utilized to treat the behavior of molten prepreg platelets, which were assumed to be incompressible and inextensible in the fiber direction during processing. A two-step simulation procedure with discretely modeled platelets and platelet-to-platelet contacts is proposed. The analysis begins with a drop simulation of the uncontrolled material deposition followed by a 3D finite element analysis of the preform compaction and bulk material flow. The proposed model provides a method for generating a composite meso-structure ab initio starting from a statistical sampling of the initial conditions for platelet deposition. The predicted platelet thickness distribution, intra-platelet fiber waviness and out-of-plane platelet undulations were found to be in good agreement with experimental measurements. An estimate of the volume of resin pockets was inferred from the compaction simulation.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Analysis of contact conditions and microstructure evolution in shear assisted processing and extrusion using smoothed particle hydrodynamics method

Shear assisted processing and extrusion (ShAPE) is a solid-phase processing technique that adds an additional shear force as compared with a conventional extrusion approach. Recently, ShAPE has demonstrated the capability of extruding high-performance aluminum alloy 7075 (AA7075) tubes at speeds up to 12.2 m/min without surface tearing. However, the relationship among the ShAPE processing parameters, thermomechanical conditions, contact conditions, heat generation, and microstructure evolution remains primarily empirical because an insightful understanding of the associated physics is still lacking. To help elucidate these relationships, this work proposes a thermomechanical meshfree model for the first time for ShAPE processing of AA7075 using the smoothed particle hydrodynamics (SPH) method. The meshfree model is first validated thoroughly by experimental data in terms of material flow, die face temperature, and extrusion force with various processing parameters. The validated model is then used to analyze the steady-state contact conditions and heat generation rates during ShAPE processing. Distributions of the average grain size of AA7075 being extruded are calculated using the SPH model output. The meshfree model results reveal that extrusions conducted at lower temperatures and higher strain rates yield more refined grains and possibly higher material strength, which is also consistent with the experimental observations.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Modeling the Size Distribution and Chemical Composition of Secondary Organic Aerosols during the Reactive Uptake of Isoprene-Derived Epoxydiols under Low-Humidity Condition

Reactive uptake of isoprene epoxydiols (IEPOX), which are isoprene oxidation products, onto acidic sulfate aerosols is recognized to be an important mechanism for the formation of isoprene-derived secondary organic aerosol (SOA). While a mechanistic understanding of IEPOX-SOA formation exists, several processes affecting their formation remain uncertain. Evaluating mechanistic IEPOX-SOA models with controlled laboratory experiments under longer atmospherically relevant timescales is critical. Here, we implement our latest understanding of IEPOX-SOA formation within a box model to simulate the measured reactive uptake of IEPOX on polydisperse ammonium bisulfate seed aerosols within an environmental Teflon chamber. The model is evaluated with single-particle measurements of size distribution, volume, density, and composition of aerosols due to IEPOX-SOA formation at timescales of hours. Additionally, we find that the model can simulate the growth of particles due to IEPOX multiphase chemistry, as reflected in increases of the mean particle size and volume concentrations, and a shift of number size distribution to larger sizes. The model also predicts the observed evolution of particle number mean diameter and total volume concentrations at the end of the experiment. We show that in addition to self-limiting effects of IEPOX-SOA coatings, accounting for the molar balance between inorganic and organic sulfate and mass accommodation coefficient of IEPOX are important parameters governing modeling of IEPOX-SOA formation. Thus, models which do not account for the molar sulfate balance and/or diffusion limitations within IEPOX-SOA coatings are likely to predict too high IEPOX-SOA formation.

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

Catalytically Active Site Mapping Realized through Energy Transfer Modeling

Abstract The demands of a sustainable chemical industry are a driving force for the development of heterogeneous catalytic platforms exhibiting facile catalyst recovery, recycling, and resilience to diverse reaction conditions. Homogeneous‐to‐heterogeneous catalyst transitions can be realized through the integration of efficient homogeneous catalysts within porous matrices. Herein, we offer a versatile approach to understanding how guest distribution and evolution impact the catalytic performance of heterogeneous host–guest catalytic platforms by implementing the resonance energy transfer (RET) concept using fluorescent model systems mimicking the steric constraints of targeted catalysts. Using the RET‐based methodology, we mapped condition‐dependent guest (re)distribution within a porous support on the example of modular matrices such as metal–organic frameworks (MOFs). Furthermore, we correlate RET results performed on the model systems with the catalytic performance of two MOF‐encapsulated catalysts used to promote CO 2 hydrogenation and ring‐closing metathesis. Guests are incorporated using aperture‐opening encapsulation, and catalyst redistribution is not observed under practical reaction conditions, showcasing a pathway to advance catalyst recyclability in the case of host–guest platforms. These studies represent the first generalizable approach for mapping the guest distribution in heterogeneous host–guest catalytic systems, providing a foundation for predicting and tailoring the performance of catalysts integrated into various porous supports.

Thompson, William J.↗
DOE OSTI2024

Catalytically Active Site Mapping Realized through Energy Transfer Modeling

Abstract The demands of a sustainable chemical industry are a driving force for the development of heterogeneous catalytic platforms exhibiting facile catalyst recovery, recycling, and resilience to diverse reaction conditions. Homogeneous‐to‐heterogeneous catalyst transitions can be realized through the integration of efficient homogeneous catalysts within porous matrices. Herein, we offer a versatile approach to understanding how guest distribution and evolution impact the catalytic performance of heterogeneous host–guest catalytic platforms by implementing the resonance energy transfer (RET) concept using fluorescent model systems mimicking the steric constraints of targeted catalysts. Using the RET‐based methodology, we mapped condition‐dependent guest (re)distribution within a porous support on the example of modular matrices such as metal–organic frameworks (MOFs). Furthermore, we correlate RET results performed on the model systems with the catalytic performance of two MOF‐encapsulated catalysts used to promote CO 2 hydrogenation and ring‐closing metathesis. Guests are incorporated using aperture‐opening encapsulation, and catalyst redistribution is not observed under practical reaction conditions, showcasing a pathway to advance catalyst recyclability in the case of host–guest platforms. These studies represent the first generalizable approach for mapping the guest distribution in heterogeneous host–guest catalytic systems, providing a foundation for predicting and tailoring the performance of catalysts integrated into various porous supports.

Thompson, William J.↗
DOE OSTI2020

Impact of boundary layer simulation on predicting radioactive pollutant dispersion: a case study for HANARO Research Reactor using the WRF-MMIF-CALPUFF modeling system

Wind plays an important role in cases of unexpected radioactive pollutant dispersion, deciding distribution and concentration of the leaked substance. The accurate prediction of wind has been challenging in numerical weather prediction models, especially near the surface because of the complex interaction between turbulent flow and topographic effect. As such, in this study, we investigated the characteristics of atmospheric dispersion of radioactive material (i.e. 137 Cs) according to the simulated boundary layer around the HANARO research nuclear reactor in Korea using the Weather Research and Forecasting (WRF)-Mesoscale Model Interface (MMIF)-California Puff (CALPUFF) model system. We examined the impacts of orographic drag on wind field, stability calculation methods, and planetary boundary layer parameterizations on the dispersion of radioactive material under a radioactive leaking scenario. We found that inclusion of the orographic drag effect in the WRF model improved the wind prediction most significantly over the complex terrain area, leading the model system to estimate the radioactive concentration near the reactor more conservatively. We also emphasized the importance of the stability calculation method and employing the skillful boundary layer parameterization to ensure more accurate low atmospheric conditions, in order to simulate more feasible spatial distribution of the radioactive dispersion in leaking scenarios.

54 ENVIRONMENTAL SCIENCES↗