Potential Approaches for Mesh Adaptation of Large Eddy Simulations with Near-Wall Treatments
Explore the source record for details and available documents.
CFD↗
SEARCH · Search NASA
Results for “CFD Modeling”
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.
Potential Approaches for Mesh Adaptation of Large Eddy Simulations with Near-Wall Treatments
Explore the source record for details and available documents.
CFD↗
Exploring the Langtry-Menter Transition Model for High Speed Applications Using FUN3D
A series of Reynolds-averaged Navier-Stokes (RANS) simulations were performed using the FUN3D flow solver to explore the capabilities of the Langtry-Menter Shear-Stress Transport (LM-SST) transition model for predicting transition for aircraft inlet applications. Two geometries were simulated: a zero-pressure-gradient flat plate and an axisymmetric cone exposed to hypersonic flow. In addition to the transition-sensitized LM-SST model investigations, simulations were run with the one-equation Spalart-Allmaras (SA) and the two-equation Menter Shear-Stress Transport (SST-V) RANS models in fully turbulent mode to identify the natural RANS model transition behavior as a function of Mach number when executed in fully turbulent mode. The flat plate simulations showed that (1) the transition model was able to predict rapid transition at a freestream Mach number of 0.2, which is expected but (2) the predicted transition location moved downstream as the freestream Mach number was increased for the simulations that used the SST-V turbulence model. The latter is significant as it is usually assumed that one- and two-equation turbulence models will produce fully turbulent flow very near the boundary layer origin. The flat plate simulation freestream Mach number trend was confirmed with simulations using the Wind-US code, which also saw a similar trend when employing the SA turbulence model. For the axisymmetric cone simulations, the transition location was highly sensitive to the inflow turbulence levels. This is significant as the prediction of the transition location is crucial when trying to predict inlet performance, especially for hypersonic vehicle applications. It was also noted that the predicted transition location for the cone when using the SST-V turbulence model agreed well with the predicted transition location from the equivalent zero-pressure-gradient flat plate case.
Transition Model↗
Unsteady CFD Simulations of a Compression Corner Geometry Using Wall-Modeled LES Methods in Loci/CHEM
As computing technology continues to improve, simulating unsteady aerodynamic environments on launch vehicle geometries becomes more viable. For NASA’s Space Launch System, the flowfield is characterized by large regions of unsteady turbulent wall-bounded flow. Simulating these unsteady wall-bounded environments is of interest to the SLS Aerodynamics team. Wall-modeled LES methods are currently the best approach for simulating these environments in engineering applications. Less-expensive computationally than fully-resolved methods while maintaining many of the benefits in simulating unsteady environments. Before applying these methods blindly to large geometries, best practices can first be applied and understood on smaller canonical problems.
Bryson Frank↗
A PIC Bootstrapping Strategy for Exascale CFD-DEM Simulation
The exascale computing era comes with the release of MFIX-Exa, a new code for studying gas-particle fluidization using CFD-DEM and PIC models on massively parallel, heterogeneous high-performance computing architectures. However, as compute resources continue to grow in scale and complexity, so too does the impetus to use them efficiently. Here, we propose a novel bootstrapping method to minimize neglected simulation time and cast the approach more broadly among a growing field of multi-fidelity methods.
Porcu, Roberto↗
Benchmarking a Pronghorn Model of the Natural Convection Shutdown Heat Removal Test Facility Against STAR-CCM+ (Poster)
Summary Verification of Pronghorn’s Net Radiation Transfer via Code-to-Code comparison with STAR-CCM+. Comparison based on Natural Convection Shutdown Heat-Removal Test Facility (NSTF). Near-perfect agreement if conjugate heat transfer inside the riser channels is replaced by fixed temperature. Difference of up to 5K if conjugate heat transfer is modeled (correlation vs. full CFD).
97 MATHEMATICS AND COMPUTING↗
Method for designing a combustion system with reduced environmentally-harmful emissions
A method for designing a combustion system which emits less of at least one environmentally-harmful emission is presented. In a describing step, an injector which introduces a fuel into a combustion chamber is described via a CFD code. In a modeling step, combustion kinetics of the fuel are modeled via a pre-processing code as the fuel mixes and reacts with an oxidizer. In a first selecting step, at least one primary scalar is derived during the modeling of the combustion kinetics. In a performing step, a table look-up is performed to obtain at least one data from a look-up database based on the primary scalar. In a second selecting step, at least one secondary scalar is selected in addition to the primary scalar(s). In a specifying step, at least one chemical pathway of formation or destruction for the secondary scalar is specified via a chemistry manager wherein the secondary scalar is representative of the environmentally-harmful emission(s) of the chemical pathway(s). In a utilizing step, the data is utilized to evaluate the chemical pathway(s) to quantify the environmentally-harmful emission(s). In an identifying step, an improvement to the combustion system is identified which reduces the environmentally-harmful emission(s).
Zambon, Andrea C.↗
Thermochemical Nonequilibrium Modeling in a Continuous-Galerkin, Finite-Element Framework
The presented work discusses the implementation, verification, and validation of Park's two-temperature model in a scalable, computational fluid dynamics (CFD) code developed at the US Department of Energy's Oak Ridge National Laboratory (ORNL). The implementation of Park's two-temperature model was verified through 0D test cases involving an adiabatic reactor and a nitrogen thermal bath. The implementation was then validated through comparisons with other validated CFD codes and experimental data on a hypersonic cylinder and double cones. These are standard benchmark test cases for thermochemical non-equilibrium (TCNE) modeling, and all data are shared publicly. The verification and validation results showed that ORNL's in-house CFD code could model complex, high-speed flow problems with and without TCNE modeling. This work is essential for future research involving 3D shock wave/boundary layer interactions (SBLIs).
Nutter, Nicole↗
Implementation of Detailed Polyethylene Pyrolysis Kinetics into CFD Simulations using Machine Learning
Municipal solid waste (MSW) and waste plastics have received significant attention due to the issues of waste generation and storage, as well as their potential as an energy resource. High-density polyethylene (HDPE) makes up a large portion of plastic waste and has been the subject of several conversion studies. However, the mechanisms associated with converting HDPE through pyrolysis and gasification are extensive and complex making them difficult to implement into high-fidelity computational fluid dynamic (CFD) simulations. For this project, a primary pyrolysis mechanism containing 42 unique species and 737 heterogeneous reactions was used to generate kinetic data over a range of operating conditions. A machine learning (ML) model was developed to replicate the results of the detailed pyrolysis mechanism while significantly increasing the computational efficiency. A deep operator network (DeepONet) architecture was adopted to train the model using time steps relevant to CFD simulations. The ML used physics-based loss functions to ensure mass conservation. The ML model has been deployed in simple MFiX CFD simulations, single particle, and an experimental drop tube reactor, and has shown promising performance compared to the original scheme.
Houston, Ross↗
Validation of NEAMS Tools Using MSRE Data
This report documents the benchmarks developed and performed for the tools and codes developed under the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program utilizing the MSRE experiment data. There are three focus areas explored in this work: fuel evolution, CFD, and system transient analyses. Fuel depletion is important for MSR reactor physics and fuel cycle analysis. Fuel depletion capability is currently being developed in the Griffin reactor physics code. To support this, a benchmark problem is developed using the MSRE data and an existing depletion code with flowing fuel analysis functionality (ADDER) to validate this new capability in Griffin. The ADDER calculation is compared with the TRITON-MSR results as well as the MSRE experimental data. This benchmark can be directly used to verify and validate the Griffin fuel depletion capability in future years. For CFD simulation, the MSRE CAD model is obtained from the Copenhagen Atomics and a mesh is generated for the MSRE inlet plenum. The Nek5000 code will be employed to find the flow distribution in the MSRE inlet plenum based on the mesh generated. For system analysis, a simple benchmark problem for three types of MSR system transients is developed. The SAM and SPECTRA codes are used to solve the simple problem, and the results from the two codes agree well. The SAM code is also employed to obtain the frequency response of the MSRE system at the 5 and 8 MWt power levels. The results agree well with the experiment measurements. These verification and validation efforts gives SAM users more confidence in MSR system transient simulations using SAM. Lastly, the SAM-Griffin Multiphysics coupling work is initiated. A simple one-way coupling of evaluating the precursor drift effect on the MSRE k-eff is demonstrated. This preliminary coupling work serves as the basis for future more complex Multiphysics simulation and modelling.
46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
Evaluation and application of the Baldwin-Lomax turbulence model in two-dimensional, unsteady, compressible boundary layers with and without separation in engine inlets
The Baldwin-Lomax model is used in many CFD codes because it is quick and easy to implement. In this paper, we discuss implementing the Baldwin-Lomax turbulence model for both steady and unsteady compressible flows. In addition, these flows may be either separated or attached. In order to apply this turbulence model to flows which may be subjected to these conditions, certain modifications should be made to the original Baldwin-Lomax model. We discuss these modifications and determine whether the Baldwin-Lomax model is a viable turbulence model that produces reasonably accurate results for high speed flows that can be found in engine inlets.
Sakowski, Barbara↗
Evaluation and application of the Baldwin-Lomax turbulence model in two-dimensional, unsteady, compressible boundary layers with and without separation in engine inlets
The Baldwin-Lomax model is used in many CFD codes because it is quick and easy to implement. In this paper, we discuss implementing the Baldwin-Lomax turbulence model for both steady and unsteady compressible flows. In addition, these flows may be either separated or attached. In order to apply this turbulence model to flows which may be subjected to these conditions, certain modifications should be made to the original Baldwin-Lomax model. We discuss these modifications and determine whether the Baldwin-Lomax model is a viable turbulence model that produces reasonably accurate results for high speed flows that can be found in engine inlets.
Sakowski, Barbara↗
Development of a Computational Fluid Dynamics Model for a High-Speed Centrifugal Compressor
Computational fluid dynamics (CFD) has become widely used in the design and analysis of turbomachinery components such as centrifugal compressors. However, CFD is only a limited representation of experimental cases and struggles to model complex flows or can lack small model details to increase computation speed. To make advancements in compressor technology, designers need tools like CFD that can help them predict flow behavior in new designs. The High Efficiency Centrifugal Compressor (HECC) was designed by United Technologies Research Center (UTRC) for NASA to investigate the difficulties behind improving centrifugal compressor technology and to provide an open case to the turbomachinery research community. CFD analysis was completed by UTRC to aid in the design process. The pre-test study significantly overpredicted the compressor’s performance, and while post-test analysis was more representative, its results have yet to be publicly available. CFD technology has also since improved in recent years. In this study, a high-fidelity computational model for the vaneless configuration of HECC operating at its design point was developed and compared to an equivalent experimental case. This model incorporated GPU versions of mesh generation and solver codes from AeroDynamic Solutions, Inc. The development of such a model is significant since there has yet to be a published numerical simulation of HECC’s vaneless configuration. Additionally, the use of GPU codes decreases computation times for the model, which has allowed for the inclusion of impeller blade fillets. Initial results showed a general overprediction of compressor performance by the model. Predictions were more accurate for the impeller compared to predictions for the full compressor stage. Spanwise analyses showed that the model tended to follow flow behavior patterns exhibited in the experiment. Beyond this study, this model will be iterated upon and experimentally validated to further examine HECC performance.
centrifugal compressor↗
Recommended DSMC Collision Model Parameters for Planetary Entry
Hybrid simulation frameworks combining Computational Fluid Dynamics (CFD) and Direct Simulation Monte Carlo (DSMC) are frequently employed to efficiently perform high-fidelity simulations of environments consisting of both continuum and rarified flow. DSMC is a stochastic, particle-based method which solves the fundamental Boltzmann equation and is therefore necessary for high-Knudsen flow where continuum-based assumptions governing CFD break down. However, the DSMC methodology is generally computationally inefficient to model the continuum regime. In a CFD/DSMC hybrid approach, obtaining an accurate, high-fidelity solution hinges on the consistent treatment of transport properties and the thermo-chemical models employed within the two solvers. In principle, in regions where CFD and DSMC are both employed, the same gas mixture under the same conditions should have the same properties, regardless of simulation type. Observed differences should be due to non-equilibrium processes, rather than differences in physical models. The goal of this work is to provide a comprehensive DSMC transport database encompassing the vast majority of known gas species encountered during Earth or other planetary atmospheric entry. This goal is largely possible due to recently performed ab initio quantum chemistry calculations. Combined with other high-fidelity data, the planned database will consist of collision integral data for over 200 neutral and ionized species and over 20000 binary collisions. From these collision integrals, Nelder-Mead optimization is used to compute collision-specific Variable Soft Sphere (VSS) collision model parameters, fitted from 300 K to 20000 K. Initial comparisons of transport properties of relevant equilibrium gas mixtures show great agreement between CFD and DSMC-derived results. The completed database can be readily applied to model binary collisions of any gas mixture containing the included species over the specified temperature range, making it a valuable tool for future planetary probe modeling efforts. An example is shown below in Fig. 1. Equilibrium mixture transport properties for a 35-species mixture composed originally of 10% air and 90% pyrolysis species of a carbon-phenolic ablator material [4] are computed using both fitted VSS Parameters and the original CFD collision integral values. Deviations in computed properties between the two approaches are less than 5% for the entire temperature range.
M. R. Gosma↗
AdapChem
AdapChem software enables high efficiency, low computational cost, and enhanced accuracy on computational fluid dynamics (CFD) numerical simulations used for combustion studies. The software dynamically allocates smaller, reduced chemical models instead of the larger, full chemistry models to evolve the calculation while ensuring the same accuracy to be obtained for steady-state CFD reacting flow simulations. The software enables detailed chemical kinetic modeling in combustion CFD simulations. AdapChem adapts the reaction mechanism used in the CFD to the local reaction conditions. Instead of a single, comprehensive reaction mechanism throughout the computation, a dynamic distribution of smaller, reduced models is used to capture accurately the chemical kinetics at a fraction of the cost of the traditional single-mechanism approach.
Oluwole, Oluwayemisi O.↗
Computational Study on the Effect of Thermal Boundary Conditions and Axial Aspect Ratio on Catalytic Oxidative Coupling of Methane
Packed-bed reactor models with a large reactor-diameter-to-particle-diameter ratio (N > 10) are commonly used in catalysis research. Their use is underpinned by strong simplifying assumptions including unidirectional plug flow, negligible radial gradients, small axial diffusion, and isothermal operation. However, the complex exothermic–endothermic nature of reactions like oxidative coupling of methane (OCM) coupled with heat transfer from the catalyst surface and convective heat losses from the reactor wall provides an interesting case study to evaluate packed-bed assumptions and understand axial and radial temperature variations that may impact conversion and selectivity. This computational study investigates the effect of thermal boundary conditions on catalyst performance in laboratory OCM reactors (N = ~13) through computational fluid dynamics (CFD) simulation. Further, the 1-D model typically used for modeling plug flow reactors assumes negligible radial temperature and species concentration gradients. This assumption is evaluated for OCM reactors over a range of geometric axial aspect ratios, defined as the reactor length divided by the diameter, comparing the 1-D with 3-D CFD results. 1-D model predictions of catalyst performance start to deviate from the 3-D model predictions as the axial aspect ratio decreases, specifically as diameter increases at constant length. In the work, the 1-D model assumption is only found to be valid for axial aspect ratios higher than 1.6, a result considered to be relevant for all exothermic catalytic reactors.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
A reduced-order model of concentration polarization in reverse osmosis systems with feed spacers
Feed spacers in reverse osmosis systems generate complex fluid flows that limit computational fluid dynamics (CFD) simulations to small length and time scales. That limits our ability to simulate mineral scaling and other membrane fouling phenomena, which occur over longer length and time scales. Thus motivated, we develop a reduced model that replaces the CFD simulation of the velocity field with an analytical model that mimics spacers. This focuses the remaining numerical effort on simulating the advection–diffusion equation governing solute transport. We motivate and validate the model with CFD simulations and bench-scale experiments of spacer filaments in three different arrangements, including cases of unsteady vortex shedding. We show that the model produces a roughly 10,000-fold speedup compared to CFD, and accurately reproduces CFD predictions of not only the average and maximum concentrations, but also the local concentration distribution along the membrane. Here, we also demonstrate the model for simulating a feed channel with a length-to-height ratio of 200. The model provides a simple testbed for exploratory studies of multispecies transport, precipitation, and membrane fouling phenomena for which simulating spacers is often prohibitive.
Computational fluid dynamics↗
Application of NEAMS Codes to Capture MSR Phenomena
This report documents work completed in FY21 under the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program’s Molten Salt Reactor (MSR) Application Drivers activity at Argonne. The common focus was on identifying the modeling and simulation functional requirements for designing, licensing, and operating MSRs and applying those capabilities already developed in NEAMS tools to example problems of interest. The four main parts of this report each focuses on a specific area of simulation physics as it relates to MSR phenomena: fuel evolution, chemistry, computational fluid dynamics (CFD), and systems analysis. In terms of fuel evolution, which includes depletion, decay, on-line separations, and transmutation, the current state of computational capabilities for modeling this behavior in liquid-fueled molten salt reactor is discussed. Some of the functional requirements to accomplish the various applications of MSR fuel depletion modeling are highlighted, followed by a summary of recent approaches and code development activities. The chemistry functional requirements were discussed in the context of several applications of high importance for MSRs, such as corrosion, salt chemistry, and species transport. Each of these types of chemistry modeling have considerable impact on various aspects of reactor applications, including informing on reactor designs, improving operational efficiencies, analyzing safety and reactivity concerns, and estimating the mechanistic source term of the reactor. A brief overview is also provided on code development activities ongoing under NEAMS relevant to chemistry modeling of MSRs. In terms of CFD applications, the state-of-the-art spectral element code, Nek5000 was used to model the fluid dynamics within a full core of the Molten Salt Fast Reactor (MSFR) concept designed as part of the Euratom EVOL project. This concept was selected as the challenge problem because of its similar features to several U.S. industry concepts. The goal was to model some of the fast MSR design challenges, including potential large internal re-circulations, the need of accurate tracking of delayed neutron precursors (DNP), and the lack of relevant thermal-hydraulics models/correlations, etc. Therefore, a series of CFD models were created for the MSFR core cavity using a k – τ model two-equation model for the turbulence. These first full core results demonstrated that any potential recirculation zones could be properly identified with the current NEAMS CFD capabilities. The development of these models will also set the stage for future testing of Nek5000’s functionalities to model other MSR thermal-fluid phenomena. Lastly, the validation of SAM against experimental data from the Molten Salt Reactor Experiment, which started in FY19, continues with the inclusion of modeling the reactivity insertion tests. This involved using SAM and its point kinetics model for flowing fuel salt to recreate the time dependent power changes and response after positive reactivity insertions at the 1, 5, and 8 MWt power levels. Through these exercises, several code modifications were suggested to the SAM development team and accommodated to enable closer agreement with solid technical and physical justifications. These include adding a moderator reactivity feedback coefficient as an available input and modifying the solution approach for the point kinetics model. Additional SAM development suggestions for flowing fuel MSRs include adding the capability to allow the moderator power change proportionally with the reactor power and enabling specification of the power and DNP distributions separately.
21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