Search NASA⌕ Search

SEARCH · Search NASA

Results for “Simulator”

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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2013

Digital Full-Scope Simulation of a Conventional Nuclear Power Plant Control Room, Phase 2: Installation of a Reconfigurable Simulator to Support Nuclear Plant Sustainability

The U.S. Department of Energy’s Light Water Reactor Sustainability program has developed a control room simulator in support of control room modernization at nuclear power plants in the U.S. This report highlights the recent completion of this reconfigurable, full-scale, full-scope control room simulator buildout at the Idaho National Laboratory. The simulator is fully reconfigurable, meaning it supports multiple plant models developed by different simulator vendors. The simulator is full-scale, using glasstop virtual panels to display the analog control boards found at current plants. The present installation features 15 glasstop panels, uniquely achieving a complete control room representation. The simulator is also full-scope, meaning it uses the same plant models used for training simulators at actual plants. Unlike in the plant training simulators, the deployment on glasstop panels allows a high degree of customization of the panels, allowing the simulator to be used for research on the design of new digital control systems for control room modernization. This report includes separate sections discussing the glasstop panels, their layout to mimic control rooms at actual plants, technical details on creating a multi-plant and multi-vendor reconfigurable simulator, and current efforts to support control room modernization at U.S. utilities. The glasstop simulator provides an ideal testbed for prototyping and validating new control room concepts. Equally importantly, it is helping create a standardized and vetted human factors engineering process that can be used across the nuclear industry to ensure control room upgrades maintain and even improve current reliability and safety.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2016

Comparison of Measured and Numerically Simulated Turbulence Statistics in a Convective Boundary Layer Over Complex Terrain

High resolution numerical simulation can provide insight into important physical processes that occur within the planetary boundary layer (PBL). The present work employs large eddy simulation (LES) nested within the Weather Forecasting and Research (WRF) framework to simulate real conditions in the convective PBL over an area of complex terrain. A multiple nesting approach has been used to downsize the grid spacing size from 12.15 km (mesoscale) to 0.03 km (LES). To minimize artifacts in theWRF-LES solutions resulting from using an inappropriate grid spacing on the WRF-Meso domain, a careful selection of grid spacing has been conducted. The WRF-LES results have been evaluated with in situ and remote sensing observations collected during the US Department of Energy supported Columbia BasinWind Energy Study (CBWES). The comparison for the first and second order moments, turbulence spectrum, and probability density function (pdf) of wind speed show good agreement between the simulation and data. Furthermore, the WRF-LES variables show a great deal of variability in space and time caused by the complex topography in the LES domain. The WRF-LES results show that the flow structures, such as roll vortices and convective cells, over the domain vary depending on both the spatial location and time of the day. In addition to basic studies related to boundary-layer meteorology, results from these simulations can be used in other applications, such as studying the wind energy resources.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2016

Intercomparison of 3D pore-scale flow and solute transport simulation methods

Multiple numerical approaches have been developed to simulate porous media fluid flow and solute transport at the pore scale. These include methods that 1) explicitly model the three-dimensional geometry of pore spaces and 2) those that conceptualize the pore space as a topologically consistent set of stylized pore bodies and pore throats. In previous work we validated a model of class 1, based on direct numerical simulation using computational fluid dynamics (CFD) codes, against magnetic resonance velocimetry (MRV) measurements of pore-scale velocities. Here we expand that validation to include additional models of class 1 based on the immersed-boundary method (IMB), lattice Boltzmann method (LBM), smoothed particle hydrodynamics (SPH), as well as a model of class 2 (a pore-network model or PNM). The PNM approach used in the current study was recently improved and demonstrated to accurately simulate solute transport in a two-dimensional experiment. While the PNM approach is computationally much less demanding than direct numerical simulation methods, the effect of conceptualizing complex three-dimensional pore geometries on solute transport in the manner of PNMs has not been fully determined. We apply all four approaches (CFD, LBM, SPH and PNM) to simulate pore-scale velocity distributions and nonreactive solute transport, and intercompare the model results with previously reported experimental observations. Experimental observations are limited to measured pore-scale velocities, so solute transport comparisons are made only among the various models. Comparisons are drawn both in terms of macroscopic variables (e.g., permeability, solute breakthrough curves) and microscopic variables (e.g., local velocities and concentrations).

