SEARCH · Engineering Papers
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.
An overview is given of compressible Navier Stokes codes and up-wind solvers. Code validation efforts, inlet analysis requirements, and computational fluid dynamics modeling uncertainties are covered. Information is given in viewgraph form.
The computational technique of automatic differentiation (AD) is applied to a three-dimensional thin-layer Navier-Stokes multigrid flow solver to assess the feasibility and computational impact of obtaining exact sensitivity derivatives typical of those needed for sensitivity analyses. Calculations are performed for an ONERA M6 wing in transonic flow with both the Baldwin-Lomax and Johnson-King turbulence models. The wing lift, drag, and pitching moment coefficients are differentiated with respect to two different groups of input parameters. The first group consists of the second- and fourth-order damping coefficients of the computational algorithm, whereas the second group consists of two parameters in the viscous turbulent flow physics modelling. Results obtained via AD are compared, for both accuracy and computational efficiency with the results obtained with divided differences (DD). The AD results are accurate, extremely simple to obtain, and show significant computational advantage over those obtained by DD for some cases.
Three examples of response surface modeling with CFD are presented for combined cycle propulsion components. The examples include a mixed-compression-inlet during hypersonic flight, a hydrogen-fueled scramjet combustor during hypersonic flight, and a ducted-rocket nozzle during all-rocket flight. Three different experimental strategies were examined, including full factorial, fractionated central-composite, and D-optimal with embedded Plackett-Burman designs. The response variables have been confined to integral data extracted from multidimensional CFD results. Careful attention to uncertainty assessment and modeling bias has been addressed. The importance of automating experimental setup and effectively communicating statistical results are emphasized.
Explore the source record for details and available documents.
The Descent Interpolated Gas-Granular Erosion Model (DIGGEM) is a hybrid computational fluid dynamics (CFD) / engineering model calibrated with flight data from Apollo. It is used for the prediction of plume induced viscous erosion and cratering during the descent of landing vehicles. The model assumes a functional relationship between local shear stress and eroded mass flux. It is used along with a suite of other tools to predict complex plume induced environments for Human Landing System (HLS) and Commercial Lunar Payload Services (CLPSS) vehicles.
Explore the source record for details and available documents.
Abstract not provided.
A Zero-D cycle simulation of the GE90-94B high bypass turbofan engine has been achieved utilizing mini-maps generated from a high-fidelity simulation. The simulation utilizes the Numerical Propulsion System Simulation (NPSS) thermodynamic cycle modeling system coupled to a high-fidelity full-engine model represented by a set of coupled 3D computational fluid dynamic (CFD) component models. Boundary conditions from the balanced, steady state cycle model are used to define component boundary conditions in the full-engine model. Operating characteristics of the 3D component models are integrated into the cycle model via partial performance maps generated from the CFD flow solutions using one-dimensional mean line turbomachinery programs. This paper highlights the generation of the high-pressure compressor, booster, and fan partial performance maps, as well as turbine maps for the high pressure and low pressure turbine. These are actually "mini-maps" in the sense that they are developed only for a narrow operating range of the component. Results are compared between actual cycle data at a take-off condition and the comparable condition utilizing these mini-maps. The mini-maps are also presented with comparison to actual component data where possible.
MFIX-Exa is a computational fluid dynamics–discrete element model (CFD-DEM) code designed to run efficiently on current and next-generation supercomputing architectures. MFIX-Exa combines the CFD-DEM expertise embodied in the MFIX code—which was developed at NETL and is used widely in academia and industry—with the modern software framework, AMReX, developed at LBNL. The fundamental physics models follow those of the original MFIX, but the combination of new algorithmic approaches and a new software infrastructure will enable MFIX-Exa to leverage future exascale machines to optimize the modeling and design of multiphase chemical reactors.
Research on the computational fluid dynamics (CFD) code evaluation with emphasis on supercomputing in reacting flows is discussed. Advantages of unstructured grids, multigrids, adaptive methods, improved flow solvers, vector processing, parallel processing, and reduction of memory requirements are discussed. As examples, researchers include applications of supercomputing to reacting flow Navier-Stokes equations including shock waves and turbulence and combustion instability problems associated with solid and liquid propellants. Evaluation of codes developed by other organizations are not included. Instead, the basic criteria for accuracy and efficiency have been established, and some applications on rocket combustion have been made. Research toward an ultimate goal, the most accurate and efficient CFD code, is in progress and will continue for years to come.
Two time-accurate Computational Fluid Dynamics (CFD) codes were used to compute several flutter points for a typical business jet model. The model consisted of a rigid fuselage with a flexible semispan wing and was tested in the Transonic Dynamics Tunnel at NASA Langley Research Center where experimental flutter data were obtained from M(sub infinity) = 0.628 to M(sub infinity) = 0.888. The computational results were computed using CFD codes based on the inviscid TSD equation (CAP-TSD) and the Euler/Navier-Stokes equations (CFL3D-AE). Comparisons are made between analytical results and with experiment where appropriate. The results presented here show that the Navier-Stokes method is required near the transonic dip due to the strong viscous effects while the TSD and Euler methods used here provide good results at the lower Mach numbers.
Full three-dimensional Navier-Stokes computational results are compared with new experimental measurements for the flowfield within a round diffusing S-duct. The present study extends previous computational and experimental results for a similar smaller scale S-duct. Predicted results are compared with the experimental static and total pressure fields, and velocity vectors. Additionally, wall pressures, velocity profiles in wall coordinates, and skin friction values are presented. The CFD results employ algebraic and k-epsilon turbulence models. The CFD computed and experimentally determined separated flowfield is carefully examined.
Full three-dimensional Navier-Stokes computational results are compared with new experimental measurements for the flowfield within a round diffusing S-duct. The present study extends previous computational and experimental results for a similar smaller scale S-duct. Predicted results are compared with the experimental static and total pressure fields, and velocity vectors. Additionally, wall pressures, velocity profiles in wall coordinates, and skin friction values are presented. The CFD results employ algebraic and k-epsilon turbulence models. The CFD computed and experimentally determined separated flowfield is carefully examined.
This report updates the high-level test plan for evaluating surface deposition on three commercial 32PTH2 spent nuclear fuel (SNF) canisters inside NUTECH Horizontal Modular Storage (NUHOMS) Advanced Horizontal Storage Modules (AHSMs) from Orano (formerly Transnuclear Inc.) and provides a summary of the surface sampling activities that have been conducted to date. The details contained in this report represent the best designs and approaches explored for testing as of this publication. Given the rapidly developing nature of this test program, some of these plans may change to accommodate new objectives or requirements. One goal of this testing is to collect defensible and detailed dust deposition measurements from the surface of dry storage canisters in a marine coastal environment to guide chloride-induced stress corrosion cracking (CISCC) research. Another goal is to provide data for the validation of computational fluid dynamics (CFD) based deposition modeling. To facilitate surface sampling, the otherwise highly prototypic dry storage systems will not contain SNF but rather will be electrically heated to mimic the decay heat and thermal hydraulic environment. Test and heater design is supported by detailed CFD modeling. Instrumentation throughout the canister, storage module, and environment will provide extensive information about the thermal-hydraulic behavior of horizontal dry cask storage systems. Manual sampling over a comprehensive portion of the canister surface at regular time intervals will offer detailed quantification and composition of the deposited particulates from a realistic storage environment. Discussions of a potential host site for the Canister Deposition Field Demonstration (CDFD) are ongoing. Until a host site is chosen, testing of key CDFD hardware components including the heater assemblies, power skid, and remote data acquisition system will continue. Functional testing of the finalized heater assemblies and test apparatus started this fiscal year. These initial heater tests have shown the assemblies are performing within design specifications. Staged surface sampling of a mockup of a canister outside the AHSM on a transfer skid was also performed. Refinements to the sampling procedures and techniques were captured from observation of these activities and lessons-learned debriefs. These updated sampling procedures and techniques are planned to be tested again in the field using the mockup in order to assure personnel are using the most accurate and repeatable methods possible prior to deployment for actual CDFD testing.
The Dragonfly Mass Spectrometer (DraMS) is an instrument on the Dragonfly mission, which will spend 7 years in deep space cruise before landing and operating on the surface of Titan. Vacuum thermal analyses are required for deep space cruise, and convection analyses are required for the Titan surface operations. Model exchanges across multiple thermal teams are needed for all phases of the mission. For DraMS, Thermal Desktop® (TD) has been the main thermal analytical tool of choice due to its capability in modeling complex thermal systems with relatively low computational power and for its availability across thermal teams. However, TD does not have computational fluid dynamics (CFD) capability and struggles to accurately capture complex convective behavior. DraMS has fans operating in tandem and gas flow behaviors are not easily predicted due to its complex flow paths. CFD software, such as Fluent, can model and predict such complex flow behaviors, but CFD models are computationally expensive, and its workflow processes are not tailored towards simulating large and complex systems. Therefore, a coupled modeling approach was chosen for DraMS: A TD model was used for simulating all the conductive, radiative, and source terms, while a Fluent CFD model was added on, as needed, to the TD model to provide the convective boundary conditions using the System Coupling software. The coupling software allows the TD and Fluent models to communicate data and arrive at a co-solved and co-converged solution. Furthermore, Thermal Iso-value Exchange (TIE) method was developed to facilitate and improve the TD-Fluent data exchange process. This paper will discuss the analytical studies that were done to verify the accuracy and usability of the coupled approach and the challenges associated, which lead to the development of the TIE approach. DraMS thermal design and co-solved analysis results will also be discussed.
The objective of this study was to demonstrate the high-fidelity numerical simulation of a modern high-bypass turbofan engine. The simulation utilizes the Numerical Propulsion System Simulation (NPSS) thermodynamic cycle modeling system coupled to a high-fidelity full-engine model represented by a set of coupled three-dimensional computational fluid dynamic (CFD) component models. Boundary conditions from the balanced, steady-state cycle model are used to define component boundary conditions in the full-engine model. Operating characteristics of the three-dimensional component models are integrated into the cycle model via partial performance maps generated automatically from the CFD flow solutions using one-dimensional meanline turbomachinery programs. This paper reports on the progress made towards the full-engine simulation of the GE90-94B engine, highlighting the generation of the high-pressure compressor partial performance map. The ongoing work will provide a system to evaluate the steady and unsteady aerodynamic and mechanical interactions between engine components at design and off-design operating conditions.
Three-dimensional (3-D) Computational Fluid Dynamics (CFD) analysis of a whole nuclear reactor core is a tremendous challenge due to the large geometric volume and complex structures. Here, this research presents a hybrid porous (HP) model to simplify a prismatic High Temperature Gas-cooled Reactor (HTGR) core, so 3-D CFD investigation can be performed on a full reactor core scale. In the HP model, the prototypic small coolant channels in the nuclear fuel blocks are lumped together to form multiple equivalent large coolant channels, and then the porous medium flow model is applied to each of them. Therefore, heat transfer in fuel blocks is computed by a hybrid combination of solid energy and porous flow energy equations. The similarity between the HP model and prototypic model is achieved by deriving the porous flow permeability, inertial resistance factor, and artificial thermophysical properties. Compared with the widely used whole porous (WP) flow model, the HP model preserves more realistic geometric structures, and therefore more accurate physical processes. The General Atomics' Modular High Temperature Gas-cooled Reactor (MHTGR) design was chosen as a prototype to demonstrate the methodology. Simulations were performed using the prototypic CFD model and HP model at steady-state forced circulation, steady-state natural circulation, and transient conditions that correspond to normal operation, extended period of pressurized cool down, and short-term transients after reactor shutdown, respectively. The comparison shows good agreement between the HP model and prototypic model in the maximum fuel temperature, average solid temperature, and helium flow rate, which demonstrates the potential applicability of the HP model for a full reactor core scale simulation in the future. As a benefit, the HP model reduces the mesh quantity by a factor of 50 from a prototypic model. Correspondingly, the computation time was reduced by a factor of at least 30.
The System Thermal Hydraulics (STH) code SAM has been coupled to the Computational Fluid Dynamics (CFD) code Simcenter STAR-CCM+ utilizing a hybrid domain overlapping method with an explicit coupling in time. The coupling aims to extend the STH code’s applicability to scenarios where local momentum and energy transfers are important yet difficult for STH codes to capture, such as three-dimensional (3D) mixing. The coupling method’s numerical stability was verified in the past against two closed-loop configurations, and it was validated against a double T-junction experiment with 3D scalar mixing. In the present work, the coupling method is validated against the TALL-3D STH/CFD coupling benchmark facility. TALL-3D is a three-legged, liquid-metal facility with a large, pool-type enclosure (test section) that exhibits 3D flow effects to be modeled by a CFD code. The rest of the system exhibits approximately 1D behavior well-predicted by an STH code. First, the present STAR-CCM+ CFD model of the 3D test section is validated against experimental data. Then, the SAM-STARCCM+ coupled model is validated against six different TALL-3D steady states, including SAM standalone model results for comparison. Lastly, the SAM-STARCCM+ coupled model is validated against two TALL-3D transients, one exhibiting flow reversal in the test section and one exhibiting nonlinear, Limit Cycle Oscillations (LCO). For the first transient, the SAM-STARCCM+ coupled model properly predicts an increase in the test section’s inlet temperature during flow reversal, and this is not predicted by the SAM standalone model. Following flow reversal, the SAM-STARCCM+ coupled model better-predicts the initial flow recovery and following oscillations as the system approaches a final natural circulation state. For the second transient, no true final steady state is observed due to LCO. Neither the SAM-STARCCM+ coupled model nor the SAM standalone model can perfectly capture the experiment’s changing oscillation frequency during the transient. However, the SAM-STARCCM+ coupled model does reproduce the oscillatory feedback observed in the system. This is a significant achievement as the SAM-STARCCM+ coupled model only uses an explicit coupling in time, as opposed to a semi-implicit coupling. In comparison, previous STH/CFD coupling efforts of the TALL-3D facility required semi-implicit coupling to obtain similar results.