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 865 records · Page 48

Trajectory Reconstruction of the Low-Earth Orbit Flight Test of an Inflatable Decelerator

The Low-Earth Orbit Flight Test of an Inflatable Decelerator (LOFTID) project conducted a flight test of a 6m inflatable aeroshell. The LOFTID test article was a secondary payload on an Atlas V launcher that carried the Joint Polar Satellite System-2 (JPSS-2) as its primary payload. The vehicle launched on November 10th, 2022. After reaching orbit, the LOFTID test article inflated the aeroshell, separated from the upper stage on an entry trajectory, and entered the atmosphere to splash down in the Pacific Ocean under parachutes. The test concept of operations is shown in Figure 1. The test article was instrumented with a variety of sensors to be used for post-flight evaluation of vehicle performance. Data from one of the key sensors for trajectory reconstruction, the Inertial Measurement Unit (IMU), was not captured in the data recorder due to a malfunction. Data from the nose cone mounted Flush Air Data Sensing (FADS) system were successfully acquired. The layout of the FADS sensors and the measured pressures during atmospheric entry are shown in Figure 2. The FADS data were combined with a Newtonian flow pressure model [1, 2] to produce estimates of the atmospheric relative trajectory. A Mach number anchoring technique given in [2] was used to stabilize estimates in high speed flight conditions. Since no IMU data were available, a trajectory simulation was used to provide the Mach number time history. The resulting estimates of the atmospheric-relative trajectory are shown in Figures 3. Given the loss of the IMU data, alternate methods for trajectory reconstruction are being explored. One approach under investigation is the use of the on-board video recorder data to be analyzed to reconstruct attitude motion. This approach is currently under investigation and will be reported on in the final paper. The Newtonian flow pressure model for the FADS analysis will also be updated with a CFD-based pressure model.

Christopher D. Karlgaard↗

Lagrangian simulations of spray interaction with a surface using a stochastic multi-regime impingement model

Modeling of spray impingement upon the piston or cylinder surface is important for hydrocarbon emission and soot predictions in high-fidelity simulation of internal combustion engines with direct injections. Due to the complexity of the physical process of spray impingement, current computational models heavily rely on the empirical studies which involve large uncertainties. A computational model is presented for spray impingement on solid walls in general computational fluid dynamics (CFD) simulations. The model adopts the Marmanis-Thoroddsen theory for determining the number of secondary droplets and a stochastic method for calculating the velocities of the secondary droplets. The numerical methods and parallel implementation of the impingement model are presented in detail. The model is implemented in Grit, which is a recently developed CPU/GPU performance-portable Lagrangian particle library, for accelerating spray simulations. The model and its implementation are tested for examples of spray impingements on flat and smooth surfaces with different impingement angles.

Ge, Wenjun↗

Complicated Geometry

Using NASA SBIR funding, CFD Research Corporation has developed CFD-GEOM, an extension of traditional computer-aided drawing (CAD) software. CFD-GEOM can provide modeling and interactivity of computational fluid dynamics (CFD) latest field, mesh-generation and allows for quick and easy updating of a grid in response to changes in the CAD model.

Source record↗

Extraction of Dynamic Inflow Models for Coaxial and Tandem Rotors from CFD Simulations

The dynamic inflow coupling with rotor/body dynamics is crucial in the analysis of stability and control law design for helicopters. Over the past several decades, finite-state inflow models for single rotor configurations in hover, forward flight, and maneuver have developed (Ref.1-3). By capturing the interference effects between rotors, the extension of pressure potential finite state inflow model has promising result for coaxial rotor configuration (Ref.4-6). Recently, the focus of the dynamic inflow modeling has shifted to tandem rotor configurations (Ref.7, 8). The development of the dynamic inflow models for tandem rotor configuration still have some limitations due to the lack of knowledge of rotor-to-rotor interference, and rotor-wake interference. Experimental methods, and computational fluid dynamics methods are commonly used to understand the rotor performance and rotor airload variations, and measure or predict inflow velocity distributions at the rotor desk. The inflow distributions are subsequently used to improve the dynamic inflow models. Tandem rotor configurations have been studied experimentally and computationally for several decades (Ref.9-12). Sweet (Ref.10) observed that a tandem rotor with 76-percent-radius overlap required 14% more induced power at hovering condition, relative to an isolated rotor of equivalent disk area. Sweet also found that, above a shaft-to-shaft distance of 1.03 diameter, the performance of the tandem rotor was nearly the same as two isolated rotors. The objective of the present study is to apply computational fluid dynamics simulations of tandem rotors for the extraction of dynamic inflow models. The extended methodology is first validated by comparing the computed induced power against test data. Subsequently inflow distributions and wake structures are analyzed.

