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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
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↗
NASA NTRS2015

Structured Overlapping Grid Simulations of Contra-rotating Open Rotor Noise

Computational simulations using structured overlapping grids with the Launch Ascent and Vehicle Aerodynamics (LAVA) solver framework are presented for predicting tonal noise generated by a contra-rotating open rotor (CROR) propulsion system. A coupled Computational Fluid Dynamics (CFD) and Computational AeroAcoustics (CAA) numerical approach is applied. Three-dimensional time-accurate hybrid Reynolds Averaged Navier-Stokes/Large Eddy Simulation (RANS/LES) CFD simulations are performed in the inertial frame, including dynamic moving grids, using a higher-order accurate finite difference discretization on structured overlapping grids. A higher-order accurate free-stream preserving metric discretization with discrete enforcement of the Geometric Conservation Law (GCL) on moving curvilinear grids is used to create an accurate, efficient, and stable numerical scheme. The aeroacoustic analysis is based on a permeable surface Ffowcs Williams-Hawkings (FW-H) approach, evaluated in the frequency domain. A time-step sensitivity study was performed using only the forward row of blades to determine an adequate time-step. The numerical approach is validated against existing wind tunnel measurements.

Grid↗
DOE OSTI2024

Fluid-structure interaction solver for transient dynamics of fracturing media

Computer implemented methods include providing a numerical simulation solid mesh comprising a plurality of solid nodes and a numerical simulation fluid mesh comprising a plurality of fluid nodes, and performing a numerical simulation with the solid mesh and fluid mesh including: determining a fluid-solid interaction between the fluid mesh and the solid mesh at a time step by (i) immersing boundaries of the solid and fluid meshes, (ii) computing current corrected fluid and solid velocities using an interaction parameter applied to an intermediate fluid-solid relative velocity, and (iii) computing a fluid-solid interaction force based on the current corrected fluid and solid velocities.

Rougier, Esteban↗
NASA NTRS1993

Experience in using SIMD and MIMD parallelism for computational fluid dynamics

One of the key objectives of the Applied Research Branch in the Numerical Aerodynamic Simulation (NAS) Systems Division at NASA Ames Research Center is the accelerated introduction of highly parallel machines into a fully operational environment. In this report we summarize some of the experiences with the parallel testbed machines at the NAS Applied Research Branch. We discuss the performance results obtained from the implementation of two computational fluid dynamics (CFD) applications, an unstructured grid solver and a particle simulation, on the Connection Machine CM-2 and the Intel iPSC/860.

Simon, Horst D.↗
NASA NTRS1993

Sloshing dynamics on rotating helium dewar tank

The generalized mathematical formulation of sloshing dynamics for partially filled liquid of cryogenic superfluid helium II in dewar containers driven by both the gravity gradient and jitter accelerations applicable to scientific spacecraft which is eligible to carry out spinning motion and/or slew motion for the purpose to perform scientific observation during the normal spacecraft operation are investigated. An example is given with Gravity Probe-B (GP-B) spacecraft which is responsible for the sloshing dynamics. The jitter accelerations include slew motion, spinning motion, atmospheric drag on the spacecraft, spacecraft attitude motions arising from machinery vibrations, thruster firing, pointing control of spacecraft, crew motion, etc. Explicit mathematical expressions to cover these forces acting on the spacecraft fluid systems are derived. The numerical computation of sloshing dynamics were based on the non-inertia frame spacecraft bound coordinate, and solve time dependent, three-dimensional formulations of partial differential equations subject to initial and boundary conditions. The explicit mathematical expressions of boundary conditions to cover capillary force effect on the liquid vapor interface in microgravity environments are also derived. The formulations of fluid moment and angular moment fluctuations in fluid profiles induced by the sloshing dynamics, together with fluid stress and moment fluctuations exerted on the spacecraft dewar containers were derived. Results were widely published in the open journals.

Hung, R. J.↗
NASA NTRS1975

Higher order accurate partial implicitization: An unconditionally stable fourth-order-accurate explicit numerical technique

The previously obtained second-order-accurate partial implicitization numerical technique used in the solution of fluid dynamic problems was modified with little complication to achieve fourth-order accuracy. The Von Neumann stability analysis demonstrated the unconditional linear stability of the technique. The order of the truncation error was deduced from the Taylor series expansions of the linearized difference equations and was verified by numerical solutions to Burger's equation. For comparison, results were also obtained for Burger's equation using a second-order-accurate partial-implicitization scheme, as well as the fourth-order scheme of Kreiss.

Graves, R. A., Jr.↗
NASA NTRS1994

Superfluid helium sloshing dynamics induced oscillations and fluctuations of angular momentum, force and moment actuated on spacecraft driven by gravity gradient or jitter acceleration associated with slew motion

