Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Particle Methods”

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 1,423 records · Page 79
All sourcesNASA NTRSDOE OSTI
NASA NTRS2015

An Experimental Investigation of the Flow Physics Associated With End Wall Losses and Large Rotor Tip Clearances as Found in the Rear Stages of a High Pressure Compressor

The focus of this work was to characterize the fundamental flow physics and the overall performance effects due to increased rotor tip clearance heights in axial compressors. Data have been collected in the three-stage axial research compressor at Purdue University with a specific focus on analyzing the multistage effects resulting from the tip leakage flow. Three separate rotor tip clearance heights were studied with nominal tip clearance heights of 1.5%, 3.0%, and 4.0% based on a constant annulus height. Overall compressor performance was investigated at four corrected speedlines (100%, 90%, 80%, and 68%) for each of the three tip clearance configurations using total pressure and total temperature rakes distributed throughout the compressor. The results have confirmed results from previous authors showing a decrease of total pressure rise, isentropic efficiency, and stall margin which is approximately linear with increasing tip clearance height. The stall inception mechanisms have also been evaluated at the same corrected speeds for each of the tip clearance configurations. Detailed flow field measurements have been collected at two loading conditions, nominal loading (NL) and high loading (HL), on the 100% corrected speedline for the smallest and largest tip clearance heights (1.5% and 4.0%). Steady detailed radial traverses of total pressure at the exit of each stator row have been supported by flow visualization techniques to identify regions of flow recirculation and separation. Furthermore, detailed radial traverses of time-resolved total pressures at the exit of each rotor row have been measured with a fast-response pressure probe. These data have helped to quantify the size of the leakage flow at the exit of each rotor. Thermal anemometry has also been implemented to evaluate the time-resolved three-dimensional components of velocity throughout the compressor and calculate blockage due to the rotor tip leakage flow throughout the compressor. These measurements have also been used to calculate streamwise vorticity. Time-resolved static pressure measurements have been collected over the rotor tips for all rotors with each of the three tip clearance configurations for up to five loading conditions along the 100% corrected speedline using fast-response piezoresistive pressure sensors. These time-resolved static pressure measurements, as well as the time-resolved total pressures and velocities have helped to reveal a profound influence of the upstream stator vane on the size and shape of the rotor tip leakage flow. Finally, a novel particle image velocimetry (PIV) technique has been developed as a proof-of- concept. In contrast to PIV methods that have been typically been utilized for turbomachinery applications in the past, the method used for this study introduced the laser light through the same access window that was also used to image the flow. This new method addresses potential concerns related to the intrusive laser-introducing techniques that have typically been utilized by other authors in the past. Ultimately, the data collected for this project represent a unique data set which contributes to build a better understanding of the tip leakage flow field and its associated loss mechanisms. These data will facilitate future engine design goals leading to small blade heights in the rear stages of high pressure compressors and aid in the development of new blade designs which are desensitized to the performance penalties attributed to rotor tip leakage flows.

turbomachinery↗
NASA NTRS2016

Model-Based Fatigue Prognosis of Fiber-Reinforced Laminates Exhibiting Concurrent Damage Mechanisms

Prognostics of large composite structures is a topic of increasing interest in the field of structural health monitoring for aerospace, civil, and mechanical systems. Along with recent advancements in real-time structural health data acquisition and processing for damage detection and characterization, model-based stochastic methods for life prediction are showing promising results in the literature. Among various model-based approaches, particle-filtering algorithms are particularly capable in coping with uncertainties associated with the process. These include uncertainties about information on the damage extent and the inherent uncertainties of the damage propagation process. Some efforts have shown successful applications of particle filtering-based frameworks for predicting the matrix crack evolution and structural stiffness degradation caused by repetitive fatigue loads. Effects of other damage modes such as delamination, however, are not incorporated in these works. It is well established that delamination and matrix cracks not only co-exist in most laminate structures during the fatigue degradation process but also affect each other's progression. Furthermore, delamination significantly alters the stress-state in the laminates and accelerates the material degradation leading to catastrophic failure. Therefore, the work presented herein proposes a particle filtering-based framework for predicting a structure's remaining useful life with consideration of multiple co-existing damage-mechanisms. The framework uses an energy-based model from the composite modeling literature. The multiple damage-mode model has been shown to suitably estimate the energy release rate of cross-ply laminates as affected by matrix cracks and delamination modes. The model is also able to estimate the reduction in stiffness of the damaged laminate. This information is then used in the algorithms for life prediction capabilities. First, a brief summary of the energy-based damage model is provided. Then, the paper describes how the model is embedded within the prognostic framework and how the prognostics performance is assessed using observations from run-to-failure experiments

