SEARCH · Engineering Papers

Results for “equation of 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 865 records · Page 48

All sources NASA NTRS DOE OSTI

NASA NTRS1985

Hybrid simulation techniques applied to the earth's bow shock

The application of a hybrid simulation model, in which the ions are treated as discrete particles and the electrons as a massless charge-neutralizing fluid, to the study of the earth's bow shock is discussed. The essentials of the numerical methods are described in detail; movement of the ions, solution of the electromagnetic fields and electron fluid equations, and imposition of appropriate boundary and initial conditions. Examples of results of calculations for perpendicular shocks are presented which demonstrate the need for a kinetic treatment of the ions to reproduce the correct ion dynamics and the corresponding shock structure. Results for oblique shocks are also presented to show how the magnetic field and ion motion differ from the perpendicular case.

Winske, D.↗

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 NTRS2003

High-Pressure Transport Properties Of Fluids: Theory And Data From Levitated Drops At Combustion-Relevant Temperatures

Although the high pressure multicomponent fluid conservation equations have already been derived and approximately validated for binary mixtures by this PI, the validation of the multicomponent theory is hampered by the lack of existing mixing rules for property calculations. Classical gas dynamics theory can provide property mixing-rules at low pressures exclusively. While thermal conductivity and viscosity high-pressure mixing rules have been documented in the literature, there is no such equivalent for the diffusion coefficients and the thermal diffusion factors. The primary goal of this investigation is to extend the low pressure mixing rule theory to high pressures and validate the new theory with experimental data from levitated single drops. The two properties that will be addressed are the diffusion coefficients and the thermal diffusion factors. To validate/determine the property calculations, ground-based experiments from levitated drops are being conducted.

Bellan, Josette↗

NASA NTRS2009

Fan Flutter Computations Using the Harmonic Balance Method

An experimental forward-swept fan encountered flutter at part-speed conditions during wind tunnel testing. A new propulsion aeroelasticity code, based on a computational fluid dynamics (CFD) approach, was used to model the aeroelastic behavior of this fan. This threedimensional code models the unsteady flowfield due to blade vibrations using a harmonic balance method to solve the Navier-Stokes equations. This paper describes the flutter calculations and compares the results to experimental measurements and previous results from a time-accurate propulsion aeroelasticity code.

Bakhle, Milind A.↗

NASA NTRS2007

Two-Dimensional Thermal Boundary Layer Corrections for Convective Heat Flux Gauges

This work presents a CFD (Computational Fluid Dynamics) study of two-dimensional thermal boundary layer correction factors for convective heat flux gauges mounted in flat plate subjected to a surface temperature discontinuity with variable properties taken into account. A two-equation k - omega turbulence model is considered. Results are obtained for a wide range of Mach numbers (1 to 5), gauge radius ratio, and wall temperature discontinuity. Comparisons are made for correction factors with constant properties and variable properties. It is shown that the variable-property effects on the heat flux correction factors become significant

Kandula, Max↗

NASA NTRS2014

Mean Flow Augmented Acoustics in Rocket Systems