The generalized mathematical formulation of sloshing dynamics for partially filled liquid of cryogenic superfluid helium II in dewar containers driven by the gravity gradient and jitter accelerations associated with slew motion for the purpose to perform scientific observation during the normal spacecraft operation are investigated. An example is given with the Advanced X-Ray Astrophysics Facility-Spectroscopy (AXAF-S) for slew motion which is responsible for the sloshing dynamics. The jitter accelerations include slew motion, spinning motion, atmospheric drag on the spacecraft, spacecraft attitude motions arising from machinery vibrations, thruster firing, pointing control of spacecraft, crew motion, etc. Explicit mathematical expressions to cover these forces acting on the spacecraft fluid systems are derived. The numerical computation of sloshing dynamics is based on the non-inertia frame spacecraft bound coordinate, and solve time-dependent, three-dimensional formulations of partial differential equations subject to initial and boundary conditions. The explicit mathematical expressions of boundary conditions to cover capillary force effect on the liquid-vapor interface in microgravity environments are also derived. The formulations of fluid moment and angular moment fluctuations in fluid profiles induced by the sloshing dynamics, together with fluid stress and moment fluctuations exerted on the spacecraft dewar containers have also been derived. Examples are also given for cases applicable to the AXAF-S spacecraft sloshing dynamics associated with slew motion.

Hung, R. J.↗
NASA NTRS2004

Dynamical Modeling of Mars' Paleoclimate

This report summarizes work undertaken under a one-year grant from the NASA Mars Fundamental Research Program. The goal of the project was to initiate studies of the response of the Martian climate to changes in planetary obliquity and orbital elements. This work was undertaken with a three-dimensional numerical climate model based on the Geophysical Fluid Dynamics Laboratory (GFDL) Skyhi General Circulation Model (GCM). The Mars GCM code was adapted to simulate various obliquity and orbital parameter states. Using a version of the model with a basic water cycle (ice caps, vapor, and clouds), we examined changes in atmospheric water abundances and in the distribution of water ice sheets on the surface. This work resulted in a paper published in the Journal of Geophysical Research - Planets. In addition, the project saw the initial incorporation of a regolith water transport and storage scheme into the model. This scheme allows for interaction between water in the pores of the near subsurface (<3m) and the atmosphere. This work was not complete by the end of the one-year grant, but is now continuing within the auspices of a three-year grant of the same title awarded by the Mars Fundamental Research Program in late 2003.

Richardson, Mark I.↗
NASA NTRS1993

Simulation of sloshing dynamics induced forces and torques actuated on dewar container driven by gravity gradient and jitter accelerations in microgravity

Some experimental spacecraft use superconducting sensors for gyro read-out and so must be maintained at a very low temperature. The boil-off from the cryogenic liquid used to cool the sensors can also be used, as the Gravity Probe B (GP-B) spacecraft does, as propellant to maintain attitude control and drag-free operation of the spacecraft. The cryogenic liquid for such spacecraft is, however, susceptible to both slosh-like motion and non-axisymmetric configurations under the influence of various kinds of gravity jitter and gravity gradient accelerations. Hence, it is important to quantify the magnitude of the liquid-induced perturbations on the spacecraft. We use the example of the GP-B to investigate such perturbations by numerical simulations. For this spacecraft disturbances can be imposed on the liquid by atmospheric drag, spacecraft attitude control maneuvers, and the earth's gravity gradient. More generally, onboard machinery vibrations and crew motion can also create disturbances. Recent studies suggest that high frequency disturbances are relatively unimportant in causing liquid motions in comparison to low frequency ones. The results presented here confirm this conclusion. After an initial calibration period, the GP-B spacecraft rotates in orbit at 0.1 rpm about the tank symmetry axis. For this rotation rate, the equilibrium liquid free surface shape is a 'doughnut' configuration for all residual gravity levels of 10(exp -6) g(sub 0) or less, as shown by experiments and by numerical simulations; furthermore, the superfluid behavior of the 1.8 K liquid helium used in GP-B eliminates temperature gradients and therefore such effects as Marangoni convection do not have to be considered. Classical fluid dynamics theory is used as the basis of the numerical simulations here, since Mason's experiments show that the theory is applicable for cryogenic liquid helium in large containers. To study liquid responses to various disturbances, we investigate and simulate three levels of gravity jitter (10(exp -6), 10(exp -7), and 10(exp -8) g(sub 0)) each at three predominant frequencies (0.1, 1.0, and 10 Hz), combined with a gravity gradient appropriate for the GP-B orbit. Dynamical evolution of sloshing dynamics excited fluid forces and torque fluctuations exerted on the dewar container driven by the combined gravity gradient and jitter accelerations are also investigated and simulated.

Hung, R. J.↗
NASA NTRS1994

Sloshing dynamics modulated angular momentum fluctuations actuated by gravity gradient acceleration associated with spacecraft slew motion