Prognostics↗
NASA NTRS2015

An Experimental Investigation of the Flow Physics Associated With End Wall Losses and Large Rotor Tip Clearances as Found in the Rear Stages of a High Pressure Compressor

The focus of this work was to characterize the fundamental flow physics and the overall performance effects due to increased rotor tip clearance heights in axial compressors. Data have been collected in the three-stage axial research compressor at Purdue University with a specific focus on analyzing the multistage effects resulting from the tip leakage flow. Three separate rotor tip clearance heights were studied with nominal tip clearance heights of 1.5%, 3.0%, and 4.0% based on a constant annulus height. Overall compressor performance was investigated at four corrected speedlines (100%, 90%, 80%, and 68%) for each of the three tip clearance configurations using total pressure and total temperature rakes distributed throughout the compressor. The results have confirmed results from previous authors showing a decrease of total pressure rise, isentropic efficiency, and stall margin which is approximately linear with increasing tip clearance height. The stall inception mechanisms have also been evaluated at the same corrected speeds for each of the tip clearance configurations. Detailed flow field measurements have been collected at two loading conditions, nominal loading (NL) and high loading (HL), on the 100% corrected speedline for the smallest and largest tip clearance heights (1.5% and 4.0%). Steady detailed radial traverses of total pressure at the exit of each stator row have been supported by flow visualization techniques to identify regions of flow recirculation and separation. Furthermore, detailed radial traverses of time-resolved total pressures at the exit of each rotor row have been measured with a fast-response pressure probe. These data have helped to quantify the size of the leakage flow at the exit of each rotor. Thermal anemometry has also been implemented to evaluate the time-resolved three-dimensional components of velocity throughout the compressor and calculate blockage due to the rotor tip leakage flow throughout the compressor. These measurements have also been used to calculate streamwise vorticity. Time-resolved static pressure measurements have been collected over the rotor tips for all rotors with each of the three tip clearance configurations for up to five loading conditions along the 100% corrected speedline using fast-response piezoresistive pressure sensors. These time-resolved static pressure measurements, as well as the time-resolved total pressures and velocities have helped to reveal a profound influence of the upstream stator vane on the size and shape of the rotor tip leakage flow. Finally, a novel particle image velocimetry (PIV) technique has been developed as a proof-of- concept. In contrast to PIV methods that have been typically been utilized for turbomachinery applications in the past, the method used for this study introduced the laser light through the same access window that was also used to image the flow. This new method addresses potential concerns related to the intrusive laser-introducing techniques that have typically been utilized by other authors in the past. Ultimately, the data collected for this project represent a unique data set which contributes to build a better understanding of the tip leakage flow field and its associated loss mechanisms. These data will facilitate future engine design goals leading to small blade heights in the rear stages of high pressure compressors and aid in the development of new blade designs which are desensitized to the performance penalties attributed to rotor tip leakage flows.

Berdanier, Reid A.↗
NASA NTRS2021

On the Cycling of 231 Pa and 230 Th in Benthic Nepheloid Layers

