Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “conditional probability”

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 505 records · Page 28
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Development and Integration of a Stochastic Clad Damage Propagation Model into PRONGHORN-SC Subchannel Analysis Code

The failure of fuel pins in nuclear reactors is intrinsically stochastic. Typically, a combination of variation in manufacturing that affects the material characteristics and the fuel assembly dimensions, variation in operating conditions, such as local power, coolant flow rate, and irradiation induced changes in material properties lead to a large uncertainty in failure margin of the fuel pins. Failure, therefore, may occur in exceptional pins with adverse combinations of these variations. Upon a metal fuel pin (U-Pu-Zr/HT9) failure, depressurization of the fuel pin takes place by release of fission gas, liquid sodium bond, and potentially solid fuel particles or molten/eutectic fuel droplets through the hole in cladding. The effect of a fission gas jet on neighbor fuel pins and possible propagation of a clad damage during normal operation was studied experimentally in 1970s and it was found that the post-failure fission gas jet insulates the jet impingement area of the target fuel pin surface and could increase the target pin’s surface temperature by as much as 100 – 200 K during the failed pin depressurization. It was concluded that the effect should not lead to fuel pin failure propagation during normal operation. In accident scenarios of sodium and lead fast reactors such as Unprotected Loss-Of-Flow (ULOF) or Unprotected Transient Over Power (UTOP), the fuel pins can be subjected to higher clad temperatures and fuel pin pressures or fuel clad mechanical/chemical interaction where thermal creep margin becomes significantly lower compared to the normal operation conditions. Therefore, possible stochastic failure and the post-failure fission gas/fuel jet impingement could be critical in order to predict fuel pin failure propagation. Pin depressurization due to fission gas release may degrade the heat transfer by formation of a gas blanket on a neighboring pin surface, which is a local phenomenon, and by causing coolant flow deceleration and starvation, which could affect a surrounding region as well. Furthermore, the potential presence of solid fuel particles or molten fuel at the time of clad failure could boost post-failure jet induced degradation even further. The present study models the U-Pu-Zr/HT9 metal fuel pin failure and stochastic clad damage propagation by biased sampling based on a Cumulative Damage Fraction (CDF) type clad failure criterion and the normal distribution of fuel failure probability density as a function of logarithm of Cumulative Damage Fraction. In addition, the effect of post-failure fission gas jet on heat transfer degradation is modeled for the target pins. This model is called stochastic Clad Damage Propagation (CDAP). The CDAP model is now fully integrated into developmental version of PRONGHORN-SC subchannel analysis code, allowing for modeling local failures and its propagation potential. Section 2 describes the components of the CDAP models. Section 3 describes the model implementation to PRONGHORN-SC and input specifications. Section 4 describes the CDAP model validation coupled to PRONGHORN-SC.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

FY24 Task 5: Leachate Disposition

