SEARCH · Search NASA

Results for “discrete phase 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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 181 records · Page 10

All sources NASA NTRS DOE OSTI

NASA NTRS1975

MTF and point-spread function for a large-area CCD imager

The MTF degradation due to lateral diffusion is calculated for a back illuminated CCD imager for typical device parameters. The discrete nature of the CCD and finite size of the photosensitive elements result in an additional MTF degradation. The Fourier transform approach is utilized to calculate the effective point spread function for these processes in the time domain. Experimental data are presented on the point spread function for a three phase, double level anodized aluminum 160 x 100 thinned and back illuminated CCD imager and compared with the theoretical results. A simple modification of the Crowell and Labuda model suggested by these results is presented.

Ando, K. J.↗

NASA NTRS1989

New measurements of the far ultraviolet scattering properties of interstellar dust

An analysis is presented for the spectra of the diffuse ultraviolet background taken during shuttle flight STS-61C (January 1986). Eight regions of the sky were observed for approximately 20 minutes each, using a spectrograph designed specifically to perform measurements of the UV background. The field of view was 3.8 deg x 8 min, with imaging along the slit to confine stellar contamination. The instrument featured a shutter mechanism to measure internal background during flight, a low-scatter holographically ruled diffraction grating, photon counting microchannel plate detectors, through baffling, and a crystal window to further attenuate stray light. The spectra covered the range 1400 to 1850 A and was binned in 50 A bands. The procedure for substracting the contribution of stars too faint to be detected as discrete sources during the observations (in general this represents a small fraction of the total intensity detected except at the longest wavelengths). A radiative transfer model used to interpret the data and set confidence intervals on the relevant parameters is described. It was found that the continuum component of the diffuse ultraviolet background arises primarily from two sources. One source is scattering of starlight by interstellar dust with an albedo of about 12 percent and a relatively isotropic phase function. A second source consists of about 150 photons/cm/sec/ster/A of extragalactic light which is attenuated by the dust in the galaxy. Although emission features possibly associated with molecular H2 are detected in one look direction, fluorescence of H2 is not a major contributor to the diffuse UV background, at least at galactic latitudes greater than -10 degrees.

Hurwitz, Mark↗

NASA NTRS1994

The reflection for dense plant canopies from the one-angle radiative transfer equation

An essential component of remote sensing of vegetation canopies from satellites is fundamental understanding. Since passive remote is driven by photons, the modeling of photon interactions with vegetation is a basic building block in that understanding. Several such photon transport models have been developed during the past two decades and continue to be developed. Different approaches have been followed including monte carlo, radiosity methods, geometric shadowing, and radiative transfer. In general, each approach has application for canopies with specific attributes. This presentation concerns the application of radiative transfer to dense vegetation canopies in which the soil does not participate. The approach taken here is novel in that a consistent theory for photon transport for non-rotationally invariant leaf scattering is developed in a canopy with a general leaf angle distribution (LAD). The theory is limited to the one-angle approximation (azimuthally averaged radiance) and is based on Chandrasekhar's analytical theory. While such a model is admittedly only approximate, it does fulfill a unique function in our search for understanding. In particular, the model is simple in its construct yet contains the essential features of canopy architecture that are mainly responsible for observed responses. Thus, this model will not only be a predictive tool but also an educational one. The mathematical setting is the radiative transfer equation in a dense (semiinfinite) canopy. The leaf scattering phase function is assumed to be Lambertian with different reflectance and transmittance. In addition, abaxial and adaxial differentiation is allowed which effectively destroys optical reciprocity. The analytical solution for the canopy BRDF is obtained by manipulation of the integral transport equation (a la Chandrasekhar) for a general LAD. With discretization of the. leaf angle, the resulting set of integral equations are solved iteratively including an acceleration procedure when the single scatter albedo is near one (in the NIR). Results will be compared to the LARS soybean canopy radiances as well as to broadleaf results from a recent Ames experiment.

Ganapol, B. D.↗

NASA NTRS2013

The W-W02 Oxygen Fugacity Buffer at High Pressures and Temperatures: Implications for f02 Buffering and Metal-silicate Partitioning