The naturally-occurring radionuclides protactinium-231 ( 231 Pa) and thorium-230 ( 230 Th) are produced at approximately uniform rates in the ocean and thought to be removed from the water column through a reversible exchange with settling particles. Recent measurements along the U.S. GEOTRACES North Atlantic transect (GA03) revealed two features which are at odds with current understanding about 231 Pa and 230 Th cycling in the ocean: (i) a sharp decrease in dissolved 231 Pa ( 231 Pa d ) and 230 Th ( 230 Th d ) activities with depth below 2000-4000 m and (ii) very high particulate 231 Pa ( 231 Pa p ) and 230 Th ( 230 Th p ) activities near the bottom, at a number of stations between the New England continental shelf and Bermuda. Concomitant measurements of light attenuation from beam transmissometry showed that both features occur in benthic nepheloid layers (BNLs), which suggests that these features may stem, at least partly, from the presence of resuspended sediment in the deep water column. Here we explore the behaviour of 231 Pa and 230 Th in BNLs by using (i) radionuclide, optical, and hydrographic data from the western segment of GA03 (west of Bermuda) and (ii) a simplified model of particle and radionuclide cycling that includes a lateral particle source. First, the BNLs observed at GA03 stations are characterized from measurements of the beam attenuation coefficient converted to particle concentrations. At all stations, particle concentrations below the clear water minimum were the highest in the bottom mixed layer, whose thickness ranged from 95 to 320 m, and decreased generally with height above the bottom. The thickness of strong BNLs varied from 482 to 1358 m and the vertical integral of particle concentration in excess to that at the clear water minimum varied from 1 x 10 4 to 2 x 10 6 mg m −2 , among different stations. Second, the particle-radionuclide model is fitted to data from stations GT11-04 (New England continental rise) and GT11-08 (Hatteras abyssal plain), where samples for radionuclide analyses were collected in the BNL. The model can reproduce simultaneously the increase of particle concentration with depth, the low 231 Pa d and 230 Th d in the BNLs, and the high 231 Pa p and 230 Th p near the bottom. According to the model, at heights less than about 300 m above the seafloor, the dissolved phase was set primarily by a balance between adsorption and desorption, with vertical turbulent mixing playing a secondary role, whilst the particulate phase behaved largely as a non-reactive constituent supplied laterally and transported vertically by particle settling and turbulent mixing. Sensitivity tests with the model suggest that lateral particle sources near continental slopes and similar reliefs can produce significant biases both in the 230 Th normalization method and in the interpretation of sediment 231 Pa/ 230 Th records. Our findings yield insights into the influence of sediment resuspension and transport on 231 Pa and 230 Th in the deep ocean and highlight the need for considering these processes in paleoceanographic applications.

Nepheloid layer↗
NASA NTRSDate not supplied

Optical Spectroscopy of Venus Aerosol Analogs and Substrate Survivability

Aerosol Rapid Analysis Combined Entry Probe/Sonde Technology (AERACEPT) is a dual descent probe and in-situ sampler of cloud and aerosol particles, used as part of the Nephele mission concept to study Venus cloud layers. Various optical spectroscopic characterization methods, including Raman spectroscopy and Surface Enhanced Raman Spectroscopy (SERS), Laser induced Breakdown Spectroscopy (LIBS), and absorbance and fluorescence spectroscopy are being assessed for viability for compatibility with AERACEPT. Specifically, the detection limits, sensitivity, specificity, and analysis cadence of these techniques for a set of Venus aerosol analogs and their mixtures are tested. We report on preliminary optical spectroscopy results of the Venus aerosol analogs, as well as optical substrate survivability in heated and acidic conditions, substrate physical robustness, and analyte wettability on substrates.

Venus↗
NASA NTRS1999

Monte Carlo mixture model of lifetime cancer incidence risk from radiation exposure on shuttle and international space station

