Search NASA⌕ Search

SEARCH · Search NASA

Results for “Numerical fluid dynamics”

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 811 records · Page 45
All sourcesNASA NTRSDOE OSTI
NASA NTRS1976

A block iterative finite element algorithm for numerical solution of the steady-state, compressible Navier-Stokes equations

An iterative method for numerically solving the time independent Navier-Stokes equations for viscous compressible flows is presented. The method is based upon partial application of the Gauss-Seidel principle in block form to the systems of nonlinear algebraic equations which arise in construction of finite element (Galerkin) models approximating solutions of fluid dynamic problems. The C deg-cubic element on triangles is employed for function approximation. Computational results for a free shear flow at Re = 1,000 indicate significant achievement of economy in iterative convergence rate over finite element and finite difference models which employ the customary time dependent equations and asymptotic time marching procedure to steady solution. Numerical results are in excellent agreement with those obtained for the same test problem employing time marching finite element and finite difference solution techniques.

Cooke, C. H.↗
DOE OSTI2025

Multiphase computational fluid dynamics modeling of reacting flows in absorption columns for carbon capture

First-principles derived computational fluid dynamics (CFD) simulations have been proposed as a fundamental tool for investigating solvent-based CO 2 absorption in packed columns due to their ability to accurately represent the underlying nonlinear, multiscale dynamics. Numerous studies have previously utilized such CFD simulations to investigate hydrodynamics of columns with structured and random packings by assessing the key hydrodynamic metrics such as the interfacial and wetted areas. While mapping such metrics for different conditions is essential to the optimization of absorption columns, it is not sufficient, as the CO 2 capture rate depends also on the coupled, nonlinear dynamics from the underlying chemical reaction kinetics, thermodynamics, and heat-transfer rates. In this work, we present detailed CFD simulation results augmented by incorporating the effects of interfacial physical mass transfer of CO 2 , heat release from chemical reaction kinetics, and thermophysical property variations from resulting temperature gradients. We demonstrate the applicability of the proposed approach in numerically assessing the performance of packed columns by evaluating key hydrodynamic quantities, CO 2 absorption rates, and temperature rise in a reference column with packings that are structurally similar to the Sulzer Mellapak™ 250.Y packing, for different solvent inflow velocities and temperatures. Predictions from simulation results are found to be consistent with the trends in experimental observations from the literature, suggesting that the predictive capabilities of the simulation framework can be leveraged to guide the future development of absorber-column designs and optimized process flowsheets.

Absorption columns↗
DOE OSTI2022

The role of breakup and coalescence in fine-scale bubble-induced turbulence. I. Dynamics

We study the effect of bubble breakup and coalescence on fine-scale dynamics of bubbly turbulent flows using direct numerical simulations. We perform two different simulations of dilute bubbly flows of void fraction 0.5%: one with bubbles breaking up and coalescing and the other without these physical processes. The volume of the fluid method is used for simulating bubbles undergoing breakup and coalescence while the bubbles are treated as rigid spheres in the immersed boundary method simulation. The energy spectrum in both types of simulation, consistent with previous studies, exhibits a -3 slope. We follow a single infinitesimal fluid element as it evolves to understand velocity gradient dynamics using conditional mean trajectories. We note finite-time divergence when the fluid element evolves under the action of inertial and pressure dynamics. The inertial, pressure, and viscous velocity gradient dynamics, when considered individually, produce the same results for bubble-induced turbulence (BIT) as with the classical homogeneous isotropic turbulence (HIT). Yet when the overall velocity gradient dynamics is considered, BIT results in non-cyclic trajectories moving toward stable node and unstable saddle while classical HIT shows cyclic behavior in their trajectories that move toward the origin. Interestingly, both the volume of fluid and immersed boundary simulations produce similar results. Therefore, there are two main takeaways from this research. First, new velocity gradient models are needed for BIT as their velocity gradients behavior is entirely different from the HIT. Second, we can neglect the bubble topology, breakup, and coalescence while studying or modeling the fine-scale dynamics of BIT.