computational fluid dynamics↗
DOE OSTI2021

Linear and Nonlinear Solvers for Simulating Multiphase Flow within Large-Scale Engineered Subsurface Systems

Simulation of multiphase flow in the subsurface is well-known to be computationally challenging. While there have been many studies that have explored approaches to overcoming these challenges, they often utilize relatively simple case studies. In this paper, we focus on the unique numerical challenges posed by modeling large-scale engineered subsurface systems, characterized by discrete features embedded in a heterogeneous natural subsurface setting. The man-made features such as shafts, tunnels, and barriers often cause multiple challenges in modeling the domain for multiphase porous media flow. This flow scenario can have a wide range of applications such as nuclear waste repositories, enhanced recovery of a petroleum reservoir, geothermal engineering, and carbon sequestration. An example of these severe numerical challenges is the case of performance assessment (PA) for Waste Isolation Pilot Plant (WIPP), the only operating deep geological repository in the US, which simulates extreme material properties of bedded salt rock formation and extreme contrast due to open excavation next to the formation. The models have extremes not only of permeability and porosity but also of the constitutive models needed for multiphase flow; additionally, they have process models like salt creep closure reducing porosity over time, fracturing in clay and anhydrite interbeds of the bedded salt, gas generation from the waste materials, and unintentional human borehole intrusions in some scenarios. Numerical simulations require the solution of coupled systems of nonlinear PDEs; in our work, we use the open-source simulator PFLOTRAN which is based on Finite Volume discretization. The solution of the nonlinear equations requires use of the Newton-Raphson iteration at each time step, which entails the solution of the linearized Jacobian system at each iteration. The effects of all the processes (i.e., large number of unknowns, highly nonlinear constitutive relations, large contrasts in material properties in short distances) lead to an ill-conditioned Jacobian matrix that severely challenges traditional linear solver, i.e., stabilized biconjugate gradient with block Jacobi incomplete LU preconditioner (BCGS-ILU) leading to non-convergence for traditional Newton-Raphson nonlinear solver causing unacceptably long computation time for each model. This paper presents linear solvers such as constrained pressure residual (CPR) two-stage preconditioner with alternate-block-factorization (ABF) and quasi- implicit pressure and explicit saturation (QIMPES) decouplers and flexible generalized residual solver (FGMRES). The new general-purpose nonlinear solver, Newton trust-region dogleg Cauchy (NTRDC), is also introduced to resolve extreme nonlinearities in the models. We demonstrate the effectiveness of each method relative to the default BCGS-Newton solver. The two best cases had nearly 50 times speed-up and achieved completion of a simulation in 14 hours that never completed due to non-convergence with the default solver. We also investigate the strong scalability of each method and discuss some of the deficiencies found for Block Jacobi preconditioner using parallel domain decomposition, and node packing effects of modern processor architecture.

Preconditioner, Nonlinear, Porous media, Multiphas↗
DOE OSTI2013

Simulating Ru L 3 -Edge X-ray Absorption Spectroscopy with Time-Dependent Density Functional Theory: Model Complexes and Electron Localization in Mixed-Valence Metal Dimers