Estimating uncertainty in lifetime cancer risk for human exposure to space radiation is a unique challenge. Conventional risk assessment with low-linear-energy-transfer (LET)-based risk from Japanese atomic bomb survivor studies may be inappropriate for relativistic protons and nuclei in space due to track structure effects. This paper develops a Monte Carlo mixture model (MCMM) for transferring additive, National Institutes of Health multiplicative, and multiplicative excess cancer incidence risks based on Japanese atomic bomb survivor data to determine excess incidence risk for various US astronaut exposure profiles. The MCMM serves as an anchor point for future risk projection methods involving biophysical models of DNA damage from space radiation. Lifetime incidence risks of radiation-induced cancer for the MCMM based on low-LET Japanese data for nonleukemia (all cancers except leukemia) were 2.77 (90% confidence limit, 0.75-11.34) for males exposed to 1 Sv at age 45 and 2.20 (90% confidence limit, 0.59-10.12) for males exposed at age 55. For females, mixture model risks for nonleukemia exposed separately to 1 Sv at ages of 45 and 55 were 2.98 (90% confidence limit, 0.90-11.70) and 2.44 (90% confidence limit, 0.70-10.30), respectively. Risks for high-LET 200 MeV protons (LET=0.45 keV/micrometer), 1 MeV alpha-particles (LET=100 keV/micrometer), and 600 MeV iron particles (LET=180 keV/micrometer) were scored on a per particle basis by determining the particle fluence required for an average of one particle per cell nucleus of area 100 micrometer(2). Lifetime risk per proton was 2.68x10(-2)% (90% confidence limit, 0.79x10(-3)%-0. 514x10(-2)%). For alpha-particles, lifetime risk was 14.2% (90% confidence limit, 2.5%-31.2%). Conversely, lifetime risk per iron particle was 23.7% (90% confidence limit, 4.5%-53.0%). Uncertainty in the DDREF for high-LET particles may be less than that for low-LET radiation because typically there is very little dose-rate dependence. Probability density functions for high-LET radiation quality and dose-rate may be preferable to conventional risk assessment approaches. Nuclear reactions and track structure effects in tissue may not be properly estimated by existing data using in vitro models for estimating RBEs. The method used here is being extended to estimate uncertainty in spacecraft shielding effectiveness in various space radiation environments.

NASA Center LaRC↗
NASA NTRS1971

A method for attitude control of a satellite in circular orbit by actively varying its inertia.

Equations are developed for a vehicle consisting of a system of particles which move in a generally prescribed fashion. These equations are then specialized to apply to a rigid vehicle with a constant amount of fluid contained in tanks and flowing through interconnecting lines. Various simplifications and specialized forms of the equations are derived, and relationships between mass and inertia that result from the fluid transfer logic and the physical constraints of the configuration are formulated. A tank arrangement and fluid transfer logic are devised which permit each vehicle product of inertia to be varied without varying the remaining elements of the inertia matrix or shifting the center of mass. Vehicle control is achieved by sensing angular rate and the local vertical and by transferring fluid as a function of these variables.

Clarke, W. D.↗
NASA NTRS1974

Geodesic synchrotron radiation in the Kerr geometry by the method of asymptotically factorized Green's functions

The scalar, electromagnetic, and gravitational geodesic-synchrotron-radiation (GSR) spectra are determined for the case of a test particle moving on a highly relativistic circular orbit about a rotating (Kerr) black hole. It is found that the spectral shape depends only weakly on the value of the angular-momentum parameter (a/M) of the black hole, but the total radiated power drops unexpectedly for a value of at least 0.95 and vanishes as the value approaches unity. A spin-dependent factor (involving the inner product of the polarization of a radiated quantum with the source) is isolated to explain the dependence of the spectral shape on the spin of the radiated field. Although the scalar wave equation is solved by separation of variables, this procedure is avoided for the vector and tensor cases by postulating a sum-over-states expansion for the Green's function similar to that found to hold in the scalar case. The terms in this sum, significant for GSR, can then be evaluated in the geometric-optics approximation without requiring the use of vector or tensor spherical harmonics.

Chrzanowski, P. L.↗
NASA NTRS2004

