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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 775 records · Page 43

All sources NASA NTRS DOE OSTI

NASA NTRS1987

Cellular automatons applied to gas dynamic problems

This paper compares the results of a relatively new computational fluid dynamics method, cellular automatons, with experimental data and analytical results. This technique has been shown to qualitatively predict fluidlike behavior; however, there have been few published comparisons with experiment or other theories. Comparisons are made for a one-dimensional supersonic piston problem, Stokes first problem, and the flow past a normal flat plate. These comparisons are used to assess the ability of the method to accurately model fluid dynamic behavior and to point out its limitations. Reasonable results were obtained for all three test cases, but the fundamental limitations of cellular automatons are numerous. It may be misleading, at this time, to say that cellular automatons are a computationally efficient technique. Other methods, based on continuum or kinetic theory, would also be very efficient if as little of the physics were included.

Long, Lyle N.↗

NASA NTRS1988

An investigation of turbulence models

The accuracy to which a turbulent boundary layer or wake can be predicted numerically depends on the validity of the turbulence closure model used. The modeling of turbulence physics is one of the most difficult problems in computational fluid dynamics (CFD). In fact, it is one of the pacing factors in the development of CFD. In general, there are three main approaches to the description of trubulence physics. First is turbulence modeling in which the Reynolds averaged Navier-Stokes equations are used and some closure approximation is made for the Reynolds stresses. A second approach to turbulence is large eddy simulation (LES) in which the computational mesh is taken to be fine enough that the large scale structure of the turbulence can be calculated directly. An empirical assumption must be made for the small scale sub-grid turbulence. The third approach is direct simulation. In this technique the Navier-Strokes equations are solved directly on a mesh which if fine enough to resolve the smallest length scale of the turbulence. The Reynolds averaged equations are not used and no closure assumption is required. These last two approaches require extensive computer resources and as such are not engineering tools. The purpose of the work was to investigate the various engineering turbulence models for accuracy and ease of programming. This involved comparison of the models with each other and with experimental data.

Doria, Michael L.↗

NASA NTRS1988

Numerical simulation methods of incompressible flows and an application to the Space Shuttle main engine

This paper discusses incompressible Navier-Stokes solution methods with an emphasis on the pseudocompressibility method. A steady-state flow solver based on the pseudocompressibility approach is then described. This flow solver code has been used to analyze the internal flow in the Space Shuttle main engine hot-gas manifold. Salient features associated with this three-dimensional realistic flow simulation are discussed. Numerical solutions relevant to the current engine analysis and the redesign effort are discussed along with experimental results. This example demonstrates the potential of computational fluid dynamics as a design tool for aerospace applications.

Chang, J. L. C.↗

NASA NTRS1989

Spectroscopic study of combustion diagnostics on hydroxyl radicals

Experimental observations of propane/air flames have been carried out. Measurements of hydroxyl (OH) radical concentration have been made using resonance line absorption techniques. A microwave-pumped low pressure discharge in argon and water vapor is employed to produce strong OH radical band radiation in the 308 nm region. This radiation is transmitted through the plume and absorption data is taken at various radial positions using an optical multichannel analyzer. This absorption data is used to compute OH number density using a model for the absorption band characteristics as a function of temperature based on an atlas of line strengths. A numerical computation of flow fields, temperature profile, OH number density and water vapor mass fraction profiles was carried out by using a technique of computational fluid dynamics (CFD). The results of CFD computation are compared with experimental observation with a good agreement.

Eskridge, Richard↗

NASA NTRS1988

Numerical simulation methods of incompressible flows and an application to the space shuttle main engine

Incompressible Navier-Stokes solution methods are discussed with an emphasis on the pseudocompressibility method. A steady-state flow solver based on the pseudocompressibility approach is then described. This flow-solver code was used to analyze the internal flow in the Space Shuttle main engine hot-gas manifold. Salient features associated with this three-dimensional realistic flow simulation are discussed. Numerical solutions relevant to the current engine analysis and the redesign effort are discussed along with experimental results. This example demonstrates the potential of computational fluid dynamics as a design tool for aerospace applications.

Chang, J. L. C.↗

NASA NTRS1991

Computational aerodynamics - The next generation

Continued advances in the various elements that comprise the field of computational fluid dynamics (CFD) are promoting a radically different approach to the aerodynamic design and analysis of aerospace vehicles and systems. The elements of CFD generally include numerical algorithm development, transition and turbulence modeling, surface modeling, and grid generation, scientific visualization and validation methodologies. This paper discusses the research progress and prospects for the future in each of these elements within NASA's CFD and Experimental Validation Program. The applicability of computational methods for the purposes of understanding complex flow phenomena, exploring aerodynamic concepts, and providing vehicle-design input is also addressed.

