Search NASA⌕ Search

SEARCH · Search NASA

Results for “Solving Fluid”

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 181 records · Page 10

High accuracy solutions of incompressible Navier-Stokes equations

In recent years, high accuracy finite difference approximations were developed for partial differential equations of elliptic type, with particular emphasis on the convection-diffusion equation. These approximations are of compact type, have a local truncation error of fourth order, and allow the use of standard iterative schemes to solve the resulting systems of algebraic equations. These high accuracy approximations are extended to the solution of Navier-Stokes equations. Solutions are obtained for the model problem of driven cavity and are compared with solutions obtained using other approximations and those obtained by other authors. It is discovered that the high order approximations do indeed produce high accuracy solutions and have a potential for use in solving important problems of viscous fluid flows.

Gupta, Murli M.↗

Performance of a parallel code for the Euler equations on hypercube computers

The performance of hypercubes were evaluated on a computational fluid dynamics problem and the parallel environment issues were considered that must be addressed, such as algorithm changes, implementation choices, programming effort, and programming environment. The evaluation focuses on a widely used fluid dynamics code, FLO52, which solves the two dimensional steady Euler equations describing flow around the airfoil. The code development experience is described, including interacting with the operating system, utilizing the message-passing communication system, and code modifications necessary to increase parallel efficiency. Results from two hypercube parallel computers (a 16-node iPSC/2, and a 512-node NCUBE/ten) are discussed and compared. In addition, a mathematical model of the execution time was developed as a function of several machine and algorithm parameters. This model accurately predicts the actual run times obtained and is used to explore the performance of the code in interesting but yet physically realizable regions of the parameter space. Based on this model, predictions about future hypercubes are made.

Barszcz, Eric↗

A finite element model of conduction, convection, and phase change near a solid/melt interface

Detailed understanding of heat transfer and fluid flow is required for many aerospace thermal systems. These systems often include phase change and operate over a range of accelerations or effective gravitational fields. An approach to analyzing such systems is presented which requires the simultaneous solution of the conservation laws of energy, momentum, and mass, as well as an equation of state. The variable property form of the governing equations are developed in two-dimensional Cartesian coordinates for a Newtonian fluid. A numerical procedure for solving the governing equations is presented and implemented in a computer program. The Galerkin form of the finite element method is used to solve the spatial variation of the field variables, along with the implicit Crank-Nicolson time marching algorithm. Quadratic Langrangian elements are used for the internal energy and the two components of velocity. Linear Lagrangian elements are used for the pressure. The location of the solid/liquid interface as well as the temperatures are determined form the calculated internal energy and pressure. This approach is quite general in that it can describe heat transfer without phase change, phase change with a sharp interface, and phase change without an interface. Analytical results from this model are compared to those of other researchers studying transient conduction, convection, and phase change and are found to be in good agreement. The numerical procedure presented requires significant computer resources, but this is not unusual when compared to similar studies by other researchers. Several methods are suggested to reduce the computational times.

Viterna, Larry A.↗

A global analysis of the electrodynamic interactions between a space station and the ionospheric plasma environment

A general analysis of the electrodynamic interactions between a space station with two biased platforms and the ionospheric plasma is presented. This problem can be separated into a far-field problem, concerned with the electromagnetic interference surrounding the entire space station, and a near-field problem, concentrated on the interactions in the vicinity of the biased platforms. The far-field problem is solved by application of plasma fluid theory. The space station will generate a radiation field composed mainly of the Alfven waves. This far-field radiation depends on the details of the near-field current collection. Computer particle simulations were performed in the near-field of the biased platform to study the plasma flow field, the sheath structure and the current collection. Approximate analytical solutions to the near-field are also obtained. The far-field and near-field solutions are coupled to provide a global description of the electrodynamic interactions.

Wang, J.↗

Finite element computation of compressible flows with the SUPG formulation

Finite element computation of compressible Euler equations is presented in the context of the streamline-upwind/Petrov-Galerkin (SUPG) formulation. The SUPG formulation, which is based on adding stabilizing terms to the Galerkin formulation, is further supplemented with a shock capturing operator which addresses the difficulty in maintaining a satisfactory solution near discontinuities in the solution field. The shock capturing operator, which has been derived from work done in entropy variables for a similar operator, is shown to lead to an appropriate level of additional stabilization near shocks, without resulting in excessive numerical diffusion. An implicit treatment of the impermeable wall boundary condition is also presented. This treatment of the no-penetration condition offers increased stability for large Courant numbers, and accelerated convergence of the computations for both implicit and explicit applications. Several examples are presented to demonstrate the ability of this method to solve the equations governing compressible fluid flow.

Le Beau, G. J.↗