Oxygen fugacity (fO2) controls multivalent phase equilibria and partitioning of redox-sensitive elements, and it is important to understand this thermodynamic parameter in experimental and natural systems. The coexistence of a metal and its oxide at equilibrium constitutes an oxygen buffer which can be used to control or calculate fO2 in high pressure experiments. Application of 1-bar buffers to high pressure conditions can lead to inaccuracies in fO2 calculations because of unconstrained pressure dependencies. Extending fO2 buffers to pressures and temperatures corresponding to the Earth's deep interior requires precise determinations of the difference in volume (Delta) V) between the buffer phases. Synchrotron x-ray diffraction data were obtained using diamond anvil cells (DAC) and a multi anvil press (MAP) to measure unit cell volumes of W and WO2 at pressures and temperatures up to 70 GPa and 2300 K. These data were fitted to Birch-Murnaghan 3rd-order thermal equations of state using a thermal pressure approach; parameters for W are KT = 306 GPa, KT' = 4.06, and αKT = 0.00417 GPa K-1. Two structural phase transitions were observed for WO2 at 4 and 32 GPa with structures in P21/c, Pnma and C2/c space groups. Equations of state were fitted for these phases over their respective pressure ranges yielding the parameters KT = 190, 213, 300 GPa, KT' = 4.24, 5.17, 4 (fixed), and αKT = 0.00506, 0.00419, 0.00467 GPa K-1 for the P21/c, Pnma and C2/c phases, respectively. The W-WO2 buffer (WWO) was extended to high pressure by inverting the W and WO2 equations of state to obtain phase volumes at discrete pressures (1-bar to 100 GPa, 1 GPa increments) along isotherms (300 to 3000K, 100 K increments). The slope of the absolute fO2 of the WWO buffer is positive with increasing temperature up to approximately 70 GPa and is negative above this pressure. The slope is positive along isotherms from 1000 to 3000K with increasing pressure up to at least 100 GPa. The WWO buffer is at a higher fO2 than the IW buffer at pressures lower than 40 GPa, and the magnitude of this difference decreases at higher pressures. This qualitatively indicates an increasingly lithophile character for W at higher pressures. The WWO buffer was quantitatively applied to W metal-silicate partitioning by using the WWO-IW buffer difference in combination with literature data on W metal-silicate partitioning to model the exchange coefficient (KD) for the Fe-W exchange reaction. This approach captures the pressure dependence of W metal-silicate partitioning using the WWO-IW buffer difference and models the activities of the components in the silicate and metallic phases using an expression of the Gibbs excess energy of mixing. Calculation of KD along a peridotite liquidus predicts a decrease in W siderophility at higher pressures that supports the qualitative behavior predicted by the WWO-IW buffer difference, and agrees with findings of others. Comparing the competing effects of temperature and pressure on W metal-silicate partitioning, our results indicate that pressure exerts a greater effect.

Shofner, G. A.↗

NASA NTRS2003

Software Estimates Costs of Testing Rocket Engines

Simulation-Based Cost Model (SiCM), a discrete event simulation developed in Extend , simulates pertinent aspects of the testing of rocket propulsion test articles for the purpose of estimating the costs of such testing during time intervals specified by its users. A user enters input data for control of simulations; information on the nature of, and activity in, a given testing project; and information on resources. Simulation objects are created on the basis of this input. Costs of the engineering-design, construction, and testing phases of a given project are estimated from numbers and labor rates of engineers and technicians employed in each phase, the duration of each phase; costs of materials used in each phase; and, for the testing phase, the rate of maintenance of the testing facility. The three main outputs of SiCM are (1) a curve, updated at each iteration of the simulation, that shows overall expenditures vs. time during the interval specified by the user; (2) a histogram of the total costs from all iterations of the simulation; and (3) table displaying means and variances of cumulative costs for each phase from all iterations. Other outputs include spending curves for each phase.

Smith, C. L.↗

NASA NTRS2012

Modifications Of Discrete Ordinate Method For Computations With High Scattering Anisotropy: Comparative Analysis

