Search NASA⌕ Search

SEARCH · Search NASA

Results for “multiphase modeling”

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.

174 records · Page 10

Visualizing turbulent mixing of gases and particles

A physical model and interactive computer graphics techniques have been developed for the visualization of the basic physical process of stochastic dispersion and mixing from steady-state CFD calculations. The mixing of massless particles and inertial particles is visualized by transforming the vector field from a traditionally Eulerian reference frame into a Lagrangian reference frame. Groups of particles are traced through the vector field for the mean path as well as their statistical dispersion about the mean position by using added scalar information about the root mean square value of the vector field and its Lagrangian time scale. In this way, clouds of particles in a turbulent environment are traced, not just mean paths. In combustion simulations of many industrial processes, good mixing is required to achieve a sufficient degree of combustion efficiency. The ability to visualize this multiphase mixing can not only help identify poor mixing but also explain the mechanism for poor mixing. The information gained from the visualization can be used to improve the overall combustion efficiency in utility boilers or propulsion devices. We have used this technique to visualize steady-state simulations of the combustion performance in several furnace designs.

Ma, Kwan-Liu↗

Numerical Predictions of Dust-Induced Heat Flux Augmentation in Hypersonic Blunt-Body Flows Using a Discontinuous Galerkin Multiphase Flow Solver

Recent interest in human-scale missions to Mars has motivated the need for high-fidelity simulations of reentry flows. During a dust storm, there can be high levels of suspended dust in the Martian atmosphere, which cannot only enhance erosion of thermal protection systems but also transfer energy and momentum to the shock layer, thereby significantly augmenting the surface heat flux. Second-order finite-volume schemes are typically employed for hypersonic flow simulations, but such schemes suffer from a number of disadvantages. An attractive alternative is discontinuous Galerkin methods, which benefit from arbitrarily high spatial order of accuracy, geometric flexibility, and other properties. To enable accurate computations of high-speed particle-laden flows, an Euler-Lagrange methodology was developed in which the Eulerian field of the carrier gas is calculated using a discontinuous Galerkin scheme while the disperse phase is treated with Lagrangian particle tracking. We discuss challenges associated with coupling these two formulations and how to handle them. Momentum and energy transfer between the carrier gas and the particle phase is considered, and the importance of accounting for interparticle collisions is assessed. In addition, we describe the physical model of the particle phase and examine effects of its uncertainties on the numerical solution. We demonstrate the performance of the Euler-Lagrange method in representative testcases, with focus on the accurate prediction of particle trajectories and heating augmentation. Quantitative comparisons with experiments are provided.

Ching, Eric J.↗

Granular Material Flows with Interstitial Fluid Effects

The research focused on experimental measurements of the rheological properties of liquid-solid and granular flows. In these flows, the viscous effects of the interstitial fluid, the inertia of the fluid and particles, and the collisional interactions of the particles may all contribute to the flow mechanics. These multiphase flows include industrial problems such as coal slurry pipelines, hydraulic fracturing processes, fluidized beds, mining and milling operation, abrasive water jet machining, and polishing and surface erosion technologies. In addition, there are a wide range of geophysical flows such as debris flows, landslides and sediment transport. In extraterrestrial applications, the study of transport of particulate materials is fundamental to the mining and processing of lunar and Martian soils and the transport of atmospheric dust (National Research Council 2000). The recent images from Mars Global Surveyor spacecraft dramatically depict the complex sand and dust flows on Mars, including dune formation and dust avalanches on the slip-face of dune surfaces. These Aeolian features involve a complex interaction of the prevailing winds and deposition or erosion of the sediment layer; these features make a good test bed for the verification of global circulation models of the Martian atmosphere.

Hunt, Melany L.↗

The Diffuse Interstellar Medium

A major component of this research is to extend existing models of thermal processes in the local interstellar medium to the inner and outer Galaxy. In completing this goal we have calculated the thermal equilibrium gas temperature of the neutral diffuse gas and constructed phase diagrams (gas pressure versus density) for gas at galactic radii between 3 and 18 kpc. An important ingredient in this computation is the far-ultraviolet (FUV) radiation field in the Galactic disk. This radiation is important since photoelectric heating via FUV radiation on dust grains is expected to be a dominant heating process in the diffuse via gas. As a check of our calculated FUV field, we compared the infrared luminosity predicted by our theory with the observations of the COBE DIRBE space satellite and found the two to be in quite good agreement. Using our phase diagrams we have predicted the thermal pressure in the Galactic plane for which a multiphase equilibrium can be maintained. In addition, by using the maximum thermal pressure allowed by the observations, we have constrained the Galactic radii over which both cold and warm gas must exist, may exist, and cannot exist.