42 ENGINEERING↗
NASA NTRS1989

Numerical Solution Of Navier-Stokes Equations

Three-dimensional, incompressible flow simulated on computer. Report discusses mathematical basis of IN3SD computer code and application of code to several test problems. Relies heavily on use of pseudocompressibility to solve numerically three-dimensional Navier-Stokes equations of viscous, incompressible flow. With further study and some modification, valuable computational tool for solution of engineering problems in fluid dynamics.

Rogers, Stuart E.↗
NASA NTRS1989

A CFD study of tilt rotor flowfields

The download on the wing produced by the rotor wake of a tilt rotor vehicle in hover is of major concern because of its severe impact on payload-carrying capability. In a concerted effort to understand the fundamental fluid dynamics that cause this download, and to help find ways to reduce it, computational fluid dynamics (CFD) is employed to study this problem. The thin-layer Navier-Stokes equations are used to describe the flow, and an implicit, finite difference numerical algorithm is the method of solution. The methodology is developed to analyze the tilt rotor flowfield. Included are discussions of computations of an airfoil and wing in freestream flows at -90 degrees, a rotor alone, and wing/rotor interaction in two and three dimensions. Preliminary results demonstrate the feasibility and great potential of the present approach. Recommendations are made for both near-term and far-term improvements to the method.

Fejtek, Ian↗
NASA NTRS1991

Computational Fluid Dynamics Symposium on Aeropropulsion

Recognizing the considerable advances that have been made in computational fluid dynamics, the Internal Fluid Mechanics Division of NASA Lewis Research Center sponsored this symposium with the objective of providing a forum for exchanging information regarding recent developments in numerical methods, physical and chemical modeling, and applications. This conference publication is a compilation of 4 invited and 34 contributed papers presented in six sessions: algorithms one and two, turbomachinery, turbulence, components application, and combustors. Topics include numerical methods, grid generation, chemically reacting flows, turbulence modeling, inlets, nozzles, and unsteady flows.

Source record↗
NASA NTRS2004

Multi-Fidelity Simulation of a Turbofan Engine With Results Zoomed Into Mini-Maps for a Zero-D Cycle Simulation

A Zero-D cycle simulation of the GE90-94B high bypass turbofan engine has been achieved utilizing mini-maps generated from a high-fidelity simulation. The simulation utilizes the Numerical Propulsion System Simulation (NPSS) thermodynamic cycle modeling system coupled to a high-fidelity full-engine model represented by a set of coupled 3D computational fluid dynamic (CFD) component models. Boundary conditions from the balanced, steady state cycle model are used to define component boundary conditions in the full-engine model. Operating characteristics of the 3D component models are integrated into the cycle model via partial performance maps generated from the CFD flow solutions using one-dimensional mean line turbomachinery programs. This paper highlights the generation of the high-pressure compressor, booster, and fan partial performance maps, as well as turbine maps for the high pressure and low pressure turbine. These are actually "mini-maps" in the sense that they are developed only for a narrow operating range of the component. Results are compared between actual cycle data at a take-off condition and the comparable condition utilizing these mini-maps. The mini-maps are also presented with comparison to actual component data where possible.

Turner, Mark G.↗
NASA NTRS2015

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↗
NASA NTRSDate not supplied

Liquefaction of Cryogenic Fluids for Production and Storage of Commodities on Extra-Terrestrial Surfaces

