Search NASA⌕ Search

SEARCH · Search NASA

Results for “Direct Numerical Simulation”

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 1,261 records · Page 70

Macrosegregation Resulting from Directional Solidification Through an Abrupt Change in Cross-Sections

Simulations of the directional solidification of two hypoeutectic alloys (Al-7Si alloy and Al-19Cu) and resulting macrosegregation patterns are presented. The casting geometries include abrupt changes in cross-section from a larger width of 9.5 mm to a narrower 3.2 mm width then through an expansion back to a width of 9.5 mm. The alloys were chosen as model alloys because they have similar solidification shrinkages, but the effect of Cu on changing the density of the liquid alloy is about an order of magnitude greater than that of Si. The simulations compare well with experimental castings that were directionally solidified in a graphite mold in a Bridgman furnace. In addition to the simulations of the directional solidification in graphite molds, some simulations were effected for solidification in an alumina mold. This study showed that the mold must be included in numerical simulations of directional solidification because of its effect on the temperature field and solidification. For the model alloys used for the study, the simulations clearly show the interaction of the convection field with the solidifying alloys to produce a macrosegregation pattern known as "steepling" in sections with a uniform width. Details of the complex convection- and segregation-patterns at both the contraction and expansion of the cross-sectional area are revealed by the computer simulations. The convection and solidification through the expansions suggest a possible mechanism for the formation of stray grains. The computer simulations and the experimental castings have been part of on-going ground-based research with the goal of providing necessary background for eventual experiments aboard the ISS. For casting practitioners, the results of the simulations demonstrate that computer simulations should be applied to reveal interactions between alloy solidification properties, solidification conditions, and mold geometries on macrosegregation. The simulations also presents the possibility of engineering the mold-material to avoid, or mitigate, the effects of thermosolutal convection and macrosegregation by selecting a mold material with suitable thermal properties, especially its thermal conductivity.

Lauer, M.↗

Direct simulation Monte Carlo and Navier-Stokes simulations of blunt body wake flows

Numerical results obtained with direct simulation Monte Carlo and Navier-Stokes methods are presented for a Mach-20 nitrogen flow about a 70-deg blunted cone. The flow conditions simulated are those that can be obtained in existing low-density hypersonic wind tunnels. Three sets of flow conditions are considered with freestream Knudsen numbers ranging from 0.03 to 0.001. The focus is on the wake structure: how the wake structure changes as a function of rare faction, what the afterbody levels of heating are, and to what limits the continuum models are realistic as rarefunction in the wake is progressively increased. Calculations are made with and without an afterbody sting. Results for the afterbody sting are emphasized in anticipation of an experimental study for the current flow conditions and model configuration. The Navier-Stokes calculations were made with and without slip boundary conditions. Comparisons of the results obtained with the two simulation methodologies are made for both flowfield structure and surface quantities.

WAKE STRUCTURE↗

Direct simulation Monte Carlo and Navier-Stokes simulations of blunt body wake flows

Numerical results obtained with direct simulation Monte Carlo and Navier-Stokes methods are presented for a Mach-20 nitrogen flow about a 70-deg blunted cone. The flow conditions simuulated are those that can be obtained in existing low-density hypersonic wind tunnels. Three sets of flow conditions are considered with freestream Knudsen numbers ranging from 0.03 to 0.001. The focus is on the wake structure: how the wake structure changes as a function of rarefaction, what the afterbody levels of heating are, and to what limits the continuum models are realistic as rarefaction in the wake is progressively increased. Calculations are made with and without an afterbody sting. Results for the after body sting are emphasizes in anticipation of an experimental study for the current flow conditions and model configuration. The Navier-Stokes calculations were made with and without slip boundary conditions. Comparisons of the results obtained with the two simulation methodologies are made for both flowfield structure and surface quantities.

Moss, James N.↗

Numerical Simulation of Lean Blowout of Alternative Fuels in 7-element Lean Direct Injector