Hessenius, Kristin A.↗

NASA NTRS2004

Shock Waves Mitigation at Blunt Bodies Using Needles and Shells Against a Supersonic Flow

The paper contains some experimental and numerical simulation test results on cylindrical blunt body drag reduction using thin spikes or shell mounted in front of a body against a supersonic flow. Experimental tests were conducted using the Aeromechanics and Gas Dynamics Laboratory facilities at the Institute of Mechanics of Moscow State University (IMMSU). Numerical simulations utilizing NASA and IM/MSU codes were conducted at the Hampton University Fluid Mechanics and Acoustics Laboratory. The main purpose of this research is to examine the efficiency of application of multiple spikes for drag reduction and flow stability at the front of a blunt body in different flight conditions, i.e. Mach number, angle of attack, etc. The principal conclusions of these test results are: multiple spike/needle application leads to decrease of drag reduction benefits by comparison with the case of one central mounted needle at the front of a blunt body, but increase lift benefits.

Gilinsky, M.↗

NASA NTRS2020

Inter-species molecular attraction effect in the development of a two-species mixing layer

Simulations of fluid mixing occurring at high-pressures mandate the use of real-fluid equations of state to capture the effect of boh intermolecular repulsive forces and attractive forces which are neglected in perfect gases. Dedicated mixing rules account for these forces by providing a generic mapping of the mixture state-space onto a pure fluid domain. However, the specific inter-species molecular forces depend on the particular characteristics of the molecular pairs under consideration, and typically require the knowledge of binary interaction coefficients k_{αβ}. The impact of k_{αβ} on a flow field is as of now unclear; further, k_{αβ} is unknown for many molecular pairs and when unknown, set to be null. This study addresses the impact of k_{αβ} on the temporal evolution of a real-fluid mixing layer and uses Direct Numerical Simulation (DNS) to explore this impact. The results show differences between augmented attraction (k_{αβ}<0) and diminished attraction (k_{αβ}>0), affecting both fluid dynamics and diffusion processes. Specifically, in a binary mixing layer, it is observed that augmented attraction leads to delayed transition and mixing layer growth, manifested in the momentum layer thickness. The fractal dimension, which is a manifestation of interface corrugation, is calculated. These results provide an hitherto undocumented mechanism affecting the development of mixing layers.

Banuti, Daniel↗

NASA NTRS1985

SPAR improved structural-fluid dynamic analysis capability

The results of a study whose objective was to improve the operation of the SPAR computer code by improving efficiency, user features, and documentation is presented. Additional capability was added to the SPAR arithmetic utility system, including trigonometric functions, numerical integration, interpolation, and matrix combinations. Improvements were made in the EIG processor. A processor was created to compute and store principal stresses in table-format data sets. An additional capability was developed and incorporated into the plot processor which permits plotting directly from table-format data sets. Documentation of all these features is provided in the form of updates to the SPAR users manual.

Pearson, M. L.↗

NASA NTRS2016

Least Squares Shadowing Sensitivity Analysis of Chaotic Flow Around a Two-Dimensional Airfoil

Gradient-based sensitivity analysis has proven to be an enabling technology for many applications, including design of aerospace vehicles. However, conventional sensitivity analysis methods break down when applied to long-time averages of chaotic systems. This breakdown is a serious limitation because many aerospace applications involve physical phenomena that exhibit chaotic dynamics, most notably high-resolution large-eddy and direct numerical simulations of turbulent aerodynamic flows. A recently proposed methodology, Least Squares Shadowing (LSS), avoids this breakdown and advances the state of the art in sensitivity analysis for chaotic flows. The first application of LSS to a chaotic flow simulated with a large-scale computational fluid dynamics solver is presented. The LSS sensitivity computed for this chaotic flow is verified and shown to be accurate, but the computational cost of the current LSS implementation is high.

Blonigan, Patrick J.↗

DOE OSTI2021

Microscale Mechanical‐Chemical Modeling of Granular Salt: Insights for Creep