Ruthenium L3-edge X-ray absorption (XA) spectroscopy probes unoccupied 4d orbitals of the metal atom and is increasingly being used to investigate the local electronic structure in ground and excited electronic states of Ru complexes. The simultaneous development of computational tools for simulating Ru L3-edge spectra is crucial for interpreting the spectral features at a molecular level. This study demonstrates that time-dependent density functional theory (TDDFT) is a viable and predictive tool for simulating ruthenium L3-edge XA spectroscopy. We systematically investigate the effects of exchange correlation functional and implicit and explicit solvent interactions on a series of RuII and RuIII complexes in their ground and electronic excited states. The TDDFT simulations reproduce all of the experimentally observed features in Ru L3-edge XA spectra within the experimental resolution (0.4 eV). Our simulations identify ligand-specific charge transfer features in complicated Ru L3-edge spectra of [Ru(CN)6]4- and RuII polypyridyl complexes illustrating the advantage of using TDDFT in complex systems. We conclude that the B3LYP functional most accurately predicts the transition energies of charge transfer features in these systems. We use our TDDFT approach to simulate experimental Ru L3-edge XA spectra of transition metal mixed-valence dimers of the form [(NC)5MII-CN-RuIII(NH3)5] (where M = Fe or Ru) dissolved in water. Our study determines the spectral signatures of electron delocalization in Ru L3-edge XA spectra. We find that the inclusion of explicit solvent molecules is necessary for reproducing the spectral features and the experimentally determined valencies in these mixed-valence complexes. This study validates the use of TDDFT for simulating Ru 2p excitations using popular quantum chemistry codes and providing a powerful interpretive tool for equilibrium and ultrafast Ru L3-edge XA spectroscopy.

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

BESTEST-GSR (Building Energy Simulation Test - Generation Simulation and Reporting) [SWR 18-23]

Building Energy Simulation Test (BESTEST) is an NREL-developed method to validate the qualitative performance of different whole building simulations engines relative to each other. https://www.nrel.gov/docs/legosti/old/6231.pdf For EnergyPlus® and OpenStudio®, generation of test cases, simulation, and reporting has been a manual or partially automated process that typically took a number of months after a new version of EnergyPlus came out. Building Energy Simulation Test - Generation Simulation and Reporting (BESTEST-GSR) was developed to fully automate this process so it is less error-prone, less time-consuming, and can be quickly run on request for any development branch of EnergyPlus and OpenStudio. The initial use case is to directly support OpenStudio and EnergyPlus development, but the software is configured to support use by external software developers that build tools on top of EnergyPlus and/or OpenStudio. Tools built on top of OpenStudio can leverage the test generation, simulation, and reporting, while tools built on EnergyPlus can leverage the test simulation and reporting. The model generation is accomplished by running OpenStudio measures for envelope, heating, or cooling test cases. At this stage the specified weather file required for a specific test case is assigned to the model. Simulation is enabled using the OpenStudio Command Line Interface (CLI) to run a series of OpenStudio and Reporting measures to a specified seed model. The specified workflow is defined with an OpenStudio OSW file. The first measure generates the test case model, while the second measure is a reporting measure that runs after the simulation to collect necessary simulation results from the EnergyPlus SQL file. The Reporting stage gathers data from individual OpenStudio OSW files and populates a series of BESTEST Microsoft Excel spreadsheets. BESTEST-GSR also gathers zip files from individual simulations for diagnostic purposes.

Goldwasser, David↗
NASA NTRS1983

Software simulator for multiple computer simulation system

A description is given of the structure and use of a computer program that simulates the operation of a parallel processor simulation system. The program is part of an investigation to determine algorithms that are suitable for simulating continous systems on a parallel processor configuration. The simulator is designed to accurately simulate the problem-solving phase of a simulation study. Care has been taken to ensure the integrity and correctness of data exchanges and to correctly sequence periods of computation and periods of data exchange. It is pointed out that the functions performed during a problem-setup phase or a reset phase are not simulated. In particular, there is no attempt to simulate the downloading process that loads object code into the local, transfer, and mapping memories of processing elements or the memories of the run control processor and the system control processor. The main program of the simulator carries out some problem-setup functions of the system control processor in that it requests the user to enter values for simulation system parameters and problem parameters. The method by which these values are transferred to the other processors, however, is not simulated.

Ogrady, E. P.↗
NASA NTRS2013