A numerical accuracy analysis of the radiative transfer equation (RTE) solution based on separation of the diffuse light field into anisotropic and smooth parts is presented. The analysis uses three different algorithms based on the discrete ordinate method (DOM). Two methods, DOMAS and DOM2+, that do not use the truncation of the phase function, are compared against the TMS-method. DOMAS and DOM2+ use the Small-Angle Modification of RTE and the single scattering term, respectively, as an anisotropic part. The TMS method uses Delta-M method for truncation of the phase function along with the single scattering correction. For reference, a standard discrete ordinate method, DOM, is also included in analysis. The obtained results for cases with high scattering anisotropy show that at low number of streams (16, 32) only DOMAS provides an accurate solution in the aureole area. Outside of the aureole, the convergence and accuracy of DOMAS, and TMS is found to be approximately similar: DOMAS was found more accurate in cases with coarse aerosol and liquid water cloud models, except low optical depth, while the TMS showed better results in case of ice cloud.

Korkin, Sergey V.↗

NASA NTRS1975

Research study on stabilization and control: Modern sampled-data control theory. Design of the large space telescope system

Conditions of self-sustained oscillations in a two-axis model of the nonlinear LST system are studied. The describing function of the CMG frictional nonlinearity of the LST system is used for the analysis, as well as continuous-data and discrete-data models of the simplified LST control system. A numerical-iterative method is described for the analysis of the two-axis system. Approximation methods and the direct plotting of the stability equation are implemented in the study. It is shown that although the dynamics of the two axes are identical, the amplitudes of self-sustained oscillations in the two axes may in principle be different. Analysis shows that the LST systems are of equal amplitudes but with 180-degree phase shift.

Kuo, B. C.↗

NASA NTRS2005

The Insoluble Carbonaceous Material of CM Chondrites as Possible Source of Discrete Organics During the Asteroidal Aqueous Phase

The larger portion of the organic carbon in carbonaceous chondrites (CC) is present as a complex and heterogeneous macromolecular material that is insoluble in acids and most solvents (IOM). So far, it has been analyzed only as a whole by microscopy (TEM) and spectroscopy (IR, NMR, EPR), which have offered and overview of its chemical nature, bonding, and functional group composition. Chemical or pyrolytic decomposition has also been used in combination with GC-MS to identify individual compounds released by these processes. Their value in the recognition of the original IOM structure resides in the ability to properly interpret the decomposition pathways for any given process. We report here a preliminary study of IOM from the Murray meteorite that combines both the analytical approaches described above, under conditions that would realistically model the IOM hydrothermal exposure in the meteorite parent body. The aim is to document the possible release of water and solvent soluble organics, determine possible changes in NMR spectral features, and ascertain, by extension, the effect of this loss on the frame of the IOM residue. Additional information is included in the original extended abstract.

Yabuta, H.↗

NASA NTRS2016

Multi-Disciplinary, Multi-Fidelity Discrete Data Transfer Using Degenerate Geometry Forms

In a typical multi-fidelity design process, different levels of geometric abstraction are used for different analysis methods, and transitioning from one phase of design to the next often requires a complete re-creation of the geometry. To maintain consistency between lower-order and higher-order analysis results, Vehicle Sketch Pad (OpenVSP) recently introduced the ability to generate and export several degenerate forms of the geometry, representing the type of abstraction required to perform low- to medium-order analysis for a range of aeronautical disciplines. In this research, the functionality of these degenerate models was extended, so that in addition to serving as repositories for the geometric information that is required as input to an analysis, the degenerate models can also store the results of that analysis mapped back onto the geometric nodes. At the same time, the results are also mapped indirectly onto the nodes of lower-order degenerate models using a process called aggregation, and onto higher-order models using a process called disaggregation. The mapped analysis results are available for use by any subsequent analysis in an integrated design and analysis process. A simple multi-fidelity analysis process for a single-aisle subsonic transport aircraft is used as an example case to demonstrate the value of the approach.

Olson, Erik D.↗

NASA NTRS2009

Magnetospheric Multiscale (MMS) Mission Commissioning Phase Orbit Determination Error Analysis

