Search NASA⌕ Search

SEARCH · Search NASA

Results for “Computational fluid dynamics code”

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,189 records · Page 66

Heating Analysis of the Lockheed Lifting Body

Computational fluid dynamic (CFD) analysis is performed on the Lockheed Lifting Body Single-Stage-to-Orbit vehicle to determine the heat transfer to the vehicle during its descent trajectory. Seven species, chemical nonequilibriurn computations using the GASP code will be completed at several trajectory points to assess the thermal protection requirements of the vehicle. Sophisticated surface boundary conditions including in-depth conduction, catalycity, and a variable temperature wall have been incorporated into the flow solver.

Palmer, Grant↗

Deploy threading in Nalu solver stack

The goal of the ExaWind project is to enable predictive simulations of wind farms composed of many MW-scale turbines situated in complex terrain. Predictive simulations will require computational fluid dynamics (CFD) simulations for which the mesh resolves the geometry of the turbines, and captures the rotation and large deflections of blades. Whereas such simulations for a single turbine are arguably petascale class, multi-turbine wind farm simulations will require exascale-class resources. The primary code in the ExaWind project is Nalu, which is an unstructured-grid solver for the acousticallyincompressible Navier-Stokes equations, and mass continuity is maintained through pressure projection. The model consists of the mass-continuity Poisson-type equation for pressure and a momentum equation for the velocity. For such modeling approaches, simulation times are dominated by linear-system setup and solution for the continuity and momentum systems. For the ExaWind challenge problem, the moving meshes greatly affect overall solver costs as re-initialization of matrices and re-computation of preconditioners is required at every time step In this Milestone, we examine the effect of threading on the solver stack performance against flat-MPI results obtained from previous milestones using Haswell performance data full-turbine simulations. Whereas the momentum equations are solved only with the Trilinos solvers, we investigate two algebraic-multigrid preconditioners for the continuity equations: Trilinos/Muelu and HYPRE/BoomerAMG. These two packages embody smoothed-aggregation and classical Ruge-Stiiben AMG methods, respectively. In our FY18 Q2 report, we described our efforts to improve setup and solve of the continuity equations under flat-MPI parallelism. While significant improvement was demonstrated in the solve phase, setup times remained larger than expected. Starting with the optimized settings described in the Q2 report, we explore here simulation performance where OpenMP threading is employed in the solver stack. For Trilinos, threading is acheived through the Kokkos abstraction where, whereas HYPRE/BoomerAMG employs straight OpenMP. We examined results for our mid-resolution baseline turbine simulation configuration (229M DOF). Simulations on 2048 Haswell cores explored the effect of decreasing the number of MPI ranks while increasing the number of threads. Both HYPRE and Trilinos exhibited similar overal solution times, and both showed dramatic increases in simulation time in the shift from MPI ranks to OpenMP threads. This increase is attributed to the large amount of work per MPI rank starting at the single-thread configuration. Decreasing MPI ranks, while increasing threads, may be increasing simulation time due to thread synchronization and start-up overhead contributing to the latency and serial time in the model. These result showed that an MPI+OpenMP parallel decomposition will be more effective as the amount per MPI rank computation per MPI rank decreases and the communication latency increases. This idea was demonstrated in a strong scaling study of our low-resolution baseline model (29M DOF) with the Trilinos-HYPRE configuration. While MPI-only results showed scaling improvement out to about 1536 cores, engaging threading carried scaling improvements out to 4128 cores — roughly 7000 DOF per core. This is an important result as improved strong scaling is needed for simulations to be executed over sufficiently long simulated durations (i.e., for many timesteps). In addition to threading work described above, the team examined solver-performance improvements by exploring communication-overhead in the HYPRE-GMRES implementation through a communicationoptimal- GMRE algorithm (CO-GMRES), and offloading compute-intensive solver actions to GPUs. To those ends, a HYPRE mini-app was allow us to easily test different solver approaches and HYPRE parameter settings without running the entire Nalu code. With GPU acceleration on the Summitdev supercomputer, a 20x speedup was achieved for the overall preconditioner and solver execution time for the mini-app. A study on Haswell processors showed that CO-GMRES provides benefits as one increases MPI ranks.