Extraction of Dynamic Inflow Models↗

Numerical Models of Human Circulatory System under Altered Gravity: Brain Circulation

A computational fluid dynamics (CFD) approach is presented to model the blood flow through the human circulatory system under altered gravity conditions. Models required for CFD simulation relevant to major hemodynamic issues are introduced such as non-Newtonian flow models governed by red blood cells, a model for arterial wall motion due to fluid-wall interactions, a vascular bed model for outflow boundary conditions, and a model for auto-regulation mechanism. The three-dimensional unsteady incompressible Navier-Stokes equations coupled with these models are solved iteratively using the pseudocompressibility method and dual time stepping. Moving wall boundary conditions from the first-order fluid-wall interaction model are used to study the influence of arterial wall distensibility on flow patterns and wall shear stresses during the heart pulse. A vascular bed modeling utilizing the analogy with electric circuits is coupled with an auto-regulation algorithm for multiple outflow boundaries. For the treatment of complex geometry, a chimera overset grid technique is adopted to obtain connectivity between arterial branches. For code validation, computed results are compared with experimental data for steady and unsteady non-Newtonian flows. Good agreement is obtained for both cases. In sin-type Gravity Benchmark Problems, gravity source terms are added to the Navier-Stokes equations to study the effect of gravitational variation on the human circulatory system. This computational approach is then applied to localized blood flows through a realistic carotid bifurcation and two Circle of Willis models, one using an idealized geometry and the other model using an anatomical data set. A three- dimensional anatomical Circle of Willis configuration is reconstructed from human-specific magnetic resonance images using an image segmentation method. The blood flow through these Circle of Willis models is simulated to provide means for studying gravitational effects on the brain circulation under auto-regulation.

Kim, Chang Sung↗

Low-dimensional modelling and prediction in fluid dynamics

Research into the development of a simplified, low-dimensional, state-variable model for a fluid dynamics system is discussed. The model would be capable of predicting the state of the system in as nearly real time as possible and would be suitable for implementation in a control loop. It is noted that state-variable modeling has an advantage for real-time health monitoring and control application over the Navier-Stokes equations since the latter require an enormous amount of real time. In this study, a dynamic procedure is developed for modeling and prediction using the CFD data generated for this model centerbody-combustor configuration. Some experiences with the NASA data available for the SSME are also reported. Application of this procedure to pressure and velocity measurements in an actual wind-tunnel experiment is described, and the applicability of the approach to health monitoring is considered.

Ghia, U.↗

An Efficient 1-D Thermal Stratification Model for Pool-Type Sodium-Cooled Fast Reactors

Investigating thermal stratification in the upper plenum of a sodium fast reactor (SFR) is presently a technology gap in SFR safety analysis. Understanding thermal stratification will promote safe operation of the SFR before its commercial deployment. Stratified layers of liquid sodium with a large vertical temperature gradient could be established in the upper plenum of an SFR during a down-power or a loss-of-flow transient. These stratified layers are unstable and could result in uncertainties for the core safety of an SFR. In order to predict the occurrence of the thermal stratification efficiently, we developed a one-dimensional (1-D) transport model to estimate the temperature profile of the ambient fluid in the upper plenum. This model demands much less computational effort than computational fluid dynamics (CFD) codes and provides calculations with higher fidelity than historical system-level codes. Two flow conditions were considered separately in the current study depending on if in-vessel components are presented in the upper plenum. For the condition where in-vessel components, specifically the upper internal structure, are presented, we assumed that the impinging sodium was evenly dispersed in the ambient fluid within the distance between the bottom of the in-vessel component and the jet inlet surface. For the condition where no in-vessel components are presented, we assumed that the impinging sodium was evenly dispersed in the ambient fluid within the jet length, which was determined through data-driven trainings. The newly developed 1-D model showed similar performance with the CFD model in both cases. However, due to the assumption of flat profiles of the impinging jet axial dispersion rate, nonnegligible discrepancies between the 1-D prediction and the measured data were observed.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Modeling DWPF Hydrogen Gas Retention and Release