Directly feeding sludge solids to the high-level waste (HLW) Waste Treatment Plant represents an alternative flowsheet seeking to initiate sludge processing as soon as possible. Key processing functions previously captured during baseline pretreatment operations include leaching and washing prior to solids concentration. These operations should be considered in the potential direct feed flowsheets to maximize waste feed loading, minimize HLW volume, and mitigate corrosion challenges associated with vitrification of high phosphate and fluoride concentrations. Additionally, single-shell tank (SST) retrievals and waste transfers to double-shell tanks (DSTs) in a direct feed flowsheet would likely also benefit from some level of leaching, washing, and solids concentration in order to reduce DST space and mission duration. These operations could occur in a new facility or potentially in available DSTs. If washing and leaching are utilized, an effective disposition pathway for the wash water and leachate solutions are needed. Three target species that benefit significantly from leaching and washing are phosphate, fluoride and aluminum. Phosphate (PO 4 3- ) and fluoride (F - ) can contribute substantially to the amount of carrier fluid needed for dissolution, and the resulting volume of liquid generated. Disposition of this retrieval solution should be evaluated in order to prevent crystallization of these anions throughout system processing. Since there is a high probability that any retrieval solutions will be at or near their PO 4 3- and F - solubility limits, evaporation or blending with a high Na supernate (>3.5 M) is not recommended for the wash water streams without a method to remove precipitants prior to solution disposal. Additionally, aluminum present in the southeast quadrant of the Hanford site represents roughly 60% of the waste solids in the initial processing tanks. These aluminum solids are in the form of gibbsite (Al(OH) 3 ) and can pose significant challenges for processing due to the fast-settling times and high solids loading associated with these materials. Easily remediated by caustic addition to the solids, these wash solutions could be processed through crystalline silicotitanate (CST) ion exchange columns to prepare the supernate solutions for disposition. The current target for feed conditions to the Low Activity Waste (LAW) melter are waste streams that contain nominally 5-6 M Na. Fractions within the tanks contain upwards of 0.2 M phosphate and fluoride in solution at 3.5 M Na. Concentrating these solutions above 5 M Na would result in an exceedance of the solubility limits, and potential for uncontrolled precipitation of the phosphate and fluoride crystal material. The resulting crystalline salt material is typically sodium fluoride phosphate, also referred to as natrophosphate (Na 7 FPO 4 ·19H 2 O). Salt phases are of importance in tank waste due to their chemical reactivity, which can result in precipitation, dissolution, or transformation, impacting any downstream processes (Bolling et al. 2020, Russell, Snow, and Peterson 2010). Salt generation and precipitation could pose challenges by causing system plugging and melter corrosion if left in the supernate stream, or limit sodium molarity of the supernate that would be accepted without incident in waste operations. To understand the impact of this salt generation, the crystallization of natrophosphate in multiple simulant feed matrices was studied to understand the implications of various tank waste supernate chemistries. Three matrices were examined: high PO 4 3- /low F - , low PO 4 3- /high F - , and an average matrix. Subsequent testing was performed with the average matrix with the inclusion of CsNO 3 , and a final run with the average matrix including CsNO 3 and a 137 Cs spike for tracer purposes.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2024

Root Cause Correlation Analysis of Software Failures via Orthogonal Defect Classification and Natural Language Processing

Systems theoretic process analysis (STPA) is becoming an increasingly popular technique to assess how complex digital software systems can fail. Rather than defining failures by their observable failure events, which may be sparse especially for safety rated nuclear digital instrumentation and control systems (DI&C), failures are defined as postulated unsafe actions under specific contextual conditions. This permits a top-down analysis of system hazards and identifies whether imposed constraints and requirements can sufficiently address undesirable hazards. However, STPA is a qualitative approach at identifying inadequacies in the development process and cannot currently be used to quantify unsafe action likelihoods for probabilistic risk assessment. Therefore, in this work, we examine the root causes of software failure and explore whether a consistent correlation can be linked to specific unsafe action classes. We implement Lbl2Vec, an unsupervised document classification and retrieval algorithm, on a database of 4,096 software defect reports acquired from various open-source software systems. By analyzing sentence structure, embedded labels, and word vectors, we show that certain defect types positively correlate to specific unsafe action classes over others. The correlations developed can be used to estimate the failure probability of safety intended DI&C systems which provides a licensing basis for nuclear plant modernization efforts.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

Root Cause Correlation Analysis of Software Failures via Orthogonal Defect Classification and Natural Language Processing