This research presents the result of numerically simulating 7-element swirl-venturi Lean Direct Injector (SV-LDI) lean blowout (LBO) experiments conducted at NASA Glenn Research Center in May of 2019. After simulating a cold flow case to confirm the pressure drop agrees well with the experiment, additional cases with two different fuels (an average jet fuel and a Gevo alcohol-to-jet fuel) from the National Jet Fuels Combustion Program (NJFCP) were computed to numerically determine the LBO condition. The procedure to approach the LBO follows the method used in the experiment where the air mass flow rate is gradually increased while the fuel supply is maintained. Transient history of global heat release rate as a function of air flow rate is presented, as well as temperature contours at different conditions to give a visual representation of the flame state. The Open National Combustion Code (OpenNCC) used in this research adopted reduced HyChem (Hybrid Chemistry) models along with k-LES turbulence model and a Lagrangian spray model that takes into account droplet internal temperature distribution affected by the shear force on the droplet surface. The transport equations and chemical reaction terms are integrated together to enhance conservation of chemical species that are especially important in the near LBO conditions. After showing the computed range of LBO agrees well with the experimental measurements, time averaged solutions of both fuels at their initial condition and at their limiting condition just before LBO are compared in detail to facilitate the understanding of LBO mechanism.

combustion↗

Effect of the volume fraction gradient on the phase interaction force model for disperse two-phase flows

In this work, the effects of the particle volume fraction gradient on fluid-particle interactions are studied. The phase interaction force is decomposed into three terms. For the first term, namely the symmetrized force density, we present theoretical reasoning and numerical evidence to assume that it is independent of the particle volume fraction gradient. The second term is the particle volume fraction gradient times a newly introduced diffusion stress. The third term is the divergence of the particle-fluid-particle (PFP) stress. If this assumption of independence of the particle volume fraction gradient for the first term can be verified, to the first order of the ratio of the mean distance between particles to the macroscopic lengthscale, all three terms can be studied and modeled in flows with uniform particle distributions. Models thus obtained are applicable to statistically inhomogeneous flows, with the second and third terms accounting for statistical inhomogeneity. To verify this assumption, numerical simulations of flows passing fixed arrays of particles are performed. Both uniform and nonuniform particle volume fractions are studied and compared for disperse multiphase flows with the particle Reynolds numbers ranging from 1 to 100, and particle volume fraction ranging from 1% to 26% in statistically steady states. It is found that the symmetrized force (first) term can be well approximated by the drag force obtained from studies of uniform flows. The diffusion stress is positive along the flow direction and negative in the directions perpendicular to the flow. In the case of moving particles, this stress could potentially cause particle aggregation in the flow direction and dispersion in the directions perpendicular to the flow. Finally, the diffusion stress is only important when there is a volume fraction gradient, while the PFP stress can be important in inhomogeneous flows with either nonuniform particle concentrations or nonuniform average relative velocities between the phases.

42 ENGINEERING↗

Chemical vapor infiltration of additively manufactured preforms: Pore-resolved simulations and experimental validation

The densification of additively manufactured porous preforms by chemical vapor infiltration (CVI) is studied using pore-resolved simulations and experiments. Experimentally, 3D printed silicon carbide (SiC) preforms are subject to CVI synthesis using methyltrichlorosilane (MTS) precursor to obtain high purity SiC/SiC composites. Optical images of the cross sections of the processed preforms are analyzed to obtain the spatial porosity distribution. The numerical method is based on a level set formulation to capture the spatial distribution and time evolution of the pore scale microstructural characteristics. The coupled transport and kinetic effects are represented using a dimensionless Thiele modulus. Simulations are initialized using representative synthetic preform geometries comprising of packed particles based on the size distribution of the powder used for 3D printing. The simulation results are validated against the experimental observations in terms of total density and the distribution of residual porosity. The densification characteristics, porosity classification, concentration profiles, and structure functions are analyzed as functions of processing temperature and Thiele modulus.

36 MATERIALS SCIENCE↗

Turbulent Simulation of the Dynamics of the Magnetotail