Hollenbach, David↗

Development of generalized pressure velocity coupling scheme for the analysis of compressible and incompressible combusting flows

The objective of this investigation has been to develop an algorithm (or algorithms) for the improvement of the accuracy and efficiency of the computer fluid dynamics (CFD) models to study the fundamental physics of combustion chamber flows, which are necessary ultimately for the design of propulsion systems such as SSME and STME. During this three year study (May 19, 1978 - May 18, 1992), a unique algorithm was developed for all speed flows. This newly developed algorithm basically consists of two pressure-based algorithms (i.e. PISOC and MFICE). This PISOC is a non-iterative scheme and the FICE is an iterative scheme where PISOC has the characteristic advantages on low and high speed flows and the modified FICE has shown its efficiency and accuracy to compute the flows in the transonic region. A new algorithm is born from a combination of these two algorithms. This newly developed algorithm has general application in both time-accurate and steady state flows, and also was tested extensively for various flow conditions, such as turbulent flows, chemically reacting flows, and multiphase flows.

Chen, C. P.↗

The Influence of the Photoionizing Radiation Spectrum on Metal-Line Ratios in Ly(alpha) Forest Clouds

Recent measurements of Si IV/C IV ratios in the high-redshift Ly(alpha) forest (Songaila & Cowie, AJ, 112, 335 (1996a); Savaglio et at., A&A (in press) (1997)) have opened a new window on chemical enrichment and the first generations of stars. However, the derivation of accurate Si/C abundances requires reliable ionization corrections, which are strongly dependent on the spectral shape of the metagalactic ionizing background and on the 'local effects' of hot stars in nearby galaxies. Recent models have assumed power-law quasar ionizing backgrounds plus a decrement at 4 Ryd to account for He II attenuation in intervening clouds. However, we show that realistic ionizing backgrounds based on cosmological radiative transfer models produce more complex ionizing spectra between 1-5 Ryd that are critical to interpreting ions of Si and C. We also make a preliminary investigation of the effects of He II ionization front nonoverlap. Because the attenuation and reemission by intervening clouds enhance Si IV relative to C the observed high Si IV/C IV ratios do not require an unrealistic Si overproduction (Si/C greater than or equal to 3 (Si/C)(solar mass)). If the ionizing spectrum is dominated by 'local effects' from massive stars, even larger Si IV/C IV ratios are possible. However, unless stellar radiation dominates quasars by more than a factor of 10, we confirm the evidence for some Si overproduction by massive stars; values Si/C approx. 2(Si/C)(solar mass) fit the measurements better than solar abundances. Ultimately, an adequate interpretation of the ratios of C IV, Si IV, and C II may require hot, collisionally ionized gas in a multiphase medium.

Giroux, Mark L.↗

In-flight surface-flow measurements on a subsonic transport high-lift flap system

As part of a multiphased program for subsonic transport high-lift flight research, flight tests were conducted on the Transport Systems Research Vehicle (B737-100 aircraft) at the NASA Langley Research Center, to obtain detailed flow characteristics of the high-lift flap system for correlation with computational and wind-tunnel investigations. Pressure distributions, skin friction, and flow-visualization measurements were made on a triple-slotted flap system for a range of flap deflections, chord Reynolds numbers (10 to 21 million), and Mach numbers (0.16 to 0.36). Experimental test results are given for representative flap settings indicating flow separation on the fore-flap element for the largest flap deflection. Comparisons of the in-flight flow measurements were made with predictions from available viscous multielement computational methods modified with simple-sweep theory. Computational results overpredicted the experimentally measured pressures, particularly in the case involving separation of the fore lap, indicating the need for better modeling of confluent boundary layers and three-dimensional sweep effects.

Yip, Long P.↗

Forming Perchlorates on Mars Through Plasma Chemistry During Dust Events

