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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Conjugate Heat Transfer Analysis in Pressurized Helium Gas Cooling Channels

Los Alamos National Laboratory (LANL) is working with NorthStar Medical Radioisotopes (NMR) on their efforts to produce Mo-99 from the irradiation of Mo-100 targets. The NMR target consists of an Inconel window that allows the electron beam to penetrate and irradiate a stack of Mo-100 disks. The irradiation process generates large amounts of heat and is cooled using pressurized helium gas flowing at 400 g/s. LANL provides both numerical analysis support and experimental support. Computational fluid dynamics (CFD) and conjugate heat transfer (CHT) analysis are performed on the NMR target to analyze the helium flow and volumetric heating through the NMR system. In this iteration of numerical analysis, we update the Mo-100 disk thickness to 0.74 mm and the gap between the Inconel window and the first disk to 0.69 mm. Estimates of the heat flux from the electron beams from Monte Carlo N-Particle code (MCNP) analysis were combined with CFD to perform the CHT analysis. The pressure drop through the targets, force applied on the first disk, beam window and target disk temperatures are some of the quantities estimated using the simulations. The pressure drop through the target was estimated as 0.198 MPa, while the force due to venturi effects across the first disk was estimated to be 100 N. The results show maximum window temperatures over 600 °C and peak disk temperatures over 1600 °C. Past work has shown that beyond 650-700 °C, the Inconel window becomes susceptible to yield and rupture. The increased temperature in these results reduces the margin of error possible on beam spot size. The disk temperatures showed a sinusoidal distribution with the 6 th – 8 th disk from each window exhibiting highest disk temperatures. The peak temperature of the stainless steel housing was close to 316 °C and that of the laminations and spacers was approximately 582 °C.

43 PARTICLE ACCELERATORS↗
NASA NTRS2002

Micro- and Macro-Fluid Dynamics and Acoustics of Resonant Liners

The objectives of this project are to perform direct numerical simulation of the micro-fluid and acoustic fields of a resonant acoustic liner and to investigate the physical processes by which incident sound waves are damped by the acoustic liner. We would like to report that our research work and results have fulfilled both objectives of the grant. The following is a summary of the important accomplishments: (1) Two dimensional direct numerical simulation of the flow and acoustic field around the cavity of resonant liner were successfully carried out; (2) The simulations of (1) were extended to include a laminar grazing flow; (3) The numerical simulations provided strong evidence that there are two principal mechanisms by which a resonant liner damps out an incident acoustic wave; (4) A validation test was performed by comparing the computed dissipation coefficients (not impedance) with impedance tube measurements done at GTRI; and (5) Some resources of this grant were used to support the development of new CAA methods. (Our work on numerical simulation of acoustic liners has benefited by the availability of these improved methods).

Tam, Christopher K. W.↗
NASA NTRS1990

Numerical methods for systems of conservation laws of mixed type using flux splitting

The essentially non-oscillatory (ENO) finite difference scheme is applied to systems of conservation laws of mixed hyperbolic-elliptic type. A flux splitting, with the corresponding Jacobi matrices having real and positive/negative eigenvalues, is used. The hyperbolic ENO operator is applied separately. The scheme is numerically tested on the van der Waals equation in fluid dynamics. Convergence was observed with good resolution to weak solutions for various Riemann problems, which are then numerically checked to be admissible as the viscosity-capillarity limits. The interesting phenomena of the shrinking of elliptic regions if they are present in the initial conditions were also observed.

Shu, Chi-Wang↗
DOE OSTI2022

Fluid Dynamics of Axial Turbomachinery: Blade- and Stage-Level Simulations and Models

The current generation of axial turbomachines is the culmination of decades of experience, and detailed understanding of the underlying flow physics has been a key factor for achieving high efficiency and reliability. Driven by advances in numerical methods and relentless growth in computing power, computational fluid dynamics has increasingly provided insights into the rich fluid dynamics involved and how it relates to loss generation. This article presents some of the complex flow phenomena occurring in bladed components of gas turbines and illustrates how simulations have contributed to their understanding and the challenges they pose for modeling. The interaction of key aerodynamic features with deterministic unsteadiness, caused by multiple blade rows, and stochastic unsteadiness, i.e., turbulence, is discussed. High-fidelity simulations of increasingly realistic configurations and models improved with help of machine learning promise to further grow turbomachinery performance and reliability and, thus, help fluid mechanics research have a greater industrial impact.

