Engineering Papers⌕ Search

SEARCH · Engineering Papers

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.

At least 973 records · Page 54

Using CFD Techniques to Predict Slosh Force Frequency and Damping Rate

Resonant effects and energy dissipation due to sloshing fuel inside propellant tanks are problems that arise in the initial design of any spacecraft or launch vehicle. A faster and more reliable method for calculating these effects during the design stages is needed. Using Computational Fluid Dynamics (CFD) techniques, a model of these fuel tanks can be created and used to predict important parameters such as resonant slosh frequency and damping rate. This initial study addresses the case of free surface slosh. Future studies will focus on creating models for tanks fitted with propellant management devices (PMD) such as diaphragms and baffles.

Marsell, Brandon↗

A Coupled Fluid-Structure Interaction Analysis of Solid Rocket Motor with Flexible Inhibitors

A capability to couple NASA production CFD code, Loci/CHEM, with CFDRC's structural finite element code, CoBi, has been developed. This paper summarizes the efforts in applying the installed coupling software to demonstrate/investigate fluid-structure interaction (FSI) between pressure wave and flexible inhibitor inside reusable solid rocket motor (RSRM). First a unified governing equation for both fluid and structure is presented, then an Eulerian-Lagrangian framework is described to satisfy the interfacial continuity requirements. The features of fluid solver, Loci/CHEM and structural solver, CoBi, are discussed before the coupling methodology of the solvers is described. The simulation uses production level CFD LES turbulence model with a grid resolution of 80 million cells. The flexible inhibitor is modeled with full 3D shell elements. Verifications against analytical solutions of structural model under steady uniform pressure condition and under dynamic condition of modal analysis show excellent agreements in terms of displacement distribution and eigen modal frequencies. The preliminary coupled result shows that due to acoustic coupling, the dynamics of one of the more flexible inhibitors shift from its first modal frequency to the first acoustic frequency of the solid rocket motor.

Yang, H. Q.↗

CFD Validation Experiment of a Mach 2.5 Axisymmetric Shock-Wave/Boundary-Layer Interaction

Experimental investigations of specific flow phenomena, e.g., Shock Wave Boundary-Layer Interactions (SWBLI), provide great insight to the flow behavior but often lack the necessary details to be useful as CFD validation experiments. Reasons include: 1.Undefined boundary conditions Inconsistent results 2.Undocumented 3D effects (CL only measurements) 3.Lack of uncertainty analysis While there are a number of good subsonic experimental investigations that are sufficiently documented to be considered test cases for CFD and turbulence model validation, the number of supersonic and hypersonic cases is much less. This was highlighted by Settles and Dodsons [1] comprehensive review of available supersonic and hypersonic experimental studies. In all, several hundred studies were considered for their database.Of these, over a hundred were subjected to rigorous acceptance criteria. Based on their criteria, only 19 (12 supersonic, 7 hypersonic) were considered of sufficient quality to be used for validation purposes. Aeschliman and Oberkampf [2] recognized the need to develop a specific methodology for experimental studies intended specifically for validation purposes.

Turbulence Modeling↗

Slosh-Induced Ullage Collapse Validation by Computational Fluid Dynamics and Application to In-space Propellant Dynamics

Ullage collapse predictions made using Computational Fluid Dynamics (CFD) were validated then made for an in-space cryogenic propellant unsettling application. Ullage collapse is the decrease in propellant tank pressure due to cooling of the ullage gas, which can be significant for cryogenic propulsion systems. Vehicles may be developed that are robust to temperature and pressure changes caused by heat leaks, e.g., propellant tanks with high maximum operation pressures, technologies yielding active thermal conditioning, and systems with ample pressurant for conditioning. These methods are robust because additional mass, cost, and complexity are used to address cryogenic fluid management (CFM). As vehicles are optimized for increased payload capacity, tracking of events such as ullage collapse and self-pressurization in tanks will become increasingly important for propellant inventory calculations. Some ullage collapse events can be studied through ground experimentation, but propellant unsettling in and around in-space coast mission phases must be studied in the absence of an acceleration field. Computational Fluid Dynamics is a powerful tool to address this need, enabling prediction of complex physics prior to flight. Validation of the chosen CFD tool and modeling methodology was done using data from a previously conducted subscale slosh-induced ullage collapse ground test performed in part by the Japan Aerospace Exploration Agency (JAXA). A small tank was partially filled with liquid nitrogen and pressurized with gaseous nitrogen at various temperatures. The tank was oscillated laterally to induce high amplitude slosh dynamics resulting in wave breakup. Successful qualitative validation of slosh wave dynamics and quantitative validation of ullage collapse promoted confidence in extending ullage collapse predictions to a flight application. Mechanisms of ullage collapse caused by main engine cutoff while in-space were identified. Both slosh dynamics leading to engine shutdown and body dynamics following shutdown were shown to contribute to propellant unsettling and ullage collapse.

