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 325 records · Page 18

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.↗

Dynamic response of a freely rotating butterfly valve in the advanced test reactor − dynamic coefficients modeling

Here, in evaluating the water hammer issue pertaining to the primary-coolant-regulating butterfly valve in the Advanced Test Reactor, the dynamic fluid body interaction (DFBI) approach was implemented in the analysis covered in Part I. Although DFBI modeling accurately and simultaneously solved the dynamic motion of the valve’s disk along with the flow field of the surrounding fluid, it shed little light on the reason behind such motion. For Part II, the reacting torque of the fluid on the disk was decomposed into representations of the dynamic coefficients in terms of stiffness, damping, and added mass. These were evaluated via simulations with steady-state static (stiffness), constant angular speed (damping), and variable angular speed (added mass) disks. Substituting the dynamic coefficients into Newton’s second law enabled the response trajectories to be obtained. Stable (by average) and unstable equilibrium positions and thrust tendencies of the valve were determined based on the stiffness coefficient (or static torque), the response amplitude was dampened or enlarged by the damping coefficient (minorly affected by added mass), and the response frequency was altered by the damping and added mass coefficients. Although the dynamic coefficient approach renders slightly different trajectories, due to the averaging effect of the torque in comparison to the DFBI method, the overall trend of the response aligns with the DFBI simulation, thus confirming the conclusion in Part I that a fix to the current butterfly valve is necessary.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗

Coupled effects of conduction in the crystal and thermo-solutal convection in a rectangular inclined enclosure

To date modeling of crystal growth of optoelectronic materials using Physical Vapor Transport has been limited to the study of the fluid phase. To achieve it, the equations of coupled heat, mass and momentum transfer in the gas have to be solved. The first objective of this study is to examine the effect of heat conduction in the crystal on the fluid flow in the neighborhood of the interface. Heat transfer boundary conditions on both interfaces were modified to take into account the additional heat flux between gas and solid. It is proved that heat conduction does not affect the fluid flow. In the presence of gravity, density gradients in the fluid phase generate convection responsible for the problem of a nonplanar growth of the interface. The second objective is to study systematically under one-g the different possible flows in order to solve this problem. Depending on the parameters, a diffusive mode and three convective modes (thermal, solutal and thermo-solutal) are observed. The competition between thermal and solutal convections leads to a mathematical condition which can be used to achieve a planar growth. It is proven that, under the physical conditions chosen, this mathematical condition cannot be thermodynamically satisfied.

Mennetrier, Christophe↗

Mass diffusion in a self-confined rotating flow

The effectiveness of fluid containment near an interior stagnation point and within a self-confined stagnation region is determined by numerically solving the species conservation equation for a bi-component mixture. The flow geometry is that of a swirling fluid stream containing a stationary eddy on the axis of rotation. The base flow is axisymmetric, and the Reynolds number is equal to 50. Schmidt numbers range from 0.1 to 10.

Torrance, K. E.↗

Progress on a Taylor weak statement finite element algorithm for high-speed aerodynamic flows

A new finite element numerical Computational Fluid Dynamics (CFD) algorithm has matured to the point of efficiently solving two-dimensional high speed real-gas compressible flow problems in generalized coordinates on modern vector computer systems. The algorithm employs a Taylor Weak Statement classical Galerkin formulation, a variably implicit Newton iteration, and a tensor matrix product factorization of the linear algebra Jacobian under a generalized coordinate transformation. Allowing for a general two-dimensional conservation law system, the algorithm has been exercised on the Euler and laminar forms of the Navier-Stokes equations. Real-gas fluid properties are admitted, and numerical results verify solution accuracy, efficiency, and stability over a range of test problem parameters.

Baker, A. J.↗

Numerical modeling of fluid and electromagnetic phenomena in an arcjet

An explicit numerical technique is used to solve the axisymmetric reduced electromagnetic field equation. The effect of an electrical arc on a viscous, axisymmetric flow is approximated using an implicit thin layer Navier-Stokes solver with additional electromagnetic source terms in conjunction with the explicit finite difference code.

Flowe, Anita C.↗

SAM: A Modern System Code for Advanced Non-LWR Safety Analysis

The System Analysis Module (SAM), developed at Argonne National Laboratory and by collaborators at other organizations, is for advanced non–light water reactor safety analysis. SAM aims to provide fast-running, modest-fidelity, whole-plant transient analysis capabilities that are essential for fast-turnaround design scoping and engineering analyses of advanced reactor concepts. To facilitate code development, SAM utilizes the MOOSE object-oriented application framework, its underlying finite element library, and linear and nonlinear solvers to leverage modern advanced software environments and numerical methods. SAM aims to solve tightly coupled physical phenomena, including fission reaction, heat transfer, fluid dynamics, and thermal-mechanical responses in advanced reactor structures, systems, and components with high accuracy and efficiency. Finally, this paper gives an overview of the SAM code development, including goals and functional requirements, physical models, current capabilities, verification and validation, software quality assurance, and examples of simulations for advanced nuclear reactor applications.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Portable detector set discloses helium leak rates