Mechanics↗
DOE OSTI2020

Free-Surface Liquid Lithium Flow Modeling and Stability Analysis for Fusion Applications

Liquid metal plasma facing components are considered an attractive design choice for fusion devices including pilot plants. Virtual prototyping of such devices includes modeling of free-surface flow of the electrically conductive liquid, which requires computational fluid dynamics (CFD) and magnetohydrodynamics (MHD) simulations. Numerical tools capable of simulating flows and heat transfer in the free-surface MHD flow were developed at PPPL based on the customized ANSYS CFX. Here, MHD is introduced using a magnetic vector potential approach. Free-surface flow capabilities are available in the code and were tested. Special stabilization procedures were derived and applied to improve convergence of the momentum equations with the source terms due to the Lorentz force and surface tension. Important characteristics of the fusion-relevant liquid metal flow is free surface smoothness and stability. Heat flux from the plasma impacts the liquid surface at a very acute angle, so any change of the free surface from axisymmetry can dramatically increase the local heat flux density and thus create excessive evaporation of liquid lithium into the plasma, which is detrimental to operations. Stability analysis of the liquid metal film flow was performed to determine applicable flow regimes. Thin film flow along horizontal wall is considered. Effects of gravity, magnetic field, and surface tension are included in the analysis.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

A numerical solver for investigating the space charge effect on the electric field in liquid argon time projection chambers

This paper reports the development of a numerical solver aimed to simulate the interaction between the space charge (i.e. ions) distribution and the electric field in liquid argon time projection chamber (LArTPC) detectors. The ion transport equation is solved by a time-accurate, cell-centered finite volume method and the electric potential equation by a continuous finite element method. The electric potential equation updates the electric field which provides the drift velocity to the ion transport equation. The ion transport equation updates the space charge density distribution which appears as the source term in the electric potential equation. The interaction between the space charge distribution and the electric field is numerically simulated within each physical time step. The convective velocity in the ion transport equation can include the background flow velocity in addition to the electric drift velocity. The numerical solver has been parallelized using the Message Passing Interface (MPI) library. Numerical tests show and verify the capability and accuracy of the current numerical solver. It is planned that the developed numerical solver, together with a Computational Fluid Dynamics (CFD) package which provides the flow velocity field, can be used to investigate the space charge effect on the electric field in large-scale particle detectors.

42 ENGINEERING↗
DOE OSTI2022

Modeling the small-scale deposition of snow onto structured Arctic sea ice during a MOSAiC storm using snowBedFoam 1.0.

The remoteness and extreme conditions of the Arctic make it a very difficult environment to investigate. In these polar regions covered by sea ice, the wind is relatively strong due to the absence of obstructions and redistributes a large part of the deposited snow mass, which complicates estimates for precipitation hardly distinguishable from blowing or drifting snow. Moreover, the snow mass balance in the sea ice system is still poorly understood, notably due to the complex structure of its surface. Quantitatively assessing the snow distribution on sea ice and its connection to the sea ice surface features is an important step to remove the snow mass balance uncertainties (i.e., snow transport contribution) in the Arctic environment. In this work we introduce snowBedFoam 1.0., a physics-based snow transport model implemented in the open-source fluid dynamics software OpenFOAM. We combine the numerical simulations with terrestrial laser scan observations of surface dynamics to simulate snow deposition in a MOSAiC (Multidisciplinary Drifting Observatory for the Study of Arctic Climate) sea ice domain with a complicated structure typical for pressure ridges. The results demonstrate that a large fraction of snow accumulates in their vicinity, which compares favorably against scanner measurements. However, the approximations imposed by the numerical framework, together with potential measurement errors (precipitation), give rise to quantitative inaccuracies, which should be addressed in future work. The modeling of snow distribution on sea ice should help to better constrain precipitation estimates and more generally assess and predict snow and ice dynamics in the Arctic.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1985

Numerical solution of strong radiation gasdynamic interactions in a hydrogen-seedant mixture