The Cryogenic Fluid In-situ Liquefaction for Landers (CryoFILL) activity has been investigating concepts for the liquefaction of cryogenic fluids produced on the surface of the Moon and Mars. CryoFILL consists of four activities: liquid nitrogen liquefaction testing that was completed in 2019, liquid oxygen liquefaction testing planned to begin in early 2022, fiber-optic sensor temperature measurement system development, and lightweight vacuum jacketed systems designed for the atmosphere of Mars. Additional activities include the development and validation of numerical modeling systems at three stages: thermodynamic models of the liquefaction process for basic energy/mass estimates, nodal models to predict liquefaction tank system level performance, and computational fluid dynamics to assess fluid phenomena occurring within the tank, specifically condensation and stratification within the ullage. These developments form a combined liquefaction and storage solution to support Lunar and Martian exploration. Key testing results from liquid nitrogen testing as well as plans for liquid oxygen testing will be discussed including key sensitivities from analytical evaluation of completed test results. Progress on the development and improvements on fiber optic sensor testing and developments will also be provided. Finally, work on the progress of the lightweight vacuum jacketed systems being co-developed by NASA and industry will be summarized.

Liquefaction↗
NASA NTRSDate not supplied

Computational Analysis of X-57 Rate Derivatives and Comparison to Wind Tunnel Data

Rate derivates of the wind tunnel model of X-57 were evaluated using hysteresis from computational fluid dynamics simulation and compared with that computed using the wind tunnel experimental data. To accomplish the task, hysteresis of the pitch, roll, and yaw oscillations were produced numerically using two different CFD solvers: STARCCM+ and Launch Ascent and Vehicle Analysis. The rate derivatives were computed using same parameter identification technique on CFD-based and wind tunnel-based hysteresis data. Results show CFD-based pitch rate derivatives have large error, up to 73%, while the lateral and direction derivates compared well with errors within 5% with respect to the that computed using the wind tunnel data. This may be due to lack of wind tunnel wall and sting model in the CFD simulation as well as lack of wall correction in the wind tunnel data which has larger impact on pitch than lateral and directional oscillations. Geometrical features such as the shape, slope of extremities, and enclosed area were compared between CFD-based and wind tunnel-based hysteresis. Results show that area enclosed by the hysteresis had largest influence on the accuracy of the computed rate derivatives given that shape of the hysteresis were roughly equivalent.

Seung Y. Yoo↗
DOE OSTI2025

Multidimensional Modeling of Mixture Formation in a Hydrogen-Fueled Heavy-Duty Optical Engine With Direct Injection

Hydrogen (H 2 ), as a carbon-free fuel, is considered as one of the most promising solutions to reduce the carbon footprint of hard-to-decarbonize energy and transportation sectors. As such, hydrogen-fueled internal combustion engines (H 2 ICEs) have recently been receiving increasing attention, particularly in applications such as on-road/off-road heavy-duty transport and combined heat and power. The direct injection (DI) of gaseous hydrogen into the combustion chamber offers great potential for achieving high power density and high engine efficiency, while mitigating the risk of backfire and reducing pre-ignition. However, the numerical simulation of H 2 DI system remains a formidable challenge associated with the high computational cost of reproducing compressible supersonic flow and shocks in narrow injector passages and in near-nozzle regions. In general, there is a lack of well-established and validated practices for the modeling of high-pressure H 2 DI in large-bore engines. Here, to this end, this study focuses on computational fluid dynamics (CFD) modeling of the mixture formation process in a heavy-duty optical engine employing a medium-pressure H 2 DI system. Both large eddy simulations (LES) and Reynolds Averaged Navier–Stokes (RANS) simulations are performed and evaluated against optical data. Gaseous hydrogen is injected into the combustion chamber via a centrally located outward opening hollow-cone injector at a pressure of 40 bar. Simulations are carried out for two injection timings, namely, −120 and −60 °CA. The numerical predictions for H 2 distribution in different horizontal and vertical planes during the compression stroke are systematically compared against optical data obtained through planar laser-induced fluorescence (PLIF) measurements. Overall, the LES approach using the Dynamic Structure model is found to have good predictive capabilities for the early jet penetration in terms of length and shape, as well as the later H 2 distributions. However, the unsteady RANS approach with the renormalization group $k - ϵ$ model, which is widely used by industry to model heavy-duty ICEs, significantly underpredicts the H 2 mixing, even at similar mesh resolution to that used in LES. These results indicate that there is a need for the improvement of mixing submodels within the RANS approach when applied to H 2 DI simulations.