Portable helium detector measuring helium leak rates makes possible the use of the inert gas helium as a tracer. This helps solve safety and contamination problems in detecting leaks in closed fluid systems.

Anderson, G. E.↗

An analytical study of reduced-gravity liquid reorientation using a simplified marker and cell technique

A computer program called HOPI was developed to predict reorientation flow dynamics, wherein liquids move from one end of a closed, partially filled, rigid container to the other end under the influence of container acceleration. The program uses the simplified marker and cell numerical technique and, using explicit finite-differencing, solves the Navier-Stokes equations for an incompressible viscous fluid. The effects of turbulence are also simulated in the program. HOPI can consider curved as well as straight walled boundaries. Both free-surface and confined flows can be calculated. The program was used to simulate five liquid reorientation cases. Three of these cases simulated actual NASA LeRC drop tower test conditions while two cases simulated full-scale Centaur tank conditions. It was concluded that while HOPI can be used to analytically determine the fluid motion in a typical settling problem, there is a current need to optimize HOPI. This includes both reducing the computer usage time and also reducing the core storage required for a given size problem.

Betts, W. S., Jr.↗

A generalized hyperbolic marching technique for three-dimensional supersonic flow with shocks

The current paper introduces a generalized, nonorthogonal coordinate system for application to supersonic flows with such large angles between streamlines of the initial Cauchy data and the marching direction that the velocity component in the marching direction is actually subsonic. Finite-difference operators are developed to solve the three-dimensional, steady, inviscid conservation equations of fluid flow referenced to this general frame. Moreover, a new method of aligning the difference mesh to the bow shock wave is presented which eliminates both the necessity for differencing the free-stream flow properties and the spurious fluctuations that often arise in conservative variables when they are differenced across the discontinuity. Results are given for the case of flow past a conical surface.

Rizzi, A. W.↗

An explicit algorithm for a fluid approach to nonlinear optics propagation using splitting and rezoning techniques

The paper presents a three-dimensional analysis of the nonlinear light matter interaction in a hydrodynamic context. It is reported that the resulting equations are a generalization of the Navier-Stokes equations subjected to an internal potential which depends solely upon the fluid density. In addition, three numerical approaches are presented to solve the governing equations using an extension of McCormack predict-corrector scheme. These are a uniform grid, a dynamic rezoned grid, and a splitting technique. It is concluded that the use of adaptive mapping and splitting techniques with McCormack two-level predictor-corrector scheme results in an efficient and reliable code whose storage requirements are modest compared with other second order methods of equal accuracy.

Mattar, F. P.↗

Vibration of a liquid with a free surface in a spinning spherical tank

The problem of forced fluid vibrations in a partially filled spinning spherical tank is solved numerically by using the finite element method. The governing equations include Coriolis acceleration and spatially homogeneous vorticity. An exponential instability is detected in the present simulation for fill ratios below 0.5 and centrifugal acceleration to thrust ratios less than 1.7. This fictitious instability appears in the model as a result of the homogeneous vortex assumption since the free slosh equations are neutrally stable in the Liapunov sense.

El-Raheb, M.↗

Collisionless electrostatic interchange instabilities

The linear Vlasov dispersion equation for electrostatic plasma instabilities driven by gravity and weak density gradients perpendicular to a uniform magnetic field is derived and solved numerically. Two interchange instabilities emerge: the well-known fluid mode at long wavelengths and a kinetic model at wavelengths short compared with the ion gyroradius. The properties of both instabilities are studied, as well as the effects of gravity on the universal and lower-hybrid density drift instabilities. The results show that the kinetic interchange generally has a larger growth rate than the fluid interchange instability, indicating that, whenever the latter is present in a collisionless plasma, the former may also be found.

Gray, S. P.↗

Studies of pressure-velocity coupling schemes for analysis of incompressible and compressible flows

Two pressure-velocity coupling schemes, both of which solve the fully implicit discretized equations governing the flow of fluids were examined, and the capability of performing large Reynolds number, low Mach number compressible flow calculations were assessed. The semi-implicit iterative SIMPLE algorithm is extended to handle transient compressible flow calculations. This extension takes into account a strong coupling between the pressure and temperature through a correction procedure, based on the equations of state. Results obtained from the extended SIMPLE algorithm are then compared to similar results obtained from the non-iterative PISO algorithm. Both time-dependent and steady state calculations were performed using an axisymmetric 2:1 pipe expansion geometry and laminar flow conditions corresponding to Reynolds number of 1000 and Mach number of 2.0. For calculations simulating a time-dependent compression/expansion process, both schemes exhibit transient features in excellent agreement with each other, and moreover, the PISO method shows a significant computational time reduction of 60 percent over the SIMPLE scheme, regardless of the time step size or grid size employed. The effects of numerical diffusion are shown to be significant in these calculations. For steady state compressible flows, however, the SIMPLE algorithm displays increasing computational efficiency over the PISO method as the time step sizes employed to reach steady state conditions are decreased.