The Magnetospheric MultiScale (MMS) mission commissioning phase starts in a 185 km altitude x 12 Earth radii (RE) injection orbit and lasts until the Phase 1 mission orbits and orientation to the Earth-Sun li ne are achieved. During a limited time period in the early part of co mmissioning, five maneuvers are performed to raise the perigee radius to 1.2 R E, with a maneuver every other apogee. The current baseline is for the Goddard Space Flight Center Flight Dynamics Facility to p rovide MMS orbit determination support during the early commissioning phase using all available two-way range and Doppler tracking from bo th the Deep Space Network and Space Network. This paper summarizes th e results from a linear covariance analysis to determine the type and amount of tracking data required to accurately estimate the spacecraf t state, plan each perigee raising maneuver, and support thruster cal ibration during this phase. The primary focus of this study is the na vigation accuracy required to plan the first and the final perigee ra ising maneuvers. Absolute and relative position and velocity error hi stories are generated for all cases and summarized in terms of the ma ximum root-sum-square consider and measurement noise error contributi ons over the definitive and predictive arcs and at discrete times inc luding the maneuver planning and execution times. Details of the meth odology, orbital characteristics, maneuver timeline, error models, and error sensitivities are provided.

Chung, Lauren R.↗

NASA NTRS1992

Laboratory for Engineering Man/Machine Systems (LEMS): System identification, model reduction and deconvolution filtering using Fourier based modulating signals and high order statistics

Several important problems in the fields of signal processing and model identification, such as system structure identification, frequency response determination, high order model reduction, high resolution frequency analysis, deconvolution filtering, and etc. Each of these topics involves a wide range of applications and has received considerable attention. Using the Fourier based sinusoidal modulating signals, it is shown that a discrete autoregressive model can be constructed for the least squares identification of continuous systems. Some identification algorithms are presented for both SISO and MIMO systems frequency response determination using only transient data. Also, several new schemes for model reduction were developed. Based upon the complex sinusoidal modulating signals, a parametric least squares algorithm for high resolution frequency estimation is proposed. Numerical examples show that the proposed algorithm gives better performance than the usual. Also, the problem was studied of deconvolution and parameter identification of a general noncausal nonminimum phase ARMA system driven by non-Gaussian stationary random processes. Algorithms are introduced for inverse cumulant estimation, both in the frequency domain via the FFT algorithms and in the domain via the least squares algorithm.

Pan, Jianqiang↗

NASA NTRS1988

Thermal effects in two-phase flow through face seals

When liquid is sealed at high temperature, it flashes inside the seal due to pressure drop and/or viscous heat dissipation. Two-phase seals generally exhibit more erratic behavior than their single phase counterparts. Thermal effects, which are often neglected in single phase seal analyses, play an important role in determining seal behavior under two-phase operation. It is necessary to consider the heat generation due to viscous shear, conduction into the seal rings and convection with the leakage flow. Analytical models developed work reasonably well at the two extremes - for low leakage rates when convection is neglected and for higher leakage rates when conduction is neglected. A preliminary model, known as the Film Coefficient Model, is presented which considers conduction and convection both, and allows continuous boiling over an extended region unlike the previous low-leakage rate model which neglects convection and always forces a discrete boiling interface. Another simplified, semi-analytical model, based on the assumption of isothermal conditions along the seal interafce, has been developed for low leakage rates. The Film Coefficient Model may be used for more accurate and realistic description.

Basu, Prithwish↗

NASA NTRS1991

High-performance multiprocessor architecture for a 3-D lattice gas model

The lattice gas method has recently emerged as a promising discrete particle simulation method in areas such as fluid dynamics. We present a very high-performance scalable multiprocessor architecture, called ALGE, proposed for the simulation of a realistic 3-D lattice gas model, Henon's 24-bit FCHC isometric model. Each of these VLSI processors is as powerful as a CRAY-2 for this application. ALGE is scalable in the sense that it achieves linear speedup for both fixed and increasing problem sizes with more processors. The core computation of a lattice gas model consists of many repetitions of two alternating phases: particle collision and propagation. Functional decomposition by symmetry group and virtual move are the respective keys to efficient implementation of collision and propagation.

Lee, F.↗

NASA NTRS2005

Lunar and Planetary Science XXXVI, Part 22