Savannah River Remediation (SRR) currently manages the risk for retained hydrogen in the Defense Waste Processing Facility (DWPF) vessels by implementing a Retained Hydrogen Program. The current program conservatively relies on bounding assumptions concerning gas release and retention to develop allowable vessel Q-times. SRR would like to establish more realistic inputs to the Retained Hydrogen Program. To address this need, SRNL conducted a literature search on work investigating the retention and release of gas from DWPF type slurries. The search revealed that the retained hydrogen is not released instantaneously upon the start of agitation, and that hydrogen gas could be released under quiescent, non-agitated, conditions. The literature search also identified a Lattice-Boltzmann computational fluid dynamics (CFD) software that could model and predict hydrogen release from DWPF vessels. SRR requested that Savannah River National Laboratory (SRNL) personnel pursue this software and focus the model for hydrogen gas retention and release to justify relaxing these conservative assumptions. This report documents the refinement of the CFD software model and application to show its ability to describe retained hydrogen release from DWPF vessels and relax the conservativisms in the Retained Hydrogen Program.

08 HYDROGEN↗

Optimization of spray breakup model parameters for predicting fuel spray and film characteristics in gasoline direct injection engines

This study investigated the behavior of gasoline direct injection (GDI) sprays using computational fluid dynamics (CFD). The authors developed an approach to identify optimal spray breakup model parameters by evaluating an error function across numerous simulations, with the goal of minimizing discrepancies from experimental data. Using the optimal setup, the simulated spray matched well with projected liquid volume distributions, liquid penetration, and spray width measured in a constant-pressure continuous-flow chamber. To further validate the approach, the same setup was tested across various fuels, injectors, and operating conditions. Subsequently, the optimal setup, along with a recently developed spray-wall interaction model, were applied to a direct-injection spark-ignited engine under late-injection conditions to predict and evaluate fuel film formation and evolution at varying engine coolant temperatures. Here, with the centrally mounted injector directing the spray toward the piston, simulations indicated that the spray tends to impinge on the piston surface. The proposed simulation framework also accurately captured the aggregate film area on the piston surface, aligning with previously published experimental results. Moreover, simulations showed that increasing the coolant temperature from cold start conditions (333 K) to warm conditions (363 K) reduced the fuel mass deposited on the piston by roughly 50%. Furthermore, for the spray-guided engine configuration studied in this work, the CFD model predicted minimal film deposition on the spark plug electrodes regardless of the coolant temperatures due to a relatively weak in-cylinder flow during the compression phase.

Computational fluid dynamics (CFD)↗

Final Plenary Session Transcript

Let s start with the report - as you know - when we are talking about flow control it is a multi-disciplinary type of work. So it involves many people and disciplines. This group first discussed the important issues associated with flow control. As you start doing flow control what are the issues to which you really have to start paying attention? That is the first part I am going to present. Then in the second part I will present some challenges - problems that we should really be looking at. So as far as the issues - if you want to control a flow, you really need to understand the flow physics, because anything that you do comes from the flow physics. The design of the controllers, your decisions on the actuators, sensors, reduced order modeling and all of that, would be helped if you understand flow physics. And you have to have a specific objective - what exactly are you controlling? Are you trying to reduce drag, eliminate separation, reduce noise, enhance mixing? So you have to have very specific control objectives. From all the talks we have seen here actuation is extremely important and it is very problem specific. It depends on what problem you are dealing with so you have to design and build actuators for that specific problem. Sensors obviously are very important, especially when you are dealing with feedback control. Consensus was that when you dealing with flow control, you must take an integrated approach; from the beginning you have to take into account every aspect of it and even maybe to modify your experiment, your geometry, to go along with the actuation, sensors and control models. Development of tools is very important in this multi-disciplinary problem. The tools include CFD, reduced order modeling, controller design, understanding and utilizing the instabilities of the flow, etc. So, in order to have success in flow control, we really need to develop these tools.

Gostelow, Paul↗

Implementation and Assessment of Menter’s Galilean-Invariant γ Transition Model in OVERFLOW

With an increased emphasis on greener air transports and sustainable aviation, the modeling of laminar-to-turbulent boundary layer transition 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. Because no single transition model is ideal for the complete spectrum of applications, it is useful to incorporate a variety of models in general-purpose CFD solvers, such as the NASA OVERFLOW Overset CFD code. While the Langtry-Menter 𝛄 − 𝑹𝒆 𝛉𝒕 model, currently available in OVERFLOW, has been widely used for CFD predictions of flows with laminar, transitional, and turbulent boundary layers, it does not possess the Galilean invariance property, a desirable attribute for rotorcraft applications. To help overcome that limitation, we have recently implemented Menter's baseline version of the SST-based γ transition model, along with a Galilean invariant stationary crossflow extension within OVERFLOW (version 2.3e). An initial assessment of the newly implemented model has been carried out using 2D benchmark cases including flat plates and the NLF-0416 airfoil, addressing several transition scenarios ranging from bypass transition due to freestream turbulence, natural transition via Tollmien-Schlichting instabilities, and transition due to a laminar separation bubble. The crossflow extension has been applied to the infinite swept NLF(2)-0415 wing and the 6:1 prolate spheroid. Wherever possible, the results were obtained on a sequence of meshes to ascertain the grid convergence behavior, which has been evaluated through global metrics such as force coefficients as well as local values of the skin-friction coefficient at selected points near and within the transition region. Overall, the model appears to be correctly implemented and the results show promise for further development using the framework of the γ transition model.

