Search NASA⌕ Search

SEARCH · Search NASA

Results for “FLUID BOUNDARY”

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 397 records · Page 22

Working Fluid Trade Study for a Two-Phase Mechanically Pumped Loop Thermal Control System

NASA Jet Propulsion Laboratory is investigating a two-phase mechanically pumped fluid loop (MPFL) technology for spacecraft thermal control for future planetary space science mission concepts. The two-phase technology combines the potential of single-phase MPFL and the passive two-phase heat pipe systems. Previous two-phase MPFL studies4,5 at JPL have identified an evaporator system with a separated flow architecture as providing a robust light mass thermal control system with very high spatial and temporal thermal stability required by the science instruments on these missions. These studies have further identified the thermophysical properties of the working fluid to be key to the performance of a two-phase fluid loop. This paper describes the methodology used in the selection of the working fluids for optimizing the performance of the two-phase MPFL. A high level model was developed that includes the constraints and boundaries driven by the system components. The performance of fluids from the REFPROP database were investigated and rated for this system. The working fluid attributes such as cost, hazardous properties, and heritage were taken into account in this selection. A typical spacecraft dissipating 1000 W and a fluid loop consisting of an evaporator, accumulator, radiator, and associated tubing components was used as an example in this study.

Daimaru, Takuro↗

Surface tension and buoyancy-driven flow in a non-isothermal liquid bridge

The Navier-Stokes-Boussinesq equations governing the transport of momentum, mass and heat in a nonisothermal liquid bridge with a temperature-dependent surface tension are solved using a vorticity-stream-function formulation together with a nonorthogonal coordinate transformation. The equations are discretized using a pseudo-unsteady semi-implicit finite difference scheme and are solved by the ADI method. A Picard-type iteration is adopted which consists of inner and outer iterative processes. The outer iteration is used to update the shape of the free surface. Two schemes have been used for the outer iteration; both use the force balance normal to the free surface as the distinguished boundary condition. The first scheme involves successive approximation by the direct solution of the distinguished boundary condition. The second scheme uses the artificial force imbalance between the fluid pressure, viscous and capillary forces at the free surface which arises when the boundary condition for force balance normal to the surface is not satisfied. This artificial imbalance is then used to change the surface shape until the distinguished boundary condition is satisfied. These schemes have been used to examine a variety of model liquid bridge situations including purely thermocapillary-driven flow situations and mixed thermocapillary- and bouyancy-driven flow.

Zhang, Yiqiang↗

Fluid physics phenomena of resistojet thrusters

This final report includes a list of publications and part of an M.S. thesis titled 'Analyses in Theoretical and Experimental Fluid Flow', by Tony G. Howell. The thesis discusses analyses of momentum and heat transfer occurring in a laminar boundary layer of a non-Newtonian power-law fluid, and experiments completed in a simulated space thruster's plume for prediction comparison.

DeWitt, Kenneth J.↗

On some problems in a theory of thermally and mechanically interacting continuous media

Using a linearized theory of thermally and mechanically interacting mixture of linear elastic solid and viscous fluid, we derive a fundamental relation in an integral form called a reciprocity relation. This reciprocity relation relates the solution of one initial-boundary value problem with a given set of initial and boundary data to the solution of a second initial-boundary value problem corresponding to a different initial and boundary data for a given interacting mixture. From this general integral relation, reciprocity relations are derived for a heat-conducting linear elastic solid, and for a heat-conducting viscous fluid. An initial-boundary value problem is posed and solved for the mixture of linear elastic solid and viscous fluid. With the aid of the Laplace transform and the contour integration, a real integral representation for the displacement of the solid constituent is obtained as one of the principal results of the analysis.

Lee, Y. M.↗

A Low Cost Thermal Control Method for Testing in a Mars Environment

Due to the unique thermal vacuum testing requirements for a Mars Rover instrument, NASA Goddard developed a low cost, high fidelity thermal control system utilizing Thermal Electric Coolers (TECs) combined with a heat rejection fluid loop to actively control 8 independent payload thermal boundary zones in a simulated Mars pressure vacuum chamber with a Carbon Dioxide atmosphere. These zones could control instrument components to a specific temperature as a function of time to simulate exact temporal flight boundary predictions.The Mars Organic Molecule Analyzer (MOMA) instrument is a dual source (pyrolysis gas chromatograph and laser desorption) mass spectrometer (MS) based package that detects and characterizes organic molecules, as part of ESA's 2020 ExoMars Rover mission to seek the signs of life on Mars.Due to the unique thermal vacuum testing requirements for a Mars Rover instrument, NASA Goddard developed a low cost, high fidelity thermal control system utilizing Thermal Electric Coolers (TECs) combined with a heat rejection fluid loop to actively control 8 independent payload thermal boundary zones in a simulated Mars pressure vacuum chamber with a Carbon Dioxide atmosphere. These zones could control instrument components to a specific temperature as a function of time to simulate exact temporal flight boundary predictions.The Mars Organic Molecule Analyzer (MOMA) instrument is a dual source (pyrolysis gas chromatograph and laser desorption) mass spectrometer (MS) based package that detects and characterizes organic molecules, as part of ESA's 2020 ExoMars Rover mission to seek the signs of life on Mars.