Burgreen, Gregory Wayne↗

Solution of the Navier-Stokes equations in three dimensions

The impact of computational fluid dynamics on the design of fluid dynamics devices is discussed. An efficient numerical procedure for solving the Navier-Stokes equations in three dimensions is discussed and flow simulations made on an inexpensive desk top work station are presented.

Maccormack, Robert W.↗

CFD applications in hypersonic flight

Design studies are underway for a variety of hypersonic flight vehicles. The National Aero-Space Plane will provide a reusable, single-stage-to-orbit capability for routine access to low earth orbit. Flight-capable satellites will dip into the atmosphere to maneuver to new orbits, while planetary probes will decelerate at their destination by atmospheric aerobraking. To supplement limited experimental capabilities in the hypersonic regime, CFD is being used to analyze the flow about these configurations. The governing equations include fluid dynamic as well as chemical species equations, which are solved with robust upwind differencing schemes. Examples of CFD applications to hypersonic vehicles suggest an important role this technology will play in the development of future aerospace systems. The computational resources needed to obtain solutions are large, but various strategies are being exploited to reduce the time required for complete vehicle simulations.

Edwards, T. A.↗

Absorption heat pump for space applications

In the first part, the performance of the Absorption Heat Pump (AHP) with water-sulfuric acid and water-magnesium chloride as two new refrigerant-absorbent fluid pairs was investigated. A model was proposed for the analysis of the new working pairs in a heat pump system, subject to different temperature lifts. Computer codes were developed to calculate the Coefficient of Performance (COP) of the system with the thermodynamic properties of the working fluids obtained from the literature. The study shows the potential of water-sulfuric acid as a satisfactory replacement for water-lithium bromide in the targeted temperature range. The performance of the AHP using water-magnesium chloride as refrigerant-absorbent pair does not compare well with those obtained using water-lithium bromide. The second part concentrated on the design and testing of a simple ElectroHydrodynamic (EHD) Pump. A theoretical design model based on continuum electromechanics was analyzed to predict the performance characteristics of the EHD pump to circulate the fluid in the absorption heat pump. A numerical method of solving the governing equations was established to predict the velocity profile, pressure - flow rate relationship and efficiency of the pump. The predicted operational characteristics of the EHD pump is comparable to that of turbomachinery hardware; however, the overall efficiency of the electromagnetic pump is much lower. An experimental investigation to verify the numerical results was conducted. The pressure - flow rate performance characteristics and overall efficiency of the pump obtained experimentally agree well with the theoretical model.

Nguyen, Tuan↗

Aerodynamic simulation on massively parallel systems

This paper briefly addresses the computational requirements for the analysis of complete configurations of aircraft and spacecraft currently under design to be used for advanced transportation in commercial applications as well as in space flight. The discussion clearly shows that massively parallel systems are the only alternative which is both cost effective and on the other hand can provide the necessary TeraFlops, needed to satisfy the narrow design margins of modern vehicles. It is assumed that the solution of the governing physical equations, i.e., the Navier-Stokes equations which may be complemented by chemistry and turbulence models, is done on multiblock grids. This technique is situated between the fully structured approach of classical boundary fitted grids and the fully unstructured tetrahedra grids. A fully structured grid best represents the flow physics, while the unstructured grid gives best geometrical flexibility. The multiblock grid employed is structured within a block, but completely unstructured on the block level. While a completely unstructured grid is not straightforward to parallelize, the above mentioned multiblock grid is inherently parallel, in particular for multiple instruction multiple datastream (MIMD) machines. In this paper guidelines are provided for setting up or modifying an existing sequential code so that a direct parallelization on a massively parallel system is possible. Results are presented for three parallel systems, namely the Intel hypercube, the Ncube hypercube, and the FPS 500 system. Some preliminary results for an 8K CM2 machine will also be mentioned. The code run is the two dimensional grid generation module of Grid, which is a general two dimensional and three dimensional grid generation code for complex geometries. A system of nonlinear Poisson equations is solved. This code is also a good testcase for complex fluid dynamics codes, since the same datastructures are used. All systems provided good speedups, but message passing MIMD systems seem to be best suited for large miltiblock applications.

Haeuser, Jochem↗