The Lunar and Planetary Science XXXVI, Part 22 is presented. The topics include: 1) Pressure Histories from Thin and Thick Shock-induced Melt Veins in Meteorites; 2) Nano-structured Minerals as Signature of Microbial Activity; 3) The Insoluble Carbonaceous Material of CM Chondrites as Possible Source of Discrete Organics During the Asteroidal Aqueous Phase; 4) Discovery of Abundant Presolar Silicates in Subgroups of Antarctic Micrometeorites; 5) Characteristics of a Seismometer for the LUNAR-A Penetrator; 6) Heating Experiments of the HaH 262 Eucrite and Implication for the Metamorphic History of Highly Metamorphosed Eucrites; 7) Measurements of Ejecta Velocity Distribution by a High-Speed Video Camera; 8) Petrological Comparison of Mongolian Jalanash Ureilite and Twelve Antarctic Ureilites; 9) Metallographic Cooling Rate of IVA Irons Revisited; 10) Inhomogeneous Temperature Distribution in Chondrules in Shock-Wave Heating Model; 11) Subsurface Weathering of Rocks and Soils at Gusev Crater; 12) Extinct Radioactivities in the Early Solar System and the Mean Age of the Galaxy; 13) Correlation of Rock Spectra with Quantitative Morphologic Indices: Evidence for a Single Rock Type at the Mars Pathfinder Landing Site; 14) Silicon Isotopic Ratios of Presolar Grains from Supernovae; 15) Current Status and Readiness on In-Situ Exploration of Asteroid Surface by MINERVA Rover in Hayabusa Mission; 16) Long Formation Period of Single CAI: Combination of O and Mg Isotope Distribution; 17) Supra-Canonical Initial 26Al/27Al Indicate a 105 Year Residence Time for CAIs in the Solar Proto-Planetary Disk; 18) Evolution of Mercury's Obliquity; 19) First Results from the Huygens Surface Science Package; 20) Polyhedral Serpentine Grains in CM Chondrites; 21) Mountainous Units in the Martian Gusev Highland Region: Volcanic, Tectonic, or Impact Related? 22) Petrography of Lunar Meteorite MET 01210, A New Basaltic Regolith Breccia; 23) Earth-Moon Impacts at ~300 Ma and ~500 Ma Ago; 24) Chemical and Isotopic Analyses of Apollo 16 Glasses: An Integrated Approach; 25) Measuring the Surfaces on Mars; 26) Melting in Martian Snowbanks; 27) The Effect of Ni on Fe-FeS Phase Relations: Implications for the Chemistry of the Martian Core; 28) Shorelines in the Western United States as Analogs for Hypothesized Shoreline Features on Mars; 29) Al and Ti Isotopic Ratios of Presolar SiC Grains of Type Z; 30) Physicochemical Conditions of Clinopyroxene Crystallization in the Parental Bodies of Ordinary Chondrites; 31) Metamorphosed CM and CI Carbonaceous Chondrites Could be from the Breakup of the Same Earth-crossing Asteroid; 32) Aqueous Alteration and Hydrogen Generation on Parent Bodies of Unequilibrated Ordinary Chondrites: Thermodynamic Modeling for the Semarkona Composition; and 33) Preliminary Analysis of Relative Abundances of Hydrothermal Alteration Products in the C1-N10, Y6-N19, and Yax-1_863.51 Impact Melt Samples, Chicxulub Structure, Mexico

Source record↗

NASA NTRSDate not supplied

Demonstration of Capability to Simulate Particle Irregular Shape and Poly-Disperse Mixtures Within Lunar Lander Plume-Surface Interaction

Plume-Surface Interaction (PSI) between lander engine plumes and regolith soil creates hazards in obscuration and contamination by particle clouds, high-energy ejecta streams, and landing area cratering damage. The MSFC Fluid Dynamics Branch is developing simulation tools to offer a predictive PSI capability to NASA customers such as the Human Lander System (HLS) and Commercial Lunar Payload Services (CLPS). The Gas-Granular Flow Solver (GGFS) is the main application tool for coupled gas-particle two-phase flow simulations to predict the range of PSI effects from onset of surface erosion to deep crater formation. GGFS features an Eulerian-Eulerian modeling approach, treating both gas and granular material as interacting continuum phases. Modeling the lunar regolith granular material fluidic characteristics poses special challenges due to complex particle shapes and mixture composition. The lunar regolith is poorly sorted with broad particle size distributions and large fines content. It has significant cohesion, due to interlocking jagged particle shapes. Eulerian granular material flow modeling requires closure formulations for the granular material constitutive models (stress, friction, collisional and kinetic energy dissipation, drag, etc.). While closure models for spherical particles are available from particle kinetic theory, closure models for realistic non-spherical particles must be extracted from unit physics Discrete Element Model (DEM) particle interaction simulations and provided in the form of tabular datasets. The effects of particle irregular shape (non-spherical shape factors, angular particle surface roughness, and interlocking features) are simulated by approximating the particle features in the form of grouped elemental spheres to form composite particles in the DEM simulations. The effects of the wide range of regolith mixture particle sizes and the strong effects of the presence of the small particle sizes results in high cohesion and low porosity of the regolith mixture. The range of particle sizes is simulated by binning the particle sizes into an appropriate finite number of particle-size species and solving the problem as a species mixture. Combining these two modeling approaches enables simulations to capture both, the contributions of the irregular particle shape and the particle size distribution. The integration and maturation of the DEM-based constitutive model database generation process and poly-disperse mixture binning approach into the GGFS simulation framework are proceeding under funding by the NASA Game Changing Development program. The status of current capabilities will be presented in comparisons of crater characteristics resulting for spherical and irregular shape particles, and for mono-, bi-, and tri-disperse mixture simulations of Apollo LM plume-surface interaction. The computational results confirm the significance of including the particle shape and mixture effects. Going forward plans for the full implementation of the general poly-disperse regolith modeling capability and maturation towards NASA project application readiness under the GCD program will be presented.

