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

An Integrated Solution for Performing Thermo-fluid Conjugate Analysis

A method has been developed which integrates a fluid flow analyzer and a thermal analyzer to produce both steady state and transient results of 1-D, 2-D, and 3-D analysis models. The Generalized Fluid System Simulation Program (GFSSP) is a one dimensional, general purpose fluid analysis code which computes pressures and flow distributions in complex fluid networks. The MSC Systems Improved Numerical Differencing Analyzer (MSC.SINDA) is a one dimensional general purpose thermal analyzer that solves network representations of thermal systems. Both GFSSP and MSC.SINDA have graphical user interfaces which are used to build the respective model and prepare it for analysis. The SINDA/GFSSP Conjugate Integrator (SGCI) is a formbase graphical integration program used to set input parameters for the conjugate analyses and run the models. The contents of this paper describes SGCI and its thermo-fluids conjugate analysis techniques and capabilities by presenting results from some example models including the cryogenic chill down of a copper pipe, a bar between two walls in a fluid stream, and a solid plate creating a phase change in a flowing fluid.

Kornberg, Oren↗

Fluids and Combustion Facility: Combustion Integrated Rack Modal Model Correlation

The Fluids and Combustion Facility (FCF) is a modular, multi-user, two-rack facility dedicated to combustion and fluids science in the US Laboratory Destiny on the International Space Station. FCF is a permanent facility that is capable of accommodating up to ten combustion and fluid science investigations per year. FCF research in combustion and fluid science supports NASA's Exploration of Space Initiative for on-orbit fire suppression, fire safety, and space system fluids management. The Combustion Integrated Rack (CIR) is one of two racks in the FCF. The CIR major structural elements include the International Standard Payload Rack (ISPR), Experiment Assembly (optics bench and combustion chamber), Air Thermal Control Unit (ATCU), Rack Door, and Lower Structure Assembly (Input/Output Processor and Electrical Power Control Unit). The load path through the rack structure is outlined. The CIR modal survey was conducted to validate the load path predicted by the CIR finite element model (FEM). The modal survey is done by experimentally measuring the CIR frequencies and mode shapes. The CIR model was test correlated by updating the model to represent the test mode shapes. The correlated CIR model delivery is required by NASA JSC at Launch-10.5 months. The test correlated CIR flight FEM is analytically integrated into the Shuttle for a coupled loads analysis of the launch configuration. The analysis frequency range of interest is 0-50 Hz. A coupled loads analysis is the analytical integration of the Shuttle with its cargo element, the Mini Payload Logistics Module (MPLM), in the Shuttle cargo bay. For each Shuttle launch configuration, a verification coupled loads analysis is performed to determine the loads in the cargo bay as part of the structural certification process.

McNelis, Mark E.↗

Time-dependent models of single- and double-layer mantle convection

Time-dependent numerical models of two-dimensional convection which simulate mantle convection are presented. Starting from a fluid initially at rest with a purely conductive temperature profile, it is found that for a Rayleigh number of ten million, the model experiences a transient period of double-layer convection which lasts on the order of hundreds of millions of years when scaled to the earth's mantle. It is suggested that transient periods of double-layer convection in the earth's mantle may have existed in the past, whether or not such a period exists today.

Boss, A. P.↗

A comparison of two plasma models

The time dependent behavior of a plasma which surrounds a highly biased conducting sphere is considered. This plasma is treated as either a cold two component fluid or as a warm plasma whose time development can be found by solving the Vlasov equation. Both models demonstrate oscillatory behavior, but the electric fields predicted by the models are quantitatively quite different in regions close to the surface of the sphere and very similar otherwise. A broadening of the electron distribution function with time is observed indicating local heating of the plasma near the surface of the sphere.

Cottam, Russell↗

The structure and evolution of Jupiter - The fluid contraction stage