Systems theoretic process analysis (STPA) is becoming an increasingly popular technique to assess how complex digital software systems can fail. Rather than defining failures by their observable failure events, which may be sparse especially for safety rated nuclear digital instrumentation and control systems (DI&C), failures are defined as postulated unsafe actions under specific contextual conditions. This permits a top-down analysis of system hazards and identifies whether imposed constraints and requirements can sufficiently address undesirable hazards. However, STPA is a qualitative approach at identifying inadequacies in the development process and cannot currently be used to quantify unsafe action likelihoods for probabilistic risk assessment. Therefore, in this work, we examine the root causes of software failure and explore whether a consistent correlation can be linked to specific unsafe action classes. We implement Lbl2Vec, an unsupervised document classification and retrieval algorithm, on a database of 4,096 software defect reports acquired from various open-source software systems. By analyzing sentence structure, embedded labels, and word vectors, we show that certain defect types positively correlate to specific unsafe action classes over others. The correlations developed can be used to estimate the failure probability of safety intended DI&C systems which provides a licensing basis for nuclear plant modernization efforts.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2020

Experimental study of vorticity-strain interactions in turbulent premixed counterflow flames

Here, the effects of heat release on interactions between vorticity (ω) and strain rate (s) in turbulent premixed CH 4 /O 2 /N 2 counterflow flames are investigated using simultaneous OH laser-induced fluorescence (LIF) and tomographic particle image velocimetry (TPIV) measurements. A comparison between the flames and a corresponding turbulent non-reacting variable density N 2 -vs-products counterflow reveals the impact of heat release on vorticity-strain rate alignment statistics. Vorticity and strain rate statistics in the flames and non-reacting flow are conditioned on distance from the local flame front and gas mixing layer interface (GMLI) contours, respectively. The magnitude, alignment, and spatial distribution of the vorticity and principal strain rates (s 1 , s 2 , s 3 ) are rather different when heat release is present. Density variations without heat release enhance the ω-s 2 alignment while significantly reducing the ω-s 3 alignment and modestly reducing the ω-s 1 alignment. In contrast, heat release at the flame front further reduces the ω-s 1 alignment but increases the ω-s 3 alignment and suppresses the preferential ω-s 2 alignment. Furthermore, increasing turbulence diminishes the effect of heat release on this preferential alignment. In regions with the largest vorticities, both the reacting and non-reacting counterflows show an increase in the probability of ω-s 2 alignment. All counterflow cases have a net positive vortex-stretching contribution to the enstrophy production with a peak production rate at the flame front or GMLI, but the peak values depend on the density variation, heat release, and turbulence level. Elucidation of the complex interplay between these factors contributes to the understanding of the dynamics of turbulence-flame interactions.

42 ENGINEERING↗
DOE OSTI2020

Probabilistic evolution of stochastic dynamical systems: A meso-scale perspective

Stochastic dynamical systems arise naturally across nearly all areas of science and engineering. Typically, a dynamical system model is based on some prior knowledge about the underlying dynamics of interest in which probabilistic features are used to quantify and propagate uncertainties associated with the initial conditions, external excitations, etc. From a probabilistic modeling standing point, two broad classes of methods exist, i.e. macro-scale methods and micro-scale methods. Classically, macro-scale methods such as statistical moments-based strategies are usually too coarse to capture the multi-mode shape or tails of a non-Gaussian distribution. Micro-scale methods such as random samples-based approaches, on the other hand, become computationally very challenging in dealing with high-dimensional stochastic systems. In view of these potential limitations, a meso-scale scheme is proposed here that utilizes a meso-scale statistical structure to describe the dynamical evolution from a probabilistic perspective. The significance of this statistical structure is twofold. First, it can be tailored to any arbitrary random space. Second, it not only maintains the probability evolution around sample trajectories but also requires fewer meso-scale components than the micro-scale samples. To demonstrate the efficacy of the proposed meso-scale scheme, a set of examples of increasing complexity are provided. Connections to the benchmark stochastic models as conservative and Markov models along with practical implementation guidelines are presented.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

MgADP Promotes Myosin Head Movement toward Actin at Low [Ca 2+ ] to Increase Force Production and Ca 2+ -Sensitivity of Contraction in Permeabilized Porcine Myocardial Strips