In situ observations indicate that the dynamical processes in the geoplasma environment generally entail localized intermittent processes and anomalous global transports. It was suggested by T. Chang that instead of considering the turbulence as a mixture of interacting waves, such type of patchy intermittency could be more easily understood in terms of the development, interaction, merging, preferential acceleration and evolution of coherent magnetic structures. In this three-year project, we have used direct numerical MHD simulations to study some aspects of the MHD dynamics in Chang's model. Our large-scale numerical calculations and simulations have been supplemented by and coordinated with theoretical studies conducted by Chang and his colleagues.

Wu, Cheng-Chin↗

Formation and evolution of protostellar accretion discs – II. From 3D simulation to a simple semi-analytic model of Class 0/I discs

Here we use a 3D radiative non-ideal magnetohydrodynamic simulation to investigate the formation and evolution of a young protostellar disc from a magnetized pre-stellar core. The simulation covers the first ${\sim }10\, {\rm kyr}$ after protostar formation and shows a massive, weakly magnetized disc with radius that initially grows and then saturates at ${\sim }30\, {\rm au}$. The disc is gravitationally unstable with prominent large-amplitude spiral arms. We use our simulation results and a series of physical arguments to construct a predictive and quantitative physical picture of Class 0/I protostellar disc evolution from several aspects, including (i) the angular-momentum redistribution in the disc, self-regulated by gravitational instability to make most of the disc marginally unstable; (ii) the thermal profile of the disc, well-approximated by a balance between radiative cooling and accretion heating; and (iii) the magnetic-field strength and magnetic-braking rate inside the disc, regulated by non-ideal magnetic diffusion. Using these physical insights, we build a simple 1D semi-analytic model of disc evolution. We show that this 1D model, when coupled to a computationally inexpensive simulation for the evolution of the surrounding pseudo-disc, can be used reliably to predict disc evolution in the Class 0/I phase. The predicted long-term evolution of disc size, which saturates at ${\sim }30\, {\rm au}$ and eventually shrinks, is consistent with a recent observational survey of Class 0/I discs. Such hierarchical modelling of disc evolution circumvents the computational difficulty of tracing disc evolution through Class 0/I phase with direct, numerically converged simulations.

79 ASTRONOMY AND ASTROPHYSICS↗

Analysis of Residual Acceleration Effects on Transport and Segregation During Directional Solidification of Tin-Bismuth in the MEPHISTO Furnace Facility

The objective of this work is to approach the problem of determining the transport conditions (and effects of residual acceleration) during the plane-front directional solidification of a tin-bismuth alloy under low gravity conditions. The work involves using a combination of 2- and 3-D numerical models, scaling analyses, ID models and the results of ground-based and low-gravity experiments. The latter are to be conducted during the MEPHISTO experiment scheduled for USMP-3 in early 1996. The models will be used to predict the response of the transport conditions and consequent solute segregation in directionally solidifying tin-bismuth melt. Real-time Seebeck voltage variations across a Sn-Bi melt during directional solidification in MEPHISTO on USMP-1 show a distinct variation which can be correlated with thruster firings. The Seebeck voltage measurement is related to the response of the instantaneous average melt composition at the melt-solid interface. This allows a direct comparison of numerical simulations with the Seebeck signals obtained on USMP-1. The effects of such accelerations on composition for a directionally solidifying Sn-Bi alloy have been simulated numerically. USMP-1 acceleration data was used to assist in our choice of acceleration magnitude and orientation. The results show good agreement with experimental observations. The USMP-3 experiments took place earlier this year (February 22 through March 6). There were several differences between the USMP-3 experiments as compared to USMP-1. Firstly a more concentrated alloy was solidified and, secondly, Primary Reaction Control System thruster burns were requested at particular times during four separate growth runs. This allowed us to monitor the response Seebeck response under well-characterized growth conditions. In addition, we carried out simulations during the experiment in order to interpret the Seebeck signal. Preliminary results are described here.

Alexander J. Iwan D.↗

Validation Data Acquisition in HTTF during PCC Events

