Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “flow physics”

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 343 records · Page 19
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Influence of the Soret effect on binary-species boundary layers at high pressure

Direct numerical simulations of binary-species temporal boundary layers at high pressure are performed. The main objective is to investigate the influence of the Soret effect on flow physics of binary-species boundary layers where the fluid has a uniform composition. The working fluid is a mixture of 25 % methane and 75 % nitrogen in mass fraction. Although the fluid composition is uniform at the initial condition, the mass fraction of methane increases near the wall when the wall temperature is hotter than the free stream temperature, whereas it decreases when the wall temperature is colder. The non-uniform mass fraction indicates that the uphill diffusion occurs near the wall. Investigation of fluctuations of the mass fraction reveals that the mass fraction fluctuates in the whole boundary layer, indicating that the uphill diffusion occurs even far from the wall. Examination of the species-mass diffusion balance for mean flow fields clarified that the Soret effect flux becomes large near the wall, and the large flux causes the non-uniform profile of the mass fraction near the wall.

Binary Diffusion Coefficient↗
NASA NTRS1998

Experimental Study of Transitional Flow Behavior in a Simulated Low Pressure Turbine

A detailed investigation of the flow physics occurring on the suction side of a simulated Low Pressure Turbine (LPT) blade was performed. A contoured upper wall was designed to simulate the pressure distribution of an actual LPT airfoil onto a flat lower plate. The experiments were carried out for the Reynolds numbers of 35,000, 70,000, 100,000 and 250,000 with four levels of freestream turbulence ranging from 1% to 4%. For the three lower Reynolds numbers, the boundary layer on the flat plate was separated and formed a bubble. The size of laminar separation bubble was measured to be inversely proportional to the freestream turbulence levels and Reynolds numbers. However, no separation was observed for the Re = 250,000 case. The transition on a separated flow was found to proceed through the formation of turbulent spots in the free shear layer as evidenced in the intermittency profiles for Re = 35,000, 70,000 and 100,000. Spectral data show no evidence of Kelvin-Helmholtz or Tollmien-Schlichting instability waves in the free shear layer over a separation bubble (bypass transition). However, the flow visualization revealed the large vortex structures just outside of the bubble and their development to turbulent flow for Re = 50,000, which is similar to that in the free shear layer (separated-flow transition). Therefore, it is fair to say that the bypass and separated-flow transition modes coexist in the transitional flows over the separation bubble for certain conditions. Transition onset and end locations and length determined from intermittency profiles decrease as Reynolds number and freestream turbulence levels increase.

Sohn, Ki Hyeon↗
NASA NTRS2007

Experimental Study of Transitional Flow Behavior in a Simulated Low Pressure Turbine

A detailed investigation of the flow physics occurring on the suction side of a simulated Low Pressure Turbine (LPT) blade was performed. A contoured upper wall was designed to simulate the pressure distribution of an actual LPT airfoil onto a flat lower plate. The experiments were carried out for the Reynolds numbers of 35,000, 70,000, 100,000, and 250,000 with four levels of freestream turbulence ranging from 1 to 4 percent. For the three lower Reynolds numbers, the boundary layer on the flat plate was separated and formed a bubble. The size of laminar separation bubble was measured to be inversely proportional to the freestream turbulence levels and Reynolds numbers. However, no separation was observed for the Re = 250,000 case. The transition on a separated flow was found to proceed through the formation of turbulent spots in the free shear layer as evidenced in the intermittency profiles for Re = 35,000, 70,000, and 100,000. Spectral data show no evidence of Kelvin-Helmholtz of Tollmien-Schlichting instability waves in the free shear layer over a separation bubble (bypass transition). However, the flow visualization revealed the large vortex structures just outside of the bubble and their development to turbulent flow for Re = 50,000, which is similar to that in the free shear layer (separated-flow transition). Therefore, it is fair to say that the bypass and separated-flow transition modes coexist in the transition flows over the separation bubble of certain conditions. Transition onset and end locations and length determined from intermittency profiles decreased as Reynolds number and freestream turbulence levels increase.

Sohn, Ki Hyeon↗
DOE OSTI2025

Concurrent two-way coupling of global and local models across internal boundaries with non-matching discretizations