Oscillatory motion in solid rocket motors and liquid engines has long been a subject of concern. Many rockets display violent fluctuations in pressure, velocity, and temperature originating from the complex interactions between the combustion process and gas dynamics. The customary approach to modeling acoustic waves inside a rocket chamber is to apply the classical inhomogeneous wave equation to the combustion gas. The assumption of a linear, non-dissipative wave in a quiescent fluid remains valid while the acoustic amplitudes are small and local gas velocities stay below Mach 0.2. The converging section of a rocket nozzle, where gradients in pressure, density, and velocity become large, is a notable region where this approach is not applicable. The expulsion of unsteady energy through the nozzle of a rocket is identified as the predominate source of acoustic damping for most rocket systems. An accurate model of the acoustic behavior within this region where acoustic modes are influenced by the presence of a steady mean flow is required for reliable stability predictions. Recently, an approach to address nozzle damping with mean flow effects was implemented by French [1]. This new approach extends the work originated by Sigman and Zinn [2] by solving the acoustic velocity potential equation (AVPE) formulated by perturbing the Euler equations [3]. The acoustic velocity potential (psi) describing the acoustic wave motion in the presence of an inhomogeneous steady high-speed flow is defined by, (del squared)(psi) − (lambda/c)(exp 2)(psi) − M(dot)[M(dot)(del)(del(psi))] − 2(lambda(M/c) + (M(dot)del(M))(dot)del(psi)−2(lambda)(psi)[M(dot)del(1/c)]=0 (1) with M as the Mach vector, c as the speed of sound, and lambda as the complex eigenvalue. French apply the finite volume method to solve the steady flow field within the combustion chamber and nozzle with inviscid walls. The complex eigenvalues and eigenvector are determined with the use of the ARPACK eigensolver. The present study employs the COMSOL Multphysics framework to solve the coupled eigenvalue problem using the finite element approach. The study requires one way coupling of the CFD High Mach Number Flow (HMNF) and mathematics module. The HMNF module evaluated the gas flow inside of a solid rocket motor using St. Robert's law modeling solid propellant burn rate, slip boundary conditions, and the supersonic outflow condition. Results from the HMNF model are used by the coefficient form of the mathematics module to determine the eigenvalues of the AVPE. The mathematics model is truncated at the nozzle sonic line, where a zero flux boundary condition is self-satisfying. The remaining boundaries are modeled with a zero flux boundary condition, assuming zero acoustic absorption on all surfaces. Pertinent results from these analyses are the complex valued eigenvalue and eigenvectors. Comparisons are made to the French results to evaluate the modeling approach. A comparison of the French results with that of the present analysis is displayed in figures 1 and 2, respectively. The graphic shows the first tangential eigenvector's real (a) and imaginary (b) values.

Fischbach, Sean R.↗

NASA NTRS2016

Modeling ARRM Xenon Tank Pressurization Using 1D Thermodynamic and Heat Transfer Equations

As a first step in understanding what ground support equipment (GSE) is required to provide external cooling during the loading of 5,000 kg of xenon into 4 aluminum lined composite overwrapped pressure vessels (COPVs), a modeling analysis was performed using Microsoft Excel. The goals of the analysis were to predict xenon temperature and pressure throughout loading at the launch facility, estimate the time required to load one tank, and to get an early estimate of what provisions for cooling xenon might be needed while the tanks are being filled. The model uses the governing thermodynamic and heat transfer equations to achieve these goals. Results indicate that a single tank can be loaded in about 15 hours with reasonable external coolant requirements. The model developed in this study was successfully validated against flight and test data. The first data set is from the Dawn mission which also utilizes solar electric propulsion with xenon propellant, and the second is test data from the rapid loading of a hydrogen cylindrical COPV. The main benefit of this type of model is that the governing physical equations using bulk fluid solid temperatures can provide a quick and accurate estimate of the state of the propellant throughout loading which is much cheaper in terms of computational time and licensing costs than a Computation Fluid Dynamics (CFD) analysis while capturing the majority of the thermodynamics and heat transfer.

Xenon storage↗

NASA NTRS2017

Modeling Xenon Tank Pressurization using One-Dimensional Thermodynamic and Heat Transfer Equations

Thermodynamic tank modeling↗

NASA NTRS2017

Implementation of a Body Force Model in OVERFLOW for Propulsor Simulations

We present an implementation of a propulsor model based on body force method into the OVERFLOW computational fluid dynamics code to model turbofan engines and propulsors of similar type. The model estimates the forces imparted on the fluid by the blade camber surfaces as a body force source terms in the momentum and energy equations over grids that represent the rotor and stator of a fan stage. We tested the implementation on three cases: 1) Source Diagnostics Test (SDT) fan with R4 rotors, which has extensive test data on aerodynamic performance and rotor wake surveys, 2) A stand-alone Aeronaut TF8000 propulsor which is an off-the-shelf propulsor used on model aircraft and 3) The D8 aircraft model with TF8000 propulsors placed in a wind tunnel. Despite missing some of the features in the velocity profiles near the endwalls (i.e. hub and casing surfaces), the OVERFLOW simulations with body force model predicted area-averaged flow speed and total pressure rise through the SDT propulsor within a few percent of the LDV (Laser Doppler Velocimetry) measurements. In the case of TF8000 propulsor on the D8 airframe, the model under-predicted mechanical power coefficient by several percent of the wind tunnel test results when the horizontal force balance condition over the airframe is targeted by tuning rotor speed. By investigating the upstream influence of the rotor swirl, it was found out that the induced swirl velocity effects upstream were relatively small and they rapidly vanished before reaching one fan diameter upstream of the fan face. The body force model provided insights on aerothermodynamics and aeromechanics of boundary layer ingesting propulsor; these insights could not be obtained by using the uniform pressure jump model.

Propulsor↗

NASA NTRSDate not supplied

Performance of Coupled Physics Solvers for Multidisciplinary Hypersonic Flow Simulations on Several Classes of Computer Architectures