A validation experimental campaign was conducted in an Integral Effect Test (IET) facility of High Temperature Gas Reactors (HTGR), the High-Temperature Test Facility (HTTF) at Oregon State University (OSU). The HTTF simulates Depressurized and Pressurized Conduction Cooldowns (DCC and PCC). This campaign required the development of a new laser spectroscopic diagnostic to measure velocity in the challenging conditions of the HTTF: low speed (~1 m/s) gas flows at HTGR prototypical temperature and 1/10th pressure. This was a collaborative effort between co-PIs at The George Washington University (GW), Oregon State University (OSU), and NASA Langley Research Center. The main accomplishments of this project include the record for dynamic range for velocimetry, highest accuracy obtained with this technique, successful deployment in an IET leading to new validation matrix for CFD. These are detailed below and in manuscript appended to this executive summary. For this project, we introduced a new class of laser spectroscopic diagnostics to Thermal-Hydraulics to measure velocity of gases; furthermore, the diagnostic was demonstrated in-situ in an IET during DCC events. In such scenarios, particles used in mainstream techniques, like Particle Image Velocimetry (PIV) are not appropriate as they settle down too rapidly and also contaminate the experimental facility. Molecular tracers stay mixed within the working gas and can seed the flow in a technique called Molecular Tagging Velocimetry (MTV). In MTV a molecular tracer is photo-dissociated by a first (write) laser into radicals or molecules. The pattern created by the write laser is then interrogated with planar laser-induced fluorescence (PLIF), the read pulse(s), which are recorded with a camera. The pattern is probed and matched at two times (interval or probe time, dt), resulting in a time-of-flight velocimetry technique. This project demonstrated a new application range for MTV in gases. MTV has been extensively used for high-speed gas dynamics: molecules do not suffer from inertia associated with particles in regions of high accelerations (such as shock waves). For high-speed flows, the probe time is on the order of 0.1-10 µs. In time-of-flight techniques, the velocity is reconstructed from the recorded displacement, dx, between time interval dt: U = dx/dt. To first order, the resolution of the pattern matching is typically a (fixed) fraction of a pixel; therefore the resolution of the velocity measurement increases with dx and, for a given flow velocity, dt. For the low-speed flows of interest here, this requires long probe time, and therefore finding tracers that are stable for a long time. After an extensive literature survey we identified 6 tracers that would be suitable. We also worked with Spectra Physics to design a custom laser system that is flexible enough to probe these various tracers. The final instrument is made of 4 lasers that are integrated into a compact and movable cart to facilitate the deployment of the instrument off-site. One unique aspect of the laser system is the tunable dye laser (needed to do PLIF of the radicals). By being pumped by two Nd:YAG lasers, this laser is able to provide two pulses with any probe time (dt) for wavelengths between 200 and 900 nm with a linewidth of 3 pm. Accommodating the dual-pulse capability required customizing the pump lasers and the dye laser, making this system unique. Additionally, a custom UV-intensified CCD camera capable of frame-straddling and exposure time as short as 10 ns was acquired. More than two-third of the instrument cost (~$200k) was covered by a cost sharing from the PI’s startup funds. Of the 6 potential tracers identified, two were tested extensively (H2O and N2O) at high-temperature and elevated pressure at GW prior to conducting the tests at OSU. For both tracers, an Ar-F laser (193 nm, 10 ns, 10 mJ/pulse) photo-dissociates trace amounts (0.1 – 5% molar fraction) of the gases into OH and NO. OH is probed at 281.905 nm and NO at 226.186 nm (4-30mJ, 10 ns pulse). NO is very stable in inert environment and we could reconstruct velocities precisely with probe delay times as long at 40 ms, compared to 2 ms for OH, resulting in increased precision. Therefore, for the tests in the HTTF, N2O was selected as the final tracer and its optimal concentration was determined to be 0.5% molar. In the GW laboratory, we obtained a velocity measurement precision of 4 mm/s at conditions mimicking the HTTF tests. We also demonstrated the benefit of molecular tracers over particles and the flexibility of our instrument by continuously measuring velocity from 1000 to 0.1 m/s in a blow down test. This represents a dynamic range on velocity of 80 dB, which is unprecedented for velocimetry! Once the diagnostic was optimized, two two-week campaigns were conducted in the HTTF. This experimental campaign required numerous practical solutions to successfully and rapidly deploy this very sensitive table-top diagnostic in-situ. This was accomplished thanks to very extensive planning from all parties involved in the project. Due to constraints on availability of the HTTF, only adiabatic DCC tests were conducted. Specifically, the velocity was measured at the exit of the hot leg into the Reactor Cavity Simulation Tank (RCST). We were able to record the gas velocity associated with the air ingress for long times (up to one hour). Depending on the initial conditions (density of the gases in the Reactor Pressure Vessel and the RCST) the lock exchange lasted 30-60 s past the valve openings. Large vortices were identified within the shear layer at the interface of the cross flow. These vortices are an indication of shear-driven instability between the counter streams of different density. After the lock-exchange, a small residual flow was still present and its strength decayed exponentially. It finds its source in a buoyancy difference between the gases in both vessels due to the mixing and diffusion of the gases in the cross-over duct and the RPV lower plenum. We did not expect to capture this part of the transient, which would have been less significant (if not negligible) in the presence of a heated core; the precision of our measurement in the HTTF was 6 mm/s, which exceeded our expectations for such a complex test! Using commercial software (Fluent) we simulated the experiment with Computational Fluid Dynamic (CFD) using Reynolds Averaged Navier-Stokes (RANS) turbulence models. We were not able to accurately reproduce the shear layer instability, mixing, and diffusion captured in our experimental data, even with large 3D meshes (>10 M points). We attributed this to the inadequacy of the RANS models for the flow of interest: our experimental Reynolds number was too low and was outside the applicability range of the RANS models. Instead, we should have employed CFD with less modeling, such as Large Eddy (LES) or even Direct Numerical (DNS) simulations. Unfortunately, we did not have the computational resources to do so. Finally, it should be noted that this experimental technique holds great promises for measuring velocity in thermalhydraulics. MTV has an application domain (highlighted in the narrative) that is complimentary to more established (and simpler) techniques, such as PIV. It has attracted the attention of several thermalhydraulics researchers and sponsors to measure velocity in harsh environments. The spectroscopic nature of this technique also enables to expand it to measure temperature and concentration simultaneously to velocity. In lieu of a narrative, 4 detailed journal publications are appended to this executive summary to describe the various development and accomplishments during this project. They include 2 published manuscripts (Experiments in Fluids and Measurement Science and Technology), and two invited contributions following the NURETH-17 conference (Nuclear Science and Technology). Our current Nuclear Engineering and Design draft is not included at this time.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