Coupling local and global models enables efficient simulation of multiscale systems, where global models capture large-scale behavior and local models, with enhanced physics, resolve finer details over a smaller region. Here, this paper presents a mathematically consistent method for coupling physics-based models of varying fidelity across adjacent, non-overlapping subdomains, even when discretizations do not match at the immersed interdomain interfaces. Incompressible Navier-Stokes equations (NSE) constitute the global model while residual-based turbulence model serves as the local high-fidelity model. In addition, a scalar advection-diffusion equation that models the convection of an active scalar field is appended to the turbulence model in the local domain. This scalar field does not have its complement in the global model, giving rise to unequal number of equations at the immersed boundary between local and global models. Interdomain coupling terms are derived via the Variational Multiscale Discontinuous Galerkin (VMDG) method with new developments in scale representation and efficient fine-scale estimation. While transient laminar flows modeled with NSE in the global domain can be resolved with relatively coarse mesh, turbulent flow calculations in the local model require much finer spatial discretizations as well as smaller time-step for appropriately resolving the turbulent flow physics. The proposed framework also accommodates non-matching meshes at the immersed boundaries. Test problems in 2D and 3D numerically showcase the concurrent two-way coupling of unknown fields across the immersed boundaries. The 3D test presents a case with an unequal number of equations, where the scalar field represents the convection of contaminant concentration. This provides more detailed physics in the local region and highlights its application in climate modeling and atmospheric sciences.