Simulation System Fidelity Assessment at the Vertical Motion Simulator

Fidelity is a word that is often used but rarely understood when talking about groundbased simulation. Assessing the cueing fidelity of a ground based flight simulator requires a comparison to actual flight data either directly or indirectly. Two experiments were conducted at the Vertical Motion Simulator using the GenHel UH-60A Black Hawk helicopter math model that was directly compared to flight data. Prior to the experiment the simulator s motion and visual system frequency responses were measured, the aircraft math model was adjusted to account for the simulator motion system delays, and the motion system gains and washouts were tuned for the individual tasks. The tuned motion system fidelity was then assessed against the modified Sinacori criteria. The first experiments showed similar handling qualities ratings (HQRs) to actual flight for a bob-up and sidestep maneuvers. The second experiment showed equivalent HQRs between flight and simulation for the ADS33 slalom maneuver for the two pilot participants. The ADS33 vertical maneuver HQRs were mixed with one pilot rating the flight and simulation the same while the second pilot rated the simulation worse. In addition to recording HQRs on the second experiment, an experimental Simulation Fidelity Rating (SFR) scale developed by the University of Liverpool was tested for applicability to engineering simulators. A discussion of the SFR scale for use on the Vertical Motion Simulator is included in this paper.

Beard, Steven D.↗
DOE OSTI2020

Toward Exascale: Overview of Large Eddy Simulations and Direct Numerical Simulations of Nuclear Reactor Flows with the Spectral Element Method in Nek5000

At the beginning of the last decade, Petascale supercomputers (i.e., computers capable of more than 1 petaFLOP) emerged. Now, at the dawn of exascale supercomputing, we provide a review of recent landmark simulations of portions of reactor components with turbulence-resolving techniques that this computational power has made possible. In fact, these simulations have provided invaluable insight into flow dynamics, which is difficult or often impossible to obtain with experiments alone. We focus on simulations performed with the spectral element method, as this method has emerged as a powerful tool to deliver massively parallel calculations at high fidelity by using large eddy simulation or direct numerical simulation. We also limit this paper to constant-property incompressible flow of a Newtonian fluid in the absence of other body or external forces, although the method is by no means limited to this class of flows. We briefly review the fundamentals of the method and the reasons it is compelling for the simulation of nuclear engineering flows. We review in detail a series of Petascale simulations, including the simulations of helical coil steam generators, fuel assemblies, and pebble beds. Even with Petascale computing, however, limitations for nuclear modeling and simulation tools remain. In particular, the size and scope of turbulence-resolving simulations are still limited by computing power and resolution requirements, which scale with the Reynolds number. In the final part of this paper, we discuss the future of the field, including recent advancements in emerging architectures such as GPUbased supercomputers, which are expected to power the next generation of high-performance computers.

computational fluid dynamics↗
DOE OSTI2015

Changes in Moisture Flux over the Tibetan Plateau during 1979-2011: Insights from a High Resolution Simulation

Net precipitation (precipitation minus evapotranspiration, P-E) changes between 1979 and 2011 from a high resolution regional climate simulation and its reanalysis forcing are analyzed over the Tibet Plateau (TP) and compared to the global land data assimilation system (GLDAS) product. The high resolution simulation better resolves precipitation changes than its coarse resolution forcing, which contributes dominantly to the improved P-E change in the regional simulation compared to the global reanalysis. Hence, the former may provide better insights about the drivers of P-E changes. The mechanism behind the P-E changes is explored by decomposing the column integrated moisture flux convergence into thermodynamic, dynamic, and transient eddy components. High-resolution climate simulation improves the spatial pattern of P-E changes over the best available global reanalysis. High-resolution climate simulation also facilitates new and substantial findings regarding the role of thermodynamics and transient eddies in P-E changes reflected in observed changes in major river basins fed by runoff from the TP. The analysis revealed the contrasting convergence/divergence changes between the northwestern and southeastern TP and feedback through latent heat release as an important mechanism leading to the mean P-E changes in the TP.