The complete evolution of a contracting star of Jovian mass consisting of a convective adiabatic homogeneous fluid is determined using stellar structure methods, improved model atmosphere calculations, and substantially improved thermodynamic properties for hydrogen and hydrogen-helium fluids. The model atmospheres are calculated in the form of time-averaged vertical temperature structures, including all relevant sources of opacity and a solar energy deposition component, and the thermodynamic properties are modified to obtain better agreement with Monte Carlo results for metallic fluids. The resultant gravitationally contracting evolutionary models are found to have two phases: an early stellar phase similar to a typical low-mass pre-main-sequence body and a later phase constituting an approach to a degenerate-dwarf cooling curve. The first phase is shown to have high luminosities and internal temperatures, while the second gives excellent agreement with the observed radius and luminosity of Jupiter. Analysis indicates that the equation of state and superadiabaticity have the strongest influence on evolution over planetary time scales.

Graboske, H. C., Jr.↗

The Hadley and Rossby regimes in a spherical atmosphere

The properties of the steady Hadley and Rossby regimes for a thermally forced rotating fluid on a sphere are studied. The two layer modified geostrophic model is employed which allows for thermal advection by the divergent wind and time dependent static stability. Heating processes are parameterized using the Newtonian approximation and Rayleigh friction is accounted for. The equations are transformed to spectral form using spherical harmonics and then truncated retaining a simple axisymmetric state and initial, one wave. A time independent Hadley circulation is obtained which is neutral to axisymmetric disturbances but unstable to wave like perturbations for intermediate values of the meridional temperature gradient, indicating the existence of both an upper and lower symmetric Hadley regime. An analytical solution for the steady Rossby circulation is determined for values of the meridional temperature gradient where the Hadley regime is unstable. Linear perturbation theory is used to show that within the steady Rossby regime two or more waves cannot exist simultaneously.

Feldstein, S. B.↗

Interfacing the Generalized Fluid System Simulation Program with the SINDA/G Thermal Program

A general purpose, one dimensional fluid flow code has been interfaced with the thermal analysis program SINDA/G. The flow code, GFSSP, is capable of analyzing steady state and transient flow in a complex network. The flow code is capable of modeling several physical phenomena including compressibility effects, phase changes, body forces (such as gravity and centrifugal) and mixture thermodynamics for multiple species. The addition of GFSSP to SINDA/G provides a significant improvement in convective heat transfer modeling for SINDA/G. The interface development was conducted in two phases. This paper describes the first (which allows for steady and quasi-steady - unsteady solid, steady fluid - conjugate heat transfer modeling). The second (full transient conjugate heat transfer modeling) phase of the interface development will be addressed in a later paper. Phase 1 development has been benchmarked to an analytical solution with excellent agreement. Additional test cases for each development phase demonstrate desired features of the interface. The results of the benchmark case, three additional test cases and a practical application are presented herein.

Schallhorn, Paul↗

Computational aeroelasticity

Recent large-scale studies of computational unsteady aerodynamics for aeroelastic applications are reviewed. The variety of fluid dynamic flow models available to address such computations are illustrated in two cases: (1) 2D Navier-Stokes (NS) equations illustrating the highest level of modeling usually employed and (2) the transonic small-disturbance potential equation representing the entry level for nonlinear flow modeling. Estimates of computer resources necessary to produce accurate converged results are given. Application of potential and NS equation codes for flows which are generally at lower angles and high speeds are addressed, including (1) the treatment of wings and configuration details using potential equation codes and (2) recent NS code calculations of complete vehicle configurations. Attention is given to high-angle conditions involving separated vortex-dominated flows. Accuracy requirements for vortex shedding over airfoils is discussed and the calculation of vorticity convected over significant distances is addressed. Steady CFD computations about delta wings at high angles are given, including a discussion of the required level of fluid dynamic flow modeling. Recent results on unsteady 'buffetlike flow' about complete vehicle models are reviewed.

Edwards, John W.↗

Computational Fluid Dynamic Analyses for the High-Lift Common Research Model Using the USM3D and FUN3D Flow Solvers

Two Navier-Stokes codes were used to compute flow over the High-Lift Common Research Model (HL-CRM) in preparation for a wind tunnel test to be performed at the NASA Langley Research Center 14-by-22-Foot Subsonic Tunnel in fiscal year 2018. Both flight and wind tunnel conditions were simulated by the two codes at set Mach numbers and Reynolds numbers over a full angle-of-attack range for three configurations: cruise, landing and takeoff. Force curves, drag polars and surface pressure contour comparisons are shown for the two codes. The lift and drag curves compare well for the cruise configuration up to 10deg angle of attack but not as well for the other two configurations. The drag polars compare reasonably well for all three configurations. The surface pressure contours compare well for some of the conditions modeled but not as well for others.