Variational Multiscale Discontinuous Galerkin (VMD↗
NASA NTRSDate not supplied

Comparison between Computational and Experimental Wall-Shear Stress and Pressure Measurements in a Supersonic Wind Tunnel

We study the performance of two selected RANS models (SST-$k$-$\omega$ and Baldwin-Lomax) to simulate the mean wall-shear and pressure data obtained from a distributed-roughness-related experimental campaign conducted at Ahmic Aerospace, LLC, located in Ohio, US. A structured compressible Navier-Stokes in-house solver, DPLR, is used to conduct the RANS simulations. As the final objective is to perform the studies under high-enthalpy flow conditions, such as inside the AHF facility at NASA-Ames Research Center (ARC), the current investigation focuses on the low-enthalpy cases using Ahmic's high-Reynolds-number supersonic wind tunnel as the initial step to simplify the flow physics. Furthermore, the wind tunnel simulations are conducted with perfectly-smooth walls to further simplify the problem and perform the initial testing of different solver-based RANS models. The Baldwin-Lomax model is observed to be more accurate than the SST model. However, when these RANS models and their roughness augmentation models are utilized to study their accuracy with respect to the experimental datasets, where the roughness effects are activated, both models perform well for a low roughness height case. However, as the roughness height increases, the SST model more accurately predicts the wall data (shear and pressure) than the Baldwin-Lomax model. Therefore, this present work has built a solid foundation for selecting the best RANS model to study the wall data under high-enthalpy flow conditions at the AHF facility soon.

Prakash Shrestha↗
NASA NTRS2016

Challenges and Progress in Aerodynamic Design of Hybrid Wingbody Aircraft with Embedded Engines

We summarize the contributions to high-fidelity capabilities for analysis and design of hybrid wingbody (HWB) configurations considered by NASA. Specifically, we focus on the embedded propulsion concepts of the N2-B and N3-X configurations, some of the future concepts seriously investigated by the NASA Fixed Wing Project. The objective is to develop the capability to compute the integrated propulsion and airframe system realistically in geometry and accurately in flow physics. In particular, the propulsion system (including the entire engine core-compressor, combustor, and turbine stages) is vastly more difficult and costly to simulate with the same level of fidelity as the external aerodynamics. Hence, we develop an accurate modeling approach that retains important physical parameters relevant to aerodynamic and propulsion analyses for evaluating the HWB concepts. Having the analytical capabilities at our disposal, concerns and issues that were considered to be critical for the HWB concepts can now be assessed reliably and systematically; assumptions invoked by previous studies were found to have serious consequences in our study. During this task, we establish firmly that aerodynamic analysis of a HWB concept without including installation of the propulsion system is far from realistic and can be misleading. Challenges in delivering the often-cited advantages that belong to the HWB are the focus of our study and are emphasized in this report. We have attempted to address these challenges and have had successes, which are summarized here. Some can have broad implications, such as the concept of flow conditioning for reducing flow distortion and the modeling of fan stages. The design optimization capability developed for improving the aerodynamic characteristics of the baseline HWB configurations is general and can be employed for other applications. Further improvement of the N3-X configuration can be expected by expanding the design space. Finally, the support of the System Analysis and Integration Element under the NASA Fixed Wing Project has enabled the development and helped deployment of the capabilities shown in this report.

design optimization↗
NASA NTRS2009

Hypersonic Viscous Flow Over Large Roughness Elements

Viscous flow over discrete or distributed surface roughness has great implications for hypersonic flight due to aerothermodynamic considerations related to laminar-turbulent transition. Current prediction capability is greatly hampered by the limited knowledge base for such flows. To help fill that gap, numerical computations are used to investigate the intricate flow physics involved. An unstructured mesh, compressible Navier-Stokes code based on the space-time conservation element, solution element (CESE) method is used to perform time-accurate Navier-Stokes calculations for two roughness shapes investigated in wind tunnel experiments at NASA Langley Research Center. It was found through 2D parametric study that at subcritical Reynolds numbers, spontaneous absolute instability accompanying by sustained vortex shedding downstream of the roughness is likely to take place at subsonic free-stream conditions. On the other hand, convective instability may be the dominant mechanism for supersonic boundary layers. Three-dimensional calculations for both a rectangular and a cylindrical roughness element at post-shock Mach numbers of 4.1 and 6.5 also confirm that no self-sustained vortex generation from the top face of the roughness is observed, despite the presence of flow unsteadiness for the smaller post-shock Mach number case.

Chang, Chau-Lyan↗
DOE OSTI2023

Frequency Recovery in Power Grids Using High-Performance Computing

Maintaining electric power system stability is paramount, especially in extreme contingencies involving unexpected outages of multiple generators or transmission lines that are typical during severe weather events. Such outages often lead to large supply-demand mismatches followed by subsequent system frequency deviations from their nominal value. The extent of frequency deviations is an important metric of system resilience, and its timely mitigation is a central goal of power system operation and control. This paper develops a novel nonlinear model predictive control (NMPC) method to minimize frequency deviations when the grid is affected by an unforeseen loss of multiple components. Our method is based on a novel multi-period alternating current optimal power flow (ACOPF) formulation that accurately models both nonlinear electric power flow physics and the primary and secondary frequency response of generator control mechanisms. We develop a distributed parallel Julia package for solving the large-scale nonlinear optimization problems that result from our NMPC method and thereby address realistic test instances on existing high-performance computing architectures. Our method demonstrates superior performance in terms of frequency recovery over existing industry practices, where generator levels are set based on the solution of single-period classical ACOPF models.

nonlinear model predictive control↗
NASA NTRS2013

Comparative Study on High-Order Positivity-preserving WENO Schemes

In gas dynamics and magnetohydrodynamics flows, physically, the density and the pressure p should both be positive. In a standard conservative numerical scheme, however, the computed internal energy is obtained by subtracting the kinetic energy from the total energy, resulting in a computed p that may be negative. Examples are problems in which the dominant energy is kinetic. Negative may often emerge in computing blast waves. In such situations the computed eigenvalues of the Jacobian will become imaginary. Consequently, the initial value problem for the linearized system will be ill posed. This explains why failure of preserving positivity of density or pressure may cause blow-ups of the numerical algorithm. The adhoc methods in numerical strategy which modify the computed negative density and/or the computed negative pressure to be positive are neither a conservative cure nor a stable solution. Conservative positivity-preserving schemes are more appropriate for such flow problems. The ideas of Zhang & Shu (2012) and Hu et al. (2012) precisely address the aforementioned issue. Zhang & Shu constructed a new conservative positivity-preserving procedure to preserve positive density and pressure for high-order WENO schemes by the Lax-Friedrichs flux (WENO/LLF). In general, WENO/LLF is too dissipative for flows such as turbulence with strong shocks computed in direct numerical simulations (DNS) and large eddy simulations (LES). The new conservative positivity-preserving procedure proposed in Hu et al. (2012) can be used with any high-order shock-capturing scheme, including high-order WENO schemes using the Roe's flux (WENO/Roe). The goal of this study is to compare the results obtained by non-positivity-preserving methods with the recently developed positivity-preserving schemes for representative test cases. In particular the more difficult 3D Noh and Sedov problems are considered. These test cases are chosen because of the negative pressure/density most often exhibited by standard high-order shock-capturing schemes. The simulation of a hypersonic nonequilibrium viscous shock tube that is related to the NASA Electric Arc Shock Tube (EAST) is also included. EAST is a high-temperature and high Mach number viscous nonequilibrium flow consisting of 13 species. In addition, as most common shock-capturing schemes have been developed for problems without source terms, when applied to problems with nonlinear and/or sti source terms these methods can result in spurious solutions, even when solving a conservative system of equations with a conservative scheme. This kind of behavior can be observed even for a scalar case (LeVeque & Yee 1990) as well as for the case consisting of two species and one reaction (Wang et al. 2012). For further information concerning this issue see (LeVeque & Yee 1990; Griffiths et al. 1992; Lafon & Yee 1996; Yee et al. 2012). This EAST example indicated that standard high-order shock-capturing methods exhibit instability of density/pressure in addition to grid-dependent discontinuity locations with insufficient grid points. The evaluation of these test cases is based on the stability of the numerical schemes together with the accuracy of the obtained solutions.

WENO↗
DOE OSTI2017

Resilient off-grid microgrids: Capacity planning and N-1 security

Over the past century the electric power industry has evolved to support the delivery of power over long distances with highly interconnected transmission systems. Despite this evolution, some remote communities are not connected to these systems. These communities rely on small, disconnected distribution systems, i.e., microgrids to deliver power. However, as microgrids often are not held to the same reliability standards as transmission grids, remote communities can be at risk for extended blackouts. To address this issue, we develop an optimization model and an algorithm for capacity planning and operations of microgrids that include N-1 security and other practical modeling features like AC power flow physics, component efficiencies and thermal limits. Lastly, we demonstrate the computational effectiveness of our approach on two test systems; a modified version of the IEEE 13 node test feeder and a model of a distribution system in a remote community in Alaska.

24 POWER TRANSMISSION AND DISTRIBUTION↗
NASA NTRS2014

Comparitive Study of High-Order Positivity-Preserving WENO Schemes

In gas dynamics and magnetohydrodynamics flows, physically, the density ρ and the pressure p should both be positive. In a standard conservative numerical scheme, however, the computed internal energy is The ideas of Zhang & Shu (2012) and Hu et al. (2012) precisely address the aforementioned issue. Zhang & Shu constructed a new conservative positivity-preserving procedure to preserve positive density and pressure for high-order Weighted Essentially Non-Oscillatory (WENO) schemes by the Lax-Friedrichs flux (WENO/LLF). In general, WENO/LLF is obtained by subtracting the kinetic energy from the total energy, resulting in a computed p that may be negative. Examples are problems in which the dominant energy is kinetic. Negative ρ may often emerge in computing blast waves. In such situations the computed eigenvalues of the Jacobian will become imaginary. Consequently, the initial value problem for the linearized system will be ill posed. This explains why failure of preserving positivity of density or pressure may cause blow-ups of the numerical algorithm. The adhoc methods in numerical strategy which modify the computed negative density and/or the computed negative pressure to be positive are neither a conservative cure nor a stable solution. Conservative positivity-preserving schemes are more appropriate for such flow problems. too dissipative for flows such as turbulence with strong shocks computed in direct numerical simulations (DNS) and large eddy simulations (LES). The new conservative positivity-preserving procedure proposed in Hu et al. (2012) can be used with any high-order shock-capturing scheme, including high-order WENO schemes using the Roe's flux (WENO/Roe). The goal of this study is to compare the results obtained by non-positivity-preserving methods with the recently developed positivity-preserving schemes for representative test cases. In particular the more di cult 3D Noh and Sedov problems are considered. These test cases are chosen because of the negative pressure/density most often exhibited by standard high-order shock-capturing schemes. The simulation of a hypersonic nonequilibrium viscous shock tube that is related to the NASA Electric Arc Shock Tube (EAST) is also included. EAST is a high-temperature and high Mach number viscous nonequilibrium ow consisting of 13 species. In addition, as most common shock-capturing schemes have been developed for problems without source terms, when applied to problems with nonlinear and/or sti source terms these methods can result in spurious solutions, even when solving a conservative system of equations with a conservative scheme. This kind of behavior can be observed even for a scalar case as well as for the case consisting of two species and one reaction.. This EAST example indicated that standard high-order shock-capturing methods exhibit instability of density/pressure in addition to grid-dependent discontinuity locations with insufficient grid points. The evaluation of these test cases is based on the stability of the numerical schemes together with the accuracy of the obtained solutions.

Schemes↗
DOE OSTI2024

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↗
DOE OSTI2022

Agent-Based Modeling for the Circular Economy: Lessons Learned From Three Case Studies

The circular economy (CE) aims at decoupling human activities from resource use, creating wealth in the process. Recently, many scholars have questioned the link between increased circularity and sustainability, resulting in many methodological approaches being developed for that purpose. This presentation summarizes the insights gained from the application of agent-based modeling (ABM) to study the techno-economic and social conditions promoting circularity and sustainability of three technologies: photovoltaic (PV) modules, hard disk drives (HDDs), and wind blades. Four main categories of agents are defined in the ABM: asset owners, service providers (e.g., refurbishers), recyclers, and manufacturers. Two main CE strategies are represented: lifetime extension (through repair or reuse) and recycling. The developed models start by projecting installed capacities and end-of-life (EOL) quantities. Then the theory of planned behavior - a social psychology model explaining behavior adoption based on attitude, peer influence, and costs - is used to model the asset owners' EOL decision (i.e., landfill, recycle or extend the lifetime of the asset). Then, the quantities of assets flowing to the recycler and service provider agents and quantities of materials flowing to manufacturers are computed. Recyclers' economies of scale are dynamically modeled, and the value generated by the CE strategies for the recyclers, service providers, and manufacturers is computed within the model. When data are available, avoided greenhouse gas emissions resulting from the CE strategies adoption are calculated exogenously from the ABM simulations. Results show that with improved used PV modules warranties, the reuse CE strategy adoption increases from 1% to 23% between 2020 and 2050. Similarly, improved standards could enhance HDDs end-users trust in data-wiping - a prerequisite to reuse - leading to a 3-fold increase in the reuse rate and avoid about 5 million tons of CO2 eq by 2050. Regarding wind blades, 5-15 years lifetime extension could reduce EOL blade quantities by 13%. High costs and logistic issues prevent blades from being recycled in greater quantities. One insight from the case studies is the necessity to have mature secondary markets for reuse to be a viable option. Interestingly, PV reuse is limited by the willingness of PV owners to purchase used modules (on the demand side), while HDDs reuse is constrained by the lack of trust toward data-wiping (limiting the supply of used HDDs). The six limits of the CE concepts presented by Korhonen et al. (2018) are finally used to interpret the results. The HDD case study is an exemplary lock-in, where the first accepted practice (shredding) retains most of the market. The PV results illustrate the technical limitation to reuse, as the growing demand cannot be supplied entirely with used PV modules. The wind case study shows the relevance of clearly defining physical flows - what type of waste should wind blades be considered, how should they be transported and landfilled? - a crucial consideration that also applies to PV modules. Finally, the three case studies highlight the relevance of studying a technology's technical, economic, and market material efficiency potentials altogether and the potential benefit of coupling ABM to life cycle assessment.

agent-based modeling↗
DOE OSTI2020

Optimal economic operation of liquid petroleum products pipeline systems

The majority of overland transport needs for crude petroleum and refined petroleum products are met using pipelines. Numerous studies have developed optimization methods for design of these systems in order to minimize construction costs while meeting capacity requirements. Here, we formulate problems to optimize the operations of existing single liquid commodity pipeline systems subject to physical flow and pump engineering constraints. The objectives are to maximize the economic value created for users of the system and to minimize operating costs. We present a general computational method for this class of continuous, non-convex nonlinear programs, and examine the use of pump operating settings and flow allocations as decision variables. Furthermore, the approach is applied to compute optimal operating regimes and perform engineering economic sensitivity analyses for a case study of a crude oil pipeline developed using publicly available data.

02 PETROLEUM↗
DOE OSTI2020

Effect of the numerical discretization scheme in Shock-Driven turbulent mixing simulations

In this work, we evaluate the effects of distinct numerical strategies in simulations of shock-driven turbulent mixing. An air-SF 6 -air gas curtain subjected to Mach 1.2 shock-waves is computed with Implicit Large-Eddy Simulation based on three numerical schemes: directional-split with Harten-Lax-van Leer (HLL) solver; directional-unsplit using HLL-Contact (HLLC) solver; and directional-unsplit with HLLC solver and a Low Mach number Correction (LMC). The results illustrate the importance of the numerical strategy to the accuracy of the predictions. Whereas both split and unsplit schemes result in a similar spatial development of the initial shock-driven instability, only the unsplit schemes can predict the turbulent mixing transition after reshock observed by the laboratory experiments. Such feature increases the mixing rate of the two fluids, this being particularly pronounced when the LMC is active due to i) the reduced flow Mach number, and ii) the larger effective Reynolds number. Since the selected mixing problem is driven by the deposition of vorticity at the fluids’ interface, the resultant flow physics is analyzed by investigating the contribution of distinct inviscid mechanisms to the production of vorticity: baroclinicity, stretching, and dilatation. As expected, it is observed that the production of vorticity is initially dominated by the baroclinicity mechanism. Yet, the relevance of the remainder mechanisms is enhanced after the reshock and may even surpass the baroclinicity term.