high resolution simulation↗
DOE OSTI2021

The AGORA High-resolution Galaxy Simulations Comparison Project. III. Cosmological Zoom-in Simulation of a Milky Way–mass Halo

In this work, we present a suite of high-resolution cosmological zoom-in simulations to z = 4 of a 10 12 M⊙ halo at z = 0, obtained using seven contemporary astrophysical simulation codes (Art-I, Enzo, Ramses, Changa, Gadget-3, Gear, and Gizmo) widely used in the numerical galaxy formation community. The physics prescriptions for gas cooling and heating and star formation are the same as the ones used in our previous Assembling Galaxies of Resolved Anatomy (AGORA) disk comparison but now account for the effects of cosmological processes such as the expansion of the universe, intergalactic gas inflow, and the cosmic ultraviolet background radiation emitted by massive stars and quasars. In this work, we introduce the most careful comparison yet of galaxy formation simulations run by different code groups, together with a series of four calibration steps each of which is designed to reduce the number of tunable simulation parameters adopted in the final run. In the first two steps, we methodically calibrate the gas physics, such as cooling and heating, in simulations without star formation. In the third step, we seek agreement on the total stellar mass produced with the common star formation prescription used in the AGORA disk comparison, in stellar-feedback-free simulations. In the last calibration step, we activate stellar feedback, where each code group is asked to set the feedback prescription to as close to the most widely used one in its code community as possible, while aiming for convergence in the stellar mass at z = 4 to the values predicted by semiempirical models. After all the participating code groups successfully complete the calibration steps, we achieve a suite of cosmological simulations with similar mass assembly histories down to z = 4. With numerical accuracy that resolves the internal structure of a target halo (≲100 physical pc at z = 4), we find that the codes overall agree well with one another, e.g., in gas and stellar properties, but also show differences, e.g., in circumgalactic medium (CGM) properties. We argue that, if adequately tested in accordance with our proposed calibration steps and common parameters, high-resolution cosmological zoom-in simulations can have robust and reproducible results. New code groups are invited to join and enrich this comparison by generating equivalent models or to test the code's compatibility on their own, by adopting the common initial conditions, the common easy-to-implement physics package, and the proposed calibration steps. Further analyses of the zoom-in simulations presented here will be presented in forthcoming reports from the AGORA Collaboration, including studies of the CGM, simulations by additional codes, and results at lower redshift.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

Monte Carlo Simulations of the Corrosion of Aluminoborosilicate Glasses

Aluminum is one of the most common components included in nuclear waste glasses. Therefore, Monte Carlo (MC) simulations were carried out to investigate the influence of aluminum on the rate and mechanism of dissolution of sodium borosilicate glasses in static conditions. The glasses studied were in the compositional range (70-2x)% SiO2 x% Al2O3 15% B2O3 (15+x)% Na2O, where 0 ≤ x ≤ 15%. The simulation results show that increasing amounts of aluminum in the pristine glasses slow down the initial rate of dissolution as determined from the rate of boron release. However, the extent of corrosion - as measured by the total amount of boron release - initially increases with addition of Al2O3, up to 5 Al2O3 mol%, but subsequently decreases with further Al2O3 addition. The MC simulations reveal that this behavior is due to the interplay between two opposing mechanisms: (1) aluminum slows down the kinetics of hydrolysis/condensation reactions that drive the reorganization of the glass surface and eventual formation of a blocking layer; and (2) aluminum strengthens the glass thereby increasing the lifetime of the upper part of its surface and allowing for more rapid formation of a blocking layer. Additional MC simulations were performed whereby a process representing the formation of a secondary aluminosilicate phase was included. Secondary phase formation draws dissolved glass components out of the aqueous solution, thereby diminishing the rate of condensation and delaying the formation of a blocking layer. As a result, the extent of corrosion is found to increase continuously with increasing Al2O3 content, as observed experimentally. For Al2O3 < 10 mol%, the MC simulations also indicate that, because the secondary phase solubility eventually controls the aluminum content in the part of the altered layer in contact with the bulk aqueous solution, the dissolved aluminum and silicon concentrations at steady state are not dependent on the Al2O3 content of the pristine aluminoborosilicate glass.