The application of hypersonic flow simulation tools to realistic flight scenarios will require the coupling of multiple physical effects to the baseline fluid dynamics. Such multiphysics effects can include the aerooelastic response of the airframe or engine components, dynamic transport of atmospheric particles, the deformation of solid-fluid interfaces that can ablate, pyrolyze, or erode, as well as a host of other processes, all of which are governed by unique sets of physical equations and models. Coupling multiple (and potentially disparate) physics solvers to a robust compressible flow solver poses additional challenges related to the stability, performance and scalability of the combined solver. The choices made during the software design process can therefore lead to a variation in simulation efficiency across different computer architectures. In this paper, we will consider two representative multiphysics hypersonic flow scenarios: the interaction of solid particulates with the flow field created by a hypersonic lifting body and the aerooelastic deformation of a model airframe under high-Mach-number flow conditions. For these simulations we explore the behavior of several hypersonic simulation tools, including Kestrel, FUN3D, US3D, and JENRE multiphysics framework, on several high performance computing systems containing various CPU and GPU architectures.

architecture↗

NASA NTRS1996

Spray Combustion Modeling with VOF and Finite-Rate Chemistry

A spray atomization and combustion model is developed based on the volume-of-fluid (VOF) transport equation with finite-rate chemistry model. The gas-liquid interface mass, momentum and energy conservation laws are modeled by continuum surface force mechanisms. A new solution method is developed such that the present VOF model can be applied for all-speed range flows. The objectives of the present study are: (1) to develop and verify the fractional volume-of-fluid (VOF) cell partitioning approach into a predictor-corrector algorithm to deal with multiphase (gas-liquid) free surface flow problems; (2) to implement the developed unified algorithm in a general purpose computational fluid dynamics (CFD) code, Finite Difference Navier-Stokes (FDNS), with droplet dynamics and finite-rate chemistry models; and (3) to demonstrate the effectiveness of the present approach by simulating benchmark problems of jet breakup/spray atomization and combustion. Modeling multiphase fluid flows poses a significant challenge because a required boundary must be applied to a transient, irregular surface that is discontinuous, and the flow regimes considered can range from incompressible to highspeed compressible flows. The flow-process modeling is further complicated by surface tension, interfacial heat and mass transfer, spray formation and turbulence, and their interactions. The major contribution of the present method is to combine the novel feature of the Volume of Fluid (VOF) method and the Eulerian/Lagrangian method into a unified algorithm for efficient noniterative, time-accurate calculations of multiphase free surface flows valid at all speeds. The proposed method reformulated the VOF equation to strongly couple two distinct phases (liquid and gas), and tracks droplets on a Lagrangian frame when spray model is required, using a unified predictor-corrector technique to account for the non-linear linkages through the convective contributions of VOF. The discontinuities within the sharp interface will be modeled as a volume force to avoid stiffness. Formations of droplets, tracking of droplet dynamics and modeling of the droplet breakup/evaporation, are handled through the same unified predictor-corrector procedure. Thus the new algorithm is non-iterative and is flexible for general geometries with arbitrarily complex topology in free surfaces. The FDNS finite-difference Navier-Stokes code is employed as the baseline of the current development. Benchmark test cases of shear coaxial LOX/H2 liquid jet with atomization/combustion and impinging jet test cases are investigated in the present work. Preliminary data comparisons show good qualitative agreement between data and the present analysis. It is indicative from these results that the present method has great potential to become a general engineering design analysis and diagnostics tool for problems involving spray combustion.

Chen, Yen-Sen↗

NASA NTRS1980

Preliminary results from a four-working space, double-acting piston, Stirling engine controls model

A four working space, double acting piston, Stirling engine simulation is being developed for controls studies. The development method is to construct two simulations, one for detailed fluid behavior, and a second model with simple fluid behaviour but containing the four working space aspects and engine inertias, validate these models separately, then upgrade the four working space model by incorporating the detailed fluid behaviour model for all four working spaces. The single working space (SWS) model contains the detailed fluid dynamics. It has seven control volumes in which continuity, energy, and pressure loss effects are simulated. Comparison of the SWS model with experimental data shows reasonable agreement in net power versus speed characteristics for various mean pressure levels in the working space. The four working space (FWS) model was built to observe the behaviour of the whole engine. The drive dynamics and vehicle inertia effects are simulated. To reduce calculation time, only three volumes are used in each working space and the gas temperature are fixed (no energy equation). Comparison of the FWS model predicted power with experimental data shows reasonable agreement. Since all four working spaces are simulated, the unique capabilities of the model are exercised to look at working fluid supply transients, short circuit transients, and piston ring leakage effects.

Daniele, C. J.↗

DOE OSTI2024

Stochastic finite volume method for uncertainty quantification of transient flow in gas pipeline networks