Theoretical and experimental results are provided for the absorption of laser energy by a gas flow. A two-dimensional numerical model described the temperature field, fluid dynamics, variable absorptivity and gas dynamics. The model was tested using a hydrogen-seedant mixture injected with 10.6 micron laser radiation. The laser mix was selected to allow reduced temperature gradients and energy source terms. It was found possible, by varying parameters such as the laser power and flash pressure, to precisely position the peak absorptance region in higher speed flows. More refined models are being tested to achieve accurate absorptance predictions in flows below Mach 0.0001.

Merkle, C. L.↗
NASA NTRS1987

Parallel computations and their impact on mechanics; Proceedings of the Symposium, ASME Winter Annual Meeting, Boston, MA, Dec. 13-18, 1987

The conference presents papers on parallel architectures and the programming environment, parallel numerical algorithms, structural mechanics applications, and fluid dynamics applications. Topics include concurrent computer architecture, the state-of-the-art in highly parallel computer systems, supercomputer programming environments, uniquely parallel algorithms, the parallel solution of nonlinear elliptic equations, highly parallel banded systems solvers, and optimal mapping of irregular finite element domains to parallel processors. Consideration is also given to parallel processing in finite element structural analysis, heirarchical parallelism in a finite element CFD algorithm, and a parallelized elliptic solver for reacting flows.

Noor, Ahmed K.↗
NASA NTRS1990

Solution of the Navier-Stokes equations in three dimensions

The impact of computational fluid dynamics on the design of fluid dynamics devices is discussed. An efficient numerical procedure for solving the Navier-Stokes equations in three dimensions is discussed and flow simulations made on an inexpensive desk top work station are presented.

Maccormack, Robert W.↗
NASA NTRS1992

Computational fluid dynamics combustion analysis evaluation

This study involves the development of numerical modelling in spray combustion. These modelling efforts are mainly motivated to improve the computational efficiency in the stochastic particle tracking method as well as to incorporate the physical submodels of turbulence, combustion, vaporization, and dense spray effects. The present mathematical formulation and numerical methodologies can be casted in any time-marching pressure correction methodologies (PCM) such as FDNS code and MAST code. A sequence of validation cases involving steady burning sprays and transient evaporating sprays will be included.

Kim, Y. M.↗
NASA NTRS1993

Tensoral: A system for post-processing turbulence simulation data

Many computer simulations in engineering and science -- and especially in computational fluid dynamics (CFD) -- produce huge quantities of numerical data. These data are often so large as to make even relatively simple post-processing of this data unwieldy. The data, once computed and quality-assured, is most likely analyzed by only a few people. As a result, much useful numerical data is under-utilized. Since future state-of-the-art simulations will produce even larger datasets, will use more complex flow geometries, and will be performed on more complex supercomputers, data management issues will become increasingly cumbersome. My goal is to provide software which will automate the present and future task of managing and post-processing large turbulence datasets. My research has focused on the development of these software tools -- specifically, through the development of a very high-level language called 'Tensoral'. The ultimate goal of Tensoral is to convert high-level mathematical expressions (tensor algebra, calculus, and statistics) into efficient low-level programs which numerically calculate these expressions given simulation datasets. This approach to the database and post-processing problem has several advantages. Using Tensoral the numerical and data management details of a simulation are shielded from the concerns of the end user. This shielding is carried out without sacrificing post-processor efficiency and robustness. Another advantage of Tensoral is that its very high-level nature lends itself to portability across a wide variety of computing (and supercomputing) platforms. This is especially important considering the rapidity of changes in supercomputing hardware.

Dresselhaus, Eliot↗
NASA NTRS1996

Sloshing of Cryogenic Helium Driven by Lateral Impulse/Gravity Gradient-Dominated/or g-Jitter-Dominated Accelerations and Orbital Dynamics

The coupling of slosh dynamics within a partially filled rotating dewar of superfluid helium 11 with spacecraft orbital dynamics is investigated in response to the environmental disturbances of (a) lateral impulses, (b) gravity gradients and (c) g-jitter forces. The purpose of this study is to investigate how the coupling of helium 11 fluid slosh dynamics driven by three cases of environmental force with spacecraft dynamics can affect the bubble deformations and their associated fluid and spacecraft mass centre fluctuations. The numerical computation of slosh dynamics is based on a rotational frame, while the spacecraft dynamics is associated with a non-rotational frame. Results show that the major contribution of orbital dynamics is driven by coupling with slosh dynamics. Neglecting the effect of slosh dynamics acting on the spacecraft may lead to the wrong results for the development of orbital and attitude control techniques.