The mathematical formulation of orbital spacecraft sloshing dynamics for a partially filled cryogenic superfluid liquid helium II in a dewar container actuated by the gravity gradient acceleration associated with slew motion is studied. The Advanced X-Ray Astrophysics Facility-Spectroscopy (AXAF-S) spacecraft is chosen as a practical example in this study. Explicit mathematical expressions which manage orbital galaxy gradient acceleration with slew motion which is acting on the spacecraft fluid systems are derived. The numerical computation of sloshing dynamics is based on the non-inertia frame spacecraft bound coordinates and the solution of time-dependent, three-dimensional formulations of partial differential equations subject to initial and boundary conditions. This study discloses the capillary effect of sloshing dynamics governed liquid-vapor interface fluctuations, angular momentum and moment fluctuations of fluid system, and also buble mass center fluctuations driven by the gravity gradient acceleration associated with slew motion which affects the stability of the orbital spacecraft fluid system in a microgravity environment.

Hung, R. J.↗
NASA NTRS1988

Numerical modeling of multidimensional flow in seals and bearings used in rotating machinery

The rotordynamic behavior of turbomachinery is critically dependent on fluid dynamic rotor forces developed by various types of seals and bearings. The occurrence of self-excited vibrations often depends on the rotor speed and load. Misalignment and rotor wobbling motion associated with differential clearance were often attributed to stability problems. In general, the rotative character of the flowfield is a complex three dimensional system with secondary flow patterns that significantly alter the average fluid circumferential velocity. A multidimensional, nonorthogonal, body-fitted-grid fluid flow model is presented that describes the fluid dynamic forces and the secondary flow pattern development in seals and bearings. Several numerical experiments were carried out to demonstrate the characteristics of this complex flowfield. Analyses were performed by solving a conservation form of the three dimensional Navier-Stokes equations transformed to those for a rotating observer and using the general-purpose computer code PHOENICS with the assumptions that the rotor orbit is circular and that static eccentricity is zero. These assumptions have enabled a precise steady-state analysis to be used. Fluid injection from ports near the seal or bearing center increased fluid-film direct dynamic stiffness and, in some cases, significantly increased quadrature dynamic stiffness. Injection angle and velocity could be used for active rotordynamic control; for example, injection, when compared with no injection, increased direct dynamic stiffness, which is an important factor for hydrostatic bearings.

Hendricks, R. C.↗
DOE OSTI2024

Numerical modeling of Joule heated ceramic melter

Computational fluid dynamics and heat transfer models of waste glass melters are being developed to support the vitrification campaign at the Waste Treatment and Immobilization Plant (WTP) as part of the United States nuclear waste stewardship and environmental protection programs. The WTP mission aims to achieve the stable and durable processing of vitrified waste in the safest, most efficient, and economical manner possible. Hanford’s tank waste is far more complex and varied than the waste treated by other vitrification facilities and both startup and long-term operational challenges are likely to be even greater than prior experience has shown. There are thousands of different waste compositions, which makes processing the waste a complex and challenging task. Models of test melters, pilot-scale, and full-scale melters have been developed to provide understanding and insight into the various physical and chemical phenomena occurring during vitrification. So the results of these modeling activities will be used to support the WTP operations to process legacy nuclear waste into a stable form for disposal.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2021

Numerical Analysis of Fuel Effects on Advanced Compression Ignition Using a Cooperative Fuel Research Engine Computational Fluid Dynamics Model

Growing environmental concerns and demand for a better fuel economy are driving forces that motivate the research for more advanced engines. Multi-mode combustion strategies have gained attention for their potential to provide high thermal efficiency and low emissions for light-duty applications. These strategies target optimizing the engine performance by correlating different combustion modes to load operating conditions. The extension from boosted spark ignition (SI) mode at high loads to advanced compression ignition (ACI) mode at low loads can be achieved by increasing the compression ratio and utilizing intake air heating. Further, in order to enable an accurate control of intake charge condition for ACI mode and rapid mode-switches, it is essential to gain fundamental insights into the autoignition process. Within the scope of ACI, homogeneous charge compression ignition (HCCI) mode is of significant interest. It is known for its potential benefits, operation at low fuel consumption, low NOx, and particulate matter (PM) emissions. In the present work, a virtual Cooperative Fuel Research (CFR) engine model is used to analyze fuel effects on ACI combustion. In particular, the effect of fuel octane sensitivity (S) (at constant Research Octane Number (RON)) on autoignition propensity is assessed under beyond-RON (BRON) and beyond-MON (BMON) ACI conditions. The three-dimensional CFR engine computational fluid dynamics (CFD) model employs a finite-rate chemistry approach with a multi-zone binning strategy to capture autoignition. Two binary blends with Research Octane Number (RON) of 90 are chosen for this study: primary reference fuel (PRF) with S = 0 and toluene–heptane (TH) blend with S = 10.8, representing paraffinic and aromatic gasoline surrogates. Reduced mechanisms for these blends are generated from a detailed gasoline surrogate kinetic mechanism. Simulation results with the reduced mechanisms are validated against experimental data from an in-house CFR engine, with respect to in-cylinder pressure, heat release rate, and combustion phasing. Thereafter, the sensitivity of combustion behavior to ACI operating condition (BRON versus BMON), air-fuel ratio (λ = 2 and 3), and engine speed (600 and 900 rpm) is analyzed for both fuels. It is shown that the sensitivity of a fuel’s autoignition characteristics to λ and engine speed significantly differs at BRON and BMON conditions. Moreover, this sensitivity is found to vary among fuels, despite the same RON. It is also observed that the presence of low-temperature heat release (LTHR) under BRON condition leads to more sequential autoignition and longer combustion duration than BMON condition. Finally, the study indicates that the octane index (OI) fails to capture the trend in the variation of autoignition propensity with S under the BMON condition.