Environmental Molecular Sciences Laboratory↗
DOE OSTI2012

Comparison of Cellulose Ib Simulations with Three Carbohydrate Force Fields

Molecular dynamics simulations of cellulose have recently become more prevalent due to increased interest in renewable energy applications, and many atomistic and coarse-grained force fields exist that can be applied to cellulose. However, to date no systematic comparison between carbohydrate force fields has been conducted for this important system. To that end, we present a molecular dynamics simulation study of hydrated, 36-chain cellulose I{beta} microfibrils at room temperature with three carbohydrate force fields (CHARMM35, GLYCAM06, and Gromos 45a4) up to the near-microsecond time scale. Our results indicate that each of these simulated microfibrils diverge from the cellulose I{beta} crystal structure to varying degrees under the conditions tested. The CHARMM35 and GLYCAM06 force fields eventually result in structures similar to those observed at 500 K with the same force fields, which are consistent with the experimentally observed high-temperature behavior of cellulose I. The third force field, Gromos 45a4, produces behavior significantly different from experiment, from the other two force fields, and from previously reported simulations with this force field using shorter simulation times and constrained periodic boundary conditions. For the GLYCAM06 force field, initial hydrogen-bond conformations and choice of electrostatic scaling factors significantly affect the rate of structural divergence. Our results suggest dramatically different time scales for convergence of properties of interest, which is important in the design of computational studies and comparisons to experimental data. This study highlights that further experimental and theoretical work is required to understand the structure of small diameter cellulose microfibrils typical of plant cellulose.

09 BIOMASS FUELS↗
DOE OSTI2021

Practical verification protocols for analog quantum simulators

Abstract Analog quantum simulation is expected to be a significant application of near-term quantum devices. Verification of these devices without comparison to known simulation results will be an important task as the system size grows beyond the regime that can be simulated classically. We introduce a set of experimentally-motivated verification protocols for analog quantum simulators, discussing their sensitivity to a variety of error sources and their scalability to larger system sizes. We demonstrate these protocols experimentally using a two-qubit trapped-ion analog quantum simulator and numerically using models of up to five qubits.

Physics↗
DOE OSTI2014

Structure of Cu{sub 64.5}Zr{sub 35.5} metallic glass by reverse Monte Carlo simulations

Reverse Monte Carlo simulations (RMC) have been widely used to generate three dimensional (3D) atomistic models for glass systems. To examine the reliability of the method for metallic glass, we use RMC to predict the atomic configurations of a “known” structure from molecular dynamics (MD) simulations, and then compare the structure obtained from the RMC with the target structure from MD. We show that when the structure factors and partial pair correlation functions from the MD simulations are used as inputs for RMC simulations, the 3D atomistic structure of the glass obtained from the RMC gives the short- and medium-range order in good agreement with those from the target structure by the MD simulation. These results suggest that 3D atomistic structure model of the metallic glass alloys can be reasonably well reproduced by RMC method with a proper choice of input constraints.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2014

Verification of gyrokinetic particle simulation of current-driven instability in fusion plasmas. I. Internal kink mode

The gyrokinetic toroidal code (GTC) capability has been extended for simulating internal kink instability with kinetic effects in toroidal geometry. The global simulation domain covers the magnetic axis, which is necessary for simulating current-driven instabilities. Here, GTC simulation in the fluid limit of the kink modes in cylindrical geometry is verified by benchmarking with a magnetohydrodynamic eigenvalue code. Gyrokinetic simulations of the kink modes in the toroidal geometry confirum that ion kinetic effects greatly reduce the growth rate even when the banana orbit width is much smaller than the radial width of the perturbed current layer at the mode rational surface.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