LES↗
NASA NTRS1993

Numerical study of base pressure characteristic curve for a four-engine clustered nozzle configuration

The objective of this study is to benchmark a four-engine clustered nozzle base flowfield with a computational fluid dynamics (CFD) model. The CFD model is a three-dimensional pressure-based, viscous flow formulation. An adaptive upwind scheme is employed for the spatial discretization. The upwind scheme is based on second and fourth order central differencing with adaptive artificial dissipation. Qualitative base flow features such as the reverse jet, wall jet, recompression shock, and plume-plume impingement have been captured. The computed quantitative flow properties such as the radial base pressure distribution, model centerline Mach number and static pressure variation, and base pressure characteristic curve agreed reasonably well with those of the measurement. Parametric study on the effect of grid resolution, turbulence model, inlet boundary condition and difference scheme on convective terms has been performed. The results showed that grid resolution had a strong influence on the accuracy of the base flowfield prediction.

Wang, Ten-See↗
DOE OSTI2023

Theory and simulations of linear and nonlinear two-dimensional Rayleigh–Taylor dynamics with variable acceleration

Interfacial Rayleigh–Taylor mixing is crucial to describing important natural and engineering processes, such as exploding supernovae, laser micromachining, hot spots in inertial confinement fusion, and optical telecommunications. These require the characterization of the time dependence of the driving acceleration. We compare our theoretical formulation based on group theory foundations with interface-capturing numerical simulations for linear and nonlinear two-dimensional Rayleigh–Taylor instabilities in a finite-sized domain with time-varying acceleration over broad ranges of Atwood numbers and acceleration exponents. Detailed corroboration between theory and simulations is provided for this foundational case. Both demonstrate the strong interfacial nature of Rayleigh–Taylor instabilities, which suggests that practical flow fields can be reconstructed from the derived fluid potential using the proposed theory. A robust agreement is also obtained for the early and late-time evolution of the amplitudes of the bubble and spike, which demonstrate that the Rayleigh–Taylor flow can transition to the mixing regime even for a single-mode initial perturbation. Corroboration with experiments of high energy density plasmas motivated by studies of supernovae is also achieved. In addition, a long-standing puzzle in Rayleigh–Taylor dynamics on the interplay between the acceleration, the shear, and the interface morphology in the theory and simulations is resolved by accounting for finite viscosity of the fluids. The characterization of Rayleigh–Taylor instabilities as a highly interfacial phenomenon provides valuable insight into its multiscale nature, which enhances the design and understanding of numerous processes of practical interest.

42 ENGINEERING↗
NASA NTRS1997

The Implicit and Explicit alpha-mu Schemes

Artificial numerical dissipation is an important issue in large Reynolds number computations. In such computations, the artificial dissipation inherent in traditional numerical schemes can overwhelm the physical dissipation and yield inaccurate results on meshes of practical size. In the present work, the space-time conservation element and solution element method is used to construct new and accurate numerical schemes such that artificial numerical dissipation will not overwhelm physical dissipation. Specifically, these schemes have the property that numerical dissipation vanishes when the physical viscosity goes to zero. These new schemes therefore accurately model the physical dissipation even when it is extremely small. The method of space-time conservation element and solution element, currently under development, is a nontraditional numerical method for solving conservation laws. The method is developed on the basis of local and global flux conservation in a space-time domain, in which space and time are treated in a unified manner. Explicit solvers for model and fluid dynamic conservation laws have previously been investigated. In this paper, we introduce a new concept in the design of implicit schemes, and use it to construct two highly accurate solvers for a convection-diffusion equation. The two schemes become identical in the pure convection case, and in the pure diffusion case. The implicit schemes are applicable over the whole Reynolds number range, from purely diffusive equations to purely inviscid (convective) equations. The stability and consistency of the schemes are analyzed, and some numerical results are presented. It is shown that, in the inviscid case, the new schemes become explicit and their amplification factors are identical to those of the Leapfrog scheme. On the other hand, in the pure diffusion case, their principal amplification factor becomes the amplification factor of the Crank-Nicolson scheme. We also construct an explicit solver with the treatment of diffusion being based on that in the implicit solvers. The explicit solver has only a CFL stability limitation on the Courant number, yet it retains the second-order spatial accuracy of the implicit schemes.

