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 307 records · Page 17

Magnetohydrodynamic and gasdynamic aspects of solar-wind flow around terrestrial planets: A critical review

Problems in space physics are discussed whose models, in simplified form, reduce to a supersonic flow scheme with a detached shock wave, namely: (1) solar wind interaction with an intrinsic planetary magnetic field; (2) solar wind interaction with the ionized component of the atmosphere of a comet; and (3) solar wind interaction with the ionosphere of a planet which does not possess its own magnetic field. The numerical study of the above problems is performed with the use of magnetogasdynamic equations for an ideal single fluid model. From the physical veiwpoint, the problems are solved in terms of as simple phenomena as possible; the principal objective is to make recently developed methods of numerical analysis of mixed flows applicable to space physics problems.

Spreiter, J. R.↗

A three-dimensional model of corotating streams in the solar wind. 1: Theoretical foundations

The theoretical and mathematical background pertinent to the study of steady, corotating solar wind structure in all three spatial dimensions (3-D) is discussed. The dynamical evolution of the plasma in interplanetary space (defined as the region beyond roughly 35 solar radii where the flow is supersonic) is approximately described by the nonlinear, single fluid, polytropic (magneto-) hydrodynamic equations. Efficient numerical techniques for solving this complex system of coupled, hyperbolic partial differential equations are outlined. The formulation is inviscid and nonmagnetic, but methods allow for the potential inclusion of both features with only modest modifications. One simple, highly idealized, hydrodynamic model stream is examined to illustrate the fundamental processes involved in the 3-D dynamics of stream evolution. Spatial variations in the rotational stream interaction mechanism were found to produce small nonradial flows on a global scale that lead to the transport of mass, energy, and momentum away from regions of relative compression and into regions of relative rarefaction.

Pizzo, V. J.↗

Conjugated heat transfer in a laminar boundary layer with heat source at the wall

This work presents the solution of the temperature field in a two-dimensional laminar incompressible flow over a conducting solid plate with a line heat source located at the fluid-solid interface perpendicular to the flow direction. A numerical scheme was used to obtain the temperature profiles as a function of the source strength, and of the properties of the fluid and the solid. The heat conduction and forced convection in the fluid and the heat conduction in the solid were solved for the case of moderate temperature rise, where the assumption of constant properties applies. The model enables the improvement of an instrument for the detection of boundary layer separation. It was found that for the actual parameters of the separation detector, a distance of 4 to 24 mm between the sensors gives an indication of 70 percent of the maximum temperature difference.

Brosh, A.↗

Solar transition region response to variations in the heating rate

The response of a numerical model for the upper chromosphere, transition region, and corona to variations in the energy input has been examined. The numerical model solves the set of one-dimensional two-fluid hydrodynamic equations in a simple vertical magnetic flux tube. The atmosphere responds to both the increase and decrease in energy deposition by smoothly readjusting the temperature gradient and the amount of material in the region of peak radiating efficiency to radiate away energy being deposited. At no time during this readjustment is a departure from a thin laminar transition region structure seen. In addition, a time-dependent description of the nonequilibrium ionization of all of the ionization stages of oxygen has been included. This calculation is coupled with the self-consistent calculations of the dynamical variables. It is found that the nonequilibrium ionization balance calculations for both heating and cooling small loops in the quiet sun predict relative ionic abundances which differ substantially from those which would be predicted by an equilibrium calculation

Mariska, J. T.↗

Programming substructure computations for elliptic problems on the CHiP system

A number of studies have been conducted with the aim to apply parallel computation to problems associated with solving finite element equations arising in structural mechanics and fluid dynamics. These studies have provided many important results. The present investigation is concerned with a set of experiments designed to test two ideas, including configurability and substructuring. The considered algorithms and tests are intended for implementation on the Configurable, Highly Parallel (CHiP) family of architecture described by Snyder (1982). The ChiP computer is composed of homogeneous processing elements (PEs) placed at regular intervals in a lattice of programmable switches. Two examples of the role of configurability and substructuring for simple iterative algorithms are considered, giving attention to conjugate gradient iterations, and tridiagonal systems of equations.

Gannon, D.↗

SINDA - Systems Improved Numerical Differencing Analyzer

Capability to model fluid flow added. SINDA is software system developed to solve physical problems governed by diffusion-type equations mathematically modeled by lumped-parameter representations. Used as general thermal analyzer with resistor-and-capacitor-network representations. Adapted to wide range of problems represented by Fourier, Possion, Laplace, or other differential equations. Written in FORTRAN and Assembler.