GPU-accelerated DNS of compressible turbulent flows

Here, this paper explores strategies to transform an existing CPU-based high-performance computational fluid dynamics solver, HyPar, for compressible flow simulations on emerging exascale heterogeneous (CPU+GPU) computing platforms. The scientific motivation for developing a GPU-enhanced version of HyPar is to simulate canonical turbulent flows at the highest resolution possible on such platforms. We show that optimizing memory operations and thread blocks results in 200x speedup of computationally intensive kernels compared with a CPU core. Using multiple GPUs and CUDA-aware MPI communication, we demonstrate both strong and weak scaling of our GPU-based HyPar implementation on the NVIDIA Volta V100 GPUs. We simulate the decay of homogeneous isotropic turbulence in a triply periodic box on grids with up to 1024 3 points (5.3 billion degrees of freedom) and on up to 1,024 GPUs. We compare the wall times for CPU-only and CPU+GPU simulations. The results presented in the paper are obtained on the Summit and Lassen supercomputers at Oak Ridge and Lawrence Livermore National Laboratories, respectively.

97 MATHEMATICS AND COMPUTING↗

Three-dimensional transient flow of spin-up in a filled cylinder with oblique gravity force

Three-dimensional transient flow profiles of spin-up in a fully liquid filled cylinder from rest with gravity acceleration at various direction are numerically simulated and studied. Particular interests are concentrated on the development of temporary reverse flow zones and Ekman layer right after the impulsive start of spin-up from rest, and decay before the flow reaching to the solid rotation. Relationship of these flow developments and differences in the Reynolds numbers of the flow and its size selection of grid points concerning the numerical instabilities of flow computations are also discussed. In addition to the gravitational acceleration along the axial direction of the cylindrical container, a series of complicated flow profiles accompanied by three-dimensional transient flows with oblique gravitational acceleration has been studies.