Myosin cross-bridges dissociate from actin following Mg 2+ -adenosine triphosphate (MgATP) binding. Myosin hydrolyses MgATP into inorganic phosphate (P i ) and Mg 2+ -adenosine diphosphate (ADP), and release of these hydrolysis products drives chemo-mechanical energy transitions within the cross-bridge cycle to power muscle contraction. Some forms of heart disease are associated with metabolic or enzymatic dysregulation of the MgATP-MgADP nucleotide pool, resulting in elevated cytosolic [MgADP] and impaired muscle relaxation. We investigated the mechanical and structural effects of increasing [MgADP] in permeabilized myocardial strips from porcine left ventricle samples. Sarcomere length was set to 2.0 µm at 28 °C, and all solutions contained 3% dextran T-500 to compress myofilament lattice spacing to near-physiological values. Under relaxing low [Ca 2+ ] conditions (pCa 8.0, where pCa = -log 10 [Ca 2+ ]), tension increased as [MgADP] increased from 0-5 mM. Complementary small-angle X-ray diffraction measurements show that the equatorial intensity ratio, I 1,1 /I 1,0 , also increased as [MgADP] increased from 0 to 5 mM, indicating myosin head movement away from the thick-filament backbone towards the thin-filament. Ca 2+ -activated force-pCa measurements show that Ca 2+ -sensitivity of contraction increased with 5 mM MgADP, compared to 0 mM MgADP. These data show that MgADP augments tension at low [Ca 2+ ] and Ca 2+ -sensitivity of contraction, suggesting that MgADP destabilizes the quasi-helically ordered myosin OFF state, thereby shifting the cross-bridge population towards the disordered myosin ON state. Together, these results indicate that MgADP enhances the probability of cross-bridge binding to actin due to enhancement of both thick and thin filament-based activation mechanisms.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2024

Optimizing Materials to Boost the Valorization of CO 2 : Tuning Cobalt–Cobalt Interactions on In 2 O 3 -Based Photothermal Catalysts

The valorization of CO 2 is an important challenge within the current panorama, since this molecule is probably the main contributor to climate change. In this study, the synthesis of materials based on a nanostructured batonnet-type indium oxide is carried out. In them, different amounts of Co are introduced, varying between 2 and 8% mol. It is verified that the most active sample in the transformation of carbon dioxide to carbon monoxide contains 6 mol %. of Co. This sample’s activity under dual excitation exceeds the thermal counterpart by more than 30%. After carrying out a complete physical and chemical characterization with the help of X-ray absorption spectroscopy and other techniques, it is shown that catalysts with amounts of cobalt equal to or below 4 mol % contain isolated single-atom species, while those with higher amounts of metal display a Co–Co interaction which triggers the evolution of the samples under reaction conditions. The optimum control of this Co–Co interaction and the nature of the final cobalt-containing species determine dual photothermal catalytic properties. In conclusion, this work establishes a structure–activity relationship to interpret the catalytic behavior of highly dispersed subnanometric cobalt species, and thus an avenue to optimize the photothermal valorization of carbon dioxide.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Multivariate Frequency Analysis of Annual Maxima Suspended Sediment Concentrations and Floods in the Jinsha River, China

Streamflow plays a critical role in sediment transport. Traditional frequency analysis of annual maximal (AM) suspended sediment concentrations (SSCs) is usually based on univariate distribution of SSCs and may provide a limited assessment of SSC state, leading to an underestimate or overestimate of the corresponding probability. To address this issue, we propose a multivariate frequency analysis framework is proposed for extreme sediment events using copula functions that account for the joint behavior of SSC and flood flow (i.e., peak discharge and volume). The uncertainties associated with copula modeling are also discussed. The Jinsha River in China is used as a case study. The results show that the copula-based multivariate method provides a more comprehensive assessment of extreme SSCs and their recurrence intervals than univariate frequency analysis. The uncertainty of copula modeling decreases with longer observation time and increases with larger return periods. The uncertainty ranges of the most likely design quantiles vary with different copulas, and the uncertainty range for the best-fit copula is not the smallest. This study is useful for the risk assessment of extreme SSC events under flood conditions.