Damico, Steve J.↗

Algebraic Generation Of Two-Dimensional Grids

Two-Boundary Grid Generation (TBGG) computer program applies interactive algebraic grid-generation technique in two dimensions. Incorporates algebraic equations that relate uniform rectangular computational domain to arbitrary physical domain. Hermite cubic interpolation performed between fixed, nonintersecting boundaries. Applied to variety of problems in which discrete techniques used to solve such partial differential equations as governing equations of fluid flow. Written in FORTRAN 77.

Smith, Robert E.↗

Mixing and non-equilibrium chemical reaction in a compressible mixing layer

The effects of compressibility, chemical reaction exothermicity, and non-equilibrium chemical modeling in a reacting plane mixing layer were investigated by means of two dimensional direct numerical simulations. The chemical reaction was irreversible and second order of the type A + B yields Products + Heat. The general governing fluid equations of a compressible reacting flow field were solved by means of high order finite difference methods. Physical effects were then determined by examining the response of the mixing layer to variation of the relevant non-dimensionalized parameters. The simulations show that increased compressibility generally results in a suppressed mixing, and consequently a reduced chemical reaction conversion rate. Reaction heat release was found to enhance mixing at the initial stages of the layer growth, but had a stabilizing effect at later times. The increased stability manifested itself in the suppression or delay of the formation of large coherent structures within the flow. Calculations were performed for a constant rate chemical kinetics model and an Arrhenius type kinetic prototype. The choice of the model was shown to have an effect on the development of the flow. The Arrhenius model caused a greater temperature increase due to reaction than the constant kinetic model. This had the same effect as increasing the exothermicity of the reaction. Localized flame quenching was also observed when the Zeldovich number was relatively large.

Steinberger, Craig J.↗

Investigation of computational and spectral analysis methods for aeroacoustic wave propagation

Most computational fluid dynamics (CFD) schemes are not adequately accurate for solving aeroacoustics problems, which have wave amplitudes several orders of magnitude smaller yet with frequencies larger than the flow field variations generating the sound. Hence, a computational aeroacoustics (CAA) algorithm should have minimal dispersion and dissipation features. A dispersion relation preserving (DRP) scheme is, therefore, applied to solve the linearized Euler equations in order to simulate the propagation of three types of waves, namely: acoustic, vorticity, and entropy waves. The scheme is derived using an optimization procedure to ensure that the numerical derivatives preserve the wave number and angular frequency of the partial differential equations being discretized. Consequently, simulated waves propagate with the correct wave speeds and exhibit their appropriate properties. A set of radiation and outflow boundary conditions, compatible with the DRP scheme and derived from the asymptotic solutions of the governing equations, are also implemented. Numerical simulations are performed to test the effectiveness of the DRP scheme and its boundary conditions. The computed solutions are shown to agree favorably with the exact solutions. The major restriction appears to be that the dispersion relations can be preserved only for waves with wave lengths longer than four or five spacings. The boundary conditions are found to be transparent to the outgoing disturbances. However, when the disturbance source is placed closer to a boundary, small acoustic reflections start appearing. CAA generates enormous amounts of temporal data which needs to be reduced to understand the physical problem being simulated. Spectral analysis is one approach that helps us in extracting information which often can not be easily interpreted in the time domain. Thus, three different methods for the spectral analysis of numerically generated aeroacoustic data are studied. First, the capabilities of two traditional methods for spectral analysis, namely, the Blackman-Tukey method and periodogram method, are compared in estimating the spectra of a simple-periodic process. The periodogram is then applied to analyze transitory-deterministic processes. Finally, these two methods are compared with a more recent method, referred as the Weighted-Overlapped-Segment-Averaging (WOSA) method, in estimating the spectra of a chaotic (random-like) process. From the demonstrative case for the spectral analyses of data generated by simple-periodic process, the periodogram method is found to give a better estimate of the steep-sloped spectra than the Blackman-Tukey method. Also, for this problem, the Hanning window is found to perform better with the periodogram method than with the Blackman-Tukey method. Finally, for the spectral analysis of data generated by the chaotic process, the periodogram method does not perform well, whereas, the WOSA and Blackman-Tukey methods give equivalently good results.

Vanel, Florence O.↗

Multi-dimensional modelling of the solar wind-LISM interaction including neutrals: A Boltzmann equation approach