Pressure updating methods for the steady-state fluid equations

We consider the steady state equations for a compressible fluid. Since we wish to solve for a range of speeds we must consider the equations in conservation form. For transonic speeds these equations are of mixed type. Hence, the usual approach is to add time derivatives to the steady state equations and then march these equations in time. One then adds a time derivative of the density to the continuity equation, a derivative of the momentum to the momentum equation and a derivative of the total energy to the energy equation. This choice is dictated by the time consistent equations. However, since we are only interested in the steady state this is not necessary. Thus we shall consider the possibility of adding a time derivative of the pressure to the continuity equation and similar modifications for the energy equation. This can then be generalized to adding combinations of time derivatives to each equation since these vanish in the steady state. When using acceleration techniques such as residual smoothing, multigrid, etc. these are applied to the pressure rather than the density. Hence, the code duplicates the behavior of the incompressible equations for low speeds.

Fiterman, A.↗

Computation of axisymmetric and ionized hypersonic flows using particle and continuum methods

Comparisons between particle and continuum simulations of hypersonic near-continuum flows are presented. The particle approach employs the direct simulation Monte Carlo (DSMC) method, and the continuum approach solves the appropriate equations of fluid flow. Both simulations have thermochemistry models for air implemented including ionization. A new axisymmetric DSMC code that is efficiently vectorized is developed for this study. In this DSMC code, particular attention is paid to matching the relaxation rates employed in the continuum approach. This investigation represents a continuum of a previous study that considered thermochemical relaxation in one-dimensional shock waves of nitrogen. Comparison of the particle and continuum methods is first made for an axisymmetric blunt-body flow of air at 7 km/s. Very good agreement is obtained for the two solutions. The two techniques also compare well for a one-dimensional shock wave in air at 10 km/s. In both applications, the results are found to be sensitive to various aspects of the chemistry model employed.

Boyd, Iain D.↗

Multistage Schemes with Multigrid for Euler and Navier-Strokes Equations: Components and Analysis

A class of explicit multistage time-stepping schemes with centered spatial differencing and multigrids are considered for the compressible Euler and Navier-Stokes equations. These schemes are the basis for a family of computer programs (flow codes with multigrid (FLOMG) series) currently used to solve a wide range of fluid dynamics problems, including internal and external flows. In this paper, the components of these multistage time-stepping schemes are defined, discussed, and in many cases analyzed to provide additional insight into their behavior. Special emphasis is given to numerical dissipation, stability of Runge-Kutta schemes, and the convergence acceleration techniques of multigrid and implicit residual smoothing. Both the Baldwin and Lomax algebraic equilibrium model and the Johnson and King one-half equation nonequilibrium model are used to establish turbulence closure. Implementation of these models is described.

Swanson, R. C.↗

On the Coupling of CDISC Design Method with FPX Rotor Code

A rotor section aerodynamics design package is developed by coupling Constrained Direct Iterative Surface Curvature (CDISC) design method with the FPX rotor code. The coupling between the CDISC design and the FPX flow analysis is fully automated. The CDISC design method employs a predictor-corrector procedure iteratively to determine a surface geometry which produces a target pressure distribution, where the target pressure distributions is either pre-defined or automatically generated through flow and geometry constraints. The FPX code is an eXtended Full-Potential rotor Computational Fluid Dynamics (CFD) code, which solves the three-dimensional unsteady full-potential equation in a strong conservative form using an implicit approximate factorization finite-difference scheme with entropy and viscosity corrections. Application of the CDISC design method coupled with the FPX rotor code is made for rotor blades in hovering motions. Several design examples are presented to demonstrate the capability of the new package in rotor section design.

Hu, Hong↗

Computational Analyses of Propulsion Aeroacoustics for Mixed Flow Nozzle Pylon Installation at Takeoff

A CFD analyses is presented for a set of baseline and noise suppression mixed flow nozzles with and without a pylon installation. The five model configurations are as follows; a baselinecore/fan dual-stream nozzle with an external plug, a chevron mixer nozzle with a peak on the symmetry plane with external plug, both of the above nozzles with an installed bifurcatingpylon and lastly a clocked chevron mixer nozzle such that a trough is aligned with the center of the pylon. The fluid flow is simulated by solving the asymptotically steady, compressible, Reynolds-averaged Navier-Stokes equations using an implicit, up-wind, flux-difference splitting finite volume scheme and standard two equation k-epsilon turbulence model with a linear stress representation. All computations are performed using the multiblock, parallel, structuredcode PAB3D. Results indicate that the clocked chevron with pylon case achieves the most optimal levels of average and peak turbulence kinetic energy and vorticity and therefore is expected to be the quietest of the five configurations tested. Further study is required to refine expressions which are indicative of noise and mate these with rigorous noise prediction models.