MOMA↗

PIV Measurements of Shock-Wave/Boundary-Layer Interactions at Mach 2.5 in a Circular Test Section

The particle image velocimetry (PIV) measurement technique has been applied to make detailed measurements in an impinging shock-wave/boundary-layer interaction (SWBLI) flowfield in a test section of circular cross-section. By virtue of this test section geometry, the flowfields studied are free from effects of the complex coupling of sidewall and corner interaction into the primary impinging shock interaction typical of a rectangular test section. The cone-cylinder shock generator hardware was arranged in both coaxial alignment with the test section to generate axisymmetric SWBLIs and in offset positions to generate swept SWBLIs. The angle of the conical section was also varied to alter the strength of the impinging shock wave. The application of the PIV technique in this test section required meticulous design to provide adequate optical access without disturbing the SWBLI flow. Details of the stereoscopic PIV apparatus and the velocity measurements obtained for a Mach 2.5 freestream are presented. Key features of the mean velocity and turbulent stress fields are highlighted. This data set should enable validation of computational fluid dynamics results for these flowfields.

Shockwave Boundary Layer Interaction↗

PIV Measurements of Shock-Wave/Boundary-Layer Interactions at Mach 2.5 in a Circular Test Section

The particle image velocimetry (PIV) measurement technique has been applied to make detailed measurements in an impinging shock-wave/boundary-layer interaction (SWBLI) flowfield in a test section of circular cross-section. By virtue of this test section geometry, the flowfields studied are free from effects of the complex coupling of sidewall and corner interaction into the primary impinging shock interaction typical of a rectangular test section. The cone-cylinder shock generator hardware was arranged in both coaxial alignment with the test section to generate axisymmetric SWBLIs and in offset positions to generate swept SWBLIs. The angle of the conical section was also varied to alter the strength of the impinging shock wave. The application of the PIV technique in this test section required meticulous design to provide adequate optical access without disturbing the SWBLI flow. Details of the stereoscopic PIV apparatus and the velocity measurements obtained for a Mach 2.5 freestream are presented. Key features of the mean velocity and turbulent stress fields are highlighted. This data set should enable validation of computational fluid dynamics results for these flowfields.

Shockwave Boundary Layer Interaction↗

Convection in Slab and Spheroidal Geometries

Three-dimensional numerical simulations of compressible turbulent thermally driven convection, in both slab and spheroidal geometries, are reviewed and analyzed in terms of velocity spectra and mixing-length theory. The same ideal gas model is used in both geometries, and resulting flows are compared. The piecewise-parabolic method (PPM), with either thermal conductivity or photospheric boundary conditions, is used to solve the fluid equations of motion. Fluid motions in both geometries exhibit a Kolmogorov-like k(sup -5/3) range in their velocity spectra. The longest wavelength modes are energetically dominant in both geometries, typically leading to one convection cell dominating the flow. In spheroidal geometry, a dipolar flow dominates the largest scale convective motions. Downflows are intensely turbulent and up drafts are relatively laminar in both geometries. In slab geometry, correlations between temperature and velocity fluctuations, which lead to the enthalpy flux, are fairly independent of depth. In spheroidal geometry this same correlation increases linearly with radius over the inner 70 percent by radius, in which the local pressure scale heights are a sizable fraction of the radius. The effects from the impenetrable boundary conditions in the slab geometry models are confused with the effects from non-local convection. In spheroidal geometry nonlocal effects, due to coherent plumes, are seen as far as several pressure scale heights from the lower boundary and are clearly distinguishable from boundary effects.

Porter, David H.↗

Flowfield-Dependent Mixed Explicit-Implicit (FDMEL) Algorithm for Computational Fluid Dynamics

