A physics-informed machine learning approach for predicting dynamic behavior of reacting flows with application to hydrogen jet flames
Not Available
Sharma, Pushan (ORCID:0000000196707147)↗
SEARCH · Search NASA
Results for “flow (dynamics)”
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.
Generative AI models for learning flow maps of stochastic dynamical systems in bounded domains
Simulating stochastic differential equations (SDEs) in bounded domains, presents significant computational challenges due to particle exit phenomena, which requires accurate modeling of interior stochastic dynamics and boundary interactions. Despite the success of machine learning-based methods in learning SDEs, existing learning methods are not applicable to SDEs in bounded domains because they cannot accurately capture the particle exit dynamics. We present a unified hybrid data-driven approach that combines a conditional diffusion model with an exit prediction neural network to capture both interior stochastic dynamics and boundary exit phenomena. Our ML model consists of two major components: a neural network that learns exit probabilities using binary cross-entropy loss with rigorous convergence guarantees, and a training-free diffusion model that generates state transitions for non-exiting particles using closed-form score functions. The two components are integrated through a probabilistic sampling algorithm that determines particle exit at each time step and generates appropriate state transitions. Here, the performance of the proposed approach is demonstrated via three test cases: a one-dimensional simplified problem for theoretical verification, a two-dimensional advection-diffusion problem in a bounded domain, and a three-dimensional problem of interest to magnetically confined fusion plasmas.
Bounded domains↗
Numerical investigation of magneto-inertial fusion targets magnetized by dynamic enforcement of helical current flow
Magnetized Liner Inertial Fusion (MagLIF) targets require premagnetization to reduce thermal conduction losses from the laser-heated fuel to the liner material during implosion. MagLIF targets are typically magnetized by external magnetic field coils which are technologically limited to providing <20 T axial magnetic fields in the fusion fuel. We present a novel target design, AutoMag-D, which employs dynamic enforcement of helical current in the liner, resulting in >20 T axial magnetic field in the fuel region prior to implosion without the need of external magnetic field coils. These self-magnetizing liners are made of helically oriented alternating electrically conductive material and electrically insulating material surrounded by a thin, conductive outer radial layer. As the liner is pulsed with a ∼ 20 MA, ∼ 100 ns rise time drive current from Z, the outer layer of conductive material is shocked and intensely Joule heated, causing it to melt, vaporize, and turn to plasma. This allows the magnetic drive field to diffuse radially inward which dynamically enforces current flow in the helical conduction paths in the liner and produces axial magnetic field in the fuel region prior to implosion of the inner liner surface. AutoMag-D liner designs do not require dielectric breakdown of electrically insulating material (as in traditional auto-magnetizing helical liners) and do not require helical return current geometries (as in dynamic screw pinches). We present results from three-dimensional radiation-magnetohydrodynamic simulations of MagLIF implosions employing AutoMag-D liner designs. Simulated AutoMag-D targets demonstrate improved fuel conditions and thermonuclear yield compared to traditional MagLIF.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Changes in soil water content and lateral flow exert large effects on soil thermal dynamics across Alaskan landscapes
Both lateral surface and subsurface water flow affect soil moisture dynamics, yet most land surface models only solve subsurface water movement vertically. Here, we use a 3D ecosystem model that considers both land surface and subsurface hydrologic processes to simulate soil moisture, which is then used to drive a 1-D vertical soil thermal model to simulate the soil moisture effects on soil thermal dynamics in central Alaska. Our coupled model improves soil temperature (ST) estimates by 43.5% in comparison with observational data. Soil moisture has little effect on ST during the wet season (-1.5%) and a substantial influence during the dry season (60%). Spatially, water lateral flow has significant impacts on both soil moisture and ST, causing model estimates for thawed areas in the transition season to increase by ~10% in the study area. Our results highlight the importance of considering dynamical soil moisture, as well as lateral flow effects, on soil thermal dynamics in permafrost regions.
54 ENVIRONMENTAL SCIENCES↗
Multimode turbulent flow measurements using magnetic resonance imaging- and laser-based techniques and computational fluid dynamics simulations
We studied the flow field characteristics of a turbulent flow over a regularized cube array with a perpendicular injection flow through the floor between the second and third cubical elements, representing the complex flow interactions of a 3D jet and the wake flows behind cubical obstacles. Four different experimental measurements were performed: two magnetic resonance imaging-based measurements for three-dimensional three-component velocity (MRV) and concentration (MRC) and two laser-based techniques, particle image velocimetry (PIV) and planar laser-induced fluorescence (PLIF), for two-dimensional two-component velocity and concentration measurement, respectively. The mainstream Reynolds number is Re = 15 000, based on the primary inlet velocity U m and channel height D H , whereas the injector Reynolds number is Re j = 3400, based on the injector velocity U j and the injector's exit width D j . Numerical simulations were performed for the studied flow configuration of turbulent flow over a regularized cube array using Reynolds-averaged Navier–Stokes (RANS) and large-eddy simulation (LES) approaches. Results obtained from experimental measurements—including MRV, MRC, PIV, and PLIF—as well as RANS and LES simulations are discussed and compared along several horizontal and vertical planes of the studied configuration. In addition, 3D turbulent flow structures, such as leading-edge vortex, horseshoe vortex, and jet shear-layer vortex, and the isosurfaces of scalar concentration successfully revealed by MRV and MRC techniques were found to be in very good agreement with those 3D features extracted from RANS and LES simulations. In conclusion, the high-resolution experimental and numerical database obtained from this study could be useful for validation and verification of numerical codes.
Computational fluid dynamics↗
Using principal component analysis to distinguish different dynamic phases in superconducting vortex matter
Vortices in type-II superconductors driven over random disorder are known to exhibit a remarkable variety of distinct nonequilibrium dynamical phases that arise owing to the competition between vortex-vortex interactions, the quenched disorder, and the drive. These include pinned states, elastic flows, plastic or disordered flows, and dynamically reordered moving crystal or moving smectic states. The plastic flow phases can be particularly difficult to characterize since the flows are strongly disordered. Here, we perform principal component analysis (PCA) on the positions and velocities of vortex matter moving over random disorder for different disorder strengths and drives. We find that PCA can distinguish the known dynamic phases as well as or better than previous measures based on transport signatures or topological defect densities. In addition, PCA recognizes distinct plastic flow regimes, a slowly changing channel flow and a moving amorphous fluid flow, that do not produce distinct signatures in the standard measurements. In conclusion, our results suggest that this position and velocity-based PCA approach could be used to characterize dynamic phases in a broader class of systems that exhibit depinning and nonequilibrium phase transitions.
75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
Development of a Method for Shape Optimization for a Gas Turbine Fuel Injector Design Using Metal-Additive Manufacturing
Adjoint shape optimization has enabled physics-based optimal designs for aerodynamic surfaces. Additive manufacturing (AM) makes it possible to manufacture complex shapes. However, there has been a gap between optimal and manufacturable surfaces due to the inherent limitations of commercial computational fluid dynamics (CFD) codes to implement geometric constraints during adjoint computation. In such cases, the design sensitivities are exported and used to perform constrained shape modifications using parametric information stored in computer aided design (CAD) files to satisfy manufacturability constraints. However, modifying the design using adjoint methods in CFD solvers and performing constrained shape modification in CAD can lead to inconsistencies due to different shape parameterization schemes. This paper describes a method to enable the simultaneous optimization of the fluid domain and impose AM manufacturability constraints, resolving one of the key issues of geometry definition for isogeometric analysis. Similar to a grid convergence study, the proposed method verifies the consistencies between shape parameterization techniques present within commercial CAD and CFD software during mesh movement as a part of the adjoint shape optimization routine. By identifying the appropriate parameters essential to a shape optimization study, the error metric between the different parameterization techniques converges to demonstrate sufficient consistencies for justifiable exchange of data between CAD and CFD. For the identified shape optimization parameters, the error metric to measure the deviation between the two parameterization schemes lies within the AM laser-powder bed fusion (L-PBF) process tolerance. Additionally, comparison for subsequent objective function calculations between iterations of the optimization loop showed acceptable differences within 1% variation between the modified geometries obtained using the two parameterization schemes. This method provides justification for the use of multiphysics guided adjoint design sensitivities computed in CFD software to perform shape modifications in CAD to incorporate AM manufacturability constraints during the shape optimization loop such that optimal designs are also additively manufacturable.
33 ADVANCED PROPULSION SYSTEMS↗
Onset of kinetic effects on Rayleigh–Taylor instability: Advective–diffusive asymmetry
In nature and engineering applications, the Rayleigh–Taylor instability (RTI) occurs over a wide range of Atwood, Reynolds, Mach, and Knudsen numbers. At low Atwood, Mach, and Knudsen numbers, the classic advective instability causes quasi-symmetric bubble and spike growth on the two sides of the interface. However, recent findings suggest that at high degrees of rarefaction, advective effects are suppressed and molecular diffusion leads to planar growth of the density fronts on either side of the interface. This study aims to investigate the flow physics of the transition from advective to diffusive behavior, focusing on the onset of the kinetic effects. Using the gas kinetic methodology, RTI is simulated over a range of Knudsen and Mach numbers in the transition regime. The simulation results reveal the various stages of transformation from advective instability to diffusive transport. For the first time, the study demonstrates the existence of a Knudsen–Mach parameter regime where the bubble side exhibits advective instability, while the other side shows a planar density front due to molecular diffusion, rather than the canonical advective spike shape. The dominance of different mechanisms on the two sides of the interface leads to the advective–diffusive asymmetry. In conclusion, the findings of this study can lead to a more comprehensive understanding of RTI over a wide range of Mach and Knudsen numbers.
42 ENGINEERING↗
Active learning of a crystal plasticity flow rule from discrete dislocation dynamics simulations
Continuum-scale material deformation models, such as crystal plasticity (CP), can significantly enhance their predictive accuracy by incorporating input from lower-scale (i.e. mesoscale) models. The procedure to generate and extract the relevant information is however typically complex and ad hoc, involving decision and intervention by domain experts, leading to long development times. In this study, we develop a principled approach for calibration of continuum-scale models using lower scale information by representing a CP flow rule as a Gaussian process model. This representation allows for efficient parameter space exploration, guided by the uncertainty embedded in the model through a process known as Bayesian optimization (BO). We demonstrate a semi-autonomous BO loop which instantiates discrete dislocation dynamics simulations whose initial conditions are automatically chosen to optimize the uncertainty of a model CP flow rule. Our self-guided computational pipeline efficiently generated a dataset and corresponding model whose error, uncertainty, and physical feature sensitivities were validated with comparison to an independent dataset four times larger, demonstrating a valuable and efficient active learning implementation readily transferable to similar material systems.
36 MATERIALS SCIENCE↗
Advanced Flexible Transformers
Advanced grid solutions are comprised of advanced transmission technologies and grid enhancing technologies. In this webinar, experts will provide participants with insights into eight advanced grid solutions. The advanced transmission technologies that will be discussed include point-to-point high voltage direct current and advanced conductoring and the grid enhancing technologies that experts will explore include topology optimization, advanced power flow control, dynamic line rating, energy storage, virtual power plants, and advanced flexible transformers
flexible power transformer, power flow controller,↗
Chemical Reactor Network Modeling of Ammonia Rich-Quench-Lean Combustion Using a Partially Stirred Reactor Approach
Ammonia is a promising alternative to hydrogen with high energy density and favorable storage and transport characteristics. However, low flammability and a propensity for high nitrogen oxide (NO x ) emissions make direct utilization challenging. Recently, two-stage rich-quench-lean (RQL) combustion strategies have shown promise in achieving low NO x emissions with ammonia. In this approach, the rich stage serves to oxidize a portion of the fuel while thermally decomposing as much of the remaining ammonia as possible, generating hydrogen. In the second (lean) stage, air is rapidly introduced, burning out the hydrogen and residual ammonia. Two-stage RQL combustion of ammonia has been investigated in the open literature both experimentally and numerically. In general, idealized chemical reactor network (CRN) models predict NO x concentrations below those of 2D/3D computational fluid dynamics models and experiments. The primary drivers of these discrepancies may be largely attributed to finite rate mixing nonadiabatic operation. The typical CRN model is comprised of a perfectly-stirred-reactor (PSR), followed by a plug-flow-reactor (PFR), meant to represent the flame, and postflame zones, respectively. In the two-stage RQL approach two PSR-PFR networks are arranged sequentially, corresponding to the rich and lean stages, with secondary air injection in between. In the authors' past work, this arrangement has demonstrated the significant sensitivity of exit NO x to the rich stage equivalence ratio, while the amount of secondary air injection was shown to be less critical. In this paper, the CRN model is extended to (1) include the impacts of heat loss and (2) utilize a partially-stirred-reactor (PaSR) approach to study the impacts of mixing on emissions performance. Varying amounts of heat loss are applied to the rich relaxation zone to understand emissions performance and changes to optimization of equivalence ratio and residence time. Premixed and nonpremixed configurations are considered in the rich stage PaSR, with varying degrees of mixing intensity to study the interaction between mixing, transport, and kinetic timescales. Critically, the impact of mixing between hot products and secondary air injection is studied to understand practical injector needs. Results show unburnt ammonia leaving the rich stage as a primary contributor to NO x emissions – driven both by increased heat loss and reduced mixing rates. Furthermore, heat losses have been shown to create conditions that are conducive to increased N 2 O formation in the lean stage. In conclusion, the results of this study will be considered in the context of developing optimized two-stage RQL combustors for ammonia.
Combustion↗
CFD simulations of Molten Salt Fast Reactor core cavity flows
Computational Fluid Dynamics (CFD) has become increasingly important in the research and development of advanced nuclear reactors. Here, in the current study, extensive CFD simulations were conducted for the coolant flow in Molten Salt Fast Reactor (MSFR) core models using the state-of-the-art spectral element flow solver Nek5000 and multiscale coarse-mesh thermal-hydraulic software Pronghorn. The underlying motivation is to seek an in-depth understanding of how the internal velocity distribution can be influenced by the MSFR core cavity shape, the Reynolds number, turbulence modeling options and the inlet boundary conditions. The CFD techniques involved in this investigation range from coarse-mesh CFD, RANS modeling, to the high-fidelity LES calculations. Specifically, a series of RANS simulations were performed for the 2-D axisymmetric core model and 3-D wedge domains to study the flow distribution inside the MSFR core. It is observed that a proper representation of the MSFR inlet channel duct is important for the prediction of internal flow distribution. It is also showcased here how researchers can leverage the Nek5000 CFD results to calibrate more efficient coarse-mesh CFD tools, like Pronghorn, for the actual MSFR design needs. Moreover, this paper highlights a 3-D LES model for an entire MSFR core using the spectral element method and demonstrates the feasibility of this modeling approach. The readiness and potential limitations of the RANS approach are examined with respect to the high-fidelity LES simulations. The present investigation lays a solid foundation as we are leveraging the high-fidelity CFD capabilities to inform MSFR design efforts.
97 MATHEMATICS AND COMPUTING↗
Fully quantum algorithm for mesoscale fluid simulations with application to partial differential equations
Fluid flow simulations marshal our most powerful computational resources. In many cases, even this is not enough. Quantum computers provide an opportunity to speed up traditional algorithms for flow simulations. We show that lattice-based mesoscale numerical methods can be executed as efficient quantum algorithms due to their statistical features. This approach revises a quantum algorithm for lattice gas automata to reduce classical computations and state preparation at every time step. For this, the algorithm approximates the qubit relative phases and subtracts them at the end of each time step. Phases are evaluated using the iterative phase estimation algorithm and subtracted using single-qubit rotation phase gates. Further, this method optimizes the quantum resource required and makes it more appropriate for near-term quantum hardware. We also demonstrate how the checkerboard deficiency that the D1Q2 scheme presents can be resolved using the D1Q3 scheme. The algorithm is validated by simulating two canonical partial differential equations: the diffusion and Burgers' equations on different quantum simulators. We find good agreement between quantum simulations and classical solutions for the presented algorithm.
97 MATHEMATICS AND COMPUTING↗
Circular economy pathways for decarbonizing aluminum and steel automotive body sheet components in the United States
Decarbonizing vehicle production is essential to reducing automotive sector emissions. This study quantifies greenhouse gas (GHG) emissions from aluminum and steel auto-body sheet components produced in the US. It evaluates the effectiveness of circular economy (CE) strategies (greater closed-loop recycling of pre-consumer scrap, post-consumer scrap, and increased manufacturing yields) to reduce supply chain emissions across different process technology and electricity grid decarbonization pathways. We combine dynamic material flow analysis (2025–2050) with cradle-to-gate life-cycle modeling to assess production emissions and the potential reductions associated with the CE strategies under frozen, moderate, and aggressive technology and grid decarbonization scenarios. Current emissions intensities are estimated at approximately 12.3 kg.CO₂eq/kg of aluminum and 4.3 kg.CO₂eq/kg of steel sheet embedded in the vehicle. Under the frozen decarbonization scenario and current levels of circularity, annual emissions attributable to US aluminum and steel auto-body sheet supply chains could rise by 54 % and 18 % respectively by 2050. Rapid deployment of the CE strategies can cut these annual emissions in 2050 by 52 % for aluminum and 44 % for steel. However, scrap quality constraints lead to saturation points, limiting these benefits unless addressed. Aggressive deployment of low-carbon production technologies and a low-carbon grid reduces the relative benefit of implementing the CE strategies; however, even under the aggressive technology and grid decarbonization scenario, the CE strategies reduce annual emissions by a further 23 %-54 % by 2050. These findings highlight the urgent need to integrate CE strategies into the sheet metal supply chain to support decarbonization efforts and help meet climate targets.
32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
Multi-modal effects on indirect noise induced by turbulent entropy fields
Planar entropy waves are commonly assumed for predicting indirect combustion noise. However, the non-planar and turbulent nature of flows found in most practical combustors challenges this assumption. In the present paper, we examine the indirect noise generated by non-planar and turbulent entropy fields through subsonic nozzles. Firstly, we introduce a new transfer function framework that accounts for the contribution of non-planar Fourier modes of the entropy field to the indirect noise spectra. When applied to a turbulent flow field, this method demonstrates a significant improvement in spectral predictions compared with a conventional approach that only considers the planar mode. Secondly, simulations show that non-planar Fourier modes become significant above a threshold frequency ƒ thresh , found in the mid- to high-frequency range. This contribution of non-planar modes is explained by two-dimensional shear effects that distort the entropy waves. A scaling relation that uses residence times along streamlines is developed for ƒ thresh , showing good agreement with simulation results. Finally, we show that the indirect noise from non-planar entropy modes found in aviation combustors can be significant at frequencies below 1 kHz, which might be relevant in situations of thermo-acoustic instabilities coupled to indirect noise.
42 ENGINEERING↗
Effect of nozzle translation on a bounded vortex flow
A bounded vortex flow is generated by a nozzle that combines azimuthally tilted downward jets in a circular array and a central suction port. When the nozzle is directed toward a bottom surface (called the impingement surface), the flow develops a strong intake vortex spanning from the impingement surface to the suction port, which acts to generate strong shear stress on the impingement surface. This strong shear stress can aid in removing particles attached to the impingement surface or in enhancement of heat or mass transfer, where in typical applications, the nozzle would be translated along the impingement surface. This paper reports on a computational study of the effect of nozzle translation on the bounded vortex flow field. The nozzle translation induces a bending of the bounded vortex in the downstream flow direction, which becomes larger with increase in the nozzle translation speed. We identify an interesting oscillatory phenomenon by which the bounded vortex can resist the bending induced by the translational flow and rotate about a mean position downstream of the suction port. At sufficiently high translation speeds, the oscillatory process can no longer resist the downwash effect from nozzle translation and the vortex eventually detaches from the impingement surface, limiting the effectiveness of this flow device.
42 ENGINEERING↗
Integrated System Planning: Emerging Software Requirements in the Power Industry
Power system planning software remains fragmented across organizational boundaries, with specialized tools for capacity expansion, production cost modeling, power flow, and dynamic analysis operating on incompatible data models and assumptions. This article argues that the fragmentation is not merely a technical problem but a predictable consequence of Conway's law: software architectures mirror the departmental structures within which they are developed. Regulatory milestones like Federal Energy Regulatory Commission (FERC) Order 888 formalized these divisions, but the roots trace back to the distinct engineering disciplines-mechanical, chemical, and electrical-that staffed generation and transmission planning departments in vertically integrated utilities. As the industry moves toward integrated system planning (ISP) that coordinates generation, transmission, and distribution investment decisions, the software ecosystem must evolve accordingly. We identify five categories of software requirements to enable this transition: coherent data inputs decoupled from individual applications, unified and extensible data schemas, modular component representations that support multiple abstraction levels, lifecycle management of planning datasets, and well-defined application programming interface (API) contracts that separate data exchange from algorithmic control. We examine how these requirements interact with three common workflow patterns-serial gate clearing, sequential multiapplication, and convergence oriented-and discuss the interface design principles each demands. We then outline a vision for platform-based planning architectures where specialized analytical services compose through standardized interfaces and where artificial intelligence (AI)/machine learning (ML) tools augment decision support within a disciplined software infrastructure. The practices proposed here offer a path from today's siloed tool collections toward collaborative planning ecosystems capable of handling the complexity of modern power system transformation.
24 POWER TRANSMISSION AND DISTRIBUTION↗
Blade Designs For Improved Multi-Phase Performance In sCO 2 Compressors; Part II - Optical Diagnostics In sCO 2 And Experimental Evaluation With Particle Image Velocimetry
This paper presents the second part of a study in which the leading-edge and suction surface of a compressor blade was modified to delay onset of phase change for sCO 2 compressors operating near the critical point. Using a first-of-its-kind apparatus for the measurement of sCO 2 flow fields, Particle Image Velocimetry (PIV) is used for local flow field measurements of two compressor blade geometries: the modified “biased-wedge,” and a conventional constant thickness blade. Utilizing the developed hardware, the feasibility of a simple, laser-based diagnostic for qualitatively measuring liquid phase regions, is also presented. The design of the optical diagnostics rig, a discussion of numerous challenges, and necessary considerations involved in performing optical-based measurements like PIV, in sCO 2 , are discussed. Velocity field measurements for the modified compressor profile show a much lower suction peak compared to a conventional blade. Furthermore, these results validate numerical results at the tested conditions, where the suction side profile of the biased wedge works to minimize the local pressure gradient.
14 SOLAR ENERGY↗