The importance of interstellar neutrals in understanding and modelling the global interaction of the solar wind with the local interstellar medium is becoming increasingly apparent. Unfortunately the self-consistent inclusion of a neutral interstellar component into time-dependent, dynamical models is formidably difficult due to the extremely large mean free paths associated with the neutrals and the creation of essentially different neutral distributions from different interaction regions of the solar wind and LISM. In full generality, one has to address the problem by treating the neutrals kinetically with the appropriate extinction and creation source terms. In this paper, a limited set of simulations will be presented in which the solar wind and interstellar plasma is described as a 2D fully compressible time-dependent fluid while the interstellar neutral distribution is derived by solving the appropriate Boltzmann equation directly.

Zank, G. P.↗

Effectiveness of Flow Control for Alleviation of Twin-Tail Buffet

Effectiveness of active flow control for twin- tail buffet alleviation is investigated. Tangen- tial leading-edge blowing (TLEB) and flow suction along the vortex cores (FSVC) of the lead- ing edges of the delta wing are used to delay the vortex breakdown flow upstream of the twin tail. The combined effect of the TLEB and FSVC is also investigated. A parametric study of the effects of the spanwise position of the suction tubes and volumetric suction flow rate on the twin-tail buffet response are also investigated. The TLEB moves the path of leading-edge vortices laterally towards the twin tail, which increases the aero- dynamic damping on the tails. The FSVC effectively delays the breakdown location at high angles of attack. The computational model consists of a sharp-edged delta wing of aspect ratio one and swept-back flexible twin tail with taper ratio of 0.23. This complex multidisciplinary problem is solved sequentially using three sets of equations for the fluid flow, aeroelastic response and grid deformation, on a dynamic multi-block grid structure. The computational model is pitched at 30 deg. angle of attack. The freestream Mach number and Reynolds number are 0.3 and 1.25 million, respectively. The model is investigated for the inboard position of the twin tails, which corresponds to a separation distance between the twin tails of 33% of the wing span.

Sheta, Essam F.↗

On Wave Processes in the Solar Atmosphere

This grant was awarded by NASA/MSFC to The University of Alabama in Huntsville (UAH) to investigate the physical processes responsible for heating and wind acceleration in the solar atmosphere, and to construct theoretical, self-consistent and time-dependent solar wind models based on the momentum deposition by finite amplitude and nonlinear Alfven waves. In summary, there are three main goals of the proposed research: (1) Calculate the wave energy spectra and wave energy fluxes carried by magnetic non- magnetic waves. (2) Find out which mechanism dominates in supplying the wave energy to different parts of the solar atmosphere. (3) Use the results obtained in (1) and (2) to construct theoretical, self-consistent and time- dependent models of the solar wind. We have completed the first goal by calculating the amount of non-radiative energy generated in the solar convection zone as acoustic waves and as magnetic tube waves. To calculate the amount of wave energy carried by acoustic waves, we have used the Lighthill-Stein theory for sound generation modified by Musielak, Rosner, Stein & Ulmschneider (1994). The acoustic wave energy fluxes for stars located in different regions of the Hertzsprung-Russell (H-R) diagram have also been computed. The wave energy fluxes carried by longitudinal and transverse waves along magnetic flux tubes have been calculated by using both analytical and numerical methods. Our analytical approach is based a theory developed by Musielak, Rosner & Ulmschnelder and Musielak, Rosner, Gall & Ulmschneider, which allows computing the wave energy fluxes for linear tube waves. A numerical approach has been developed by Huang, Musielak & Ulmschneider and Ulmschneider & Musielak to compute the energy fluxes for nonlinear tube waves. Both methods have been used to calculate the wave energy fluxes for stars located in different regions of the HR diagram (Musielak, Rosner & Ulmschneider 1998; Ulmschneider, Musielak & Fawzy 1998). Having obtained the wave energy fluxes for acoustic and magnetic tube waves, we have investigated the behavior of these waves in the solar and stellar atmospheres. The results of our extensive studies have been published in many papers and presented at numerous scientific meetings. In these studies we have investigated different aspects of propagation of acoustic and magnetic waves, the efficiency of energy transfer along magnetic structures in the solar atmosphere, and behavior of Alfven waves in stgeady and expanding solar and stellar atmospheres. Recently, we have used some of these results to construct first purely theoretical, two component and time-dependent models of solar and stellar chromospheres. Finally, to address the third goal, we have constructed first fully theoretical, self-consistent and time dependent wind models based on the momentum deposition by non-linear Alfven waves. The full set of single-fluid MHD equations with the background flow has been solved by using a modified version of the ZEUS MHD code. The constructed wind models are radially symmetric with the magnetic field decreasing radially and the initial outflow is described by the standard Parker wind solution. In contrast to previous studies, no assumptions regarding wave linearity, wave damping, and wave-flow interaction are made; the models thus naturally account for the backreaction of the wind on the waves as well as for the nonlinear interaction between different types of MHD waves. The models have been used to explain the origin of fast speed streams in solar coronal holes. The obtained results clearly demonstrate that the momentum deposition by Alfven waves in the solar wind can be sufficient to explain the origin of fast stream components of the solar wind. The range of wave amplitudes required to obtain the desired results seems to be in good agreement with recent observations.

