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 775 records · Page 43
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Design Optimization of A-Type Heat Exchangers based on High Performance, Non-Round Tubes

Air-to-refrigerant heat exchangers (HXs) have been the topic of considerable research as they are the fundamental heat transfer components of HVAC&R systems. For residential and commercial applications, indoor heat exchangers are often in A-type configurations to reduce air-conditioner system footprint while still delivering the required cooling or heating. Classical HX design practices typically assume a uniform frontal air velocity profile to facilitate the use of conventional tube-fin airside heat transfer and pressure drop correlations for performance prediction. However, it is well-reported throughout the literature that A-type HXs experience significant airflow maldistribution which can severely degrade HX performance. Additionally, recent advancements in simulation software such as Computational Fluid Dynamics (CFD), Finite Element Analysis (FEA) and optimization algorithms have led researchers to consider primary heat transfer surface optimization to achieve highly compact HXs which do not require fins. In this paper, an A-type HX which utilizes shape-optimized non-round tubes is optimized for minimum duct size and airside pressure drop. This is achieved through coupling of an experimentally validated finite volume HX model with automated 2D CFD simulations of airflow through the HX for airside thermal-hydraulic performance evaluation. Preliminary results show that the optimal designs deliver similar capacity while achieving 20% reductions in airside pressure drop and duct cross-section area. Additionally, the optimal design apex angles in this study were similar to those from previous A-Type HX optimization studies for conventional tube-fin HX geometries, suggesting that the apex angle for optimal A-Type configuration performance may be independent of tube geometry.

42 ENGINEERING↗
NASA NTRS2007

Numerical Simulations of Single Flow Element in a Nuclear Thermal Thrust Chamber

The objective of this effort is to develop an efficient and accurate computational methodology to predict both detailed and global thermo-fluid environments of a single now element in a hypothetical solid-core nuclear thermal thrust chamber assembly, Several numerical and multi-physics thermo-fluid models, such as chemical reactions, turbulence, conjugate heat transfer, porosity, and power generation, were incorporated into an unstructured-grid, pressure-based computational fluid dynamics solver. The numerical simulations of a single now element provide a detailed thermo-fluid environment for thermal stress estimation and insight for possible occurrence of mid-section corrosion. In addition, detailed conjugate heat transfer simulations were employed to develop the porosity models for efficient pressure drop and thermal load calculations.

Cheng, Gary↗
NASA NTRS2008

Feasibility of Conducting J-2X Engine Testing at the Glenn Research Center Plum Brook Station B-2 Facility