Rivers, Melissa↗

Retrieving the global distribution of the threshold of wind erosion from satellite data and implementing it into the Geophysical Fluid Dynamics Laboratory land–atmosphere model (GFDL AM4.0/LM4.0)

Dust emission is initiated when surface wind velocities exceed the threshold of wind erosion. Many dust models used constant threshold values globally. Here we use satellite products to characterize the frequency of dust events and land surface properties. By matching this frequency derived from Moderate Resolution Imaging Spectroradiometer (MODIS) Deep Blue aerosol products with surface winds, we are able to retrieve a climatological monthly global distribution of the wind erosion threshold (V(threshold)) over dry and sparsely vegetated surfaces. This monthly two-dimensional threshold velocity is then implemented into the Geophysical Fluid Dynamics Laboratory coupled land–atmosphere model (AM4.0/LM4.0). It is found that the climatology of dust optical depth (DOD) and total aerosol optical depth, surface PM10 dust concentrations, and the seasonal cycle of DOD are better captured over the “dust belt” (i.e., northern Africa and the Middle East) by simulations with the new wind erosion threshold than those using the default globally constant threshold. The most significant improvement is the frequency distribution of dust events, which is generally ignored in model evaluation. By using monthly rather than annual mean V(threshold), all comparisons with observations are further improved. The monthly global threshold of wind erosion can be retrieved under different spatial resolutions to match the resolution of dust models and thus can help improve the simulations of dust climatology and seasonal cycles as well as dust forecasting.

Bing Pu↗

Effect of fluid forces on rotor stability of centrifugal compressors and pumps

A simple two dimensional model for calculating the rotordynamic effects of the impeller force in centrifugal compressors and pumps is presented. It is based on potential flow theory with singularities. Equivalent stiffness and damping coefficients are calculated for a machine with a vaneless volute formed as a logarithmic spiral. It is shown that for certain operating conditions, the impeller force has a destablizing effect on the rotor.

Colding-Jorgensen, J.↗

A Modified Wall Matching Treatment to Account for Local Solid to Fluid Thermal Coupling

The wall-matching methodology of Wilcox is modified to include a solid-wall, thermal-conduction model. This coupled fluid-thermal-structure model is derived assuming that the wall thermal-structure behavior is locally one-dimensional and that structural deformations, due to thermally induced stresses, are not significant. The one-dimensional coupled fluid-thermal-structure model is derived such that the wall temperature is removed as an independent boundary condition variable. The one-dimensional coupled fluid-thermal-structure model is also derived for the general case of an arbitrary mixture of thermally prefect gases and a wall of arbitrary thickness and conductivity by using a compressible, streamwise-pressure-gradient-corrected, wall-matching function and Fourier's law of heat conduction. The resulting model was implemented in the VULCAN CFD code as a new boundary condition type. VULCAN was then used to simulate a two-dimensional Mach 6 wind tunnel facility nozzle flow to demonstrate/validate the one-dimensional coupled fluid-thermal-structure model. The nozzle internal-wall surface temperature and heat transfer distributions computed using the one-dimensional coupled fluid-thermal-structure model are compared to wall temperature and heat transfer distributions from an iterative multi-dimensional analysis obtained by coupling the VULCAN CFD code and the MSC/NASTRAN-thermal code. The one-dimensional coupled fluid-thermal-structure model analysis is shown to be very robust and in excellent agreement with the multi-dimensional iteratively coupled analysis. It is also shown that the one-dimensional analysis can be used as an initial guess for the multi-dimensional iteratively coupled analysis.

White, Jeffery A.↗

A Modified Wall Matching Treatment to Account for Local Solid to Fluid Thermal Coupling