Massey, Steven J.↗

Survey of Aerothermodynamics Facilities Useful for the Design of Hypersonic Vehicles Using Air-Breathing Propulsion

The dream of producing an air-breathing, hydrogen fueled, hypervelocity aircraft has been before the aerospace community for decades. However, such a craft has not yet been realized, even in an experimental form. Despite the simplicity and beauty of the concept, many formidable problems must be overcome to make this dream a reality. This paper summarizes the aero/aerothermodynamic issues that must be addressed to make the dream a reality and discusses how aerothermodynamics facilities and their modem companion, real-gas computational fluid dynamics (CFD), can help solve the problems blocking the way to realizing the dream. The approach of the paper is first to outline the concept of an air-breathing hypersonic vehicle and then discuss the nose-to-tail aerothermodynamics issues and special aerodynamic problems that arise with such a craft. Then the utility of aerothermodynamic facilities and companion CFD analysis is illustrated by reviewing results from recent United States publications wherein these problems have been addressed. Papers selected for the discussion have k e n chosen such that the review will serve to survey important U.S. aero/aerothermodynamic real gas and conventional wind tunnel facilities that are useful in the study of hypersonic, hydrogen propelled hypervelocity vehicles.

Arnold, James O.↗

Workload Characterization of CFD Applications Using Partial Differential Equation Solvers

Workload characterization is used for modeling and evaluating of computing systems at different levels of detail. We present workload characterization for a class of Computational Fluid Dynamics (CFD) applications that solve Partial Differential Equations (PDEs). This workload characterization focuses on three high performance computing platforms: SGI Origin2000, EBM SP-2, a cluster of Intel Pentium Pro bases PCs. We execute extensive measurement-based experiments on these platforms to gather statistics of system resource usage, which results in workload characterization. Our workload characterization approach yields a coarse-grain resource utilization behavior that is being applied for performance modeling and evaluation of distributed high performance metacomputing systems. In addition, this study enhances our understanding of interactions between PDE solver workloads and high performance computing platforms and is useful for tuning these applications.

Waheed, Abdul↗

Inflow/Outflow Boundary Conditions with Application to FUN3D

Several boundary conditions that allow subsonic and supersonic flow into and out of the computational domain are discussed. These boundary conditions are demonstrated in the FUN3D computational fluid dynamics (CFD) code which solves the three-dimensional Navier-Stokes equations on unstructured computational meshes. The boundary conditions are enforced through determination of the flux contribution at the boundary to the solution residual. The boundary conditions are implemented in an implicit form where the Jacobian contribution of the boundary condition is included and is exact. All of the flows are governed by the calorically perfect gas thermodynamic equations. Three problems are used to assess these boundary conditions. Solution residual convergence to machine zero precision occurred for all cases. The converged solution boundary state is compared with the requested boundary state for several levels of mesh densities. The boundary values converged to the requested boundary condition with approximately second-order accuracy for all of the cases.

Carlson, Jan-Renee↗

Bioreactors Drive Advances in Tissue Engineering

It was an unlikely moment for inspiration. Engineers David Wolf and Ray Schwarz stopped by their lab around midday. Wolf, of Johnson Space Center, and Schwarz, with NASA contractor Krug Life Sciences (now Wyle Laboratories Inc.), were part of a team tasked with developing a unique technology with the potential to enhance medical research. But that wasn t the focus at the moment: The pair was rounding up colleagues interested in grabbing some lunch. One of the lab s other Krug engineers, Tinh Trinh, was doing something that made Wolf forget about food. Trinh was toying with an electric drill. He had stuck the barrel of a syringe on the bit; it spun with a high-pitched whirr when he squeezed the drill s trigger. At the time, a multidisciplinary team of engineers and biologists including Wolf, Schwarz, Trinh, and project manager Charles D. Anderson, who formerly led the recovery of the Apollo capsules after splashdown and now worked for Krug was pursuing the development of a technology called a bioreactor, a cylindrical device used to culture human cells. The team s immediate goal was to grow human kidney cells to produce erythropoietin, a hormone that regulates red blood cell production and can be used to treat anemia. But there was a major barrier to the technology s success: Moving the liquid growth media to keep it from stagnating resulted in turbulent conditions that damaged the delicate cells, causing them to quickly die. The team was looking forward to testing the bioreactor in space, hoping the device would perform more effectively in microgravity. But on January 28, 1986, the Space Shuttle Challenger broke apart shortly after launch, killing its seven crewmembers. The subsequent grounding of the shuttle fleet had left researchers with no access to space, and thus no way to study the effects of microgravity on human cells. As Wolf looked from Trinh s syringe-capped drill to where the bioreactor sat on a workbench, he suddenly saw a possible solution to both problems. It dawned on me that rotating the wall of the reactor would solve one of our fundamental fluid mechanical problems, specifically by removing the velocity gradient of the tissue culture fluid media near the reactor s walls, says Wolf. It looked as though it would allow us to suspend the growing cells within the reactor without introducing turbulent fluid mechanical conditions.