We report experimental evidences to support a new formation mechanism, multiphase redox plasma chemistry, for perchlorate on Mars observed during the Phoenix mission, whose high concentrations and high ClO4/Cl ratio cannot be fully interpreted by photochemistry. This chemical reaction occurs between Cl-bearing minerals on the Mars surface and free radicals generated by electrostatic discharge (ESD) during Mars dust events (dust storms, dust devils, and grain saltation). We conducted simulated ESD experiments in a Mars chamber with pure CO2, CO2+H2O(g), and Mars Simulate Gas Mixture at Martian atmospheric pressure. We directly observed (1) the instantaneous generation of atmospheric free radicals CO2+, CO+, OI, HIII, HII, OH, ArI, N2, and N2+in normal glow discharge (NGD), detected by in situ plasma emission spectroscopy, and O3by UV and Mid-IR spectroscopy; (2) the fast transformation of NaCl to NaClO3and NaClO4detected by laser Raman spectroscopy, with oxychlorine enrichment at the sample surfaces confirmed by ion chromatography. Through two sets of experimental comparison, we found that the oxidation power of ESD-electron is three orders of magnitude higher than that of UVC-photon. We scaled our experimental results to the modeled ESD in Mars dust events and Mars surface UV radiation level, and concluded that plasma chemistry occurred during Mars dust events can be an additional important formation mechanism for the large amounts of perchlorates observed during various missions to Mars.

Wu, Zhongchen↗

A Study into Validating a Coupled Method of Characteristics and Direct Simulation Monte Carlo Method Against Empirical Data

The following will outline the methodology and results of validating a coupled Method of Characteristics (MOC) and Direct Simulation Monte Carlo (DSMC) method. This research focused specifically on modeling plume impingement, induced by Reaction Control System (RCS) thrusters that flew on the National Aeronautics and Space Administration’s (NASA’s) space shuttle Discovery. For each simulation, the continuum portion of the RCS thruster was simulated using MOC for solving hyperbolic Partial Differential Equations (PDEs) and computed with the NASA code, Reacting and Multi-phase Program (RAMP). The solution was then implemented as a starting condition into the NASA DSMC code, Direct Simulation and Monte Carlo Analysis Code (DAC). Typically, DSMC models rely on code-to-code validation for fidelity. The significance of this research is in its ability to validate its models against empirical data. Prior to computing solutions for these simulations, the mesh size and structure were optimized and many variants of DSMC input parameters were iterated on in order to acquire a reliable, mesh-independent, fully optimized numerical solution. This research will discuss the mathematical formulation of MOC for nozzle flow and DSMC for rarefied gases. Additionally, it will provide an explanation of how to implement these mathematical concepts into the two solvers: RAMP and DAC. Ultimately, this research will demonstrate that the overall process illustrated produces results in good agreement with empirical data. As a consequence, the methodology presented is granted an increased level of confidence and will greatly contribute to the aerospace industry and its effort in understanding and predicting rarefied flow fields.

Direct Simulation Monte Carlo Analysis Code↗

On Lateral Viscosity Contrast in the Mantle and the Rheology of Low-Frequency Geodynamics

Mantle-wide heterogeneity is largely controlled by deeply penetrating thermal convective currents. These thermal currents are likely to produce significant lateral variation in rheology, and this can profoundly influence overall material behaviour. How thermally related lateral viscosity variations impact models of glacio-isostatic and tidal deformation is largely unknown. An important step towards model improvement is to quantify, or bound, the actual viscosity variations that characterize the mantle. Simple scaling of viscosity to shear-wave velocity fluctuations yields map-views of long- wavelength viscosity variation. These give a general quantitative description and aid in estimating the depth dependence of rheological heterogeneity throughout the mantle. The upper mantle is probably characterized by two to four orders of magnitude variation (peak-to-peak). Discrepant time-scales for rebounding Holocene shorelines of Hudson Bay and southern Iceland are consistent with this characterization. Results are given in terms of a local average viscosity ratio, (Delta)eta(bar)(sub i), of volumetric concentration, phi(sub i). For the upper mantle deeper than 340 km the following reasonable limits are estimated for (delta)eta(bar) approx. equal 10(exp -2): 0.01 less than or equal to phi less than or equal to 0.15. A spectrum of ratios (Delta)eta(bar)(sub i) less than 0.1 at concentration level eta(sub i) approx. equal 10(exp -6) - 10(exp -1) in the lower mantle implies a spectrum of shorter time-scale deformational response modes for second-degree spherical harmonic deformations of the Earth. Although highly uncertain, this spectrum of spatial variation allows a purely Maxwellian viscoelastic rheology simultaneously to explain all solid tidal dispersion phenomena and long-term rebound-related mantle viscosity. Composite theory of multiphase viscoelastic media is used to demonstrate this effect.