Development of a Method for Local Electron Temperature and Density Measurements in the Divertor of the JET Tokamak

Plasma volume recombination in the divertor, a process in which charged particles recombine to neutral atoms, contributes to plasma detachment and hence cooling at the divertor target region. Detachment has been observed at JET and other tokamaks and is known to occur at low electron temperatures (T(sub e)<1 eV) and at high electron density (n(sub e)>10(exp 20)/m(exp 3)). The ability to measure such low temperatures is therefore of interest for modelling the divertor. In present work we report development of a new spectroscopic technique for investigation of local electron density (n(sub e)) and temperature (T,) in the outer divertor at JET.

Jupen, C.↗
NASA NTRS1989

Diamondlike carbon protective coatings for optical windows

Diamondlike carbon (DLC) films were deposited on infrared transmitting optical windows and were evaluated as protective coatings for these windows exposed to particle and rain erosion. The DLC films were deposited on zinc selenide (ZnSe) and zinc sulfide (ZnS) by three different ion beam methods: (1) sputter deposition from a carbon target using an 8-cm argon ion source; (2) direct deposition by a 30-cm hollow cathode ion source with hydrocarbon gas in argon; and (3) dual beam direct deposition by the 30-cm hollow cathode ion source and an 8-cm argon ion source. In an attempt to improve the adherence of the DLC films on ZnSc and ZnS, ion beam cleaning, ion implantation with helium and neon ions, or sputter deposition of a thin, ion beam intermediate coating was employed prior to deposition of the DLC film. The protection that the DLC films afforded the windows from particle and rain erosion was evaluated, along with the hydrogen content, adherence, intrinsic stress, and infrared transmittance of the films. Because of the elevated stress levels in the ion beam sputtered DLC films and in those ion beam deposited with butane, films thicker than 0.1 micron and with good adherence on ZnS and ZnSe could not be generated. An intermediate coating of germanium successfully allowed the DLC films to remain adherent to the optical windows and caused only negligible reduction in the specular transmittance of the ZnS and ZnSe at 10 microns.

Swec, Diane M.↗
NASA NTRS1989

Diamondlike carbon protective coatings for optical windows

Diamondlike carbon (DLC) films were deposited on infrared transmitting optical windows and were evaluated as protective coatings for these windows exposed to particle and rain erosion. The DLC films were deposited on zinc selenide (ZnSe) and zinc sulfide (ZnS) by three different ion beam methods: (1) sputter deposition from a carbon target using an 8-cm argon ion source; (2) direct deposition by a 30-cm hollow cathode ion source with hydrocarbon gas in argon; and (3) dual beam direct deposition by the 30-cm hollow cathode ion source and an 8-cm argon ion source. In an attempt to improve the adherence of the DLC films on ZnSc and ZnS, ion beam cleaning, ion implantation with helium and neon ions, or sputter deposition of a thin, ion beam intermediate coating was employed prior to deposition of the DLC film. The protection that the DLC films afforded the windows from particle and rain erosion was evaluated, along with the hydrogen content, adherence, intrinsic stress, and infrared transmittance of the films. Because of the elevated stress levels in the ion beam sputtered DLC films and in those ion beam deposited with butane, films thicker than 0.1 micron and with good adherence on ZnS and ZnSe could not be generated. An intermediate coating of germanium successfully allowed the DLC films to remain adherent to the optical windows and caused only negligible reduction in the specular transmittance of the ZnS and ZnSe at 10 microns.

Swec, Diane M.↗
NASA NTRS1992

Dual-polarization characteristics of the radar ocean return in the presence of rain

