Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “computational fluid dynamics simulation”

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,153 records · Page 64

Longitudinal Aerodynamic Characteristics and Wing Pressure Distributions of a Blended-Wing-Body Configuration at Low and High Reynolds Numbers

Force balance and wing pressure data were obtained on a 0.017-Scale Model of a blended-wing-body configuration (without a simulated propulsion system installation) to validate the capability of computational fluid dynamic codes to predict the performance of such thick sectioned subsonic transport configurations. The tests were conducted in the National Transonic Facility of the Langley Research Center at Reynolds numbers from 3.5 to 25.0 million at Mach numbers from 0.25 to 0.86. Data were obtained in the pitch plane only at angles of attack from -1 to 8 deg at Mach numbers greater than 0.25. A configuration with winglets was tested at a Reynolds number of 25.0 million at Mach numbers from 0.83 to 0.86.

Re, Richard J.↗

CFD Modeling of Helium Pressurant Effects on Cryogenic Tank Pressure Rise Rates in Normal Gravity

A recently developed computational fluid dynamics modeling capability for cryogenic tanks is used to simulate both self-pressurization from external heating and also depressurization from thermodynamic vent operation. Axisymmetric models using a modified version of the commercially available FLOW-3D software are used to simulate actual physical tests. The models assume an incompressible liquid phase with density that is a function of temperature only. A fully compressible formulation is used for the ullage gas mixture that contains both condensable vapor and a noncondensable gas component. The tests, conducted at the NASA Marshall Space Flight Center, include both liquid hydrogen and nitrogen in tanks with ullage gas mixtures of each liquid's vapor and helium. Pressure and temperature predictions from the model are compared to sensor measurements from the tests and a good agreement is achieved. This further establishes the accuracy of the developed FLOW-3D based modeling approach for cryogenic systems.

Grayson, Gary↗

Equivalent Source Method Applied to Launch Acoustic Simulations

Aeroacoustic simulations of the launch environment are described. A hybrid computational fluid dynamics (CFD)/computational aeroacoustic (CAA) approach is developed in order to accurately and efficiently predict the sound pressure level spectrum on the launch vehicle and surrounding structures. The high-fidelity CFD code LAVA (Launch Ascent and Vehicle Analysis), is used to generate pressure time history at select locations in the flow field. A 3D exterior Helmholtz solver is then used to iteratively determine a set of monopole sources which mimic the noise generating mechanisms identified by the CFD solver. The acoustic pressure field generated from the Helmholtz solver is then used to evaluate the sound pressure levels.

Housman, Jeffrey A.↗

Mitigation of Autoignition Due to Premixing in a Hypervelocity Flow Using Active Wall Cooling

Preinjection of fuel on the forebody of an airbreathing vehicle is a proposed method to gain access to hypervelocity flight Mach numbers. However, this creates the possibility of autoignition either near the wall or in the core of the flow, thereby consuming fuel prematurely as well as increasing the amount of pressure drag on the vehicle. The computational fluid dynamics code VULCAN was used to conduct three dimensional simulations of the reacting flow in the vicinity of hydrogen injectors on a flat plate at conditions relevant to a Mach 12 notional flight vehicle forebody to determine the location where autoignition occurs. Active wall cooling strategies were formulated and simulated in response to regions of autoignition. It was found that tangential film cooling using hydrogen or helium were both able to nearly or completely eliminate wall autoignition in the flow domain of interest.

Axdahl, Erik↗

Design of a Lunar Plume-Surface Interaction Measurement System

Plume-surface interactions between a rocket plume and the lunar surface will be studied in-situ during two of NASA’s upcoming Commercial Lunar Payload Services Program missions. The payload, Stereo Cameras for Lunar Plume-Surface Studies (SCALPSS), will employ a multiple-camera photogrammetry system to obtain stereo images of the lunar regolith before, during, and after significant erosion and subsequent crater formation. The evaluation of measurement system capabilities and design process of the payload is informed by computational fluid dynamics predictions, accuracy modeling based on experimental data, camera simulation, and lander design, all of which are combined in the Virtual Diagnostic Interface. The second SCALPSS mission, traveling to the moon on Firefly Aerospace’s Blue Ghost lander in early 2023, aims to build on the design and complexity of the first payload, which is launching in early 2022 on Intuitive Machines’ Nova-C lander. The Blue Ghost SCALPSS mission will include two additional cameras and a total of three different lens focal lengths which will target specific points in the vehicle’s descent to obtain quantitative and accurate 3D reconstruction of the lunar surface both prior to and after crater formation.