Chang, Sin-Chung↗
NASA NTRS2005

A New Unsteady Model for Dense Cloud Cavitation in Cryogenic Fluids

A new unsteady, cavitation model is presented wherein the phase change process (bubble growth/collapse) is coupled to the acoustic field in a cryogenic fluid. It predicts the number density and radius of bubbles in vapor clouds by tracking both the aggregate surface area and volume fraction of the cloud. Hence, formulations for the dynamics of individual bubbles (e.g. Rayleigh-Plesset equation) may be integrated within the macroscopic context of a dense vapor cloud i.e. a cloud that occupies a significant fraction of available volume and contains numerous bubbles. This formulation has been implemented within the CRUNCH CFD, which has a compressible real fluid formulation, a multi-element, unstructured grid framework, and has been validated extensively for liquid rocket turbopump inducers. Detailed unsteady simulations of a cavitating ogive in liquid nitrogen are presented where time-averaged mean cavity pressure and temperature depressions due to cavitation are compared with experimental data. The model also provides the spatial and temporal history of the bubble size distribution in the vapor clouds that are shed, an important physical parameter that is difficult to measure experimentally and is a significant advancement in the modeling of dense cloud cavitation.

Hosangadi, A.↗
NASA NTRS1988

Computational fluid dynamics drag prediction: Results from the Viscous Transonic Airfoil Workshop

Results from the Viscous Transonic Airfoil Workshop are compared with each other and with experimental data. Test cases used include attached and separated transonic flows for the NACA 0012 airfoil. A total of 23 sets of numerical results from 15 different author groups are included. The numerical method used vary widely and include: 16 Navier-Stokes methods, 2 Euler boundary layer methods, and 5 potential boundary layer methods. The results indicate a high degree of sophistication among the numerical methods with generally good agreement between the various computed and experimental results for attached or moderately separated cases. The agreement for cases with larger separation is only fair and suggests additional work is required in this area.

Holst, Terry L.↗
NASA NTRSDate not supplied

Numerical Modeling of No Vent Filling of a Cryogenic Tank with Thermo-dynamic Vent System Assisted Injector

This paper presents a multi-node finite volume model of No Vent Filling (NVF) of a cryogenic tank with Thermo-dynamic Vent System (TVS) assisted injector, using the Generalized Fluid System Simulation Program (GFSSP), a general purpose flow network code. NVF tests were conducted in a CRYOgenic Orbital Testbed (CRYOTE) tank which was filled by liquid nitrogen using a TVS assisted injector with vent valve closed during the entire filling process. In TVS assisted injector, the liquid flow splits into two streams: one stream is routed to a Joule-Thomson (J-T) orifice where the flow immediately flashes from liquid to vapor or two phase mixture because the downstream to J-T leg is maintained at vacuum level.; the other stream is injected into the tank after being cooled by cold vapor of J-T leg in a heat exchanger. The flow through J-T leg is also used to cool the outer metal matrix of the injector which in turn cools the vapor in the tank ullage. The cooling of vapor in the ullage by cold injector surface reduces ullage pressure that allows liquid to enter and fill the tank. An integrated numerical model of the test set up was developed. The model included the dual lines of TVS assisted injector, tank and tank wall. The tank was discretized into multiple fluid nodes and branches to represent the ullage and liquid nitrogen and multiple solid nodes to represent the tank wall and structure. The heat transfer between solid to fluid was calculated from pool boiling correlations which include film, transition, and nucleate boiling, as well as natural convection during pre- and post-boiling. The model also accounts for the condensation of vapor at ullage-injector interface and when it comes in contact with the liquid spray. The predicted pressure, resident mass, wall temperature in the tank were compared with the test data.

Alok Majumdar↗