Search NASA⌕ Search

SEARCH · Search NASA

Results for “two fluid model”

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 451 records · Page 25
All sourcesNASA NTRSDOE OSTI
NASA NTRS2010

Computational Analysis of Arc-Jet Wedge Tests Including Ablation and Shape Change

Coupled fluid-material response analyses of arc-jet wedge ablation tests conducted in a NASA Ames arc-jet facility are considered. These tests were conducted using blunt wedge models placed in a free jet downstream of the 6-inch diameter conical nozzle in the Ames 60-MW Interaction Heating Facility. The fluid analysis includes computational Navier-Stokes simulations of the nonequilibrium flowfield in the facility nozzle and test box as well as the flowfield over the models. The material response analysis includes simulation of two-dimensional surface ablation and internal heat conduction, thermal decomposition, and pyrolysis gas flow. For ablating test articles undergoing shape change, the material response and fluid analyses are coupled in order to calculate the time dependent surface heating and pressure distributions that result from shape change. The ablating material used in these arc-jet tests was Phenolic Impregnated Carbon Ablator. Effects of the test article shape change on fluid and material response simulations are demonstrated, and computational predictions of surface recession, shape change, and in-depth temperatures are compared with the experimental measurements.

Goekcen, Tahir↗
NASA NTRS2010

Hybrid Model of Inhomogeneous Solar Wind Plasma Heating by Alfven Wave Spectrum: Parametric Studies