42 ENGINEERING↗
DOE OSTI2024

Time-lapse seismic inversion for CO 2 saturation with SeisCO2Net: An application to Frio-II site

Seismic monitoring of geological CO 2 storage (GCS) involves highly nonlinear seismic inversion and petrophysical inversion, making it challenging to estimate CO 2 volume efficiently and detect possible early CO 2 leakages. Deep learning (DL) using convolutional neural networks (CNNs) has shown promise in solving highly nonlinear seismic inversion problems. However, direct estimation of CO 2 plume extent/saturation from time-lapse seismic gathers using DL is still underexplored, with no reported field applications to date. The investigation of field data is primarily hindered by scarcity of field data for neural network training. Other obstacles include highly nonlinear seismic-petrophysics inverse relationship, and presence of noise in field seismic data. We introduce SeisCO2Net, a deep CNN that predicts CO 2 saturation maps directly from time-lapse full waveform shot gathers. For training, we use site-specific geological information, fluid flow physics, rock physics, and seismic modeling to generate synthetic datasets that closely resemble the CO 2 storage site. Synthetic tests show promising results, inspiring us to apply SeisCO2Net's trained weights on field data collected at Frio-II GCS site by leveraging transfer learning principles. As reference, we compare SeisCO2Net's predicted CO 2 saturation maps with results obtained from physics-based inversion. Our analyses show both methods display similar CO 2 plume shapes, reasonable CO 2 plume characteristics, and comparable saturation values. Our results suggest pre-training CNNs on physics-informed synthetic datasets and then applying the learned weights to field data is a viable approach to estimating field CO 2 saturation. This method effectively addresses the scarcity of field training data, thus encouraging the feasibility of long-term GCS monitoring.