CFD modeling↗

Evaluation of Transport-Equations-Based Transition Models for High-Speed Boundary Layers Using OVERFLOW

Accurate modeling of laminar-turbulent transition is crucial for the design of hypersonic flight systems. However, the current transition models used in production CFD codes are insufficient for high-speed flows. Many extensions to low-speed models have been suggested; however, a thorough verification and validation effort is needed before these models can be used in design settings. Challenges include potentially missing details of the model implementation requirements and/or a complete specification of the input parameters needed to replicate the test findings. A meaningful assessment of the generalization capability of these models is also hindered by a lack of information regarding the specific flow configurations and associated grids employed for model calibration. As a key first step toward model verification, we present an independent assessment of two recently proposed models for high-speed transition, namely, a model within the SST-𝛾 framework and a model based on the SST-𝛾 – 𝜈 𝐿 equations. These models are implemented in the NASA OVERFLOW 2.3e solver and their performance in predicting first mode, second mode, and crossflow transition has been evaluated for several test cases in the supersonic and hypersonic regimes. Besides the test cases employed by the model developers, which could have also been used for model calibration, the present assessment includes supplementary configurations that contribute to an unbiased assessment of the models. The outcomes presented in this study indicate the potential for the models to be applied to high-speed flight configurations. Key steps toward future improvements to these models are also outlined.

CFD↗

Generation of Mars Helicopter Rotor Model for Comprehensive Analyses

The present research is aimed at providing a performance model for the Mars Helicopter (MH), to understand the complexity of the flow, and identify future regions of improvement. The low density of the Martian atmosphere and the relatively small MH rotor, result in very low chord-based Reynolds number flows. The low density and Reynolds numbers reduce the lifting force and lifting efficiency, respectively. The high drag coefficients in subcritical flow, especially for thicker sections, are attributed to laminar separation from the rear of the airfoil. In the absence of test data, efforts have been made to explore these effects using prior very low Reynolds number research efforts. The rotor chord-based Reynolds number range is observed to be subcritical, which makes boundary layer transition unlikely to occur. The state of the two-dimensional rotor boundary layer in hover is approximated by calculating the instability point, laminar separation point, and the transition location to provide understanding of the flow state in the high Mach-low Reynolds number regime. The results are used to investigate the need for turbulence modeling in Computational Fluid Dynamics (CFD) calculations afterwards. The goal is to generate a performance model for the MH rotor for a free wake analysis, because the computational budget for a complete Navier-Stokes solution for a rotating body-fitted rotor is substantial. In this study, a Reynolds-Averaged Navier-Stokes (RANS) based approach is used to generate the airfoil deck using C81Gen with stitched experimental data for very high angles of attack. A full Grid Resolution Study is performed and over 4,500 cases are completed to create the full airfoil deck. The laminar separation locations are predicted within the accuracy of the approximate method when compared with the CFD calculations. The model is presented through airfoil data tables (c81 files) that are used by comprehensive rotor analysis codes such as CAMRADII, or the mid-fidelity CFD solver RotCFD. Finally, the rotor performance is compared with experimental data from the 25ft Space Simulator at the NASA Jet Propulsion Laboratory (JPL) and shows good correlation for the rotor Figure of Merit over the available thrust range.

Model↗

Automated Boundary Conditions for Wind Tunnel Simulations

Computational fluid dynamic (CFD) simulations of models tested in wind tunnels require a high level of fidelity and accuracy particularly for the purposes of CFD validation efforts. Considerable effort is required to ensure the proper characterization of both the physical geometry of the wind tunnel and recreating the correct flow conditions inside the wind tunnel. The typical trial-and-error effort used for determining the boundary condition values for a particular tunnel configuration are time and computer resource intensive. This paper describes a method for calculating and updating the back pressure boundary condition in wind tunnel simulations by using a proportional-integral-derivative controller. The controller methodology and equations are discussed, and simulations using the controller to set a tunnel Mach number in the NASA Langley 14- by 22-Foot Subsonic Tunnel are demonstrated.

Carlson, Jan-Renee↗

Development & Experimental Validation of a Generalized Resistance-Capacitance Model for Numerical Simulation of Phase-Change Material Embedded Heat Exchangers