17 WIND ENERGY↗

Decrease time-to-solution through improved linear-system setup and solve

The goal of the ExaWind project is to enable predictive simulations of wind farms composed of many MW-scale turbines situated in complex terrain. Predictive simulations will require computational fluid dynamics (CFD) simulations for which the mesh resolves the geometry of the turbines, and captures the rotation and large deflections of blades. Whereas such simulations for a single turbine are arguably petascale class, multi-turbine wind farm simulations will require exascale-class resources. The primary code in the ExaWind project is Nalu, which is an unstructured-grid solver for the acoustically-incompressible Navier-Stokes equations, and mass continuity is maintained through pressure projection. The model consists of the mass-continuity Poisson-type equation for pressure and a momentum equation for the velocity. For such modeling approaches, simulation times are dominated by linear-system setup and solution for the continuity and momentum systems. For the ExaWind challenge problem, the moving meshes greatly affect overall solver costs as re-initialization of matrices and re-computation of preconditioners is required at every time step We describe in this report our efforts to decrease the setup and solution time for the mass-continuity Poisson system with respect to the benchmark timing results reported in FY18 Q1. In particular, we investigate improving and evaluating two types of algebraic multigrid (AMG) preconditioners: Classical Ruge-Stfiben AMG (C-AMG) and smoothed-aggregation AMG (SA-AMG), which are implemented in the Hypre and Trilinos/MueLu software stacks, respectively. Preconditioner performance was optimized through existing capabilities and settings.

17 WIND ENERGY↗

Deploy threading in Nalu solver stack

The goal of the ExaWind project is to enable predictive simulations of wind farms composed of many MW-scale turbines situated in complex terrain. Predictive simulations will require computational fluid dynamics (CFD) simulations for which the mesh resolves the geometry of the turbines, and captures the rotation and large deflections of blades. Whereas such simulations for a single turbine are arguably petascale class, multi-turbine wind farm simulations will require exascale-class resources. The primary code in the ExaWind project is Nalu, which is an unstructured-grid solver for the acoustically-incompressible Navier-Stokes equations, and mass continuity is maintained through pressure projection. The model consists of the mass-continuity Poisson-type equation for pressure and a momentum equation for the velocity. For such modeling approaches, simulation times are dominated by linear-system setup and solution for the continuity and momentum systems. For the ExaWind challenge problem, the moving meshes greatly affect overall solver costs as re-initialization of matrices and re-computation of preconditioners is required at every time step.

17 WIND ENERGY↗

Fuselage-Induced Velocity Model for LRTA, RTA, and HART Fuselages

In order to make the aerodynamic fuselage-rotor interference effects available to comprehensive rotor codes, a simple analytical model of the fuselage-induced velocities within the volume of rotor blade operation above the fuselage is developed here for the following bodies used in wind tunnel experiments: the Large Rotor Test Apparatus (LRTA), the Rotor Test Apparatus (RTA), and the Higher Harmonic Control Aeroacoustic Rotor Test (HART).While the first two are used in the National Full-Scale Aerodynamics Complex (NFAC) at NASA Ames, California, the third one is used by DLR in the large low-speed facility of the German-Dutch wind tunnel in the Netherlands. The fuselage-induced velocity model is based on parameter identification of isolated fuselage-induced velocity data (computed by means of computational fluid dynamics, CFD) and is intended to be generic enough to be used for real helicopter fuselages as well. The accuracies obtained in reproducing the CFD data show a remaining average error of less or equal 5 of the peak-to-peak induced velocity range, which is considered sufficient for comprehensive code analysis.

HART↗

Critical analysis of velocimetry methods for particulate flows from synthetic data