photogrammetry↗

Advancing xEMU Lunar Dust Mitigation Devices

Fine, electrically charged, glass like dust particles caused significant damage to the Apollo EMU9 during lunar EVAs, identified as one of the greatest challenges to future exploration. Passive Lunar Dust Mitigation Devices (LDMD) were developed, within a SBIR Phase II, to prohibit this dust from interrupting venting space suit component operation. A Computational Fluid-Dynamics and Discrete Element Method (CFD-DEM) Simulation Tool was developed at the University of Colorado, Boulder to predict venting gas flow ability to self-clean adhered dust particles from LDMD surfaces. Lunar dust properties (i.e., adhesion and cohesion strengths) required to complete Simulation Tool analysis are relatively unknown due to considerable differences between the Earth and the Moon (i.e., gravity, humidity) and due to an absence of dust particles in their native state. Analytically determining gas flow velocity, density and direction within the fluid Boundary Layer, microns from LDMD surfaces presented a second challenge. Dusty Plasma Chamber testing is being performed at Auburn University, Auburn to observe electrostatically charged dust behavior as it adheres to LDMD prototypes and specific geometry and is then blown away by metered gas flow. Test articles were developed to offer insight into the impact of different flow geometries, surface roughness and dust removal within the gas Boundary-Layer. Observed dust behavior is currently being developed to support the CFD-DEM analysis. Many Simulation Tool analytical cases have been processed to support the intention of completing sensitivity studies to assess how different dust adherence values and Boundary Layer fluid properties impact LDMD self-cleaning effectivity.

Thomas J Stapleton↗

Numerical Simulation of a High-Lift Configuration Embedded with High Momentum Fluidic Actuators

Numerical simulations have been performed for a vertical tail configuration with deflected rudder. The suction surface of the main element of this configuration, just upstream of the hinge line, is embedded with an array of 32 fluidic actuators that produce oscillating sweeping jets. Such oscillating jets have been found to be very effective for flow control applications in the past. In the current paper, a high-fidelity computational fluid dynamics (CFD) code known as the PowerFLOW R code is used to simulate the entire flow field associated with this configuration, including the flow inside the actuators. A fully compressible version of the PowerFLOW R code valid for high speed flows is used for the present simulations to accurately represent the transonic flow regimes encountered in the flow field due to the actuators operating at higher mass flow (momentum) rates required to mitigate reverse flow regions on a highly-deflected rudder surface. The computed results for the surface pressure and integrated forces compare favorably with measured data. In addition, numerical solutions predict the correct trends in forces with active flow control compared to the no control case. The effect of varying the rudder deflection angle on integrated forces and surface pressures is also presented.

Vatsa, Veer N.↗

Integrated Heat Exchanger-Phase Change Material Thermal Energy Storage System