Musielak, Z. E.↗

Adaptive Suction and Blowing for Twin-Tail Buffet Control

Adaptive active flow control for twin-tail buffet alleviation is investigated. The concept behind this technique is to place control ports on the tail outer and inner surfaces with flow suction or blowing applied through these ports in order to minimize the pressure difference across the tail. The suction or blowing volume flow rate from each port is proportional to the pressure difference across the tail at this location. A parametric study of the effects of the number and location of these ports on the buffet response is carried out. The computational model consists of a sharp-edged delta wing of aspect ratio one and swept-back flexible twin tail with taper ratio of 0.23. This complex multidisciplinary problem is solved sequentially using three sets of equations for the fluid flow, aeroelastic response and grid deformation, using a dynamic multi-block grid structure. The computational model is pitched at 30 deg angle of attack. The freestream Mach number and Reynolds number are 0.3 and 1.25 million, respectively. The model is investigated for the inboard position of the twin tails, which corresponds to a separation distance between the twin tails of 33% of the wing span. Comparison of the time history and power spectral density responses of the tails for various distributions of the control ports are presented and discussed.

Kandil, Osama A.↗

Overcoming Geometry-Induced Stiffness with IMplicit-Explicit (IMEX) Runge-Kutta Algorithms on Unstructured Grids with Applications to CEM, CFD, and CAA

My goal is to develop and implement efficient, accurate, and robust Implicit-Explicit Runge-Kutta (IMEX RK) methods [9] for overcoming geometry-induced stiffness with applications to computational electromagnetics (CEM), computational fluid dynamics (CFD) and computational aeroacoustics (CAA). IMEX algorithms solve the non-stiff portions of the domain using explicit methods, and isolate and solve the more expensive stiff portions using implicit methods. Current algorithms in CEM can only simulate purely harmonic (up to lOGHz plane wave) EM scattering by fighter aircraft, which are assumed to be pure metallic shells, and cannot handle the inclusion of coatings, penetration into and radiation out of the aircraft. Efficient MEX RK methods could potentially increase current CEM capabilities by 1-2 orders of magnitude, allowing scientists and engineers to attack more challenging and realistic problems.

Kanevsky, Alex↗

An Optimized Multicolor Point-Implicit Solver for Unstructured Grid Applications on Graphics Processing Units

In the field of computational fluid dynamics, the Navier-Stokes equations are often solved using an unstructuredgrid approach to accommodate geometric complexity. Implicit solution methodologies for such spatial discretizations generally require frequent solution of large tightly-coupled systems of block-sparse linear equations. The multicolor point-implicit solver used in the current work typically requires a significant fraction of the overall application run time. In this work, an efficient implementation of the solver for graphics processing units is proposed. Several factors present unique challenges to achieving an efficient implementation in this environment. These include the variable amount of parallelism available in different kernel calls, indirect memory access patterns, low arithmetic intensity, and the requirement to support variable block sizes. In this work, the solver is reformulated to use standard sparse and dense Basic Linear Algebra Subprograms (BLAS) functions. However, numerical experiments show that the performance of the BLAS functions available in existing CUDA libraries is suboptimal for matrices representative of those encountered in actual simulations. Instead, optimized versions of these functions are developed. Depending on block size, the new implementations show performance gains of up to 7x over the existing CUDA library functions.

Zubair, Mohammad↗

Introduction and Overview of the Object-Oriented Turbomachinery Analysis Code (OTAC)