J M Brodnick↗

CFD simulations of Molten Salt Fast Reactor core cavity flows

Computational Fluid Dynamics (CFD) has become increasingly important in the research and development of advanced nuclear reactors. Here, in the current study, extensive CFD simulations were conducted for the coolant flow in Molten Salt Fast Reactor (MSFR) core models using the state-of-the-art spectral element flow solver Nek5000 and multiscale coarse-mesh thermal-hydraulic software Pronghorn. The underlying motivation is to seek an in-depth understanding of how the internal velocity distribution can be influenced by the MSFR core cavity shape, the Reynolds number, turbulence modeling options and the inlet boundary conditions. The CFD techniques involved in this investigation range from coarse-mesh CFD, RANS modeling, to the high-fidelity LES calculations. Specifically, a series of RANS simulations were performed for the 2-D axisymmetric core model and 3-D wedge domains to study the flow distribution inside the MSFR core. It is observed that a proper representation of the MSFR inlet channel duct is important for the prediction of internal flow distribution. It is also showcased here how researchers can leverage the Nek5000 CFD results to calibrate more efficient coarse-mesh CFD tools, like Pronghorn, for the actual MSFR design needs. Moreover, this paper highlights a 3-D LES model for an entire MSFR core using the spectral element method and demonstrates the feasibility of this modeling approach. The readiness and potential limitations of the RANS approach are examined with respect to the high-fidelity LES simulations. The present investigation lays a solid foundation as we are leveraging the high-fidelity CFD capabilities to inform MSFR design efforts.

97 MATHEMATICS AND COMPUTING↗

Airfoil Computational Fluid Dynamics - 2k shapes, 25 AoA's, 3 Re numbers

This dataset contains aerodynamic quantities - including flow field values (momentum, energy, and vorticity) and summary values (coefficients of lift, drag, and momentum) - for 1,830 airfoil shapes computed using the HAM2D CFD (computational fluid dynamics) model. The airfoil shapes were designed using the separable shape tensor parameterization that encodes two-dimensional shapes as elements of the Grassmann manifold. This data-driven approach learns two independent spaces of parameter from a collection of sample airfoils. The first captures large-scale, linear perturbations, and the second defines small-scale, higher-order perturbations. For this dataset, we used the G2Aero database of over 19,000 airfoil shapes to learn a parameter space that captured a wide array of shape characteristics. We sampled airfoil designs over both parameter spaces to explore the full range of possible shape variations. The aerodynamic quantities for the generated airfoil were obtained using the HAM2D code, which is a finite-volume Reynolds-averaged Navier-Stokes (RANS) flow solver. We employ a fifth-order WENO scheme for spatial reconstruction with Roe's flux difference scheme for inviscid flux and second-order central differencing for viscous flux. A preconditioned GMRES method is applied for implicit integration. The Spalart-Allmaras 1-eq turbulence model is used for the turbulence closure, and the Medida-Baeder 2-eq transition model is applied to account for the effects of laminar turbulent transition. The airfoil grid is generated with a total of 400 points on the airfoil surface, the initial wall-normal spacing of y+ = 1, and an outer boundary located at 300 chord lengths away from the wall. The CFD simulations are performed at a freestream Mach number of 0.1, for or three different Reynolds' numbers (3M, 6M, and 9M), and for 25 angles of attack from -4 deg. to 20 deg. with 1 degree increments. Across all these various parameters, this dataset includes the results from over 250,000 CFD simulations. The simulations were performed using the Bridges-2 system at the Pittsburgh Supercomputing Center in February 2023 as part of the INTEGRATE project funded by the Advanced Research Projects Agency - Energy, in the U.S. Department of Energy. The data was collected, reformatted, and preprocessed for this OEDI submission in July 2023 under the Foundational AI for Wind Energy project funded by the U.S. Department of Energy Wind Energy Technologies Office. This dataset is intended to serve as a benchmark against which new artificial intelligence (AI) or machine learning (ML) tools may be tested. Baseline AI/ML methods for analyzing this dataset have been implemented, and a link to their repository containing those models has been provided. The .h5 data file structure can be found in the GitHub Repository resource under explore_airfoil_2k_data.ipynb.