Nonlinear gyrokinetic simulations of the I-mode high confinement regime and comparisons with experiment

For the first time, nonlinear gyrokinetic simulations of I-mode plasmas are performed and compared with experiment. I-mode is a high confinement regime, featuring energy confinement similar to H-mode, but without enhanced particle and impurity particle confinement [D. G. Whyte et al., Nucl. Fusion 50, 105005 (2010)]. As a consequence of the separation between heat and particle transport, I-mode exhibits several favorable characteristics compared to H-mode. The nonlinear gyrokinetic code GYRO [J. Candy and R. E. Waltz, J Comput. Phys. 186, 545 (2003)] is used to explore the effects of E × B shear and profile stiffness in I-mode and compare with L-mode. The nonlinear GYRO simulations show that I-mode core ion temperature and electron temperature profiles are more stiff than L-mode core plasmas. Scans of the input E × B shear in GYRO simulations show that E × B shearing of turbulence is a stronger effect in the core of I-mode than L-mode. The nonlinear simulations match the observed reductions in long wavelength density fluctuation levels across the L-I transition but underestimate the reduction of long wavelength electron temperature fluctuation levels. The comparisons between experiment and gyrokinetic simulations for I-mode suggest that increased E × B shearing of turbulence combined with increased profile stiffness are responsible for the reductions in core turbulence observed in the experiment, and that I-mode resembles H-mode plasmas more than L-mode plasmas with regards to marginal stability and temperature profile stiffness.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

High fidelity simulation and analysis of liquid jet atomization in a gaseous crossflow at intermediate Weber numbers

Recent advances in numerical methods coupled with the substantial enhancements in computing power and the advent of high performance computing have presented first principle, high fidelity simulation as a viable tool in the prediction and analysis of spray atomization processes. The credibility and potential impact of such simulations, however, has been hampered by the relative absence of detailed validation against experimental evidence. The numerical stability and accuracy challenges arising from the need to simulate the high liquid-gas density ratio across the sharp interfaces encountered in these flows are key reasons for this. In this work we challenge this status quo by presenting a numerical model able to deal with these challenges, employing it in simulations of liquid jet in crossflow atomization and performing extensive validation of its results against a carefully executed experiment with detailed measurements in the atomization region. We then proceed to the detailed analysis of the flow physics. The computational model employs the coupled level set and volume of fluid approach to directly capture the spatiotemporal evolution of the liquid-gas interface and the sharp-interface ghost fluid method to stably handle high liquid-air density ratio. Adaptive mesh refinement and Lagrangian droplet models are shown to be viable options for computational cost reduction. Moreover, high performance computing is leveraged to manage the computational cost. The experiment selected for validation eliminates the impact of inlet liquid and gas turbulence and focuses on the impact of the crossflow aerodynamic forces on the atomization physics. Validation is demonstrated by comparing column surface wavelengths, deformation, breakup locations, column trajectories and droplet sizes, velocities, and mass rates for a range of intermediate Weber numbers. Analysis of the physics is performed in terms of the instability and breakup characteristics and the features of downstream flow recirculation, and vortex shedding. Formation of “Λ” shape windward column waves is observed and explained by the combined upward and lateral surface motion. The existence of Rayleigh-Taylor instability as the primary mechanism for the windward column waves is verified for this case by comparing wavelengths from the simulations to those predicted by linear stability analyses. Physical arguments are employed to postulate that the type of instability manifested may be related to conditions such as the gas Weber number and the inlet turbulence level. The decreased column wavelength with increasing Weber number is found to cause enhanced surface stripping and early depletion of liquid core at higher Weber number. A peculiar “three-streak-two-membrane” liquid structure is identified at the lowest Weber number and explained as the consequence of the symmetric recirculation zones behind the jet column. It is found that the vortical flow downstream of the liquid column resembles a von Karman vortex street and that the coupling between the gas flow and droplet transport is weak for the conditions explored.

42 ENGINEERING↗