Abstract Microscale numerical modeling is potentially an effective approach for understanding salt creep, but it is quite challenging because natural salt grain boundaries can have complex and evolving geometry, and because contacts between these deformable salt mineral grains are dynamic as a result of compaction, chemical reaction, fluid flow, and heat transfer. In this study, we have overcome these challenges and developed a new microscale mechanical‐chemical (MC) model to analyze creep of salt at the microscale, accounting for coupled deformation, dynamic contacts, and chemical reaction in granular systems with realistic geometric representations. The MC model was realized by linking a new microscale mechanical code based on the numerical manifold method (NMM) to a reactive transport code named Crunch. We simulated the processes of reorganization of the salt grains, microfracturing, and pressure solution that contribute to creep at larger scales. Based on this first quantitative microscale model, we found that sharp corners of mineral grains can dominate the contact dynamics, microfracturing, and pressure solution when salt is compacted, thus governing the structural changes and porosity loss of the system. We found that pressure solution, which preferentially dissolves sharp corners and edges, can lead to relatively high porosity loss in the system, thus playing an important role in the creep of salt. Our analysis shows that the dynamic changes of salt granular systems involving grain relocation and pressure solution can occur repeatedly and continuously, thus contributing to longer‐term creep of salt at larger scales.

58 GEOSCIENCES↗

NASA NTRS2001

Novel Shapes of Miscible Interfaces Observed

The dynamics of miscible displacements in a cylindrical tube are being investigated experimentally and numerically, with a view to understand the complex processes that occur, for example, in enhanced oil recovery, hydrology, and filtration. We have observed complex shapes of the interface between two liquids that mix with each other when the less viscous liquid is displaced by the more viscous one in a tube. A less viscous fluid that displaces a more viscous fluid is known to propagate in the form of a "finger," and a flight experiment proposed by Maxworthy et al. to investigate the miscible-interface dynamics is currently being developed by NASA. From the current theory of miscible displacements, which was developed for a porous medium satisfying Darcy's law, it can be shown that in the absence of gravity the interface between the fluids is destabilized and thus susceptible to fingering only when a more viscous fluid is displaced by a less viscous one. Therefore, if the interface is initially flat and the more viscous fluid displaces the less viscous fluid, the interface ought to be stable and remain flat. However, numerical simulations by Chen and Meiburg for such displacement in a cylindrical tube show that the interface is unstable and a finger of the more viscous fluid is indeed formed. Preliminary experiments performed at the NASA Glenn Research Center show that not only can fingering occur when the more viscous fluid displaces a less viscous one in a cylindrical tube, but also that under certain conditions the advancing finger achieves a sinuous or snakelike shape. These experiments were performed using silicone oils in a vertical pipette of small diameter. In the initial configuration, the more viscous fluid rested on top of the less viscous one, and the interface was nominally flat. A dye was added to the upper liquid for ease of observation of the interface between the fluids. The flow was initiated by draining the lower fluid from the bottom of the pipette, at speeds less than 0.1 mm/sec.

Balasubramaniam, Ramaswamy↗

NASA NTRS1998

1-D Numerical Analysis of RBCC Engine Performance

An RBCC engine combines air breathing and rocket engines into a single engine to increase the specific impulse over an entire flight trajectory. Considerable research pertaining to RBCC propulsion was performed during the 1960's and these engines were revisited recently as a candidate propulsion system for either a single-stage-to-orbit (SSTO) or two-stage-to-orbit (TSTO) launch vehicle. There are a variety of RBCC configurations that had been evaluated and new designs are currently under development. However, the basic configuration of all RBCC systems is built around the ejector scramjet engine originally developed for the hypersonic airplane. In this configuration, a rocket engine plays as an ejector in the air-augmented initial acceleration mode, as a fuel injector in scramjet mode and the rocket in all rocket mode for orbital insertion. Computational fluid dynamics (CFD) is a useful tool for the analysis of complex transport processes in various components in RBCC propulsion systems. The objective of the present research was to develop a transient 1-D numerical model that could be used to predict flow behavior throughout a generic RBCC engine following a flight path.

Han, Samuel S.↗

NASA NTRS1994

Multidisciplinary propulsion simulation using the numerical propulsion system simulator (NPSS)

Implementing new technology in aerospace propulsion systems is becoming prohibitively expensive. One of the major contributions to the high cost is the need to perform many large scale system tests. The traditional design analysis procedure decomposes the engine into isolated components and focuses attention on each single physical discipline (e.g., fluid for structural dynamics). Consequently, the interactions that naturally occur between components and disciplines can be masked by the limited interactions that occur between individuals or teams doing the design and must be uncovered during expensive engine testing. This overview will discuss a cooperative effort of NASA, industry, and universities to integrate disciplines, components, and high performance computing into a Numerical propulsion System Simulator (NPSS).

Claus, Russel W.↗

DOE OSTI2023

N-particle irreducible actions for stochastic fluids