Observations of the solar wind plasma at 0.3 AU and beyond show that a turbulent spectrum of magnetic fluctuations is present. Remote sensing observations of the corona indicate that heavy ions are hotter than protons and their temperature is anisotropic (T(sub perpindicular / T(sub parallel) >> 1). We study the heating and the acceleration of multi-ion plasma in the solar wind by a turbulent spectrum of Alfvenic fluctuations using a 2-D hybrid numerical model. In the hybrid model the protons and heavy ions are treated kinetically as particles, while the electrons are included as neutralizing background fluid. This is the first two-dimensional hybrid parametric study of the solar wind plasma that includes an input turbulent wave spectrum guided by observation with inhomogeneous background density. We also investigate the effects of He++ ion beams in the inhomogeneous background plasma density on the heating of the solar wind plasma. The 2-D hybrid model treats parallel and oblique waves, together with cross-field inhomogeneity, self-consistently. We investigate the parametric dependence of the perpendicular heating, and the temperature anisotropy in the H+-He++ solar wind plasma. It was found that the scaling of the magnetic fluctuations power spectrum steepens in the higher-density regions, and the heating is channeled to these regions from the surrounding lower-density plasma due to wave refraction. The model parameters are applicable to the expected solar wind conditions at about 10 solar radii.

Ofman, L.↗
NASA NTRS2011

DeCoM Validation

Deepak Condenser Model (DeCoM) was validated in a two-step process. The first step validated DeCoM/TTH/FloCAD against a simplified 1D (1D fluid flow and 1D heat flow) scenario. Through this simplification the fluid calculations were simple enough that the results were comparable to hand calculations. Once a method was compared against 1D case, the next step was the verification against a 2D test case. 2D is only for heat flow, while the fluid flow is still 1D. The test data was obtained from the project Geoscience Laser Altimeter System (GLAS) which performed a Looped Heat Pipe (LHP) test and contained necessary condenser information; such as a test point and condenser configuration. From the GLAS LHP test report, the construction of the LHP condenser and correlation was roughly possible. The assumptions and limitations are discussed, within this paper, of the 2D correlation. The success of DeCoM is illustrated by the results after correlation to the test data. Results show the differences between the model and the test setup, as well as explaining the probable cause for the differences; another set of results were produced after taking into account these differences, such as parasitic heat leak. Substantial effort has been made in listing the possibilities to further develop DECOM, thus to make it more reliable, user-friendly, and better integrated with generic SINDA models. The comparison of DeCoM to FloCAD and TTH LHP should not be taken as a definite answer to choosing DeCoM as a condenser modeler, rather as an alternate to obtaining preliminary results without the complexity of TTH LHP and license of FloCAD.

Patel, Deepak↗
NASA NTRS1998

A Method for Generating Reduced-Order Linear Models of Multidimensional Supersonic Inlets

Simulation of high speed propulsion systems may be divided into two categories, nonlinear and linear. The nonlinear simulations are usually based on multidimensional computational fluid dynamics (CFD) methodologies and tend to provide high resolution results that show the fine detail of the flow. Consequently, these simulations are large, numerically intensive, and run much slower than real-time. ne linear simulations are usually based on large lumping techniques that are linearized about a steady-state operating condition. These simplistic models often run at or near real-time but do not always capture the detailed dynamics of the plant. Under a grant sponsored by the NASA Lewis Research Center, Cleveland, Ohio, a new method has been developed that can be used to generate improved linear models for control design from multidimensional steady-state CFD results. This CFD-based linear modeling technique provides a small perturbation model that can be used for control applications and real-time simulations. It is important to note the utility of the modeling procedure; all that is needed to obtain a linear model of the propulsion system is the geometry and steady-state operating conditions from a multidimensional CFD simulation or experiment. This research represents a beginning step in establishing a bridge between the controls discipline and the CFD discipline so that the control engineer is able to effectively use multidimensional CFD results in control system design and analysis.

Chicatelli, Amy↗
NASA NTRS1996

Modeling Orifice Pulse Tube Coolers

We have developed a calculational model that treats all the components of an orifice pulse tube cooler. We base our analysis on 1-dimensional thermodynamic equations for the regenerator and we assume that all mass flows, pressure oscillations and temperature oscillations are small and sinusoidal. Non-linear pressure drop effects are included in the regenerator to account for finite pressure amplitude effects. The resulting mass flows and pressures are matched at the boundaries with the other components of the cooler: compressor, aftercooler, cold heat exchanger, pulse tube, hot heat exchanger, orifice and reservoir. The results of the calculation are oscillating pressures, mass flows and enthalpy flows in the main components of the cooler. By comparing with the calculations of other available models, we show that our model is very similar to REGEN 3 from NIST and DeltaE from Los Alamos National Lab. Our model is much easier to use than other available models because of its simple graphical interface and the fact that no guesses are required for the operating pressures or mass flows. In addition, the model only requires a few minutes of running time allowing many parameters to be optimized in a reasonable time. A version of the model is available for use over the World Wide Web at http://irtek.arc.nasa.gov. Future enhancements include adding a bypass orifice and including second order terms in steady mass streaming and steady heat transfer. A two-dimensional anelastic approximation of the fluid equations will be used as the basis for the latter analysis. Preliminary results are given in dimensionless numbers appropriate for oscillating compressible flows. The model shows how transverse heat transfer reduces enthalpy flow, particularly for small pulse tubes. The model also clearly shows mass recirculation in the open tube on the order of the tube length. They result from the higher order Reynolds stresses. An interesting result of the linearized approach is that the steady mass streaming does not affect the enthalpy flow at second order. The major effect of recirculating mass streaming is to increase transverse temperature gradients, which leads to higher entropy production and reduced efficiency.

Roach, Kittel P.↗
NASA NTRS1977

Finite element analysis of low speed viscous and inviscid aerodynamic flows

A weak interaction solution algorithm was established for aerodynamic flow about an isolated airfoil. Finite element numerical methodology was applied to solution of each of differential equations governing potential flow, and viscous and turbulent boundary layer and wake flow downstream of the sharp trailing edge. The algorithm accounts for computed viscous displacement effects on the potential flow. Closure for turbulence was accomplished using both first and second order models. The COMOC finite element fluid mechanics computer program was modified to solve the identified equation systems for two dimensional flows. A numerical program was completed to determine factors affecting solution accuracy, convergence and stability for the combined potential, boundary layer, and parabolic Navier-Stokes equation systems. Good accuracy and convergence are demonstrated. Each solution is obtained within the identical finite element framework of COMOC.

Baker, A. J.↗
NASA NTRS1985

Numerical solution of the full-potential equation for rotors and oblique wings using a new wake model

A three-dimensional, full-potential, quasi-steady code TFAR1 is proposed for calculating the transonic flow past a lifting helicopter rotor blade and oblique wing. The TFAR1 uses a two-dimensional nonlinear wake-model that allows a jump in velocity potential to propagate with the local fluid flow in the wake. Rotor calculations were made for a single blade at an advance ratio of 0.3, a rotational tip Mach number of 0.7, and at 0-degree incidence. A 1/7-scale model of the Cobra Operational Load Survey (OLS) rotor blade is calculated, and the pressure distributions are compared to the measurements for azimuth angles 0, 30, 60, 90, 120, and 150 degrees at the 95 percent spanwise station of the OLS blade. Furthermore, an oblique wing with Korn airfoil was calculated at the high transonic free-stream Mach number of 0.9791, zero incidence, and yaw angle of 40 degrees. The TFAR1, coupled with a helicopter performance code CAMRAD (Johnson, 1981), provides a full-potential code for calculating the entire flow field for a multiple-bladed rotor in transonic lifting forward flight.

Chang, I.-C.↗
NASA NTRS1989

The transient behavior of weldpools with a deformed free surface

Transient heat transfer, fluid flow and phase change in a weldpool with a deformed free surface are modeled using a two-dimensional coordinate transfornmation. The theoretical approach makes it possible to accurately model high Marangoni number flows even when the free surface is not horizontal. The method is used to examine the role of surface deformation and surface tension forces in controlling weld penetration during cathode spot welding of steel and aluminum. The calculated flow fields, temperature distributions and weldpool shapes are consistent with experimental findings, and are in reasonable agreement with previously published modeling efforts.

Thompson, M. E.↗
NASA NTRS1990

Unstructured mesh methods for CFD

Mesh generation methods for Computational Fluid Dynamics (CFD) are outlined. Geometric modeling is discussed. An advancing front method is described. Flow past a two engine Falcon aeroplane is studied. An algorithm and associated data structure called the alternating digital tree, which efficiently solves the geometric searching problem is described. The computation of an initial approximation to the steady state solution of a given poblem is described. Mesh generation for transient flows is described.

Peraire, J.↗
NASA NTRS2011

Six-Tube Freezable Radiator Testing and Model Correlation

Freezable radiators offer an attractive solution to the issue of thermal control system scalability. As thermal environments change, a freezable radiator will effectively scale the total heat rejection it is capable of as a function of the thermal environment and flow rate through the radiator. Scalable thermal control systems are a critical technology for spacecraft that will endure missions with widely varying thermal requirements. These changing requirements are a result of the spacecraft s surroundings and because of different thermal loads rejected during different mission phases. However, freezing and thawing (recovering) a freezable radiator is a process that has historically proven very difficult to predict through modeling, resulting in highly inaccurate predictions of recovery time. These predictions are a critical step in gaining the capability to quickly design and produce optimized freezable radiators for a range of mission requirements. This paper builds upon previous efforts made to correlate a Thermal Desktop(TradeMark) model with empirical testing data from two test articles, with additional model modifications and empirical data from a sub-component radiator for a full scale design. Two working fluids were tested, namely MultiTherm WB-58 and a 50-50 mixture of DI water and Amsoil ANT.

Lillibridge, Sean↗
NASA NTRS2012

Six-Tube Freezable Radiator Testing and Model Correlation

Freezable Radiators offer an attractive solution to the issue of thermal control system scalability. As thermal environments change, a freezable radiator will effectively scale the total heat rejection it is capable of as a function of the thermal environment and flow rate through the radiator. Scalable thermal control systems are a critical technology for spacecraft that will endure missions with widely varying thermal requirements. These changing requirements are a result of the spacecraft?s surroundings and because of different thermal loads rejected during different mission phases. However, freezing and thawing (recov ering) a freezable radiator is a process that has historically proven very difficult to predict through modeling, resulting in highly inaccurate predictions of recovery time. These predictions are a critical step in gaining the capability to quickly design and produce optimized freezable radiators for a range of mission requirements. This paper builds upon previous efforts made to correlate a Thermal Desktop(TM) model with empirical testing data from two test articles, with additional model modifications and empirical data from a sub-component radiator for a full scale design. Two working fluids were tested: MultiTherm WB-58 and a 50-50 mixture of DI water and Amsoil ANT.

Lilibridge, Sean T.↗
NASA NTRS2024

Tail-Cone Thruster Fan Rotor Flutter Analysis

The aeroelastic stability of a tail cone thruster fan was analyzed for three structural modes and various nodal diameter patterns. Computational fluid dynamics modeling with blade vibrations was used in the analyses performed at design (peak-efficiency) and two near-stall conditions at the design rotational speed. Aerodynamic damping was used to assess flutter since this integrally bladed rotor design will have minimal structural damping. Aerodynamic damping was found to be low at the design condition. For the first mode, aerodynamic damping was seen to decrease towards stall, vanishing for a zero nodal diameter pattern. Further decrease in mass flow led to negative aerodynamic damping (flutter) for the first and second nodal diameters in the near stall condition. The flowfield at this lowest mass flow rate exhibited some additional unsteadiness at higher frequency. The sensitivity of results to time step and amplitudes of vibration was examined and found to be very satisfactory. The largest uncertainty in the results comes from the observed sensitivity of the aerodynamic damping to details of the variation of grid deformations in the blade passage mesh from the blade passage to the boundaries. For the parametric variations analyzed, the calculated aerodynamic damping is lower if the grid deformations near the blade surface are higher. Further work is needed to understand the cause of this sensitivity and find possible solutions to reduce the effects on flutter predictions.

Fan↗
NASA NTRS1993

Cavitation Modeling in Euler and Navier-Stokes Codes

Many previous researchers have modeled sheet cavitation by means of a constant pressure solution in the cavity region coupled with a velocity potential formulation for the outer flow. The present paper discusses the issues involved in extending these cavitation models to Euler or Navier-Stokes codes. The approach taken is to start from a velocity potential model to ensure our results are compatible with those of previous researchers and available experimental data, and then to implement this model in both Euler and Navier-Stokes codes. The model is then augmented in the Navier-Stokes code by the inclusion of the energy equation which allows the effect of subcooling in the vicinity of the cavity interface to be modeled to take into account the experimentally observed reduction in cavity pressures that occurs in cryogenic fluids such as liquid hydrogen. Although our goal is to assess the practicality of implementing these cavitation models in existing three-dimensional, turbomachinery codes, the emphasis in the present paper will center on two-dimensional computations, most specifically isolated airfoils and cascades. Comparisons between velocity potential, Euler and Navier-Stokes implementations indicate they all produce consistent predictions. Comparisons with experimental results also indicate that the predictions are qualitatively correct and give a reasonable first estimate of sheet cavitation effects in both cryogenic and non-cryogenic fluids. The impact on CPU time and the code modifications required suggests that these models are appropriate for incorporation in current generation turbomachinery codes.

Deshpande, Manish↗
NASA NTRS1987

Principles for problem aggregation and assignment in medium scale multiprocessors

One of the most important issues in parallel processing is the mapping of workload to processors. This paper considers a large class of problems having a high degree of potential fine grained parallelism, and execution requirements that are either not predictable, or are too costly to predict. The main issues in mapping such a problem onto medium scale multiprocessors are those of aggregation and assignment. We study a method of parameterized aggregation that makes few assumptions about the workload. The mapping of aggregate units of work onto processors is uniform, and exploits locality of workload intensity to balance the unknown workload. In general, a finer aggregate granularity leads to a better balance at the price of increased communication/synchronization costs; the aggregation parameters can be adjusted to find a reasonable granularity. The effectiveness of this scheme is demonstrated on three model problems: an adaptive one-dimensional fluid dynamics problem with message passing, a sparse triangular linear system solver on both a shared memory and a message-passing machine, and a two-dimensional time-driven battlefield simulation employing message passing. Using the model problems, the tradeoffs are studied between balanced workload and the communication/synchronization costs. Finally, an analytical model is used to explain why the method balances workload and minimizes the variance in system behavior.

Nicol, David M.↗
NASA NTRS2004

Transient Two-Dimensional Analysis of Side Load in Liquid Rocket Engine Nozzles

Two-dimensional planar and axisymmetric numerical investigations on the nozzle start-up side load physics were performed. The objective of this study is to develop a computational methodology to identify nozzle side load physics using simplified two-dimensional geometries, in order to come up with a computational strategy to eventually predict the three-dimensional side loads. The computational methodology is based on a multidimensional, finite-volume, viscous, chemically reacting, unstructured-grid, and pressure-based computational fluid dynamics formulation, and a transient inlet condition based on an engine system modeling. The side load physics captured in the low aspect-ratio, two-dimensional planar nozzle include the Coanda effect, afterburning wave, and the associated lip free-shock oscillation. Results of parametric studies indicate that equivalence ratio, combustion and ramp rate affect the side load physics. The side load physics inferred in the high aspect-ratio, axisymmetric nozzle study include the afterburning wave; transition from free-shock to restricted-shock separation, reverting back to free-shock separation, and transforming to restricted-shock separation again; and lip restricted-shock oscillation. The Mach disk loci and wall pressure history studies reconfirm that combustion and the associated thermodynamic properties affect the formation and duration of the asymmetric flow.

Wang, Ten-See↗
NASA NTRS2019

Pilot Sensitivity to Simulator Flight Dynamics Model Formulation for Stall Training

A piloted simulation study was performed in the Cockpit Motion Facility at the National Aeronautics and Space Administration Langley Research Center. The research was motivated by the desire to reduce the commercial transport airplane fatal accident rate due to in-flight loss of control. The purpose of this study, which focused on a generic T-tail transport airplane, was to assess pilot sensitivity to flight dynamics model formulation used during a simulator stall recognition and recovery training/demonstration profile. To accomplish this, the flight dynamics model was designed with many configuration options. The model options were based on recently acquired static and dynamic stability and control data from sources that included wind tunnel, water tunnel, and computational fluid dynamics. The results, which are specific to a transport airplane stall recognition and recovery guided demonstration scenario, showed the two most important aerodynamic effects (other than stick pusher) to model were stall roll- off and the longitudinal static stability characteristic associated with the pitch break.

Cunningham, Kevin↗
NASA NTRS1993

Liquid hydrogen and liquid oxygen feedline passive recirculation analysis

The primary goal of the National Launch System (NLS) program was to design an operationally efficient, highly reliable vehicle with minimal recurring launch costs. To achieve this goal, trade studies of key main propulsion subsystems were performed to specify vehicle design requirements. These requirements include the use of passive recirculation to thermally condition the liquid hydrogen (LH2) and liquid oxygen (LO2) propellant feed systems and Space Transportation Main Engine (STME) fuel pumps. Rockwell International (RI) proposed a joint independent research and development (JIRAD) program with Marshall Space Flight Center (MSFC) to study the LH2 feed system passive recirculation concept. The testing was started in July 1992 and completed in November 1992. Vertical and sloped feedline designs were used. An engine simulator was attached at the bottom of the feedline. This simulator had strip heaters that were set to equal the corresponding heat input from different engines. A computer program is currently being used to analyze the passive recirculation concept in the LH2 vertical feedline tests. Four tests, where the heater setting is the independent variable, were chosen. While the JIRAD with RI was underway, General Dynamics Space Systems (GDSS) proposed a JIRAD with MSFC to explore passive recirculation in the LO2 feed system. Liquid nitrogen (LN2) is being used instead of LO2 for safety and economic concerns. To date, three sets of calibration tests have been completed on the sloped LN2 test article. The environmental heat was calculated from the calibration tests in which the strip heaters were turned off. During the LH2 testing, the environmental heat was assumed to be constant. Therefore, the total heat was equal to the environmental heat flux plus the heater input. However, the first two sets of LN2 calibration tests have shown that the environmental heat flux varies with heater input. A Systems Improved Numerical Differencing Analyzer and Fluid Integrator (SINDA/FLUINT) model is currently being built to determine if this variation in environmental heat is due to a change in the wall temperature. During the third set of calibration tests, a faulty reference junction was found. Based on this anomaly with the reference junction, the heat flux calculations from the first two calibration sets are now considered questionable.

Holt, Kimberly Ann↗
NASA NTRS1991

A Local Condensation Analysis Representing Two-phase Annular Flow in Condenser/radiator Capillary Tubes

NASA's effort for the thermal environmental control of the Space Station Freedom is directed towards the design, analysis, and development of an Active Thermal Control System (ATCS). A two phase, flow through condenser/radiator concept was baselined, as a part of the ATCS, for the radiation of space station thermal load into space. The proposed condenser rejects heat through direct condensation of ATCS working fluid (ammonia) in the small diameter radiator tubes. Analysis of the condensation process and design of condenser tubes are based on the available two phase flow models for the prediction of flow regimes, heat transfer, and pressure drops. The prediction formulas use the existing empirical relationships of friction factor at gas-liquid interface. An attempt is made to study the stability of interfacial waves in two phase annular flow. The formulation is presented of a stability problem in cylindrical coordinates. The contribution of fluid viscosity, surface tension, and transverse radius of curvature to the interfacial surface is included. A solution is obtained for Kelvin-Helmholtz instability problem which can be used to determine the critical and most dangerous wavelengths for interfacial waves.

Karimi, Amir↗