33 ADVANCED PROPULSION SYSTEMS↗
NASA NTRS2004

Magnetohydrodynamic Convection in the Outer Core and its Geodynamic Consequences

The Earth's fluid outer core is in vigorous convection through much of the Earth's history. In addition to generating and maintaining Earth s time-varying magnetic field (geodynamo), the core convection also generates mass redistribution in the core and a dynamical pressure field on the core-mantle boundary (CMB). All these shall result in various core-mantle interactions, and contribute to surface geodynamic observables. For example, electromagnetic core-mantle coupling arises from finite electrically conducting lower mantle; gravitational interaction occurs between the cores and the heterogeneous mantle; mechanical coupling may also occur when the CMB topography is aspherical. Besides changing the mantle rotation via the coupling torques, the mass-redistribution in the core shall produce a spatial-temporal gravity anomaly. Numerical modeling of the core dynamical processes contributes in several geophysical disciplines. It helps explain the physical causes of surface geodynamic observables via space geodetic techniques and other means, e.g. Earth's rotation variation on decadal time scales, and secular time-variable gravity. Conversely, identification of the sources of the observables can provide additional insights on the dynamics of the fluid core, leading to better constraints on the physics in the numerical modeling. In the past few years, our core dynamics modeling efforts, with respect to our MoSST model, have made significant progress in understanding individual geophysical consequences. However, integrated studies are desirable, not only because of more mature numerical core dynamics models, but also because of inter-correlation among the geophysical phenomena, e.g. mass redistribution in the outer core produces not only time-variable gravity, but also gravitational core-mantle coupling and thus the Earth's rotation variation. They are expected to further facilitate multidisciplinary studies of core dynamics and interactions of the core with other components of the Earth.

Kuang, Weijia↗
NASA NTRS1989

Assessment of numerical techniques for unsteady flow calculations

The characteristics of unsteady flow motions have long been a serious concern in the study of various fluid dynamic and combustion problems. With the advancement of computer resources, numerical approaches to these problems appear to be feasible. The objective of this paper is to assess the accuracy of several numerical schemes for unsteady flow calculations. In the present study, Fourier error analysis is performed for various numerical schemes based on a two-dimensional wave equation. Four methods sieved from the error analysis are then adopted for further assessment. Model problems include unsteady quasi-one-dimensional inviscid flows, two-dimensional wave propagations, and unsteady two-dimensional inviscid flows. According to the comparison between numerical and exact solutions, although second-order upwind scheme captures the unsteady flow and wave motions quite well, it is relatively more dissipative than sixth-order central difference scheme. Among various numerical approaches tested in this paper, the best performed one is Runge-Kutta method for time integration and six-order central difference for spatial discretization.

Hsieh, Kwang-Chung↗
NASA NTRS2016

Experimental, Numerical, and Analytical Slosh Dynamics of Water and Liquid Nitrogen in a Spherical Tank

Understanding, predicting, and controlling fluid slosh dynamics is critical to safety and improving performance of space missions when a significant percentage of the spacecraft's mass is a liquid. Computational fluid dynamics simulations can be used to predict the dynamics of slosh, but these programs require extensive validation. Many experimental and numerical studies of water slosh have been conducted. However, slosh data for cryogenic liquids is lacking. Water and cryogenic liquid nitrogen are used in various ground-based tests with a spherical tank to characterize damping, slosh mode frequencies, and slosh forces. A single ring baffle is installed in the tank for some of the tests. Analytical models for slosh modes, slosh forces, and baffle damping are constructed based on prior work. Select experiments are simulated using a commercial CFD software, and the numerical results are compared to the analytical and experimental results for the purposes of validation and methodology-improvement.

slosh↗