42 ENGINEERING↗
DOE OSTI2024

Pseudorapidity dependence of anisotropic flow and its decorrelations using long-range multiparticle correlations in Pb–Pb and Xe–Xe collisions

The pseudorapidity dependence of elliptic (v 2 ), triangular (v 3 ), and quadrangular (v 4 ) flow coefficients of charged particles measured in Pb–Pb collisions at a centre-of-mass energy per nucleon pair of $\sqrt{S_{NN}}$ = 5.02 TeV and in Xe–Xe collisions at $\sqrt{S_{NN}}$ = 5.44 TeV with ALICE at the LHC are presented. The measurements are performed in the pseudorapidity range − 3.5 < η < 5 for various centrality intervals using two- and multi-particle cumulants with the subevent method. The flow probability density function (p.d.f.) is studied with the ratio of flow coefficient v 2 calculated with four- and two-particle cumulant, and suggests that the variance of flow p.d.f. is independent of pseudorapidity. The decorrelation of the flow vector in the longitudinal direction is probed using two-particle correlations. The results measured with respect to different reference regions in pseudorapidity exhibit differences, argued to be a result of saturating decorrelation effect above a certain pseudorapidity separation, in contrast to previous publications which assign this observation to non-flow effects. The results are compared to 3 + 1 dimensional hydrodynamic and the AMPT transport model calculations. Neither of the models is able to simultaneously describe the pseudorapidity dependence of measurements of anisotropic flow and its fluctuations. The results presented in this work highlight shortcomings in our current understanding of initial conditions and subsequent system expansion in the longitudinal direction. Therefore, they provide input for its improvement.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Representing surface heterogeneity in land–atmosphere coupling in E3SMv1 single-column model over ARM SGP during summertime

Abstract. The Earth's land surface features spatial and temporal heterogeneity over a wide range of scales below those resolved by current Earth system models (ESMs). State-of-the-art land and atmosphere models employ parameterizations to represent their subgrid heterogeneity, but the land–atmosphere coupling in ESMs typically operates on the grid scale. Communicating the information on the land surface heterogeneity with the overlying atmospheric boundary layer (ABL) remains a challenge in modeling land–atmosphere interactions. In order to account for the subgrid-scale heterogeneity in land–atmosphere coupling, we implement a new coupling scheme in the Energy Exascale Earth system model version 1 (E3SMv1) that uses adjusted surface variances and covariance of potential temperature and specific water content as the lower boundary condition for the atmosphere model. The new lower boundary condition accounts for both the variability of individual subgrid land surface patches and the inter-patch variability. The E3SMv1 single-column model (SCM) simulations over the Atmospheric Radiation Measurement (ARM) Southern Great Plain (SGP) site were performed to assess the impacts. We find that the new coupling parameterization increases the magnitude and diurnal cycle of the temperature variance and humidity variance in the lower ABL on non-precipitating days. The impacts are primarily attributed to subgrid inter-patch variability rather than the variability of individual patches. These effects extend vertically from the surface to several levels in the lower ABL on clear days. We also find that accounting for surface heterogeneity increases low cloud cover and liquid water path (LWP). These cloud changes are associated with the change in cloud regime indicated by the skewness of the probability density function (PDF) of the subgrid vertical velocity. In precipitating days, the inter-patch variability reduces significantly so that the impact of accounting for surface heterogeneity vanishes. These results highlight the importance of accounting for subgrid heterogeneity in land–atmosphere coupling in next-generation ESMs.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Equatorial waves triggering extreme rainfall and floods in southwest Sulawesi, Indonesia