Hung, R. J.↗
NASA NTRS2001

Velocity Measurement in a Dual-Mode Supersonic Combustor using Particle Image Velocimetry

Temporally and spatially-resolved, two-component measurements of velocity in a supersonic hydrogen-air combustor are reported. The combustor had a single unswept ramp fuel injector and operated with an inlet Mach number of 2 and a flow total temperature approaching 1200 K. The experiment simulated the mixing and combustion processes of a dual-mode scramjet operating at a flight Mach number near 5. The velocity measurements were obtained by seeding the fuel with alumina particles and performing Particle Image Velocimetry on the mixing and combustion wake of the ramp injector. To assess the effects of combustion on the fuel air-mixing process, the distribution of time-averaged velocity and relative turbulence intensity was determined for the cases of fuel-air mixing and fuel-air reacting. Relative to the mixing case, the near field core velocity of the reacting fuel jet had a slower streamwise decay. In the far field, downstream of 4 to 6 ramp heights from the ramp base, the heat release of combustion resulted in decreased flow velocity and increased turbulence levels. The reacting measurements were also compared with a computational fluid dynamics solution of the flow field. Numerically predicted velocity magnitudes were higher than that measured and the jet penetration was lower.

Goyne, C. P.↗
NASA NTRS2000

The fluid dynamic and shear environment in the NASA/JSC rotating-wall perfused-vessel bioreactor

The rotating-wall perfused-vessel (RWPV) bioreactor, used for both microgravity and Earth-based cell science experiments, is characterized in terms of the fluid dynamic and fluid shear stress environment. A numerical model of the flow field is developed and verified with laser Doppler velocimeter measurements. The effects of changes in operating conditions, including rotation rates and fluid perfusion rates, are investigated with the numerical model. The operating conditions typically used for ground-based experiments (equal rotation of the inner and outer cylinders) leads to flow patterns with relatively poor mass distribution characteristics. Approximately 50% of the inlet-perfused fluid bypasses the bulk of the fluid volume and flows to the perfusion exit. For operating conditions typical in microgravity, small differential rotation rates between the inner and outer cylinders lead to greatly improved flow distribution patterns and very low fluid shear stress levels over a large percentage of the fluid volume. Differences in flow patterns for the different operating conditions are explored. Large differences in the hydrodynamic environments for operating conditions typical of true microgravity and ground-based "microgravity simulations" are demonstrated.

Bioreactors↗
NASA NTRS2005

CFD - Mature Technology?

Over the past 30 years, numerical methods and simulation tools for fluid dynamic problems have advanced as a new discipline, namely, computational fluid dynamics (CFD). Although a wide spectrum of flow regimes are encountered in many areas of science and engineering, simulation of compressible flow has been the major driver for developing computational algorithms and tools. This is probably due to a large demand for predicting the aerodynamic performance characteristics of flight vehicles, such as commercial, military, and space vehicles. As flow analysis is required to be more accurate and computationally efficient for both commercial and mission-oriented applications (such as those encountered in meteorology, aerospace vehicle development, general fluid engineering and biofluid analysis) CFD tools for engineering become increasingly important for predicting safety, performance and cost. This paper presents the author's perspective on the maturity of CFD, especially from an aerospace engineering point of view.

Kwak, Dochan↗
NASA NTRS2015

2nd-Order CESE Results For C1.4: Vortex Transport by Uniform Flow

The Conservation Element and Solution Element (CESE) method was used as implemented in the NASA research code ez4d. The CESE method is a time accurate formulation with flux-conservation in both space and time. The method treats the discretized derivatives of space and time identically and while the 2nd-order accurate version was used, high-order versions exist, the 2nd-order accurate version was used. In regards to the ez4d code, it is an unstructured Navier-Stokes solver coded in C++ with serial and parallel versions available. As part of its architecture, ez4d has the capability to utilize multi-thread and Messaging Passage Interface (MPI) for parallel runs.

Space-time CE/SE method↗