Experimental data are presented on the polarimetric and dual-wavelength characteristics of the ocean surface in the presence of rain. To explain a portion of the variability observed in scatter plots under rain conditions, a storm model is used that incorporates measured drop size distributions. The fairly large variability indicates that effects of drop size distribution and the presence of partially melted particles can introduce a significant error in the estimate of attenuation. This effect is especially significant in the case of a 10-GHz radar under high rain rates. A surface reference method at this frequency will tend to overestimate the rain attenuation unless melting layer attenuation is properly taken into account. Observations of the cross-polarization return in stratiform rain over an ocean surface show three distinct components. Two of these correspond to aspherical, nonaligned particles in the melting layer seen in the direct and mirror-image returns. The remaining part depends both on the off-nadir depolarization by the surface and on the rain medium. A possible mechanism for this latter effect is the bistatic scattering from the rain to the surface.

Meneghini, R.↗
NASA NTRS1972

Solar system applications of Mie theory and of radiative transfer of polarized light

A theory of the multiple scattering of polarized light is discussed using the doubling method of van de Hulst. The concept of the Stokes parameters is derived and used to develop the form of the scattering phase matrix of a single particle. The diffuse reflection and transmission matrices of a single scattering plane parallel atmosphere are expressed as a function of the phase matrix, and the symmetry properties of these matrices are examined. Four matrices are required to describe scattering and transmission. The scattering matrix that results from the addition of two identical layers is derived. Using the doubling method, the scattering and transmission matrices of layers of arbitrary optical thickness can be derived. The doubling equations are then rewritten in terms of their Fourier components. Computation time is reduced since each Fourier component doubles independently. Computation time is also reduced through the use of symmetry properties.

Whitehill, L. P.↗
NASA NTRS1988

Growth of diamond by RF plasma-assisted chemical vapor deposition

A system has been designed and constructed to produce diamond particles by inductively coupled radio-frequency, plasma-assisted chemical vapor deposition. This is a low-pressure, low-temperature process used in an attempt to deposit diamond on substrates of glass, quartz, silicon, nickel, and boron nitride. Several deposition parameters have been varied including substrate temperature, gas concentration, gas pressure, total gas flow rate, RF input power, and deposition time. Analytical methods employed to determine composition and structure of the deposits include scanning electron microscopy, absorption spectroscopy, scanning Auger microprobe spectroscopy, and Raman spectroscopy. Analysis indicates that particles having a thin graphite surface, as well as diamond particles with no surface coatings, have been deposited. Deposits on quartz have exhibited optical bandgaps as high as 4.5 eV. Scanning electron microscopy analysis shows that particles are deposited on a pedestal which Auger spectroscopy indicates to be graphite. This is a phenomenon that has not been previously reported in the literature.

Meyer, Duane E.↗
NASA NTRS1991

Annular denuder systems: Precision and applicability to stratospheric sampling

Annular denuder systems (ADS) are designed to sample both gases and particles from ambient air without cross interference. The ADS's were used to detect chemical gradients over rural forest canopies. Accurate and precise methods for the sampling and analysis of low concentrations of gaseous and particulate species over short sampling periods were developed. Typical analytical precision data for nitrate and sulfate are shown. Both sampling and analytical precision data for a 2 hour test using four collocated ADS's are presented. On this particular day, the overall precision for HNO3 was unusually poor. The overall precision for the major constituents of interest is generally about 3 to 5 percent for typical ambient concentrations. The ADS's could potentially be used to determine the concentrations of HNO3 and fine-particle nitrate in the stratosphere, and to simultaneously characterize the chemical composition of stratospheric aerosols. Assuming 50 mb pressure, a 2 hour sampling time, a 50 Lpm sampling rate, and a low blank value, an overal precision of about 15 percent should be possible for 1 ppbv of HNO3 or nitrate, which would include nitric-acid trihydrate (NAT); and pump-driven systems could accommodate more complex systems involving denuders and filter packs in sequence. Alternate coatings could possibly be developed to selectively sample species such as Cl2 and ClNO3. The ADS's might also be useful in laboratory studies of heterogeneous chemical reactions involving NAT aerosols.

Randtke, Stephen J.↗
NASA NTRS2005

Effective Radius of Ice Cloud Particle Populations Derived from Aircraft Probes