2k↗

Airfoil Computational Fluid Dynamics - 9k shapes, 2 AoA's

This dataset contains aerodynamic quantities - including flow field values (momentum, energy, and vorticity) and summary values (coefficients of lift, drag, and momentum) - for 8,996 airfoil shapes, computed using the HAM2D CFD (computational fluid dynamics) model. The airfoil shapes were designed using the separable shape tensor parameterization that encodes two-dimensional shapes as elements of the Grassmann manifold. This data-driven approach learns two independent spaces of parameter from a collection of sample airfoils. The first captures large-scale, linear perturbations, and the second defines small-scale, higher-order perturbations. For this data, we used the G2Aero database of over 19,000 airfoil shapes to learn a parameter space that captured a wide array of shape characteristics. We fixed the linear deformations to be the mean over the database and sampled new shapes over a four-dimensional parameter space of higher-order perturbation. This sampling approaches allows for isolated analysis of non-linear airfoil shape deformations while holding other aspects (e.g., airfoil thickness) approximately constant. The aerodynamic quantities for the generated airfoil were obtained using the HAM2D code, which is a finite-volume Reynolds-averaged Navier-Stokes (RANS) flow solver. We employ a fifth-order WENO scheme for spatial reconstruction with Roe's flux difference scheme for inviscid flux and second-order central differencing for viscous flux. A preconditioned GMRES method is applied for implicit integration. The Spalart-Allmaras 1-eq turbulence model is used for the turbulence closure, and the Medida-Baeder 2-eq transition model is applied to account for the effects of laminar turbulent transition. The airfoil grid is generated with a total of 400 points on the airfoil surface, the initial wall-normal spacing of y+ = 1, and an outer boundary located at 300 chord lengths away from the wall. The CFD simulations are performed at a freestream Mach number of 0.1, Reynolds number of 9M, and at two angles of attack, 4 deg. and 12 deg. The simulations were performed using the Bridges-2 system at the Pittsburgh Supercomputing Center in February 2023 as part of the INTEGRATE project funded by the Advanced Research Projects Agency - Energy in the U.S. Department of Energy. The data was collected, reformatted, and preprocessed for this OEDI submission in July 2023 under the Foundational AI for Wind Energy project funded by the U.S. Department of Energy Wind Energy Technologies Office. This dataset is intended to serve as a benchmark against which new artificial intelligence (AI) or machine learning (ML) tools may be tested. Baseline AI/ML methods for analyzing this dataset have been implemented, and a link to their repository containing those models has been provided. The .h5 data file structure can be found in the GitHub Repository resource under explore_airfoil_9k_data.ipynb.

9k↗

Implementation of Two Local Correlation-Based Transition Models in OVERFLOW 2.3e