On the basis of detailed analysis of a case study and long-term climatology, it is shown that equatorial waves and their interactions serve as precursors for extreme rain and flood events in the central Maritime Continent region of southwest Sulawesi, Indonesia. Meteorological conditions on January 22, 2019, leading to heavy rainfall and devastating flooding in this area were studied. It is shown that a convectively coupled Kelvin wave (CCKW) and a convectively coupled Rossby wave (CCERW) embedded within the larger-scale envelope of the Madden-Julian Oscillation (MJO) enhanced convective phase, contributed to the onset of a mesoscale convective system which developed over the Java Sea. Low-Level convergence from the CCKW forced mesoscale convective organization and orographic ascent of moist air over the slopes of southwest Sulawesi. Climatological analysis shows that 92 % of December-January-February floods and 76% of extreme rain events in this region were immediately preceded by positive low-level westerly wind anomalies. It is estimated that both CCKWs and CCERWs propagating over Sulawesi double the chance of floods and extreme rain event development, which the probability of such hazardous events occurring during their combined activity is eight times greater than on a random day. While the MJO is a key component shaping tropical atmospheric variability, it is shown that its usefulness as a single factor for extreme weather-driven hazard prediction is limited.

Latos, Beata↗
DOE OSTI2025

Constraining the 3He + 3He Gamow energy probed in high energy density plasmas at the National Ignition Facility

Polar-direct-drive implosions at the National Ignition Facility generated large plasma volumes to study the 3He + 3He fusion reaction. The ion temperature, which determines the Gamow peak energy, was constrained by isolating the thermal contribution to the D3He-proton spectral width in a 3He plasma doped with deuterium. X-ray penumbral imaging was used to measure electron temperature, density, and hotspot volume, which was subsequently used to model the spectral broadening from plasma stopping power. Results showed 30% of the D3He-proton spectral width was due to stopping power, with residual flows contributing ≈10%. The 3He temperature was determined as T3He = 12.4 ± 3.2 keV, corresponding to a Gamow energy of 95 ± 14 keV. These experiments achieved the lowest Gamow energy to date for studying 3He + 3He fusion in high energy density plasma, approaching conditions in the Sun.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

A combined ensemble-volume average homogenization method for lattice structures with defects under dynamic and static loading

In the study of lattices structures, both experiments and numerical simulations are often conducted with small samples. Using combined ensemble and volume averaging, this work introduces a method to extract a macroscopic constitutive response of a lattice material from numerical simulations performed in periodic domains. The domain size needed to obtain statistically accurate results is investigated. Similar to molecular dynamics, the concept of the virial stress is introduced after homogenized equations are derived using the ensemble averaging method. Under static conditions, the virial stress is shown to agree with the volume averaged solid stress. Using the homogenization method, constitutive relations for this stress can be obtained from systems with uniform strains. Application of such obtained constitutive relations to more general cases results in an error proportional to the square of the ratio between the lattice length scale and the macroscopic length scale. Taking advantage of this property, numerical simulations are performed in systems with a uniform gradient of the average velocity. The volume average method is then used to accelerate convergence when studying lattices with defects. To avoid the artificial numerical time scale from the size of a representative volume element divided by the wave speed, a numerical scheme is developed to enforce a spatially uniform velocity gradient within the computational domain while allowing fluctuations of the velocity or displacement to develop naturally. To account for probability distribution of lattice defects, the stress is calculated as the ensemble-volume averaged value. For dynamic systems, energy dissipation properties are also studied.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Direct observation of the superallowed α-decay of 104 Te

The radioactivity of the α particle is among the most compelling evidence for the existence of cluster structures in atomic nuclei. During the decay process, a pre-existing α particle tunnels through the potential barrier formed by the residual nucleus1,2. The degree of preformation of the α particle, a strongly bound system of two protons and two neutrons, is extracted from the data by dividing the α-decay probability by the barrier penetrability for a given particle energy. The preformation probability changes rapidly near nuclear shell closures, which is direct evidence that clustering is connected to nuclear structure3. Enhanced preformation was observed in the lightest α-particle emitters, spherical tellurium and xenon isotopes decaying to magic isotopes of tin. Here we show the most extreme case of α-particle preformation from the measurement of the decay of tellurium-104 (104Te). With a half-life of , 104Te is the fastest ground-state α-emitting nucleus known so far. The deduced preformation demonstrates that the enhancement is greater for 104Te than for any other nucleus. One nuclear model that can explain our observation postulates that the α particle can exist only in the low-nuclear-matter-density regions on the surface of the nucleus. The uniquely high preformation for 104Te is attributed to its relation to doubly magic tin-100 (100Sn), creating conditions conducive to form an α particle.