The effective radius(r(sub e)) is a crucial variable in representing the radiative properties of cloud layers in general circulation models. This parameter is proportional to the condensed water content (CWC) divided by the extinction (sigma). For ice cloud layers, parameterizations for r(sub e), have been developed from aircraft in-situ measurements 1) indirectly, using data obtained from particle spectrometer probes and assumptions or observations about particle shape and mass to get the ice water content (IWC) and area to get sigma, and recently 2) from probes that measure IWC and sigma directly. This study compares [IWC/sigma] derived from the two methods using data sets acquired from comparable instruments on two aircraft, one sampling clouds at mid-levels and the other at upper-levels during the CRYSTAL-FACE field program in Florida in 2002. The sigma and IWC derived by each method are compared and evaluated in different ways for each aircraft data set. Direct measurements of sigma exceed those derived indirectly by a factor of two to two and a half. The IWC probes, relying on ice sublimation, appear to measure accurately except when the IWC is high or the particles too large to sublimate completely during the short transit time through the probe. The IWC estimated from the particle probes are accurate when direct measurements are available to provide constraints and useful information in high IWC/large particle situations. Because of the discrepancy in sigma estimates between the direct and indirect approaches, there is a factor of 2 to 3 difference in [IWC/sigma] between them. Although there are significant uncertainties involved in its use, comparisons with several independent data sources suggest that the indirect method is the more accurate of the two approaches. However, experiments are needed to resolve the source of the discrepancy in sigma.

Heymsfield, Andrew J.↗
NASA NTRS1996

Numerical Modeling of Spray Combustion with an Unstructured-Grid Method

The present unstructured-grid method follows strictly the basic finite volume forms of the conservation laws of the governing equations for the entire flow domain. High-order spatially accurate formulation has been employed for the numerical solutions of the Navier-Stokes equations. A two-equation k-epsilon turbulence model is also incorporated in the unstructured-grid solver. The convergence of the resulted linear algebraic equation is accelerated with preconditioned Conjugate Gradient method. A statistical spray combustion model has been incorporated into the present unstructured-grid solver. In this model, spray is represented by discrete particles, rather than by continuous distributions. A finite number of computational particles are used to predict a sample of total population of particles. Particle trajectories are integrated using their momentum and motion equations and particles exchange mass, momentum and energy with the gas within the computational cell in which they are located. The interaction calculations are performed simultaneously and eliminate global iteration for the two-phase momentum exchange. A transient spray flame in a high pressure combustion chamber is predicted and then the solution of liquid-fuel combusting flow with a rotating cup atomizer is presented and compared with the experimental data. The major conclusion of this investigation is that the unstructured-grid method can be employed to study very complicated flow fields of turbulent spray combustion. Grid adaptation can be easily achieved in any flow domain such as droplet evaporation and combustion zone. Future applications of the present model can be found in the full three-dimensional study of flow fields of gas turbine and liquid propulsion engine combustion chambers with multi-injectors.

Shang, H. M.↗
NASA NTRS1998

Gas-Kinetic Theory Based Flux Splitting Method for Ideal Magnetohydrodynamics

A gas-kinetic solver is developed for the ideal magnetohydrodynamics (MHD) equations. The new scheme is based on the direct splitting of the flux function of the MHD equations with the inclusion of "particle" collisions in the transport process. Consequently, the artificial dissipation in the new scheme is much reduced in comparison with the MHD Flux Vector Splitting Scheme. At the same time, the new scheme is compared with the well-developed Roe-type MHD solver. It is concluded that the kinetic MHD scheme is more robust and efficient than the Roe- type method, and the accuracy is competitive. In this paper the general principle of splitting the macroscopic flux function based on the gas-kinetic theory is presented. The flux construction strategy may shed some light on the possible modification of AUSM- and CUSP-type schemes for the compressible Euler equations, as well as to the development of new schemes for a non-strictly hyperbolic system.

Xu, Kun↗