The wall-matching methodology of Wilcox is modified to include a solid-wall, thermal-conduction model. This coupled fluid-thermal-structure model is derived assuming that the wall thermal-structure behavior is locally one-dimensional and that structural deformations, due to thermally induced stresses, are not significant. The one-dimensional coupled fluid-thermal-structure model is derived such that the wall temperature is removed as an independent boundary condition variable. The one-dimensional coupled fluid-thermal-structure model is also derived for the general case of an arbitrary mixture of thermally prefect gases and a wall of arbitrary thickness and conductivity by using a compressible, streamwise-pressure-gradient-corrected, wall-matching function and Fourier's law of heat conduction. The resulting model was implemented in the VULCAN CFD code as a new boundary condition type. VULCAN was then used to simulate a two-dimensional Mach 6 wind tunnel facility nozzle flow to demonstrate/validate the one-dimensional coupled fluid-thermal-structure model. The nozzle internal-wall surface temperature and heat transfer distributions computed using the one-dimensional coupled fluid-thermal-structure model are compared to wall temperature and heat transfer distributions from an iterative multi-dimensional analysis obtained by coupling the VULCAN CFD code and the MSC/NASTRAN-thermal code. The one-dimensional coupled fluid-thermal-structure model analysis is shown to be very robust and in excellent agreement with the multi-dimensional iteratively coupled analysis. It is also shown that the one-dimensional analysis can be used as an initial guess for the multi-dimensional iteratively coupled analysis.

White, Jeffery A.↗

An Analog Computer Study of the Low-Frequency Flow Dynamics of Two Nuclear-Rocket Cold-Flow Engine Systems

An analog computer study is made of the low-frequency flow dynamics observed with liquid hydrogen flowing in two different nuclear-rocket cold-flow (unfueled) engine systems. The prediction of the low-frequency response depends largely upon the generation of the proper mass-flow storage in the system. This is achieved by relating, in the model lumping, the average values of fluid enthalpy and a certain thermodynamic fluid property. A two-phase state equation is derived for thermodynamic equilibrium between the gas and liquid phases. Parametric studies are made which demonstrate the effect of the pump characteristics on the flow damping.

Fox, Harry W., Jr.↗

Freezable Radiator Model Correlation and Full Scale Design

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 space craft s surroundings and because of different thermal loads during different mission phases. However, freezing and thawing (recovering) a radiator is a process that has historically proven very difficult to predict through modeling, resulting in highly inaccurate predictions of recovery time. This paper summarizes efforts made to correlate a Thermal Desktop (TM) model with empirical testing data from two test articles. A 50-50 mixture of DowFrost HD and water is used as the working fluid. Efforts to scale this model to a full scale design, as well as efforts to characterize various thermal control fluids at low temperatures are also discussed.

Lillibridge, Sean T.↗

Characterization of the Ignition Over-Pressure/Sound Suppression Water in the Space Launch System Mobile Launcher Using Volume of Fluid Modeling

The Space Launch System (SLS) Vehicle consists of a Core Stage with four RS-25 engines and two Solid Rocket Boosters (SRBs). This vehicle is launched from the Launchpad using a Mobile Launcher (ML) which supports the SLS vehicle until its liftoff from the ML under its own power. The combination of the four RS-25 engines and two SRBs generate a significant Ignition Over-Pressure (IOP) and Acoustic Sound environment. One of the mitigations of these environments is the Ignition Over-Pressure/Sound Suppression (IOP/SS) subsystem installed on the ML. This system consists of six water nozzles located parallel to and 24 inches downstream of each SRB nozzle exit plane as well as 16 water nozzles located parallel to and 53 inches downstream of the RS-25 nozzle exit plane. During launch of the SLS vehicle, water is ejected through each water nozzle to reduce the intensity of the transient pressure environment imposed upon the SLS vehicle. While required for the mitigation of the transient pressure environment on the SLS vehicle, the IOP/SS subsystem interacts (possibly adversely) with other systems located on the Launch Pad. One of the other systems that the IOP/SS water is anticipated to interact with is the Hydrogen Burn-Off Igniter System (HBOI). The HBOI system's purpose is to ignite the unburned hydrogen/air mixture that develops in and around the nozzle of the RS-25 engines during engine start. Due to the close proximity of the water system to the HBOI system, the presence of the IOP/SS may degrade the effectiveness of the HBOI system. Another system that the IOP/SS water may interact with adversely is the RS-25 engine nozzles and the SRB nozzles. The adverse interaction anticipated is the wetting, to a significant degree, of the RS-25 nozzles resulting in substantial weight of ice forming and water present to a significant degree upstream of the SRB nozzle exit plane inside the nozzle itself, posing significant additional blockage of the effluent that exits the nozzle upon motor start leading to detrimental effects. The purpose of the CFD simulations were to i) characterize the location of the IOP/SS water after it is ejected from the IOP/SS nozzles, ii) characterize the interaction of the IOP/SS system with the HBOI system and iii) characterize the interaction of the IOP/SS water with the RS-25 nozzles and the SRB nozzles.