Cox, Ian↗
DOE OSTI2025

Voltage Probability Density Function Shaping Control Strategy Considering Grid Operational Uncertainties

It is well-known that power systems operation always affected by various uncertainties which make the bus voltage a random process that can be characterized by its probability density function (PDF) at any time instant. In this context, this paper presents a novel PDF-based voltage control framework for power systems. By modeling voltage as a stochastic process, we formulate a stochastic differential equationthat captures grid uncertainties. The associated Fokker-Planck-Kolmogorov equation is derived to describe the evolution of the voltage PDF, which enables the formulation of a PDF-shaping control strategy. To simplify the PDF control formulation, a B-spline neural network is introduced for real-time estimation and regulation of the voltage distribution. The proposed PDF control law updates voltage references for energy storage systems and synchronous generators using real-time PDF measurements and feedback signals. The proposed method is validated on a modified Kundur’s two-area system. Simulation results demonstrate that the controller can significantly improve the voltage stability under stochastic conditions, highlighting its effectiveness in modern inverter-rich grids.

Gui, Yonghao [ORNL] (ORCID:0000000250435534)↗
DOE OSTI2023

COMPUTATIONAL MODELING OF IGNITION AND PREMIXED FLAME PROPAGATION INITIATED BY A PRE-CHAMBER TURBULENT JET

Addressing the pressing need for reduced carbon emissions, Turbulent Jet Ignition (TJI) emerges as a promising technology for ultra-lean combustion, offering enhanced thermal efficiencies and minimized cyclic variability in spark-ignited engines. To facilitate rapid testing and integration of this technology, a robust computational modeling framework is crucial. This study delves into the predictive capabilities of computational models for main-chamber ignition and premixed flame propagation using a single-cycle TJI rig measured by Biswas et al. (Applied Thermal Engineering, vol 106, 2016). Employing an open-source compressible flow simulation solver with Large Eddy Simulation (LES) for turbulence modeling, the investigation integrates the conventional Laminar Finite Rate Chemistry (LFRC) model alongside the transported Probability Density Method (PDF) for turbulence-chemistry interaction. A fully-consistent Eulerian Monte-Carlo Fields (EMCF) method is utilized to approximate the transported PDF, while Interaction by Exchange with Mean is employed to close micro-mixing terms in stochastic differential equations. A reduced chemical reaction mechanism with 21 species and 84 reactions (DRM-19) is used for solving chemical kinetics, and a double Gaussian energy deposition model is used to approximate the spark ignition in the pre-chamber. An unstructured O-grid mesh with 0.3 million cells in the prechamber and 1 million cells in the main chamber is employed. Results are divided into two phases: pre-chamber initialization and full TJI simulations. Validation of the predicted pre-chamber flame propagation and the lean ignition in the main-chamber is carried out by using available experimental data. Under quiescent conditions, both the LFRC and transported PDF methods largely underestimate the flame speed and subsequent pressure growth in the pre-chamber. A linear momentum forcing technique is applied to investigate the impact of initial turbulence in the pre-chamber, demonstrating a notable influence on flame propagation. Fine-tuning of the forcing coefficient reproduces the sudden pressure growth observed in the experiment. The experimentally validated pre-chamber simulation serves as the initial condition for the full TJI simulations. It is found that the LFRC model fails to predict lean-ignition in the main-chamber, resulting in a misfiring event. Incorporation of turbulence-chemistry interaction using the transported PDF method substantially improves the prediction of the ignition event in the main-chamber, achieving fair qualitative agreement and quantitative validation of combustion parameters within ±10% of the reported experimental data. The rich simulation results consisting of a full set of statistical description of the thermo-chemical states enable us to gain deep insights into the ignition mechanisms in the main chamber, which is limited when done experimentally. A novel dual ignition phenomenon is revealed in the TJI rig for the first time. Initially, a primary ignition kernel is formed at a downstream location which eventually detaches from the main jet. As the jet momentum decreases, a secondary ignition event follows, this time at a more upstream location which eventually combines with the primary ignition kernel to form a single connected flame front. Investigation of these ignition sequences in chemical composition space reveal distinct differences between the two. The primary ignition event in the main-chamber is followed by a large concentration of active radicals from the pre-chamber jet, accelerating the chain-branching steps, characterizing what has been referred to as flame ignition. In contrast, the secondary ignition occurs in the absence of active radicals in the pre-chamber jet, hence characterized as jet ignition. Further analysis of the effect of pre-chamber jet characteristics on lean ignition in the main-chamber is conducted by setting up cases with different initial pressure ratios (por) between the two chambers, a non-dimensional parameter, ranging from 1.2 to 3.2. As the initial pressure ratio increases, jet momentum increases, with dual ignition observed in cases above por= 2.2. Case with por= 3.2 lead to misfiring. The effect of ignition sequence on global combustion characteristics of TJI is analyzed. Dual ignition events lead to non-monotonicity in combustion characteristics such as global reaction progress variable, flame penetration, and global heat release rate. In dual ignition events, although the rate of fuel consumption and global heat release rate is initially lower, the secondary ignition leads to a sudden increase in flame surface area, resulting in a sudden jump and promoting the overall performance of the TJI system.