Ivins, Erik R.↗

Properties and processing characteristics of low density carbon cloth phenolic composites

Ply-lift and pocketing are two critical anomalies of carbon cloth phenolic composites (CCPC) in rocket nozzle applications. Ply lift occurs at low temperatures when the A/P and in-plane permeabilities of the composite materials are still very low and in-plane porous paths are blocked. Pocketing occurs at elevated temperatures when in-plane permeability is reduced by the A/P compressive stress. The thermostructural response of CCPC in a rapid heating environment involves simultaneous heat, mass, and momentum transfer along with the degradation of phenolic resin in a multiphase system with temperature- and time-dependent material properties as well as dynamic processing conditions. Three temperature regions represent the consequent chemical reactions, material transformations, and property transitions, and provide a quick qualitative method for characterizing the thermostructural behavior of a CCPC. In order to optimize the FM5939 LDCCP (low density carbon cloth phenolic) for the nozzle performance required in the Advanced Solid Rocket Motor (ASRM) program, a fundamental study on LDCCP materials was conducted. The cured composite has a density of 1.0 +/- 0.5 gm/cc which includes 10 to 25 percent void volume. The weight percent of carbon microballoon is low (7-15 percent). However, they account for approximately one third of the volume and historically their percentages have not been controlled very tightly. In addition, the composite properties show no correlation with microballoon weight percent or fiber properties (e.g. fiber density or fiber moisture adsorption capacity). Test results concerning the ply-lift anomaly in the MNASA motor firings were: (1) Steeper ply angle (shorter path lenght) designs minimized/eliminated by lifting, (2) material with higher void volume ply lifted less frequently, (3) materials with high (greater than 9 percent) microballoon content had a higher rate of ply lifting, and (4) LDCCP materials failed at microballoon-resin interfaces. The objectives of this project are: (1) to investigate the effects of carbon microballoon and cabosil fillers as well as fiber heat treatment on plylift-related mechanical properties, (2) to develop a science-based thermostructural process model for the carbon phenolics. The model can be used in the future for the selection of the improved ASRM materials, (3) to develop the micro-failure mechanisms for the ply-lift initiation and propagation processes during the thermoelastic region of phenolic degradation, i.e. postcuring and devolatilization.

Wang, C. Jeff↗

In-flight pressure distributions and skin-friction measurements on a subsonic transport high-lift wing section

Flight experiments are being conducted as part of a multiphased subsonic transport high-lift research program for correlation with wind-tunnel and computational results. The NASA Langley Transport Systems Research Vehicle (B737-100 aircraft) is used to obtain in-flight flow characteristics at full-scale Reynolds numbers to contribute to the understanding of 3-D high-lift, multi-element flows including attachment-line transition and relaminarization, confluent boundary-layer development, and flow separation characteristics. Flight test results of pressure distributions and skin friction measurements were obtained for a full-chord wing section including the slat, main-wing, and triple-slotted, Fowler flap elements. Test conditions included a range of flap deflections, chord Reynolds numbers (10 to 21 million), and Mach numbers (0.16 to 0.40). Pressure distributions were obtained at 144 chordwise locations of a wing section (53-percent wing span) using thin pressure belts over the slat, main-wing, and flap elements. Flow characteristics observed in the chordwise pressure distributions included leading-edge regions of high subsonic flows, leading-edge attachment-line locations, slat and main-wing cove-flow separation and reattachment, and trailing-edge flap separation. In addition to the pressure distributions, limited skin-friction measurements were made using Preston-tube probes. Preston-tube measurements on the slat upper surface suggested relaminarization of the turbulent flow introduced by the pressure belt on the slat leading-edge surface when the slat attachment line was laminar. Computational analysis of the in-flight pressure measurements using two-dimensional, viscous multielement methods modified with simple-sweep theory showed reasonable agreement. However, overprediction of the pressures on the flap elements suggests a need for better detailed measurements and improved modeling of confluent boundary layers as well as inclusion of three-dimensional viscous effects in the analysis.

Yip, Long P.↗