We develop a weakly intrusive framework to simulate the propagation of uncertainty in solutions of generic hyperbolic partial differential equation systems on graph-connected domains with nodal coupling and boundary conditions. The method is based on the Stochastic Finite Volume (SFV) approach and can be applied for uncertainty quantification (UQ) of the dynamical state of fluid flow over actuated transport networks. The numerical scheme has specific advantages for modeling intertemporal uncertainty in time-varying boundary parameters, which cannot be characterized by strict upper and lower (interval) bounds. We describe the scheme for a single pipe, and then formulate the controlled junction Riemann problem (JRP) that enables the extension to general network structures. In conclusion, we demonstrate the method's capabilities and performance characteristics using a standard benchmark test network.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2023

Application of a Distributed Element Roughness Model to Additively Manufactured Internal Cooling Channels

Abstract Design for cooling effectiveness in turbine blades relies on accurate models for dynamic losses and heat transfer of internal cooling passages. Metal additive manufacturing (AM) has expanded the design space for these configurations, but can give rise to large-scale roughness features. The range of roughness length scales in these systems makes morphology resolved computational fluid dynamics (CFD) impractical. However, volumetric roughness models can be leveraged, as they have computational costs orders of magnitude lower. In this work, a discrete element roughness model (DERM), based on the double-averaged Navier–Stokes equations, is presented and applied to additively manufactured rough channels, representative of gas turbine blade cooling passages. Unique to this formulation of DERM is a generalized sheltering-based treatment of drag, a two-layer model for spatially averaged Reynolds stresses, and explicit treatment of dispersion. Six different AM rough surface channel configurations are studied, with roughness trough to peak sizes ranging from 15% to 60% nominal channel passage half-width, and the roughness Reynolds number ranges from Rek = 60 to 300. DERM predictions for spatially and temporally averaged mean flow quantities are compared to previously reported direct numerical simulation results. Good agreement in the mean velocity profiles, stress balances, and drag partitions are observed. While DERM models are typically calibrated to specific deterministic roughness morphologies at comparatively small roughness Reynolds numbers, the present more generalized DERM formulation has wider applicability. Here, it is demonstrated that the model can accommodate random roughness of large scale, typical of AM.

Engineering↗

NASA NTRS1976

On the flow in an annulus surrounding a whirling cylinder

When fluid in an annulus between two cylinders is set in motion by whirling movements of one or both of the cylinders, dynamic forces are imposed by the fluid on the cylinders. Knowledge of these forces is frequently important, indeed often critical, to the engineer designing rotor systems or journal bearings. Quite general solutions of the Navier-Stokes equations are presented for this problem and are limited only by restrictions on the amplitude of the whirl motion. From these solutions, the forces are derived under a wide variety of circumstances, including large and small annular widths, high and low Reynolds numbers, and the presence and absence of a mean flow created by additional net rotation of one or both of the cylinders.

Brennen, C.↗

NASA NTRS1982

Fluid dynamical description of relativistic nuclear collisions

On the basis of both a conventional relativistic nuclear fluid dynamic model and a two fluid generalization that takes into account the interpenetration of the target and projectile upon contact, collisions between heavy nuclei moving at relativistic speeds are calculated. This is done by solving the relevant equations of motion numerically in three spatial dimensions by use of particle in cell finite difference computing techniques. The effect of incorporating a density isomer, or quasistable state, in the nuclear equation of state at three times normal nuclear density, and the effect of doubling the nuclear compressibility coefficient are studied. For the reaction 20Ne + 238U at a laboratory bombarding energy per nucleon of 393 MeV, the calculated distributions in energy and angle of outgoing charged particles are compared with recent experimental data both integrated over all impact parameters and for nearly central collisions.

Nix, J. R.↗

NASA NTRS1984

Alfven waves in a thermally stratified fluid

The properties of Alfven waves propagating along a uniform horizontal field in a highly conducting incompressible medium in the presence of strong convective instability are examined in the Boussinesq approximation. In particular, it is sought to determine whether there are exact solutions to the dynamical equations in the presence of convective forces. It is shown that a class of exact solutions of arbitrary amplitude, but of limited form, which may be of some physical interest, does exist. For large amplitudes, any mixtures of polarization states are shown to cause scattering into new modes.

Parker, E. N.↗

NASA NTRS1989

Computer Simulation of Aircraft Aerodynamics

The role of Ames Research Center in conducting basic aerodynamics research through computer simulations is described. The computer facilities, including supercomputers and peripheral equipment that represent the state of the art, are described. The methodology of computational fluid dynamics is explained briefly. Fundamental studies of turbulence and transition are being pursued to understand these phenomena and to develop models that can be used in the solution of the Reynolds-averaged Navier-Stokes equations. Four applications of computer simulations for aerodynamics problems are described: subsonic flow around a fuselage at high angle of attack, subsonic flow through a turbine stator-rotor stage, transonic flow around a flexible swept wing, and transonic flow around a wing-body configuration that includes an inlet and a tail.

Inouye, Mamoru↗