Despite significant achievements in computational fluid dynamics, there still remain many fluid flow phenomena not well understood. For example, the prediction of temperature distributions is inaccurate when temperature gradients are high, particularly in shock wave turbulent boundary layer interactions close to the wall. Complexities of fluid flow phenomena include transition to turbulence, relaminarization separated flows, transition between viscous and inviscid incompressible and compressible flows, among others, in all speed regimes. The purpose of this paper is to introduce a new approach, called the Flowfield-Dependent Mixed Explicit-Implicit (FDMEI) method, in an attempt to resolve these difficult issues in Computational Fluid Dynamics (CFD). In this process, a total of six implicitness parameters characteristic of the current flowfield are introduced. They are calculated from the current flowfield or changes of Mach numbers, Reynolds numbers, Peclet numbers, and Damkoehler numbers (if reacting) at each nodal point and time step. This implies that every nodal point or element is provided with different or unique numerical scheme according to their current flowfield situations, whether compressible, incompressible, viscous, inviscid, laminar, turbulent, reacting, or nonreacting. In this procedure, discontinuities or fluctuations of an variables between adjacent nodal points are determined accurately. If these implicitness parameters are fixed to certain numbers instead of being calculated from the flowfield information, then practically all currently available schemes of finite differences or finite elements arise as special cases. Some benchmark problems to be presented in this paper will show the validity, accuracy, and efficiency of the proposed methodology.

Garcia, S. M.↗

An Application of CFD to Guide Forced Boundary-Layer Transition for Low-Speed Tests of a Hybrid Wing-Body Configuration

A hybrid transition trip-dot sizing and placement test technique was developed in support of recent experimental research on a hybrid wing-body configuration under study for the NASA Environmentally Responsible Aviation project. The approach combines traditional methods with Computational Fluid Dynamics. The application had three-dimensional boundary layers that were simulated with either fully turbulent or transitional flow models using established Reynolds-Averaged Navier-Stokes methods. Trip strip effectiveness was verified experimentally using infrared thermography during a low-speed wind tunnel test. Although the work was performed on one specific configuration, the process was based on fundamental flow physics and could be applicable to other configurations.

Luckring, James M.↗

Convection equation modeling: A non-iterative direct matrix solution algorithm for use with SINDA

The determination of the boundary conditions for a component-level analysis, applying discrete finite element and finite difference modeling techniques often requires an analysis of complex coupled phenomenon that cannot be described algebraically. For example, an analysis of the temperature field of a coldplate surface with an integral fluid loop requires a solution to the parabolic heat equation and also requires the boundary conditions that describe the local fluid temperature. However, the local fluid temperature is described by a convection equation that can only be solved with the knowledge of the locally-coupled coldplate temperatures. Generally speaking, it is not computationally efficient, and sometimes, not even possible to perform a direct, coupled phenomenon analysis of the component-level and boundary condition models within a single analysis code. An alternative is to perform a disjoint analysis, but transmit the necessary information between models during the simulation to provide an indirect coupling. For this approach to be effective, the component-level model retains full detail while the boundary condition model is simplified to provide a fast, first-order prediction of the phenomenon in question. Specifically for the present study, the coldplate structure is analyzed with a discrete, numerical model (SINDA) while the fluid loop convection equation is analyzed with a discrete, analytical model (direct matrix solution). This indirect coupling allows a satisfactory prediction of the boundary condition, while not subjugating the overall computational efficiency of the component-level analysis. In the present study a discussion of the complete analysis of the derivation and direct matrix solution algorithm of the convection equation is presented. Discretization is analyzed and discussed to extend of solution accuracy, stability and computation speed. Case studies considering a pulsed and harmonic inlet disturbance to the fluid loop are analyzed to assist in the discussion of numerical dissipation and accuracy. In addition, the issues of code melding or integration with standard class solvers such as SINDA are discussed to advise the user of the potential problems to be encountered.

Schrage, Dean S.↗

Adaptive Immersed Boundary Simulations for the Launch Environment

A high-fidelity computational fluid dynamics simulation of a next generation heavy lift space vehicle during launch is presented. The purpose of the simulation is to evaluate the acoustic overpressures during ignition to permit re-design of the launch site to safely handle heavy lift vehicles. The simulation is performed using the Launch, Ascent, and Vehicle Aerodynamics (LAVA) code, an immersed boundary block-structured Cartesian adaptive mesh refinement based solver. A verification and validation study of LAVA in the launch environment context is also performed, comparing to flight data and previous simulations of a Space Shuttle launch

Barad, Michael F.↗

Analytic corrections to CFD heating predictions accounting for changes in surface catalysis