42 ENGINEERING↗
DOE OSTI2023

Communication-Constrained Robust Control and Learning of Grid-Connected

The electric grid of things (EGoT) promises great potential for innovative grid services by tapping into vast load flexibility. However, the unique characteristics of EGoT, being a part of the cyber-physical electric power system, present both opportunities and challenges, especially concerning supply-demand balancing, stability, and communication constraints. Traditionally, centralized control was employed to ensure balance and stability in power systems. However, with the massive influx of EGoT devices, new strategies are needed to efficiently coordinate and control these distributed devices for optimal grid operations. While some studies have explored efficiency and economic models, there remains a gap in ensuring reliability under everyday operations and resilience during extreme conditions. Addressing this gap, this project develops the technology for an Energy Service Interface (ESI) that includes novel pricing, control, learning, and distributed optimization algorithms, which will enable utilities to recruit EGoT assets for crucial grid services such as load flexibility, voltage regulation, and situation-awareness. The key novelty of the proposed technology is the careful distribution of learning and control functions across utility and EGoT asset owners such that provably efficient and resilient grid operations are attained while respecting communication and information-exchange constraints. Specifically, the project team develops machine-learning enhanced load modeling methods to allow EGoT asset owners to learn their load capability and flexibility, and develops pricing-based and decentralized learning-based control so that asset owners can coordinate to meet system-wide demand-supply balance and reliability goals. For extreme situations involving high-impact, low-probability catastrophic events (termed the “black-sky” events), the team also develops (1) a “Feeder-Operating Center-on-a-Laptop” (FOCAL) software that can assist utility personnel in leveraging EGoT assets to accelerate the service recovery of damaged feeders, and (2) distributed optimization algorithms that can coordinate the operation points of EGoT devices under severe communication constraints. The proposed technology has been extensively tested and evaluated through simulations and on a testbed. In summary, as we transition into a more interconnected and digital power grid era, our project’s findings and developments offer a pivotal step toward guaranteeing both efficiency and resilience in the face of both everyday operations and rare “black-sky” events.

24 POWER TRANSMISSION AND DISTRIBUTION↗