The purpose of this study is to experimentally investigate the thermal performance of an innovative thermal energy storage (TES) system that combines the advantages of the phase-change material (PCM)/graphite foam latent heat TES medium developed at Argonne National Laboratory (Argonne) and the internally supported plate-fin (ISPE) cell architecture heat transfer fluid (HTF) flow channels developed at Brayton Energy (Brayton). Several essential tasks were accomplished: (1) Thermal property characterization. Thermal properties of the graphite foam were characterized, providing necessary data for experimental result analysis and numerical simulation. (2) Design and optimization of lab-scale test module. Based on Brayton’s full-scale heat exchanger (HX)-TES system, the experimental test module was designed, optimized, and fabricated. (3) Thermal performance testing and data analysis. Five cycle tests were successfully conducted—including one with approximately 3.5 psig of pressure applied to the diaphragms—to investigate the thermal performance of the experimental test module for charging and discharging. Temperature profiles were generated for each charging test and discharging test as a function of time. The temperature profiles clearly show three TES stages: sensible heat (temperature increase), latent heat (melting), and sensible heat (temperature increase) for the charging process. Similarly, the temperature profiles clearly show three thermal energy release stages: sensible heat (temperature decrease), latent heat (solidification), and sensible heat (temperature decrease). Melting and solidification of the PCM generally occurred in relatively narrow temperature ranges, indicated by the flattened temperature regions in the temperature profiles. These phase changes ranged approximately 3°C for melting and 3.5°C for solidification. The charging and discharging temperature profiles were similar for similar experimental parameter tests whether or not pressure was applied to the diaphragm to eliminate the gap between the HX surface and the TES subsystem. This indicates that the effect of a small gap between the HX surface and the TES subsystem is insignificant for charging and discharging. (4) Comparison of experimental data and simulation results. We compared the experimental data to the numerical simulation results. Numerical simulations were conducted by using the ANSYS FLUENT 2019 R3 commercial computational fluid dynamics software. The predicted phase-change times agreed reasonably well with those from the experimental data. In most cases, the estimated time differences between the relative phase changes were within 16%. The predicted start and end times for the charging process agreed well with those from the experimental data. However, the simulation results showed earlier start and end times than the experimental data for the discharging process. Overall, the experimental data and its comparison with the simulation predictions verified the technical viability of the integrated ISPF HX-PCM/graphite foam latent-heat TES system.

25 ENERGY STORAGE↗

Experience with Aero- and Fluid-Dynamic Testing for Engineering and CFD Validation

Ever since computations have been used to simulate aerodynamics the need to ensure that the computations adequately represent real life has followed. Many experiments have been performed specifically for validation and as computational methods have improved, so have the validation experiments. Validation is also a moving target because computational methods improve requiring validation for the new aspect of flow physics that the computations aim to capture. Concurrently, new measurement techniques are being developed that can help capture more detailed flow features pressure sensitive paint (PSP) and particle image velocimetry (PIV) come to mind. This paper will present various wind-tunnel tests the author has been involved with and how they were used for validation of various kinds of CFD. A particular focus is the application of advanced measurement techniques to flow fields (and geometries) that had proven to be difficult to predict computationally. Many of these difficult flow problems arose from engineering and development problems that needed to be solved for a particular vehicle or research program. In some cases the experiments required to solve the engineering problems were refined to provide valuable CFD validation data in addition to the primary engineering data. All of these experiments have provided physical insight and validation data for a wide range of aerodynamic and acoustic phenomena for vehicles ranging from tractor-trailers to crewed spacecraft.

experiment↗

Turbulence Simulations of Transonic Flows over NACA-0012 and OAT15A Airfoils

Three different simulation approaches, namely unsteady Reynolds-averaged Navier-Stokes (URANS), delayed detached-eddy simulation (DDES), and wall-modeled large-eddy simulation (WMLES) are employed to simulate transonic flow over an NACA-0012, OAT15A airfoils at different angles of attack covering pre- and post-buffet-onset regimes. The freestream Mach number is 0.75, and the Reynolds number based on the chord length is 10 and 3million, respectively. The NASA FUN3D solver is used for the simulations, which is an unstructured, compressible flow solver. The URANS simulations are performed using the Spalart-Allmaras (SA) model with the compressibility correction, the DDES predictions are based on the SA model, and the WMLES are performed using an equilibrium wall-model. The unsteady RANS simulations, only with the compressibility correction, predict the pre- and post- buffet characteristics, which compare well with the experimental results. DDES results predicted a lower buffet onset angle compared to experiment. The predicted shock locations are upstream of the locations predicted by URANS. Using a fine grid in the spanwise direction, WMLES predictions show buffeting consistent with the experiment.

Computational Fluid Dynamics↗