Integral boundary-layer solution techniques applicable to the problem of determining aerodynamic heating rates of hypersonic vehicles in the vicinity of stagnation points and windward centerlines are briefly summarized. A new approach for combining the insight afforded by integral boundary-layer analysis with comprehensive (but time intensive) computational fluid dynamic (CFD) flowfield solutions of the thin-layer Navier-Stokes equations is described. The approach extracts CFD derived quantities at the wall and at the boundary layer edge for inclusion in a post-processing boundary-layer analysis. It allows a designer at a workstation to address two questions, given a single CFD solution. (1) How much does the heating change for a thermal protection system with different catalytic properties than was used in the original CFD solution? (2) How does the heating change at the interface of two different TPS materials with an abrupt change in catalytic efficiency? The answer to the second question is particularly important, because abrupt changes from low to high catalytic efficiency can lead to localized increase in heating which exceeds the usually conservative estimate provided by a fully catalytic wall assumption.

Gnoffo, Peter A.↗

Analytic Corrections to CFD Heating Predictions Accounting for Changes in Surface Catalysis

A new approach for combining the insight afforded by integral boundary-layer analysis with comprehensive (but time intensive) computational fluid dynamic (CFD) flowfield solutions of the thin-layer Navier-Stokes equations is described. The approach extracts CFD derived quantities at the wall and at the boundary layer edge for inclusion in a post-processing boundary-layer analysis. It allows a designer at a work-station to address two questions, given a single CFD solution. (1) How much does the heating change for a thermal protection system (TPS) with different catalytic properties than was used in the original CFD solution? (2) How does the heating change at the interface of two different TPS materials with an abrupt change in catalytic efficiency? The answer to the second question is particularly important, because abrupt changes from low to high catalytic efficiency can lead to localized increase in heating which exceeds the usually conservative estimate provided by a fully catalytic wall assumption. Capabilities of this approach for application to Reusable Launch Vehicle (RLV) design are demonstrated. If the definition of surface catalysis is uncertain early in the design process, results show that fully catalytic wall boundary conditions provide the best baseline for CFD design points.

Gnoffo, Peter A.↗

Two-Dimensional Dynamic Simulation of a Continuous Foil Bearing

In this paper, the two dimensional(radial and circumferential) transient Navier-Stokes equations are used to solve the hydrodynamic problem in conjunction with the time dependent motion of the journal, and the deformable, spring supported foil. The elastic deformation of the foil and its supports are simulated by a finite element model. The time-dependent Navier-Stokes formulation is used to solve for the interaction between the fluid lubricant, the motion of the journal and the deformable foil boundary. The steady state, the quasi-transient and the full transient dynamic simulation of the foil-fluid journal interaction are examined on a comparative basis. For the steady state simulation, the fluid lubricant pressures are evaluated for a particular journal position, by means of an iterative scheme until convergence is achieved in both the fluid pressures and the corresponding foil deformation. For the quasi-transient case, the transient motion of the journal is calculated using a numerical integration scheme for the velocity and displacement of the journal. The deformation of the foil is evaluated through numerical iteration in feedback mode with the fluid film pressure generated by the journal motion until convergence at every time step is achieved. For the full transient simulation, a parallel real-time integration scheme is used to evaluate simultaneously the new journal position and the new deformed shape of the foil at each time step. The pressure of the fluid lubricant is iterated jointly with the corresponding journal position and the deformed foil geometry until convergence is achieved. A variable time-stepping Newmark-Beta integration procedure is used to evaluate the transient dynamics at each time step of the bearing.

Braun, M. Jack↗

An Application of CFD to Guide Forced Boundary-Layer Transition at Subsonic and Transonic Speeds for a Lambda-Wing Configuration

A hybrid transition trip-dot sizing and placement test technique was extended to support recent experimental research on a moderately-swept hybrid wing-body Uninhabited Combat Air Vehicle configuration under study in a North Atlantic Treaty Organization Applied Vehicle Technology task group. The approach combines traditional test techniques with Computational Fluid Dynamics. The application had three-dimensional boundary layers that were simulated with either fully turbulent or mixed laminar-turbulent flow models using established Reynolds-averaged Navier-Stokes methods. Forced transition effectiveness was verified experimentally using sublimating chemicals during low-speed and transonic wind tunnel tests. Assessments were performed between tripped low Reynolds number data and un-tripped high Reynolds number data. The hybrid test technique is based on fundamental flow physics and as such could be applicable to other configurations.

UCAV↗