58 GEOSCIENCES↗
DOE OSTI2022

Programming approaches for scalability, performance, and portability of combustion physics codes

Here, this paper presents the process, strategy, and results associated with porting a typical combustion physics flow solver to current state-of-the-art and future massively-parallel computer architectures. Major focus is placed on the distinct algorithmic structure of these types of codes and how it can be integrated with modern programming paradigms for heterogeneous platforms (i.e., distributed many-core systems with accelerators). An end-to-end case study is presented that exemplifies the process in a generic manner, which then serves as a clear guide with respect to the strategy and best practices leading to a robust and adaptable framework that performs well, is durable over time, is portable, and requires minimal human-effort. This end is accomplished beginning with the use of a mature, validated, structured, multiblock code framework optimized for application of both Large Eddy Simulation (LES) and Direct Numerical Simulation (DNS). This code has been ported to a variety of platforms over the past decade, including most recently the Oak Ridge Leadership Computing Facility’s “Summit” Platform. The experience gained on these multiple platforms provides general insights and thus the results presented are not specific to any one code or platform other than the overarching trend toward distributed many-core systems with accelerators in order to move toward exascale performance. The resultant performance and scalability of the ported code is demonstrated on a real-world application; a state-of-the-art rotating detonation rocket engine simulation that matches the complex geometry and boundary conditions imposed as part of a companion experimental campaign.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

Nonlinear saturation of the slab ITG instability and zonal flow generation with fully kinetic ions

Fully kinetic turbulence models are of interest for their potential to validate or replace gyrokinetic models in plasma regimes where the gyrokinetic expansion parameters are marginal. Here, we demonstrate fully kinetic ion capability by simulating the growth and nonlinear saturation of the ion-temperature-gradient instability in shearless slab geometry assuming adiabatic electrons and including zonal flow dynamics. The ion trajectories are integrated using the Lorentz force, and the cyclotron motion is fully resolved. Linear growth and nonlinear saturation characteristics show excellent agreement with analogous gyrokinetic simulations across a wide range of parameters. The fully kinetic simulation accurately reproduces the nonlinearly generated zonal flow. In conclusion, this work demonstrates nonlinear capability, resolution of weak gradient drive, and zonal flow physics, which are critical aspects of modeling plasma turbulence with full ion dynamics.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