Particle tracking methods that extract high-fidelity particle velocity data from high speed video of particle laden flows is a common experimental technique applied to chemical processes. These measurements are used to better understand the motion of particles and fluids in complex systems and create data against which computational models are validated. However, the methods, codes, and experimental setups all have limitations. It is imperative that practitioners verify the methods and their implementation as well as understand the limitations of experimental setups. This work focuses on quantifying the visible depth of field in a high particle concentration fluidized bed. Following a precedent set by the particle imaging velocimetry community, a particle velocity field is manufactured using a computational fluid dynamics and discrete element method simulation. Photo realistic high-speed videos are rendered based on the simulated data using the three-dimensional creation software Blender. Particle velocities are extracted from the synthetic high-speed videos using three variants of Particle Tracking Velocimetry and Optical Flow Velocimetry methodologies. Here, the tracked results are then compared to the known solution, quantifying the error associated with the assumed visible depth. The results indicate that at depth of one particle diameter, all three particle tracking codes give accurate measurements, largely within 5%. However, the error increases when the full bed video measurements are compared to the known solution at one particle diameter, i.e., mimicking a validation study. Finally, for some statistics the constant depth assumption only increases the error slightly, for others significantly.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Design of a wing shape for study of hypersonic crossflow transition in flight

Computational fluid dynamics methods were used in the design of a wing shape for study of hypersonic crossflow transition in flight. The flight experiment is to be performed on the delta wing of the first stage of a Pegasus launch vehicle as a piggy-back experiment to support boundary-layer stability code development and validation. The design goal is to obtain crossflow-induced transition at 20-40 percent of the chord for a flight Mach number of approximately six. The present paper describes the design and analysis process utilized to obtain desired glove shape. A variety of schemes were used in the design, ranging from simple empirical crossflow correlations to three-dimensional Navier-Stokes codes in conjunction with linear stability/N-factor computations. The sensitivity to various parameters, such as trajectory variations, allowable wing thickness, leading-edge radius and surface temperature, is also discussed.

Godil, A.↗

Three-dimensional flow model development for thermal mixing and stratification modeling in reactor system transients analyses

Mixing, thermal-stratification, and mass transport phenomena in large pools or enclosures play major roles for the safety of reactor systems. Depending on the fidelity requirement and computational resources, various modeling methods, from the 0-D perfect mixing model to 3-D Computational Fluid Dynamics (CFD) models, are available. Each is associated with its own advantages and shortcomings. It is very desirable to develop an advanced and efficient thermal mixing and stratification modeling capability embedded in a modern system analysis code to improve the accuracy of reactor safety analyses and to reduce modeling uncertainties. An advanced system analysis tool, SAM, is being developed at Argonne National Laboratory for advanced non-LWR reactor safety analysis. While SAM is being developed as a system-level modeling and simulation tool, a reduced-order three-dimensional module is under development to model the mull-dimensional flow and thermal mixing and stratification in large enclosures of reactor systems. This paper provides an overview of the three-dimensional finite element flow model in SAM, including the governing equations, stabilization scheme, and solution methods. Additionally, several verification and validation tests are presented, including lid-driven cavity flow, natural convection inside a cavity, laminar flow in a channel of parallel plates. Here, based on the comparisons with the analytical solutions and experimental results, it is demonstrated that the developed 3-D fluid model can perform very well for a wide range of flow problems.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Development of Preliminary VERA-CS Crud Induced Localized Corrosion Modeling Capability