West, Jeff↗

Numerical simulations of loops heated to solar flare temperatures. III - Asymmetrical heating

A numerical model is defined for asymmetric full solar flare loop heating and comparisons are made with observational data. The Dynamic Flux Tube Model is used to describe the heating process in terms of one-dimensional, two fluid conservation equations of mass, energy and momentum. An adaptive grid allows for the downward movement of the transition region caused by an advancing conduction front. A loop 20,000 km long is considered, along with a flare heating system and the hydrodynamic evolution of the loop. The model was applied to generating line profiles and spatial X-ray and UV line distributions, which were compared with SMM, P78-1 and Hintori data for Fe, Ca and Mg spectra. Little agreement was obtained, and it is suggested that flares be treated as multi-loop phenomena. Finally, it is concluded that chromospheric evaporation is not an effective mechanism for generating the soft X-ray bursts associated with flares.

Cheng, C.-C.↗

Numerical Simulation of No Vent Chill and Fill of a Large Liquid Hydrogen Tank

On orbit transfer of cryogenic liquids is necessary for many future space missions. The techniques required to transfer cryogens in micro gravity are quite different from ground system applications. In normal gravity, cryogenic tanks are usually filled from the bottom at nearly atmospheric pressure. The vapor generated due to heat transfer from warm tank walls, is allowed to vent from the top of the tank while the tank is being filled. If the normal gravity technique is used on orbit, the uncertainty of liquid and vapor distributions in unsettled condition may result in dumping of large amounts of liquid propellant overboard. The Vented Chill and No Vent Fill (VCNVF) process is a methodology [1] used to reduce propellant loss by holding the vent valve open during initial wall chilldown, preventing over-pressurization of the tank and stalling of the incoming flow, but closing the tank vent valve after chilling the tank wall sufficiently. An accurate simulation of no vent tank chill and fill is necessary to estimate the amount of propellant necessary to chill the tank, determine the maximum sustainable flowrate to fill the tank without stalling, and establish a timeline for the filling process. Numerical modeling of filling a cryogenic tank is complex because the flow is two-phase and the process involves both boiling and condensation with heat transfer between solid and fluid. In recent years, progress has been made to develop a numerical modeling technique [2,3,4] that has been verified by comparing with the test data from a small flight tank using liquid nitrogen [5]. The purpose of the present paper is to apply this technique in a larger tank with liquid hydrogen and compare with the test data of Chato [1]. The K-site Test Facility [1], shown in Figure 1, consists of a Test Tank, spray system, instrumentation, and the vacuum chamber. The test tank selected was ellipsoidal with an 87 inch major diameter and a 1.2 to 1 major to minor axis ratio. The two ends are joined by a short 1.5 inch cylindrical section. The tank is made of 2219 aluminum chemically milled to a nominal thickness of 0.087 inches. The tank weighs 329.25 pounds, and the tank volume is 175 ft 3 . The GFSSP model of the Test Tank is shown in Figure 2. Node 13 is a boundary node that represents the supply tank which is supplying liquid hydrogen at a given pressure and temperature. The total measured flowrate is evenly distributed through Branches 131 through 139. Nodes 1 through 9 represent the test tank and they are connected with metal solid nodes 14 through 22 through fluid-to-solid conductors that allow convective heat transfer between the solid and fluid nodes. The predicted tank pressure, filling rate, wall temperature and fluid temperature will be compared with test data.

Alok Majumdar↗