A trade study of the feasibility of conducting J-2X testing in the Glenn Research Center (GRC) Plum Brook Station (PBS) B-2 facility was initiated in May 2006 with results available in October 2006. The Propulsion Test Integration Group (PTIG) led the study with support from Marshall Space Flight Center (MSFC) and Jacobs Sverdrup Engineering. The primary focus of the trade study was on facility design concepts and their capability to satisfy the J-2X altitude simulation test requirements. The propulsion systems tested in the B-2 facility were in the 30,000-pound (30K) thrust class. The J-2X thrust is approximately 10 times larger. Therefore, concepts significantly different from the current configuration are necessary for the diffuser, spray chamber subsystems, and cooling water. Steam exhaust condensation in the spray chamber is judged to be the key risk consideration relative to acceptable spray chamber pressure. Further assessment via computational fluid dynamics (CFD) and other simulation capabilities (e.g. methodology for anchoring predictions with actual test data and subscale testing to support investigation.

Schafer, Charles F.↗
DOE OSTI2020

Optimization of a Cyclone Using Multiphase Flow Computational Fluid Dynamics

The U.S. Department of Energy National Energy Technology Laboratory's (NETL) 50 kW th chemical looping reactor (CLR) has an underperforming cyclone, which was designed using empirical correlations. To improve the performance of this cyclone using computational fluid dynamics (CFD)-based modeling simulations, four critical design parameters including the vortex tube radius and length, barrel radius, and the inlet width and height were optimized. Here, NETL's open source multiphase flow with interphase exchange (MFiX) CFD code has been used to model a series of cyclones by systematically varying the geometric design parameters. To perform the optimization process, the surrogate modeling and sensitivity analysis followed by the optimization capability in nodeworks was used. The basic methodology for the process is to employ a statistical design of experiments (DOE) method to generate sampling simulations that fill the design space. Corresponding CFD models are then created, executed, and postprocessed. A response surface is created to characterize the relationship between input parameters and the quantities of interest (QoI). Finally, the CFD-surrogate is used by an optimization method to find the optimal design condition based on the objective and constraints prescribed. The resulting optimal cyclone has a larger diameter and longer vortex tube, a larger diameter barrel, and a taller and narrower solids inlet. The improved design has a predicted pressure drop 11 times lower than the original design while reducing the mass loss by a factor of 2.3.

42 ENGINEERING↗
NASA NTRS2008

An Experimental and Numerical Study of a Supersonic Burner for CFD Model Development

A laboratory scale supersonic burner has been developed for validation of computational fluid dynamics models. Detailed numerical simulations were performed for the flow inside the combustor, and coupled with finite element thermal analysis to obtain more accurate outflow conditions. A database of nozzle exit profiles for a wide range of conditions of interest was generated to be used as boundary conditions for simulation of the external jet, or for validation of non-intrusive measurement techniques. A set of experiments was performed to validate the numerical results. In particular, temperature measurements obtained by using an infrared camera show that the computed heat transfer was larger than the measured value. Relaminarization in the convergent part of the nozzle was found to be responsible for this discrepancy, and further numerical simulations sustained this conclusion.

Magnotti, G.↗
DOE OSTI2022

An Investigation into the Effects of Swirl on the Performance and Emissions of an Opposed-Piston Two-Stroke Engine using Large Eddy Simulations

Opposed-piston two-stroke (OP-2S) engines have the potential to achieve higher thermal efficiency than a conventional four-stroke diesel engine. However, the uniflow scavenging process is difficult to control over a wider range of speed and loads due to its sensitivity to pressure dynamics, port timings, and port design. Specifically, the angle of the intake ports can be used to generate swirl which has implications for open and closed cycle effects. This study proposes an analysis of the effects of port angle on the in-cylinder flow distribution and combustion performance of an OP-2S using computational fluid dynamics engine. Large Eddy Simulation (LES) was used to model turbulence given its ability to predict in-cylinder mixing and cyclic variability. A three-cylinder model was validated to experimental data collected by Achates Power and the grid was verified using an LES quality approach from the literature. The model was used to simulate port angles from 12 to 29 degrees at constant pressure and temperature boundary conditions. Results indicated that the higher bulk swirl ratio generated by larger port angles tends to trap more internal residuals. This effect on the scavenging performance, combined with the larger trapped swirl ratio, also has a significant impact on the combustion performance in a two-stroke engine. It was concluded that there exists a tradeoff of efficiency and emissions that must be considered when increasing the port angle in a uniflow two-stroke engine. 1. Introduction Two-stroke engines were developed in the 1800s and are commonly used in practice today for lightweight power applications such as motorcycles or handheld power tools, as well as large-bore low speed engines for stationary power and marine applications. Two-stroke engines provide the distinct advantage of higher torque/power output but present challenges due to its coupled intake and exhaust process. Four-stroke engines have the benefit of using two full piston strokes to induct intake charge and to expel exhaust gases. Controlling the trapped charge composition in a twostroke engine via a scavenging process is crucial for reliable and efficient engine operation.

33 ADVANCED PROPULSION SYSTEMS↗
NASA NTRS1992

Multi-phase simulations of coaxial injector combustion

A multiphase computational fluid dynamics code (ARICC-3D) is presented and results of two simulations are discussed. The numerical framework of the CFD code is reviewed as well as some of the two-phase physical submodels. The simulations performed include a single coaxial element injector and a multielement injector using LOX/Hydrogen reactants. The single element injector simulation verified the interaction among the code's submodels. The multielement injector simulation transient results include the chamber response to a transverse pressure wave with and without a chamber baffle. The results of these simulations demonstrate the current capabilities and their limitations to model complex two-phase combustion phenomena. Possible ways to exceed these limitations are suggested.

Liang, P. Y.↗
NASA NTRSDate not supplied

High-Fidelity CFD Verification Workshop 2024 Summary: Spalart-Allmaras QCR2000-R Turbulence Model

This paper summarizes solutions submitted for the Reynolds-averaged Navier-Stokes (RANS) test suite of the High-Fidelity CFD Verification Workshop. The goal of the workshop is to establish standards for verification of computational fluid dynamics (CFD) approaches to simulations of steady and unsteady turbulent flows. The RANS verification studies focus on a one-equation Spalart-Allmaras model with quadratic constitutive relation and rotation correction, SA-neg-QCR2000-R. The verification test cases are a two-dimensional subsonic flow around a Joukowski airfoil, a three-dimensional subsonic flow around an extruded NACA 0012 wing in a tunnel, and a subsonic flow around a wing-body configuration developed for verification of solvers participating in the 5 𝑡 ℎ High-Lift Prediction Workshop. The turbulencemodel formulation, geometry, flow conditions, grids, and reference solutions are described in detail. Solutions for the test cases are computed by seven established CFD solvers on adaptedand fixed-grid families using different discretization approaches. While some noticeable differences between solutions remain, the results achieved by contributing solvers show that different solutions computed for the same RANS model on different grid families can converge to a common limit with grid refinement. The apparent requirements for grid convergence are a well designed family of grids that provide sufficient resolution in important areas and a strong solver capable of deep iterative convergence on each grid. For each test case in the study, the variation between aerodynamic forces computed by different solvers on the finest grids of different families is less than 2%.

Boris Diskin↗
NASA NTRSDate not supplied

Parametric Mechanism Design Through Numerical Optimization and Physics Simulation

Design-Build-Test approaches for developing spaceflight hardware are prohibitively time and cost intensive and often lead to suboptimal mechanism designs. Approaches that couple machine learning and high-fidelity physics simulation could eliminate the need for hardware prototyping and dramatically accelerate the engineering design cycle, ultimately reducing cost. This work presents a modular NASA-developed toolchain to optimize hardware mechanisms in a virtual environment using numerical optimization and multi-body physics simulation. The toolchain enables multi-objective optimization, generates parametric CAD files that can be further post-processed by an end user, and can be expanded to optimize full systems and non-mechanical parameters such as feedback control variables. We demonstrate the toolchain through an independently verifiable design problem that optimizes wheel radius to achieve a desired linear velocity in a rigid-body physics environment when the wheel rotates at a constant angular speed, and then post-process the parametric CAD file of the optimal design generated by the tool before ultimately manufacturing it via 3D printing. We end with a discussion of how the toolchain can incorporate other analysis tools, including finite element analysis, computational fluid dynamics, and granular media simulations.

Optimization↗
NASA NTRSDate not supplied

An Optimization-Based Toolchain for Parametric Mechanism Design

Design-Build-Test approaches for developing spaceflight hardware are prohibitively time and cost intensive and often lead to suboptimal mechanism designs. Approaches that couple machine learning and high-fidelity physics simulation could eliminate the need for hardware prototyping and dramatically accelerate the engineering design cycle, ultimately reducing cost. This work presents a modular NASA-developed toolchain to optimize hardware mechanisms in a virtual environment using numerical optimization and multi-body physics simulation. The toolchain enables multi-objective optimization, generates parametric CAD files that can be further post-processed by an end user, and can be expanded to optimize full systems and non-mechanical parameters such as feedback control variables. We demonstrate the toolchain through an independently verifiable design problem that optimizes wheel radius to achieve a desired linear velocity in a rigid-body physics environment when the wheel rotates at a constant angular speed, and then post-process the parametric CAD file of the optimal design generated by the tool before ultimately manufacturing it via 3D printing. We end with a discussion of how the toolchain can incorporate other analysis tools, including finite element analysis, computational fluid dynamics, and granular media simulations.

Optimization↗
NASA NTRSDate not supplied

Modeling Self-Pressurization and Spray Bar Pressure Control of a Cryogenic Storage Tank in Normal Gravity

This paper presents Computational Fluid Dynamics (CFD) models for simulating self-pressurization and spray bar cooling processes in a large-scale liquid hydrogen storage tank under normal gravity. The self-pressurization model utilizes the kinetics-based Schrage equation with the Volume-Of-Fluid (VOF) model to account for the evaporative and condensing interfacial mass fluxes. It uses a laminar approach to model natural convection and compares the pressures and temperatures predicted by the CFD model with experimental data from the Multipurpose Hydrogen Test Bed (MHTB) self-pressurization experiment. Additionally, the predicted interfacial mass flux is presented. A CFD model for simulating pressure control using a spray bar is also presented. An Eulerian-Lagrangian approach models the interactions between the discrete droplets and the continuous ullage phase. The spray model is coupled with the VOF model by tracking particles in the ullage, removing particles from the ullage when they reach the interface, and then adding their contributions to the liquid. The T-sat model is used to calculate the droplet-ullage heat and mass transfer. In this model, droplets warm up to the saturation temperature corresponding to the ullage vapor pressure thereafter evaporating while remaining at the saturation temperature. The tank pressure, vapor, and liquid temperature evolutions are validated against data from the MHTB spray bar mixing experiment.

Self-Pressurization↗
DOE OSTI2023

Optimal gas cooling with minimal aero-opticaldistortion for next generation high-average-power lasers (Progress report, Year 2)

Next-generation lasers that combine high-peak-power ultrashort pulses with high repetition rates have the potential to efficiently deliver the high-average-power outputs required for inertial fusion energy generation, laser-based particle acceleration, and to generate compact secondary x- and 𝛾-ray radiation sources, among others. The combination of long gain length with high surface-to-volume ratio to remove residual heat makes gas-cooled multi-slab amplifier heads an inherently aperture- and average power scalable technology that enables this new class of solid-state lasers, and their thermal management becomes the most pressing technical challenge. As the index-of refraction depends on the gas density, heat removal from the slab surfaces unavoidably entails scattering of propagated light from regions of high temperature fluctuations. The requirement to maximize convective heat extraction from the gain medium under the constraint of minimum optical aberration of the beam renders amplifier slab cooling a technological challenge at the intersection of aero-optics, fluid-mechanics and heat transfer. The goal of this research is to leverage computational fluid dynamics and aero optics simulations of different levels of fidelity to devise a numerical modeling and optimization framework for gas-cooled multi-slab amplifiers. The project will be conducted in close collaboration with the Advanced Photon Technologies Group at Lawrence Livermore National Laboratory, where a gas-cooled test stand for laser materials and thermal management structures will be developed and installed. The final goal of the collaborative effort is to manufacture and benchmark a physical prototype of an numerically optimized amplifier head design. Measurement data on solid-state and thermo-optical distortions will be used for the validation of the numerical models. The key technical challenge in this multi-objective design problem is that of balancing optimal heat removal with minimal optical distortion. To facilitate computationally efficient design optimization and, at the same time, deepen our basic understanding of the interaction of light with turbulence, a complementary study that combines high and low-fidelity flow simulations and aero-optical models is proposed. Along the lines of these two objectives and solution strategies, the project is organized into two thrusts.

42 ENGINEERING↗
DOE OSTI2021

Direct numerical simulation of supercritical oxy-methane mixing layers with CO 2 substituted counterparts

Direct numerical simulations (DNS) of temporally developing, three dimensional, CH 4 /CO 2 , CH 4 /O 2 and CO 2 /O 2 mixing layers, are conducted at a supercritical pressure of 300 atm. To effectively model the supercritical regime, the employed formulation includes the compressible form of the governing equations, the cubic Peng-Robinson equation of state and a generalized formulation for heat and mass flux vectors derived from non-equilibrium thermodynamics and fluctuation theory. A linear inviscid stability analysis is also performed for each case, to determine its most unstable wavelength. Flow visualizations reveal the presence of high density gradient magnitude regions for all three mixing layers, with conditional averages indicating increased presence of heavier fluid species within these regions. No significant departures are observed from perfect gas behavior, with compressibility factors very close to unity for all three mixing cases. Applicability of presumed probability density function methods (PDF) is examined for the three supercritical mixing layers. An a priori analysis is also conducted to investigate various simplifying assumptions employed in modeling various subgrid scale (SGS) flux models. Two additional terms are identified in the large eddy simulations (LES) equations, the gradient of SGS contribution of pressure in the momentum equation and the gradient of SGS contribution of heat flux in energy equation, whose magnitudes are similar and comparable with their respective resolved terms. The performance of the scale similarity model to represent these additional terms is investigated. Lastly, the performance of Smagorinsky, gradient and scale similarity models are also investigated.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS2015

Modeling of Non-Isothermal Cryogenic Fluid Sloshing

A computational fluid dynamic model was used to simulate the thermal destratification in an upright self-pressurized cryostat approximately half-filled with liquid nitrogen and subjected to forced sinusoidal lateral shaking. A full three-dimensional computational grid was used to model the tank dynamics, fluid flow and thermodynamics using the ANSYS Fluent code. A non-inertial grid was used which required the addition of momentum and energy source terms to account for the inertial forces, energy transfer and wall reaction forces produced by the shaken tank. The kinetics-based Schrage mass transfer model provided the interfacial mass transfer due to evaporation and condensation at the sloshing interface. The dynamic behavior of the sloshing interface, its amplitude and transition to different wave modes, provided insight into the fluid process at the interface. The tank pressure evolution and temperature profiles compared relatively well with the shaken cryostat experimental test data provided by the Centre National D'Etudes Spatiales.

Cryogenics↗
DOE OSTI2020

A Comparison of Interface Growth Models Applied to Rayleigh–Taylor and Richtmyer–Meshkov Instabilities

Sophisticated numerical models that contain fluid interfaces rely upon interface evolution models to approximate the transition to turbulence near interfaces, in the presence of Rayleigh–Taylor (RTI) or Richtmyer–Meshkov (RMI) instability. Semi-analytical models have been developed in recent decades to predict the interface growth from an initial state into the nonlinear regime. Two of these models are considered in this study. They are the Goncharov and the z-models. Both of these interface models have strengths and weaknesses, which are examined here. Both of them have been implemented in the xRAGE compressible flow solver, which models fluid interfaces. The flow solver provides the fluids acceleration as a body force to the interface model. The purpose of such interface model is to evolve the early times of interface position as a subgrid model within a compressible flow simulation in order to then initialize a turbulence model. In this work, the interface models are assessed and compared for their evolution of RTI and RMI. The z-model performed better than the Goncharov model for all cases that were explored.

42 ENGINEERING↗
NASA NTRSDate not supplied

OVERFLOW Analysis of Supersonic Retropropulsion Testing on a Blunt Mars Entry Vehicle Concept

Supersonic retropropulsion (SRP) flow over a Hypersonic Inflatable Aerodynamic Decelerator (HIAD) blunt-body vehicle was simulated using the Overflow Computational Fluid Dynamics (CFD) solver. Two nozzle configurations were tested (1E and 1F) as a subset of the seven total configurations in the Descent System Study (DSS) testing campaign. A generalized, Adaptive Mesh Refinement (AMR) shock capturing and plume refinement technique was developed and calibrated for producing automatic and unique grid systems optimized for any given testing condition. Solution independence from grid resolution was determined by successively increasing grid refinement until asymptotic convergence of the mean aerodynamic loads was observed. Overflow Reynolds-averaged Navier-Stokes (RANS) solutions produced realistic SRP flow phenomena, including Mach disk normal shocks contained in high-thrust plumes and bow shock triple-points. Solutions for the 1E nozzle configuration were steady. Conversely, a subset of the 1F test conditions were unsteady and exhibited periodic, non-sinusoidal expansions and contractions of the streamwise-oriented plume, resulting in unsteady aerodynamic loads on the vehicle. Initial investigations into the effect of turbulence modeling fidelity demonstrated significant differences between the aerodynamic loads simulated with RANS and Detached Eddy Simulation (DES). Multiple flow mechanisms were identified as root causes of these differences, including bow shock shape augmentation and reduced entrainment due to reduced turbulence in the DES simulations. Initial comparisons of pitching moment control authority and aerodynamic drag performance between the tested nozzle configurations demonstrated that the 1E configuration may be more advantageous for the Mars Entry, Descent, and Landing (EDL) task.

Supersonic Retropropulsion↗
NASA NTRSDate not supplied

OVERFLOW Analysis of Supersonic Retropropulsion Testing on a Blunt Mars Entry Vehicle Concept

Supersonic retropropulsion (SRP) flow over a Hypersonic Inflatable Aerodynamic Decelerator (HIAD) blunt-body vehicle was simulated using the Overflow Computational Fluid Dynamics (CFD) solver. Two nozzle configurations were tested (1E and 1F) as a subset of the seven total configurations in the Descent System Study (DSS) testing campaign. A generalized, Adaptive Mesh Refinement (AMR) shock capturing and plume refinement technique was developed and calibrated for producing automatic and unique grid systems optimized for any given testing condition. Solution independence from grid resolution was determined by successively increasing grid refinement until asymptotic convergence of the mean aerodynamic loads was observed. Overflow Reynolds-averaged Navier-Stokes (RANS) solutions produced realistic SRP flow phenomena, including Mach disk normal shocks contained in high-thrust plumes and bow shock triple-points. Solutions for the 1E nozzle configuration were steady. Conversely, a subset of the 1F test conditions were unsteady and exhibited periodic, non-sinusoidal expansions and contractions of the streamwise-oriented plume, resulting in unsteady aerodynamic loads on the vehicle. Initial investigations into the effect of turbulence modeling fidelity demonstrated significant differences between the aerodynamic loads simulated with RANS and Detached Eddy Simulation (DES). Multiple flow mechanisms were identified as root causes of these differences, including bow shock shape augmentation and reduced entrainment due to reduced turbulence in the DES simulations. Initial comparisons of pitching moment control authority and aerodynamic drag performance between the tested nozzle configurations demonstrated that the 1E configuration may be more advantageous for the Mars Entry, Descent, and Landing (EDL) task.

Supersonic Retropropulsion↗
DOE OSTI2024

LES simulations of a vacuum membrane distillation channel with geometric alterations

3D LES simulations were carried out to study the performance of a vacuum membrane distillation module. A wiggly wall profile with/without embedded stiffeners was considered to alleviate axial and radial temperature polarization, the cause of performance decrease in membrane distillation. Results of the flow field show the wiggles create unsteady vortex shedding inducing intense mixing in the channel. Vortex shedding intensity increases as the Reynolds number increases or stiffeners are added. 3D results show a 56% improvement in flux when moving from a flat sheet membrane to a wiggly membrane with stiffeners at a constant Reynolds number, corresponding to an improvement from 11.7 to 41.7 in the Nusselt number, showing that the alterations improved the flux performance by enhancing the heat transfer along the membrane surface and therefore alleviating temperature polarization, a critical bottleneck in membrane distillation systems. Further, a merit criterion was defined based on the Nusselt number and friction factor, and a 40% increase in merit was shown switching from a flat to a wiggly channel, while a 97% merit increase was seen going from flat to wiggly with stiffeners. A 42% enhancement in the flux was also seen moving from a straight channel to a wiggly channel at a higher Reynolds number which highlights the importance of strategically choosing the mass flow rate, as well as inducing flow separation and vortex shedding in the channel to promote mixing and dissipate the thermal boundary layer. The variation in results between the wiggly module with stiffeners and the wiggly module at high Reynolds numbers suggests enhancing mixing structures can be more impactful on flux for this geometry than increasing the flow rate to a turbulent/transitional regime. However, both are preferable for peak system flux performance. Furthermore, a 2D approximation was used to perform simulations on more extended channels to examine the length degradation. The modules with wiggly channels performed at the same flux level with a doubling of the length. On the other hand, the flat sheet modules experienced length degradation by temperature polarization and dropped in flux yield by around 13% of the short-channel value. This work illustrates that modeling the system and understanding how the performance decreases as the membrane surface area increases are critical for a larger module (scaling up from a lab to a prototype module) to maintain high flux performance.

42 ENGINEERING↗