The objective of this milestone was to develop a crud-induced localized corrosion (CILC) screening capability in VERA-CS using existing codes and capabilities that had been developed in the Consortium for Advanced Simulation of Light Water Reactors (CASL) as part of the efforts to resolve the crud challenge problem. The milestone work focused on improving existing CASL software and implementing clad corrosion modeling into VERA-CS. Clad corrosion modeling was implemented by creating an application interface (API) in the Cicada software so that its clad heat transfer and corrosion models could be used in VERA-CS. Cicada was coupled to the subchannel code CTF for solving crud and corrosion growth. Additionally, a clad corrosion model was added directly to CTF to serve as an alternate approach. The rod thermal hydraulic reconstruction (ROTHCON) process, which is used for capturing computational fluid dynamics (CFD)-predicted rod-surface heat and turbulence behavior that can impact crud-induced power shift (CIPS) and CILC prediction, was also implemented into VERA-CS. Usability improvements were made to this feature, which was developed over several previous milestones, so that the VERA-CS user can link spacer grids in a VERA-CS model to an applicable pregenerated spacer grid data file if available. Several assessments were performed to demonstrate the proper functioning of the new capability. Sensitivity studies were performed to demonstrate the effect of crud lithium content, oxide thermal conductivity, Cicada mesh refinement, and ROTHCON coupling mesh refinement on predicted amounts of crud and oxide growth. A full-scale demonstration was performed by running Seabrook Cycle 5 with the CILC capability. Some general observations were made from the assessments. It was observed that predicted corrosion thickness is less than anticipated for CILC cases, so further work is needed to better understand and improve the corrosion modeling. Refining the coupling mesh in the ROTHCON process leads to increased corrosion and crud thickness, which can affect where crud and corrosion grow, how thick it gets, and how many rods experience crud and corrosion growth. Results indicate that applying the ROTHCON process and refining the coupling mesh leads to only a slight increase in total core maximum crud and corrosion thickness, but it leads to a significant increase in the number of rods with a maximum thickness that surpasses a set threshold value (i.e., more rods grow thicker corrosion and crud). Computational performance of the simulation is degraded by a factor of 5{6 with high coupling mesh refinement. The impact of the ROTHCON coupling mesh refinement on total core crud mass is less significant; however, using the ROTHCON process with CFD-generated data over the Yao{Hochreiter{Leech (YHL) model in CTF leads to significantly different amounts of total core crud mass. Future work on the CILC screening capability will include model improvement and validation, development and testing of a CILC-risk screening criteria, and performance improvements.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

High Fidelity CFD Simulations Supporting the KP-FHR

Kairos Power, LLC, is developing its version of the Fluoride-cooled High-temperature Reactor, the KP-FHR. The design uses a pebble bed core with fluoride salt as a coolant. The pebbles used in the KP-FHR have a diameter of 4 cm, with a shell fuel region where TRISO particles are embedded. A Pebble bed core design is adopted by several Gen IV reactors, They boast many benefits, such as fuel integrity, highly efficient heat transfer, and passive safety. However, it is challenging to accurately predict temperature and flow inside a pebble bed. Traditional approaches use the porous media model, which regards the pebble bed as a continuous medium, but with different temperature fields representing different levels, such as the fluid temperature, pebble surface temperature, and pebble center temperature. Empirical heat transfer correlations are adopted to calculate the heat transfer coefficient between different phases. However, empirical correlations are usually validated with experimental data, which usually lacks detail inside the pebble bed. The available experimental data is also generally at a high Reynolds number, which falls outside of the conditions of KP-FHR. Explicit computational fluid dynamics (CFD) simulations of randomly packed pebble beds have only become feasible recently. This is thanks to the rapid development of computational power and scalable algorithms. In this work, we used the Spectral Element Method (SEM) CFD code NekRS to simulate the randomly packed pebble bed in a cylindrical container. NekRS, which is the GPU variant of Nek5000, but refactored to utilize the computational power of GPUs using the OCCA library to run on hybrid architecture high performance computing systems. It was initially developed with the libParamunal library, but truncated and tuned for large-scale turbulence simulation. As a result, the SEM reaches higher precision with the same degrees of freedom by using a high-order Lagrange polynomial basis distributed on Gauss-Lobatto-Legendre quadrature inside each element, compared to lower-order methods, such the Finite Volume Method and Finite Element Method. The report is divided into five parts. We start with a general discussion of the pebble bed reactor, along with a specific investigation into the KP-FHR. The second part presents the numerical methodology. In the third part, we study a modular pebble bed with 1741 pebbles in a container of 7 pebble-diameter radius. Beyond LES simulations done by NekRS, we also leveraged the thermal radiation model in OpenFOAM to study heat transfer under no-forced-flow scenarios. Then, in the fourth part we simulated a pebble bed similar to the size of the Hermes Test Reactor. The total number of pebbles is in these simulations is 34,374. The container radius is 14 pebble-diameters. Finally, the report concludes in part five, with a discussion of future work.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

A Measurement and Simulation Based Methodology for Cache Performance Modeling and Tuning

We present a cache performance modeling methodology that facilitates the tuning of uniprocessor cache performance for applications executing on shared memory multiprocessors by accurately predicting the effects of source code level modifications. Measurements on a single processor are initially used for identifying parts of code where cache utilization improvements may significantly impact the overall performance. Cache simulation based on trace-driven techniques can be carried out without gathering detailed address traces. Minimal runtime information for modeling cache performance of a selected code block includes: base virtual addresses of arrays, virtual addresses of variables, and loop bounds for that code block. Rest of the information is obtained from the source code. We show that the cache performance predictions are as reliable as those obtained through trace-driven simulations. This technique is particularly helpful to the exploration of various "what-if' scenarios regarding the cache performance impact for alternative code structures. We explain and validate this methodology using a simple matrix-matrix multiplication program. We then apply this methodology to predict and tune the cache performance of two realistic scientific applications taken from the Computational Fluid Dynamics (CFD) domain.

Waheed, Abdul↗

Recent Enhancements to USM3D Unstructured Flow Solver for Unsteady Flows

The NASA USM3D unstructured flow solver is undergoing extensions to address dynamic flow problems in support of NASA and NAVAIR efforts to study the applicability of Computational Fluid Dynamics tools for the prediction of aircraft stability and control characteristics. The initial extensions reported herein include two second-order time stepping schemes, Detached-Eddy Simulation, and grid motion. This paper reports the initial code verification and validation assessment of the dynamic flow capabilities of USM3D. The cases considered are the classic inviscid shock-tube problem, low Reynolds number wake shedding from a NACA 0012 airfoil, high Reynolds number DES-based wake shedding from a 4-to-1 length-to-diameter cylinder, and forced pitch oscillation of a NACA 0012 airfoil with inviscid and turbulent flow.

Pandya, Mohagna J.↗

Development of a Reduced-Order Three-Dimensional Flow Model for Thermal Mixing and Stratification Simulation during Reactor Transients

Mixing, thermal-stratification, and mass transport phenomena in large pools or enclosures play major roles for the safety of reactor systems. Depending on the fidelity requirement and computational resources, various modeling methods, from the 0-D perfect mixing model to 3-D Computational Fluid Dynamics (CFD) models, are available. Each is associated with its own advantages and shortcomings. It is very desirable to develop an advanced and efficient thermal mixing and stratification modeling capability embedded in a modern system analysis code to improve the accuracy of reactor safety analyses and to reduce modeling uncertainties. An advanced system analysis tool, SAM, is being developed at Argonne National Laboratory for advanced non-LWR reactor safety analysis. While SAM is being developed as a system-level modeling and simulation tool, a reduced-order three-dimensional module is under development to model the multi-dimensional flow and thermal mixing and stratification in large enclosures of reactor systems. This paper provides an overview of the three-dimensional finite element flow model in SAM, including the governing equations, stabilization scheme, and solution methods. Additionally, several verification and validation tests are presented, including lid-driven cavity flow, natural convection inside a cavity, laminar flow in a channel of parallel plates. Based on the comparisons with the analytical solutions and experimental results, it is demonstrated that the developed 3-D fluid model can perform very well for a wide range of flow problems.

Hu, Rui↗

External-Compression Supersonic Inlet Design Code

A computer code named SUPIN has been developed to perform aerodynamic design and analysis of external-compression, supersonic inlets. The baseline set of inlets include axisymmetric pitot, two-dimensional single-duct, axisymmetric outward-turning, and two-dimensional bifurcated-duct inlets. The aerodynamic methods are based on low-fidelity analytical and numerical procedures. The geometric methods are based on planar geometry elements. SUPIN has three modes of operation: 1) generate the inlet geometry from a explicit set of geometry information, 2) size and design the inlet geometry and analyze the aerodynamic performance, and 3) compute the aerodynamic performance of a specified inlet geometry. The aerodynamic performance quantities includes inlet flow rates, total pressure recovery, and drag. The geometry output from SUPIN includes inlet dimensions, cross-sectional areas, coordinates of planar profiles, and surface grids suitable for input to grid generators for analysis by computational fluid dynamics (CFD) methods. The input data file for SUPIN and the output file from SUPIN are text (ASCII) files. The surface grid files are output as formatted Plot3D or stereolithography (STL) files. SUPIN executes in batch mode and is available as a Microsoft Windows executable and Fortran95 source code with a makefile for Linux.