Accurate prediction of the laminar-to-turbulent boundary layer transition is highlighted as a vital component of NASA's CFD Vision 2030 Study. With an increased emphasis on greener air transports and sustainable aviation, transition modeling is anticipated to have an added significance, particularly in the applications related to laminar flow technology. However, unmanned aerial vehicles, crewed reentry vehicles, and ground-to-flight extrapolation all benefit from transition models. It is useful to incorporate a variety of models in CFD solvers because there isn't a single transition model that is ideal for the complete spectrum of applications. While the Langtry-Menter γ–Reθt model has been widely used for CFD predictions of flows with laminar, transitional, and turbulent boundary layers, it does not meet the criteria for Galilean invariance, a desirable attribute for rotorcraft applications. To help overcome that limitation, we have implemented Menter's SST-based γ transition model within NASA’s OVERFLOW Overset CFD code (version 2.3e). In addition, an SA-based implementation of the Langtry-Menter γ–Reθt transition model has also been implemented in view of the increased robustness of the SA-turbulence model for external aerodynamic applications. An initial assessment of the implemented models has been carried out using benchmark 2D test cases involving flat plates and the NLF-0416 airfoil, focusing on bypass, Tollmien-Schlichting, and laminar separation bubble-induced transition scenarios. In comparison to the Langtry-Menter γ–Reθt model, both models were found to produce increased computation efficiency, principally as a result of the decrease from four to three auxiliary transport equations. Results from the initial analysis will be presented, along with additional details pertaining to the deployment of these models.

CFD modeling↗

Requirements and Challenges for CFD Validation within the High-Lift Common Research Model Ecosystem

The High-Lift Common Research Model (CRM-HL) ecosystem is envisioned to become an industry standard set of test cases for high-lift aerodynamics prediction, with data generated providing a strong foundation for CFD validation. Since its original development, preliminary wind tunnel tests have been conducted to further understand and configure the geometry, and several more are in the planning phases. It is anticipated that the data generated in these tests will be used extensively by the CFD community for research, development, and validation of new techniques that yield more accurate results. To derive the maximum benefit from planned wind tunnel testing, data acquisition will be focused in several key areas to address specific shortcomings in CFD predictive capabilities. Significant challenges currently exist in areas including, but not limited to, characterizing high-lift flow phenomena, quantifying effects due to configuration changes, understanding wind tunnel and model installation effects, and establishing data uncertainty. This paper will discuss the current state-of-the-art in high-lift CFD prediction, key limitations of current CFD methods and tools, establishing a validation dialog between test and CFD practitioners, and the flow phenomena of highest interest, in an effort to begin to address these challenges. To this end, it is expected that the data generated from CRM-HL tests will form a strong foundation for future predictive capability.

Adam M Clark↗

Three-Dimensional CFD Simulations in the Film Boiling Regime During Liquid Nitrogen Chilldown Process

The paper presents a three-dimensional CFD simulation for predicting the liquid nitrogen chilldown process through a SS-304 tube. A dispersed mixture model has been employed to predict the two-phase flow behavior whereas the Lee phase change model has been used to predict the boiling process. Cryogenic chilldown process starts with a pure vapor region and proceeds to film, transition and nucleate boiling regimes before reaching the single phase liquid flow and complete chilldown of the transfer line. The present numerical model predicted the film boiling regime accurately and the slope of the chilldown curve deviates slightly from the experimental curve towards the transition and nucleate boiling regime, which needs further investigation. The contours of liquid volume fraction and temperature profiles has been presented to understand the inverted annular film flows in the film boiling regime in detail. In future, the present CFD model will be extended to a wide range of chilldown conditions and different cryogenic fluids to develop a robust methodology.

CFD↗

The Effect of Turbulence Modeling on the Mixing Characteristics of Several Fuel Injectors at Hypervelocity Flow Conditions

CFD analysis is presented on the effects of turbulence modeling choices on the mixing characteristics and performance of three fuel injectors at hypervelocity flow conditions. The analyses were carried out with the VULCAN-CFD solver using Reynolds-Averaged Simulations (RAS). The hypervelocity flow conditions match the high Mach number flow of the experiments conducted as a part of the Enhanced Injection and Mixing Project (EIMP) at the NASA Langley Research Center. The three injectors are the baseline configurations used in the experiments and represent three categories of injectors typically considered individually or in combination for fueling high-speed propulsive devices. The current work discusses the impact of the turbulence model and the turbulent Schmidt number on the mixing flow field behavior and the mixing performance as described by the one-dimensional values of the Mach number, total pressure recovery, and the mixing efficiency. Because planar laser induced fluorescence (PLIF) images are available from the EIMP experiments, the sensitivity of the synthetic LIF signal to turbulence modeling choices is also examined to determine whether PLIF can be extended beyond its intended qualitative visualization purpose and used to guide CFD turbulence model and parameter selections. It is found that the mixing performance, as quantified using mixing efficiency, exhibits a strong sensitivity to both turbulence model choice and turbulent Schmidt number value. However, the synthetic LIF signal only demonstrates a modest level of sensitivity, which suggests that PLIF is of limited use for guiding CFD turbulence model and parameter selections.