This is an overview of the Object-Oriented Turbomachinery Analysis Code (OTAC). OTAC is an object-oriented code developed to perform 2-D meridional flowthrough design and analysis of turbomachines. OTAC uses an implicit formulation of the governing equations relevant to meanline and streamline simulations to solve for the thermodynamic and velocity state of the fluid at the exit of each blade row. OTAC is written using the interpreted language capability of the Numerical Propulsion System Simulation (NPSS) code; as such, OTAC requires a version of NPSS to run. A series of models and test cases have been developed to demonstrate the capability and features of OTAC.

Turbomachinery↗

Hall2De Simulations of a Magnetic Nozzle

The two-dimensional, fluid-based Hall thruster code Hall2De is adapted to simulate a low temperature propulsive magnetic nozzle. Electron cyclotron resonance (ECR) thrusters developed at the Office National d’Etudes et de Recherches Aérospatiales (ONERA) in France and the University of Michigan are investigated. For each device, the simulated domain is defined from the nozzle throat to a location 10 thruster radii downstream. A numerical mesh is defined based on the measured magnetic field, and the governing fluid equations for singlycharge ions, electrons, and neutrals are solved. The boundary conditions for the throat of the magnetic nozzle are based on measured values of electron temperature and ion velocity. Multiple boundary conditions for the electron temperature in the downstream regions of the domain are considered. The simulated ion velocity and potential agrees well with experiment, while the simulated divergence efficiencies and thrusts are 25% lower and 45-55% higher than experiment. When a Neumann boundary condition for electron temperature is applied to the far-plume boundaries, the nozzle was seen to be purely isothermal. This suggests that the thermal conductivity along field lines in real devices is much lower than classical simulations predict, and anomalous resistivity may be needed to reproduce the temperature distribution observed in experiment. The larger beam divergence observed in our simulation may also be related to the lack of anomalous resistivity.

Jorns, Benjamin A.↗

Lumped Parameter Models of the Central Nervous System for VIIP Research

INTRODUCTION: Current long-duration missions to the International Space Station and future exploration-class missions beyond low-Earth orbit, such as to Mars and asteroids, expose astronauts to increased risk of Visual Impairment and Intracranial Pressure (VIIP) syndrome [1]. It has been hypothesized that the headward shift of cerebral spinal fluid (CSF) and blood in microgravity may cause significant elevation of intracranial pressure (ICP), which in turn induces VIIP syndrome through biomechanical pathways [1, 2]. However, there is insufficient evidence to confirm this hypothesis. In this light, we are developing lumped-parameter models of fluid transport in the central nervous system (CNS) as a means to simulate the influence of microgravity on ICP. The CNS models will also be used in concert with the lumped parameter and finite element models of the eye described in the realted IWS abstracts submitted by Nelson et al., Feola et al. and Ethier et al. METHODS: We have developed a nine compartment CNS model (Figure 1) capable of both time-dependent and steady state fluid transport simulations, based on the works of Stevens et al. [3]. The breakdown of compartments within the model includes: vascular (3), CSF (2), brain (1) and extracranial (3). The boundary pressure in the Central Arteries [A] node is prescribed using an oscillating pressure function PA(t) simulating the carotid pulsatile pressure wave as developed by Linninger et al. [4]. For each time step, pressures are integrated through time using an adaptive-timestep 4th and 5th order Runga-Kutta solver. Once pressures are found, constitutive equations are used to solve for flowrates (Q) between each compartment. In addition to fluid flow between the different compartments, compliance (C) interactions between neighboring compartments are represented. We are also developing a second CNS model based on the works of Linninger et al. [4] which takes a more granular approach to represent the interactions of the intracranial and spinal compartments with the inclusion of arteries, arterioles, capillaries, venules, veins, venous sinus, and ventricles. The flow through the arteries, veins and CSF compartments are governed by continuity, momentum and distensibility balance equations. Furthermore, unlike the Stevens et al. approach, the Monro-Kellie doctrine of constant cranial volume and the bi-phasic nature of the brain parenchyma are implemented. These features appear to be more consistent with the physiologic and anatomical behavior of the CNS, and follow a modeling philosophy similar to the lumped parameter eye model that is intended to be integrated with the CNS model. However, Linninger’s approach has never been implemented to include hydrostatic gradient and microgravity simulation capabilities. Therefore, we aim at implement this modeling approach for spaceflight simulations and assess its overall applicability to VIIP research. OBJECTIVES: We will present verification and validation test results for both models, as well as head-to-head comparison to explore their strengths and limitations with respect to mathematical implementation and physiological significance for VIIP research. In doing so, we hope to provide some guidance to the HRP research community on how to appropriately leverage lumped parameter models for space biomedical research.

Vera, J.↗