Turbulence Simulations of Transonic Flows over NACA-0012 and OAT15A Airfoils

Three different simulation approaches, namely unsteady Reynolds-averaged Navier-Stokes (URANS), delayed detached-eddy simulation (DDES), and wall-modeled large-eddy simulation (WMLES) are employed to simulate transonic flow over an NACA-0012, OAT15A airfoils at different angles of attack covering pre- and post-buffet-onset regimes. The freestream Mach number is 0.75, and the Reynolds number based on the chord length is 10 and 3million, respectively. The NASA FUN3D solver is used for the simulations, which is an unstructured, compressible flow solver. The URANS simulations are performed using the Spalart-Allmaras (SA) model with the compressibility correction, the DDES predictions are based on the SA model, and the WMLES are performed using an equilibrium wall-model. The unsteady RANS simulations, only with the compressibility correction, predict the pre- and post- buffet characteristics, which compare well with the experimental results. DDES results predicted a lower buffet onset angle compared to experiment. The predicted shock locations are upstream of the locations predicted by URANS. Using a fine grid in the spanwise direction, WMLES predictions show buffeting consistent with the experiment.

Computational Fluid Dynamics↗

CFD Simulation of Floating Body Motion with Mooring Dynamics: Coupling MoorDyn with OpenFOAM

It is increasingly popular to use Computational Fluid Dynamics (CFD) models to study floating structures subjected to ocean waves, especially when it comes to applications of floating offshore wind turbines and Wave Energy Converters (WECs). Mooring dynamics are currently lacking in most of these applications. This paper presents a coupled simulation study of moored floating body motion by coupling two open-source libraries: a finite volume CFD toolbox, OpenFOAM, and a lumped-mass mooring model, MoorDyn. The instantaneous floating body position and velocity are passed from the body motion solver in the CFD model to the mooring model to calculate the fairlead kinematics. The mooring reaction forces, which are calculated by MoorDyn after updating the mooring system states, are then returned to the body motion solver to update the floating body motion. Both mesh deformation and overset mesh methods are used as the mesh motion solver in the CFD model to account for the floating body motion. The coupled model was validated against experimental measurements for a floating box moored with four catenary lines under the action of regular waves, which came from a preliminary test campaign for WECs. Apart from the lumped-mass mooring model, the present work also coupled a quasi-static mooring model and a finite element model with the floating body motion solver in OpenFOAM. The mooring line tensions predicted by these models were compared. The coupled model equipped with three mooring line codes may be further used to carry out survivability studies of FOWTs and WECs subject to severe sea states.

floating body↗

Simulation of Atmospheric-Entry Capsules in the Subsonic Regime

The accuracy of Computational Fluid Dynamics predictions of subsonic capsule aerodynamics is examined by comparison against recent NASA wind-tunnel data at high-Reynolds-number flight conditions. Several aspects of numerical and physical modeling are considered, including inviscid numerical scheme, mesh adaptation, rough-wall modeling, rotation and curvature corrections for eddy-viscosity models, and Detached-Eddy Simulations of the unsteady wake. All of these are considered in isolation against relevant data where possible. The results indicate that an improved predictive capability is developed by considering physics-based approaches and validating the results against flight-relevant experimental data.

Simulation↗

Fluid and gyrokinetic turbulence in open field-line, helical plasmas

Two-fluid Braginskii codes have simulated open-field line turbulence for over a decade, and only recently, it has become possible to study these systems with continuum gyrokinetic codes. This work presents a first-of-its-kind comparison between fluid and (long-wavelength) gyrokinetic models in open field-lines, using the GDB and Gkeyll codes to simulate interchange turbulence in the Helimak device at the University of Texas [T. N. Bernard et al., Phys. Plasmas 26, 042301 (2019)]. Partial agreement is attained in a number of diagnostic channels when the GDB sources and sheath boundary conditions (BCs) are selected carefully, especially the heat-flux BCs that can drastically alter the temperature. The radial profile of the fluctuation levels is qualitatively similar and quantitatively comparable on the low-field side, although statistics such as moments of the probability density function and the high-frequency spectrum show greater differences. This comparison indicates areas for future improvement in both simulations, such as sheath BCs, and improvements in GDB like particle conservation and spatially varying thermal conductivity, in order to achieve better fluid-gyrokinetic agreement and increase fidelity when simulating experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Computational fluid dynamics modeling and analysis of silica nanoparticle synthesis in a flame spray pyrolysis reactor