Drozda, Tomasz G.↗

Infrared imaging - A validation technique for computational fluid dynamics codes used in STOVL applications

The need for a validation technique for computational fluid dynamics (CFD) codes in STOVL applications has led to research efforts to apply infrared thermal imaging techniques to visualize gaseous flow fields. Specifically, a heated, free-jet test facility was constructed. The gaseous flow field of the jet exhaust was characterized using an infrared imaging technique in the 2 to 5.6 micron wavelength band as well as conventional pitot tube and thermocouple methods. These infrared images are compared to computer-generated images using the equations of radiative exchange based on the temperature distribution in the jet exhaust measured with the thermocouple traverses. Temperature and velocity measurement techniques, infrared imaging, and the computer model of the infrared imaging technique are presented and discussed. From the study, it is concluded that infrared imaging techniques coupled with the radiative exchange equations applied to CFD models are a valid method to qualitatively verify CFD codes used in STOVL applications.

Hardman, R. R.↗

Validation of Slosh Modeling Approach Using STAR-CCM+

Without an adequate understanding of propellant slosh, the spacecraft attitude control system may be inadequate to control the spacecraft or there may be an unexpected loss of science observation time due to higher slosh settling times. Computational fluid dynamics (CFD) is used to model propellant slosh. STAR-CCM+ is a commercially available CFD code. This paper presents a comparison between STAR-CCM+ liquid slosh modeling results and experimental, empirically, and analytically derived results. The geometries examined are a bare right cylinder tank and a right cylinder

Benson, David↗

Aerodynamic and Structural Evaluation of an SMA Slat-Cove Filler Using Computational and Experimental Tools at Model Scale

Transport class aircraft produce a significant amount of airframe noise during approach and landing due to exposed geometric discontinuities that are hidden during cruise. The leading-edge slat is a primary contributor to this noise. In previous work, use of a slat-cove filler (SCF) has proven to reduce airframe noise by filling the cove aft of the slat, eliminating the circulation region within the cove. The goal of this work is to extend and improve upon past experimental and computational efforts on the evaluation of a scaled high-lift wing with a super elastic shape memory alloy (SMA) SCF. Recent turbulence measurements of the Texas A&M University 3ft-by-4ft wind tunnel allow for more accurate representation of the flow through the test section in computational fluid dynamics (CFD) analysis. The finite volume models used in CFD analysis are coupled to structural finite element models using a framework compatible with an SMA constitutive model and significant deformation, enabling fluid-structure interaction (FSI) analysis of the SCF. Both fully-deployed and retraction/deployment cases are considered. The displacement of the SCF on the experimental model is measured at various stages of retraction/deployment using a laser displacement sensor and digital image correlation system. Due to a lack of structural stiffness in the 3D-printed plastic slat during retraction and SCF stowage, a rigid steel slat is incorporated into the physical model and preliminary wind tunnel tests are conducted at multiple angles of attack in both retracted and deployed configurations.

Scholten, William↗

Evaluation of semi-empirical soot models for nonpremixed flames with increased stoichiometric mixture fraction and strain

Due to the size and complexity of industrial-scale systems and computing limitations, semi-empirical soot models are often employed in CFD simulations rather than computationally-expensive detailed models. Many of these models were developed for specific applications with unique characteristic timescales and/or validated only under fuel-air combustion conditions. Hence, their use in different contexts, such as oxy-combustion, could lead to inaccurate predictions. In this study, twelve semi-empirical models (1-step or 2-step) are evaluated on their ability to respond to changes in stoichiometric mixture fraction (Z st ) and strain in a series of ethylene counterflow flames which span the experimental sooting-to-non-sooting (yellow to blue) transition. A unique approach of plotting soot formation rate in normalized, local equivalence ratio space is introduced to aid in the analysis. Results show that no existing model is able to predict a blue flame when Z st is increased beyond the experimentally-measured sooting limit. Many models give the counter-intuitive result of increasing peak and/or integrated soot volume fraction (svf) as Z st increased, contrary to experimental observations. Two models are able to successfully predict a blue flame upon increase of strain to the sooting limit value. For two-step models, the strong dependence of growth rate on surface area results in high sensitivity of svf to changing flame boundary conditions, but with mixed accuracy. There is still a significant need for a robust semi-empirical model which can accurately predict soot fraction in systems with oxygen-enrichment or variable strain.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Computational Fluid Dynamics (CFD) applications in rocket propulsion analysis and design