Slater, John W.↗

Methodology for Thermodynamic Analysis Coupled with Computational Fluid Dynamics Modeling for Casting a Novel Aluminum–Cerium Alloy

In this paper, we present the results of computational fluid dynamics (CFD) analysis to assess castability, porosity, molten metal fluidity and other technological properties of a new Al-Ce alloy. A thermodynamic analysis of a new Al-Ce alloy is done to obtain the physical properties as a function of temperature across both solid and liquid phases. These properties are then used to build the CFD of model the full casting process from initial pouring through to final solidification of the part, here an example of a heavy duty mor mount is used. While several excellent commercial CFD codes exists this process shows that CFD can be used to assess the capability of the alloy to properly fill the mold, as well as give predictions where scattered porosity or large-scale defects may occur in the casting. Further, unlike the commercially available software (e.g., ProCAST, SOLIDCast , MAGMASOFT®) the complex 3D-analysis of the stress / strain fields in cast parts is not performed at this time. However, the availability of free software for assessing the required thermodynamic / thermo-physical properties of new alloys (OpenCALPHAD) and CFD codes such as OpenFOAM® makes the developed option attractive and economical, especially for the analysis of new Al-Ce alloys, for which the available data does not exist.

36 MATERIALS SCIENCE↗

Implicit solution methods in computational fluid dynamics