Flame Spray Pyrolysis (FSP) is a method for large-scale production of nanoparticles and nanoscale powders employed in a wide range of industrial applications. Particle size and morphology are complex functions of the physicochemical phenomena occurring in the FSP reactor. An extensive study of FSP-related phenomena can be utilized to develop effective strategies for achieving desired particle size/morphology and scaling up the overall yield of an FSP system. In this work, a computational fluid dynamics (CFD) model of an FSP reactor is developed to simulate the coupling of key phenomena involved in the particle synthesis process: liquid spray breakup and evaporation, mixing, combustion, and particle formation/growth of silica nanoparticles. Herein, the particle sizes and their distributions from the CFD simulations are validated against experimental data. Subsequently, the simulations are utilized to investigate the impact of process parameters on the resultant flame dynamics and particle growth. Firstly, the CFD results show that the particle sizes are strongly correlated with the precursor concentration in the solvent. At lower precursor concentrations, the spread of the distribution is relatively insensitive to the value of the concentration. At higher concentrations, the spread is higher as the collision probability between particles is higher. Secondly, increasing the pilot flow rate increases the length of the pilot flames impacting the local ignition location of the spray flame. Lastly, it is shown that the dispersion gas flow rate strongly influences the spray flame shape. This shape can be used for control of particle growth as it helps determine the regions of high temperature and the residence time of the particles in the high temperature region enabling the design and process optimization of the FSP reactor.

42 ENGINEERING↗

Hardness Prediction by Incorporating Heat Transfer and Molten Pool Fluid Flow in a Multi-pass, Multi-layer Weld for Onsite Repair of Grade 91 steel

Introduction: In the current fleet of fossil-fired power plants, creep strength enhanced ferritic steels (CSEF) are used to sustain the harsh service conditions. Enhanced properties of Grade 91 steel result from tempered martensite with a fine distribution of MX and M23C6 carbides. Grade 91 steel is subjected to onsite welding repair to remedy their degradation due to extreme service condition. Knowledge of weld repairability of these steels, such as as-welded hardness distribution, is essential to establishing sound repair procedures. Experimental trial and error tests can consume a lot of time as many welding variables need to be studied. For numerical modelling, most of the multi-pass multi-layer models are based on finite element method, which are limited to solve the heat conduction equation and ignore convective heat transfer due to melt flow. Moreover, the mesh has to be pre-built based on a known or assumed weld cross-section geometry. These finite element based models thus have limited predictive capability as defects are not considered and nugget size are pre-assumed. This research aims at developing a thermal and microstructure evolution model incorporating molten pool dynamics in a multi-pass multi-layer material deposition to predict the as-welded hardness distribution. Technical Approach: All the thermal, physical, and metallurgical properties of Grade 91 as a function of temperature are collected from the literature and inputted into the thermo-fluid model based on Flow-3D, a computational fluid dynamics software. A multi-pass, multi-layer material deposition is simulated where the melting of filler wire into the molten pool is directly considered based on the volume of fluid (VOF) method. The flow behaviour of the molten pool is used to understand the formation of deposition geometry and defects. The temperature profiles during the multi-pass, multi-layer welding are calculated. The results computed using the new model are compared against the experimental data of fusion zone geometry and thermal cycles. Hardness prediction in the heat-affected zone (HAZ) are made using Johnson-Mehl-Avrami (JMA) equation for solid-state phase transformation kinetics. The JMA parameters are extracted from the experimental data available in the literature. For comparison, a standard finite element heat conduction model is also developed to predict the thermal cycles and hardness distribution in the multi-pass, multi-layer weld. Expected Result: Results obtained using the molten pool dynamic simulation versus the finite element heat conduction model are compared. Specifically, the effects of convective heat transfer on the accuracy of the calculated thermal history, bead shape and size, and HAZ hardness distribution are examined.