Computational Fluid Dynamics (CFD) has been used in recent applications to affect subcomponent designs in liquid propulsion rocket engines. This paper elucidates three such applications for turbine stage, pump stage, and combustor chamber geometries. Details of these applications include the development of a high turning airfoil for a gas generator (GG) powered, liquid oxygen (LOX) turbopump, single-stage turbine using CFD as an integral part of the design process. CFD application to pump stage design has emphasized analysis of inducers, impellers, and diffuser/volute sections. Improvements in pump stage impeller discharge flow uniformity have been seen through CFD optimization on coarse grid models. In the area of combustor design, recent CFD analysis of a film cooled ablating combustion chamber has been used to quantify the interaction between film cooling rate, chamber wall contraction angle, and geometry and their effects of these quantities on local wall temperature. The results are currently guiding combustion chamber design and coolant flow rate for an upcoming subcomponent test. Critical aspects of successful integration of CFD into the design cycle includes a close-coupling of CFD and design organizations, quick turnaround of parametric analyses once a baseline CFD benchmark has been established, and the use of CFD methodology and approaches that address pertinent design issues. In this latter area, some problem details can be simplified while retaining key physical aspects to maintain analytical integrity.

Mcconnaughey, P. K.↗

Investigating the Effects of Chemical Mechanism on Soot Formation Under High-Pressure Fuel Pyrolysis

We performed Computational Fluid Dynamics (CFD) simulations using a Reynolds-Averaged Navier-Stokes (RANS) turbulence model of high-pressure spray pyrolysis with a detailed chemical kinetic mechanism encompassing pyrolysis of n-dodecane and formation of polycyclic aromatic hydrocarbons. We compare the results using the detailed mechanism and those found using several different reduced chemical mechanisms to experiments carried out in an optically accessible, high-pressure, constant-volume combustion chamber. Three different soot models implemented in the CONVERGE CFD software are used: an empirical soot model, a method of moments, and a discrete sectional method. There is a large variation in the prediction of the soot between different combinations of chemical mechanisms and soot model. Furthermore, the amount of soot produced from all models is substantially less than experimental measurements. All of this indicates that there is still substantial work that needs to be done to arrive at simulations that can be relied on to accurately predict soot formation.

02 PETROLEUM↗

Nonequilibrium radiative heating prediction method for aeroassist flowfields with coupling to flowfield solvers

A method for predicting radiation adsorption and emission coefficients in thermochemical nonequilibrium flows is developed. The method is called the Langley optimized radiative nonequilibrium code (LORAN). It applies the smeared band approximation for molecular radiation to produce moderately detailed results and is intended to fill the gap between detailed but costly prediction methods and very fast but highly approximate methods. The optimization of the method to provide efficient solutions allowing coupling to flowfield solvers is discussed. Representative results are obtained and compared to previous nonequilibrium radiation methods, as well as to ground- and flight-measured data. Reasonable agreement is found in all cases. A multidimensional radiative transport method is also developed for axisymmetric flows. Its predictions for wall radiative flux are 20 to 25 percent lower than those of the tangent slab transport method, as expected, though additional investigation of the symmetry and outflow boundary conditions is indicated. The method was applied to the peak heating condition of the aeroassist flight experiment (AFE) trajectory, with results comparable to predictions from other methods. The LORAN method was also applied in conjunction with the computational fluid dynamics (CFD) code LAURA to study the sensitivity of the radiative heating prediction to various models used in nonequilibrium CFD. This study suggests that radiation measurements can provide diagnostic information about the detailed processes occurring in a nonequilibrium flowfield because radiation phenomena are very sensitive to these processes.

Hartung, Lin C.↗