Peter A Liever↗

NASA NTRSDate not supplied

Demonstration of Capability to Simulate Particle Irregular Shape and Poly-Disperse Mixtures Within Lunar Lander Plume-Surface Interaction

Peter A. Liever↗

NASA NTRS1989

A wall interference assessment/correction interface measurement system for the NASA/ARC 12-ft PWT

Development of complex air vehicle configurations is placing increasing demands on wind tunnel testing capabilities. A major area of concern is wall induced interference. Recent developments in wall interference technology provide a means for assessing and correcting for the wall induced interference using information contained in the distribution of flow variables measured at, or near, the wall. The restoration of the NASA-ARC 12-ft pressure wind tunnel (PWT) provides an opportunity to incorporate a measurement system with which wall interference assessment/correction (WIAC) technology can be applied. In this first phase of the development of a WIAC system for the PWT, the design criteria for the placement and the geometry of wall static pressure orifices were determined with a three step approach. First, the operational environment of the PWT was analyzed as to the requirements for the WIAC system. Second, appropriate wall interference theories were evaluated against the requirements determined from the operational environment. Third, the flow about representative models in the PWT was calculated and, specifically, the pressure signatures at the location of the test section wall were obtained. The number of discrete pressure measurements and their locations were determined by curve fitting the pressure distribution through the discrete measurements and evaluating the resulting error.

Source record↗

NASA NTRS2003

Finite-Rate Ablation Boundary Conditions for Carbon-Phenolic Heat-Shield

A formulation of finite-rate ablation surface boundary conditions, including oxidation, nitridation, and sublimation of carbonaceous material with pyrolysis gas injection, has been developed based on surface species mass conservation. These surface boundary conditions are discretized and integrated with a Navier-Stokes solver. This numerical procedure can predict aerothermal heating, chemical species concentration, and carbonaceous material ablation rate over the heatshield surface of re-entry space vehicles. In this study, the gas-gas and gas-surface interactions are established for air flow over a carbon-phenolic heatshield. Two finite-rate gas-surface interaction models are considered in the present study. The first model is based on the work of Park, and the second model includes the kinetics suggested by Zhluktov and Abe. Nineteen gas phase chemical reactions and four gas-surface interactions are considered in the present model. There is a total of fourteen gas phase chemical species, including five species for air and nine species for ablation products. Three test cases are studied in this paper. The first case is a graphite test model in the arc-jet stream; the second is a light weight Phenolic Impregnated Carbon Ablator at the Stardust re-entry peak heating conditions, and the third is a fully dense carbon-phenolic heatshield at the peak heating point of a proposed Mars Sample Return Earth Entry Vehicle. Predictions based on both finite-rate gas- surface interaction models are compared with those obtained using B' tables, which were created based on the chemical equilibrium assumption. Stagnation point convective heat fluxes predicted using Park's finite-rate model are far below those obtained from chemical equilibrium B' tables and Zhluktov's model. Recession predictions from Zhluktov's model are generally lower than those obtained from Park's model and chemical equilibrium B' tables. The effect of species mass diffusion on predicted ablation rate is also examined.

Chen, Y.-K.↗