Modeling, Grade 91 steel, Multi-pass multi-layer d↗

Verification methods for drift–diffusion reaction models for plasma simulations

Abstract Compared to other computational physics areas such as codes for general computational fluid dynamics, the documentation of verification methods for plasma fluid codes remains under developed. Current analytical solutions for plasma are often highly limited in terms of testing highly coupled physics, due to the harsh assumptions needed to derive even simple plasma equations. This work highlights these limitations, suggesting the method of manufactured solutions (MMSs) as a potential option for future verification efforts. To demonstrate the flexibility of MMS in verifying these highly coupled systems, the Multiphysics Object-Oriented Simulation Environment (MOOSE) framework was utilized. Thanks to the MOOSE framework’s robustness and modularity, as well as to its physics module capabilities and ecosystem applications (i.e. Zapdos and the chemical reaction network) developed for plasma physics modeling and simulation, this report lays the groundwork for a structured method of conducting plasma fluid code verification.

DeChant, Corey (ORCID:0000000293779732)↗

TRUST Sensors in Environments: Thermocouples (SE-TC), Release FY25

The Delivery Environments Testbeds to Reduce Uncertainty in Simulations and Tests (TRUST) project is a broad project intended to analyze simplified problems experimentally and with modeling and simulation. The purpose of analyzing these simplified problems is to extend solution methods to more complex problems, as well as understand deficiencies and gaps in knowledge of methods currently used in more complex analyses. The TRUST project encompasses several smaller testbeds intended to isolate individual phenomena. The testbed under consideration in this report is the Sensors in Environments: Thermocouples testbed. In previous years, the purpose of this testbed was to quantify uncertainty of thermocouple sensors. To accomplish this, an aluminum plate was placed in a thermal chamber and subject to various types of thermal loading. Thermocouples were placed in various locations on the aluminum plate in various configurations (e.g., embedded in the plate, placed under Kapton tape), and an effort was made to quantify uncertainty in these measurements. Finite element simulations were performed to investigate how sensitive these measurements were to parameters such as the boundary conditions on the plate and material properties. However, a fundamental source of uncertainty in this analysis was the convective heat transfer from the plate. Convective heat transfer is a complex physical phenomenon comprised of a number of interacting sub-processes and is difficult to predict accurately a priori. As such, the main purpose of this testbed in FY25 was to better understand, both experimentally and numerically, the convective heat transfer from the plate. This is a highly applicable problem to several more complex problems, as convective heat transfer occurs in nearly all problems where a body is moving through air. Numerically, this required a two-step approach. First, the air flow in the thermal chamber was in vestigated using computational fluid dynamics. The commercial solver Fluent was used to perform these simulations. From these simulations, a heat transfer coefficient over the surface of the plate was calculated. This heat transfer was then used as boundary conditions for finite element heat transfer simulations within the plate, which were performed using Abaqus. Significant effort was devoted to automating the handoff between these two solvers. Experimentally, previous thermocouple results in the plate were used to validate the time-dependent thermal profiles produced from Abaqus. Further experimental efforts were performed both to help validate the Fluent simulations and to inform its boundary conditions. For example, hot-wire anemometers were used to measure the velocity in the chamber, which would be particularly useful in understanding the chamber inlet velocity. Thermocouple measurements were also taken in the chamber, instead of only on the plate, to serve as validation evidence for the Fluent simulations. Numerical results showed that the Fluent to Abaqus workflow matched previous plate thermocouple measurements well. This type of handoff is useful for more complex experiments, or those that are not able to be examined in as great of detail as this testbed, as it was performed without any experimental input. Experimental results, however, were more mixed. The anemometers proved unreliable, with inconsistent measurements across all anemometers, even at locations that were nearly identical. On the other hand, the thermocouples provided a relatively rich view of the temperature field in the chamber.

42 ENGINEERING↗