Hung, R. J.↗

Shape Optimization of Vortex Generators to Control Mack Mode Amplification

This paper demonstrates the potential to use shape optimization for the design of vortex generators in an axisymmetric boundary layer. This shape optimization increases the amplitude of stationary streaks created by streamwise vorticity, with the goal of reducing the amplification of Mack mode instabilities. The test case under consideration matches a trajectory point during the ascent phase of the HIFiRE-1 flight experiment. Wall-mounted vortex generators are added to the cone and their shape is optimized to control the amplification of Mack mode instabilities that are known to initiate laminar-turbulent transition in this flow. An empirical objective function is developed based on the previous studies of optimal streaks. A constraint on the maximum streak amplitude is added to avoid the initiation of bypass transition. The shape optimization is conducted using SU2, an open-source suite for multiphysics simulation and design. A significant improvement is observed in an integral metric of the streak amplitude, while the maximum streak amplitude is maintained close to the baseline level. The qualitative features of the optimized geometry are discussed.

Vortex generators↗

3D micromechanical simulation of PBX composites

Previous research studied the constitutive response and interface strength of a bonded bicrystal system containing two HMX explosive grains bound with a HTPB polymer binder by performing FE simulations on real geometry imaged with micro Computed Tomography (µCT). The parameters generated from this past study were successfully applied to 2D mesoscale simulations of µCT imaged HMX-HTPB polycrystal samples. Presently, the mechanical response of PBXs containing a Nitro-Plasticized Estane (NPE) binder with similar microstructural geometry to the HMX-HTPB polycrystal samples were studied. Imaging of a HMX-NPE bicrystal sample allowed for simulation of the new material pairing. However, imaging using µCT techniques on the new polycrystal samples was difficult due to the decrease in contrast between the HMX grains and the NPE binder. Mesoscale polycrystal simulations utilizing a single µCT imaged microstructure were performed on a 3D representative volume element (RVE) comparing the delamination behavior and mechanical response of both HMX-HTPB and HMX-NPE PBXs. Furthermore, this technique shows the ability to virtually explore the mechanical response of a range of hypothetical materials that share common microstructural geometric characteristics in addition to studying numerically altered geometries.

36 MATERIALS SCIENCE↗

Reproducing curvature effects due to differential diffusion in tabulated chemistry for premixed flames

Tabulated chemistry reduces the cost of modeling reacting flows by transporting fewer scalars. It also has the benefit of introducing fewer terms that must be closed for Reynolds-averaged Navier-Stokes and large eddy simulation modeling. In this study of premixed lean hydrogen–air flames, an extension to an existing tabulated chemistry model was developed by including differential diffusion and thermal diffusion of multiple species, using a mixture-averaged diffusion model. The transport equations for the tabulation variables (a progress variable and a mixture fraction-like variable) have been formulated, and curvature effects are intrinsic to the transport equations derived from first principles. The mathematical derivation is general, such that other fuel/air mixtures could be examined. The new tabulated chemistry model was then validated using a range of configurations, including flames with curvature, extinction, and hot spots. The model validation indicated excellent agreement with detailed chemistry simulations for lean premixed hydrogen–air mixtures. Conditional means of the progress variable source term were shown to agree well with detailed chemistry simulations of two-dimensional freely propagating flames and three-dimensional turbulent flames. Finally, progress variable source terms three times higher than one-dimensional flat flame results were identified and validated with detailed chemistry, and super-adiabatic hot spots were also predicted in the turbulent flame configuration.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Tools for Design and Scale-Up of Solar Thermochemical Reactors: Cooperative Research and Development Final Report, CRADA Number CRD-13-00530