Source record↗

Coupled CFD-PBE Predictions of Renal Stone Size Distributions in the Nephron in Microgravity

In this paper, a deterministic model is developed to assess the risk of critical renal stone formation for astronauts during space travel. A Population Balance Equation (PBE) model is used to compute the size distribution of a population of nucleating, growing and agglomerating renal calculi as they are transported through different sections of the nephron. The PBE model is coupled to a Computational Fluid Dynamics (CFD) model that solves for steady state flow of urine and transport of renal calculi along with the concentrations of ionic species, calcium and oxalate, in the nephron using an Eulerian two-phase mathematical framework. Parametric simulation are performed to study stone size enhancement and steady state volume fraction distributions in the four main sections of the nephron under weightlessness conditions. Contribution of agglomeration to the stone size distribution and effect of wall friction on the stone volume fraction distributions are carefully examined. Case studies using measured astronaut urinary calcium and oxalate concentrations in microgravity as input indicate that under nominal conditions the largest stone sizes developed in Space will be still considerably below the critical range for problematic stone development. However, results also indicate that the highest stone volume fraction occurs next to the tubule and duct walls. This suggests that there is an increased potential for wall adhesion with the possibility of evolution towards critical stone sizes.

Microgravity↗

Coupled CFD-PBE Predictions of Renal Stone Size Distributions in the Nephron in Microgravity

In this paper, a deterministic model is developed to assess the risk of critical renal stone formation for astronauts during space travel. A Population Balance Equation (PBE) model is used to compute the size distribution of a population of nucleating, growing and agglomerating renal calculi as they are transported through different sections of the nephron. The PBE model is coupled to a Computational Fluid Dynamics (CFD) model that solves for steady state flow of urine and transport of renal calculi along with the concentrations of ionic species, calcium and oxalate, in the nephron using an Eulerian two-phase mathematical framework. Parametric simulation are performed to study stone size enhancement and steady state volume fraction distributions in the four main sections of the nephron under weightlessness conditions. Contribution of agglomeration to the stone size distribution and effect of wall friction on the stone volume fraction distributions are carefully examined. Case studies using measured astronaut urinary calcium and oxalate concentrations in microgravity as input indicate that under nominal conditions the largest stone sizes developed in Space will be still considerably below the critical range for problematic stone development. However, results also indicate that the highest stone volume fraction occurs next to the tubule and duct walls. This suggests that there is an increased potential for wall adhesion with the possibility of evolution towards critical stone sizes.

Crystal Growth↗

Galactic cosmic-ray mediation of a spherical solar wind flow. 1: The steady state cold gas hydrodynamical approximation

Realistic models of the outer heliosphere should consider that the interstellar cosmic-ray pressure becomes comparable to pressures in the solar wind at distances more than 100 AU from the Sun. The cosmic-ray pressure dynamically affects solar wind flow through deceleration. This effect, which occurs over a scale length of the order of the effective diffusion length at large radial distances, has important implications for cosmic-ray modulation and acceleration. As a first step toward solution of this nonlinear problem, a steady state numerical model was developed for a relatively cold spherical solar wind flow which encounters the confining isotropic pressure of the surrounding Galactic medium. This pressure is assumed to be dominated by energetic particles (Galactic cosmic rays). The system of equations, which are solved self-consistently, includes the relevant hydrodynamical equations for the solar wind flow and the spherical cosmic-ray transport equation. To avoid the closure parameter problem of the two-fluid model, the latter equation is solved for the energy-dependent cosmic-ray distribution function.

Le Roux, J. A.↗

A coupled fluid-structure dynamic analysis for water impact loads

In this paper a nonlinear transient, hydroelastic method is developed for response analysis of the Space Shuttle recoverable booster on water impact. At each time step, coupled equations of fluid flow and structure dynamics are solved with an iterative, self-correcting process. Rigid body motions are represented by large angle equations. Local deformations are represented by a NASTRAN-generated model of the three-dimensional structure. The fluid is represented by series solutions of the potential flow equation which include the effects of the local structure motions. The resulting water loads and structural stresses obtained will aid in qualification of the final structural design. Comparisons are made to experimental impact data to validate the method.

Herting, D. N.↗