We construct one- and two-particle irreducible (1PI and 2PI) effective actions for the stochastic fluid dynamics of a conserved density undergoing diffusive motion. We compute the 1PI action in one-loop order and the 2PI action in two-loop approximation. We derive a set of Schwinger-Dyson equations and regularize the resulting equations using Pauli-Villars fields. We numerically solve the Schwinger-Dyson equations for a non-critical fluid. We find that higher-loop effects summed by the Schwinger-Dyson renormalize the non-linear coupling. We also find indications of a diffuson-cascade, the appearance of n-loop correction with smaller and smaller exponential suppression.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2023

Rheology of concentrated fiber suspensions with a load-dependent friction coefficient

Here we investigate the effects of fiber aspect ratio, roughness, flexibility, and volume fraction on the rheology of concentrated suspensions in a steady shear flow using direct numerical simulations. We model the fibers as inextensible continuous flexible slender bodies with the Euler-Bernoulli beam equation governing their dynamics suspended in an incompressible Newtonian fluid. The fiber dynamics and fluid flow coupling is achieved using the immersed boundary method. In addition, the fiber surface roughness might lead to interfiber contacts, resulting in normal and tangential forces between the fibers, which follow Coulomb’s law of friction. The surface roughness is modeled as hemispherical protrusions on the fiber surfaces. Their deformation results in a normal load-dependent friction coefficient. Our simulations accurately predict the experimentally observed shear thinning in fiber suspensions. Furthermore, we find that the suspension viscosity η increases with increasing the volume fraction, roughness, fiber rigidity, and aspect ratio. The increase in η is the macroscopic manifestation of a similar increase in the microscopic contact contribution to the total stress with these parameters. In addition, we observe positive and negative first N 1 and second N 2 normal stress differences, respectively, with |N 2 | < |N 1 |, in agreement with previous experiments. Last, we propose a modified Maron-Pierce law to quantify the reduction in the jamming volume fraction by increasing the fiber aspect ratio and roughness. Our results and analysis establish the use of fiber surface tribology to tune the suspension flow behavior.

42 ENGINEERING↗

DOE OSTI2024

Neutrino Trapping and Out-of-Equilibrium Effects in Binary Neutron-Star Merger Remnants

We study out-of-thermodynamic-equilibrium effects in neutron-star mergers with 3D general-relativistic neutrino-radiation large-eddy simulations. During mergers, the cores of the neutron stars remain cold (T ~ a few MeV) and out of thermodynamic equilibrium with trapped neutrinos originating from the hot collisional interface between the stars. However, within ~2 to 3 ms matter and neutrinos reach equilibrium everywhere in the remnant massive neutron star. Furthermore, our results show that dissipative effects, such as bulk viscosity, if present, are only active for a short window of time after the merger.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2025

LDRD Abbreviated report: High-Order General-Discrete-Ordinates Method Enabling Efficient Deterministic Transport in Hydrodynamic Simulations

Deterministic transport simulations for national-security and energy applications often operate in high-dimensional phase-space, where accuracy and cost both become major challenges. A common numerical artifact in such problems is the “ray-effect,” which appears as unphysical streaks. Beyond misinterpretation, these artifacts can contaminate tightly coupled physics, such as fluid dynamics, radiation-hydrodynamics, and laser-plasma interactions, eroding the predictive capability of entire multiphysics workflows. Our objective was to make high-dimension studies practical on modern hardware while mitigating the ray-effect without relying on prohibitively expensive sampling approaches such as Monte Carlo methods. We developed the Generic Discretization Library (GenDiL), a Graphics Processing Unit (GPU)-first framework that uses high-order Discontinuous Galerkin (DG) methods and matrix-free algorithms to reduce memory usage and improve computational efficiency, critical for phase-space simulations. GenDiL supports phase-space adaptivity in both mesh size and polynomial order (hp-adaptivity) to place resolution only where it is needed. A central capability is Local Dimensional Refinement (LDR), which couples lower-dimension continuum models to higher-dimension kinetic models through stable and conservative interfaces, so that high-fidelity physics is applied only in regions where it is essential. Building on the GenDiL framework, we developed the General SN (GSN) family of algorithms as a true generalization of the polar SN approach (discrete ordinates, often denoted SN). Rather than tying discrete ordinates to a specific polar change of coordinates, GSN formulates transport on an arbitrary change of coordinates chosen to reduce ray-effect. We studied two complementary variants: an analytic variant, where the coordinate map is prescribed in advance by a closed-form function; and a data-driven variant, where a quantity of interest, such as the net flux, guides the coordinate system. GenDiL provides the library infrastructure for efficient GPU execution, but the GSN concept is algorithmic and independent of any one library. Across representative high-dimension tests, including non-symmetric solutions, both variants delivered strong ray-effect mitigation at practical cost, moving four- to six-dimensional analysis toward repeatable, routine studies.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