Latent heat thermal energy storage (LHTES) using phase change material (PCM) has attracted increased attention as a viable solution for overcoming the mismatch between energy supply and demand for renewable energy-based systems. PCM-embedded heat exchangers (PCM-HX) have the potential to significantly improve thermal performance due to high storage capacity and low temperature variation during the phase change process. Most models for simulating LHTES heat transfer use Computational Fluid Dynamics (CFD) simulations, which have high computational costs resulting from considering the complex and time-dependent physics relevant to PCM-HXs. In this paper, a Generalized Resistance Capacitance-based Model (GRCM) was developed to predict the thermal performance of arbitrary PCM-HXs in a computationally efficient manner without compromising modeling accuracy. The GRCM is exercised for three case studies: (i) verification for a single-slabbed finned PCM-HX, (ii) verification and validation for a copper foam/paraffin composite PCM-HX, and (iii) validation for a straight tube annular finned PCM-HX. The copper foam PCM-HX uses an electric heater at the top of HX, while the other two configurations utilize water as heat transfer fluid. For the single-slabbed finned PCM-HX melting case, the mean deviation in average PCM temperature predicted by the GRCM compared to the CFD model was between 0.56 – 0.73 K, with maximum temperature deviation of 2.68 K. For the HTF outlet temperature, the validation results showed that GRCM prediction matches very well with experimental data, with mean temperature deviation of 0.24 K during melting case, while for solidification case was 0.34 K. These results showcase the GRCM’s capability for accurately reproducing the thermal characteristics of PCM-HXs with considerably lower computational effort.

42 ENGINEERING↗

Evaluation and application of the Baldwin-Lomax turbulence model in two-dimensional, unsteady, compressible boundary layers with and without separation in engine inlets

There is a practical need to model high speed flows that exist in jet engine inlets. The boundary layers that form in these inlets may be turbulent or laminar and either separated or attached. Also, unsteady supersonic inlets may be subject to frequent changes in operating conditions. Some changes in the operating conditions of the inlets may include varying the inlet geometry, bleeds and bypasses, and rotating or translating the centerbody. In addition, the inlet may be either started or unstarted. Therefore, a CFD code, used to model these inlets, may have to run for several different cases. Also, since the flow conditions through an unsteady inlet may be continually fluctuating, the CFD code which models these flows may have to be run over many time steps. Therefore, it would be beneficial that the code run quickly. Many turbulence models, however, are cumbersome to implement and require a lot of computer time to run, since they add to the number of differential equations to be solved to model a flow. The Baldwin-Lomax turbulence model is a popular model. It is an algebraic, eddy viscosity model. The Baldwin-Lomax model is used in many CFD codes because it is quick and easy to implement. In this paper, we will 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 will 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↗

Comparison of CFD and DSMC Using Calibrated Transport Parameters

Hypersonic re-entry flows span a wide range of length scales where regions of both rarefied and continuum flow exist. Traditional computational fluid dynamics (CFD) techniques do not provide an accurate solution for the rarefied regions of such ‘mixed’ flow fields. Although direct simulation Monte Carlo (DSMC) can be used to accurately capture both the continuum and rarefied features of ‘mixed’ flow fields, they are computationally expensive when employed to simulate the low Knudsen number continuum regimes. Thus, a hybrid framework for seamlessly combining the two methodologies, CFD and DSMC, continues to be a topic of significant research effort. Ensuring consistency in the reaction kinetics and transport models employed within CFD and DSMC is a crucial requirement for obtaining a reliable solution from a hybrid framework for combined continuum/rarefied high speed flows. This paper represents one of the first studies to utilize the calibrated transport parameters developed to ensure consistency between CFD and DSMC solvers. The new variable soft sphere (VSS) parameters are compared to both previous “standard” variable hard sphere (VHS) parameters and also to solutions from the CFD transport properties that the new parameters were developed to reproduce.

Liechty, D. S.↗

Turbulence program for propulsion systems

Program goals at the Center for Modeling of Turbulence and Transition (CMOTT), NASA Lewis Research Center, are (1) to develop reliable turbulence (including bypass transition) and combustion models for complex flows in propulsion systems and (2) to integrate developed models into deliverable CFD tools for propulsion systems in collaboration with industry. This viewgraph presentation covers the following topics: development of turbulence and combustion models; collaboration with industry and technology transfer; isotropic eddy viscosity models; algebraic Reynolds stress models; scalar turbulence models; second order closure models; multiple scale k-epsilon models; and PDF modeling of turbulent reacting flows.

Shih, Tsan-Hsing↗