A class of implicit finite difference schemes of the Beam and Warming approximate factorization type will be addressed. The development and analysis of various aspects of this class of schemes will be given along with the motivations behind many of the choices. Various acceleration and efficiency modifications such as matrix reduction, diagonalization and flux split schemes will be presented. The methods are demonstrated in fully vectorized codes for a CRAY type architecture. The emphasis will be on the Euler equations in generalized coordinates.

Pulliam, Thomas H.↗

Data-driven model for cross ventilation potential in high-density cities based on coupled CFD simulation and machine learning

Effective urban ventilation through decent urban planning and building design can alleviate the deterioration of the urban built environment. Furthermore, natural ventilation requirements and guidelines in current building codes and standards are either qualitative or quantitative but subject to an absolute indoor airspeed threshold without considering the outdoor wind environment. To fill this gap, this paper develops an urban-scale coupled indoor and outdoor computational fluid dynamics (CFD) model and defines a novel ventilation index to assess natural ventilation potential. The index considers wind environments of both indoor and outdoor spaces.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

MHD simulation on magnetic compression of field reversed configurations with NIMROD

Abstract Magnetic compression has long been proposed a promising method for plasma heating in a field reversed configuration (FRC). However, it remains a challenge to fully understand the physical mechanisms underlying the compression process, due to its highly dynamic nature beyond the one-dimensional (1D) adiabatic theory model (Spencer et al 1983 Phys. Fluids 26 1564). In this work, magnetohydrodynamics simulations on the magnetic compression of FRCs using the NIMROD code (Sovinec et al 2004 J. Comput. Phys. 195 355) and their comparisons with the 1D theory have been performed. The effects of the assumptions of the theory on the compression process have been explored, and the detailed profiles of the FRC during compression have been investigated. The pressure evolution agrees with the theoretical prediction under various initial conditions. The axial contraction of the FRC can be affected by the initial density profile and the ramping rate of the compression magnetic field, but the theoretical predictions on the FRC’s length in general and the relation r s = 2 r o in particular hold approximately well during the whole compression process, where r s is the major radius of FRC separatrix and r o is that of the magnetic axis. The evolutions of the density and temperature can be affected significantly by the initial equilibrium profile and the ramping rate of the compression magnetic field. During the compression, the major radius of the FRC is another parameter that is susceptible to the ramping rate of the compression field. Basically, for the same magnetic compression ratio, the peak density is higher and the FRC’s radius r s is smaller than the theoretical predictions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