NREL will be collaborating with the Participant on a United States - Australia Solar Energy Collaboration (USASEC) Project Number 1-USO034 "Tools for design and scale-up of solar thermochemical reactors." The grant funds for the Participant's 3.5 year project number 1-US034 commencing on 1 February 2013 have been awarded to the Participant by the Australian Renewable Energy Agency and NREL will be collaborating with the Participant during the final 28 months of this project. This project seeks to provide basic knowledge required to design solar thermochemical reactors able to perform the required energy conversions. In several proposed and demonstrated reactors, concentrated sunlight directly irradiates small solid particles suspended in fluid, enabling very high heat transfer rates to the particles which are the sites of chemical reaction. The reactors, therefore, involve the complex and couple dynamics of turbulent, chemically reacting, particle-laden flows and their interaction with concentrating solar radiation. A strong understanding of these coupled interactions will be crucial important in predicting and optimizing the performance of prototype reactors, but this understanding does not yet exist, since they have never been studied in any fundamental way. The project has a assembled an internationally leading team from The University of New South Wales (UNSW) and the University of Adelaide in Australia and the NREL in the United States to address this key gap in available know-how. The project will use U.S. Dept. of Energy (DOE) supercomputers, among the most powerful available worldwide, with cutting-edge software tools to perform first-principles simulations of the relevant interactions. These studies will be combined with detailed laser-based measurements in Australia to provide the first comprehensive databases concerning the governing phenomena in directly irradiated solar-thermochemical reactors. The outcomes will be the basic scientific knowledge, engineering knowhow and modeling tools necessary to design new reactor concepts and then scale up from the laboratory bench to practical size systems.

14 SOLAR ENERGY↗

Low-Reynolds-number oscillating boundary layers on adiabatic slopes

We investigate the instabilities and transition mechanisms of Boussinesq stratified boundary layers on sloping boundaries when subjected to oscillatory body forcing parallel to the slope. We examine idealized forms of boundary layers on hydraulically smooth abyssal slopes in tranquil mid- to low-latitude regions, where low-wavenumber internal tides gently heave isopycnals up and down adiabatic slopes in the absence of mean flows, high-wavenumber internal tides, shelf breaks, resonant tide–bathymetry interactions (critical slopes) and other phenomena associated with turbulence ‘hot spots’. In non-rotating low-Reynolds-number flow, increased stratification on the downslope phase has a relaminarizing effect, while on the upslope phase we find transition-to-turbulence pathways arise from shear production triggered by gravitational instabilities. When rotation is significant (low slope Burger numbers) we find that boundary layer turbulence is sustained throughout the oscillation period, resembling stratified Stokes–Ekman layer turbulence. Simulation results suggest that oscillating boundary layers on smooth slopes at low Reynolds number ($Re\leqslant 840$), unity Prandtl number and slope Burger numbers greater than unity do not cause significant irreversible turbulent buoyancy flux (mixing), and that flat-bottom dissipation rate models derived from the tide amplitude are accurate within an order of magnitude.

42 ENGINEERING↗

Interface-Resolved Simulations of Reactor Flows

The presented review paper outlines the recent progress of the high-resolution single- and two-phase flow simulations of reactor-relevant flows. Rapid development of high-performance computing capabilities is creating exciting opportunities to study complex reactor thermal hydraulic phenomena. Today’s advances in thermal hydraulic analysis, machine learning techniques and interface resolved simulations will help pave the way to the next level of understanding of two-phase flow behavior in complex geometries. This paper consists of two major parts: (i) history and review of interface tracking simulations for nuclear thermal hydraulics in recent years and (ii) several opportunities to apply those advanced tools in the future. First part will discuss typical computational methods used for those simulations, provide some examples of the past work, as well as computational cost estimates and affordability of such simulations for research and industrial applications. In the second part some specific examples are discussed which could be analyzed using exascale supercomputers being designed and projected to be online in the early 2020s. New generation methodologies are required in order to take full advantage of those capabilities to greatly enhance the scientific understanding of complex two-phase flow phenomena in nuclear reactors under normal operation and postulated accident conditions.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