Search NASA⌕ Search

SEARCH · Search NASA

Results for “diffusion 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 217 records · Page 12

Mathematical model of diffusion-limited evolution of multiple gas bubbles in tissue

Models of gas bubble dynamics employed in probabilistic analyses of decompression sickness incidence in man must be theoretically consistent and simple, if they are to yield useful results without requiring excessive computations. They are generally formulated in terms of ordinary differential equations that describe diffusion-limited gas exchange between a gas bubble and the extravascular tissue surrounding it. In our previous model (Ann. Biomed. Eng. 30: 232-246, 2002), we showed that with appropriate representation of sink pressures to account for gas loss or gain due to heterogeneous blood perfusion in the unstirred diffusion region around the bubble, diffusion-limited bubble growth in a tissue of finite volume can be simulated without postulating a boundary layer across which gas flux is discontinuous. However, interactions between two or more bubbles caused by competition for available gas cannot be considered in this model, because the diffusion region has a fixed volume with zero gas flux at its outer boundary. The present work extends the previous model to accommodate interactions among multiple bubbles by allowing the diffusion region volume of each bubble to vary during bubble evolution. For given decompression and tissue volume, bubble growth is sustained only if the bubble number density is below a certain maximum.

Evaluation Studies↗

Some basic mathematical methods of diffusion theory

An introductory treatment of the fundamentals of diffusion theory is presented, starting with molecular diffusion and leading up to the statistical methods of turbulent diffusion. A multilayer diffusion model, designed to permit concentration and dosage calculations downwind of toxic clouds from rocket vehicles, is described. The concepts and equations of diffusion are developed on an elementary level, with emphasis on atmospheric applications.

Giere, A. C.↗

Second order closure modeling of turbulent buoyant wall plumes

Non-intrusive measurements of scalar and momentum transport in turbulent wall plumes, using a combined technique of laser Doppler anemometry and laser-induced fluorescence, has shown some interesting features not present in the free jet or plumes. First, buoyancy-generation of turbulence is shown to be important throughout the flow field. Combined with low-Reynolds-number turbulence and near-wall effect, this may raise the anisotropic turbulence structure beyond the prediction of eddy-viscosity models. Second, the transverse scalar fluxes do not correspond only to the mean scalar gradients, as would be expected from gradient-diffusion modeling. Third, higher-order velocity-scalar correlations which describe turbulent transport phenomena could not be predicted using simple turbulence models. A second-order closure simulation of turbulent adiabatic wall plumes, taking into account the recent progress in scalar transport, near-wall effect and buoyancy, is reported in the current study to compare with the non-intrusive measurements. In spite of the small velocity scale of the wall plumes, the results showed that low-Reynolds-number correction is not critically important to predict the adiabatic cases tested and cannot be applied beyond the maximum velocity location. The mean and turbulent velocity profiles are very closely predicted by the second-order closure models. but the scalar field is less satisfactory, with the scalar fluctuation level underpredicted. Strong intermittency of the low-Reynolds-number flow field is suspected of these discrepancies. The trends in second- and third-order velocity-scalar correlations, which describe turbulent transport phenomena, are also predicted in general, with the cross-streamwise correlations better than the streamwise one. Buoyancy terms modeling the pressure-correlation are shown to improve the prediction slightly. The effects of equilibrium time-scale ratio and boundary condition are also discussed.

Zhu, Gang↗

A numerical analysis of high-temperature heat pipe startup from the frozen state

Continuum and rarefied vapor flows co-exist along the heat pipe length for most of the startup period. A two-region model is proposed in which the vapor flow in the continuum region is modeled by the compressible Navier-Stokes equations, and the vapor flow in the rarefied region is simulated by a self-diffusion model. The two vapor regions are linked with appropriate boundary conditions, and heat pipe wail, wick, and vapor flow are solved as a conjugate problem. The numerical solutions for the entire heat pipe startup process from the frozen state are compared with the corresponding experimental data with good agreement.

Cao, Y.↗

The ionosphere of Triton

A model of the atmospheric temperature structure and composition inferred from the Voyager (UVS) solar occultations was used together with a one-dimensional chemical diffusive model to interpret the Voyager Radio Science Spectrometer (RSS) ingress measurements of Triton's electron density. Although N2(+) is the major ion created, N(+) produced by dissociative ionization is the dominant ion. Reaction of thermospheric H2, produced by Lyman-alpha dissociation of CH4 in the lower atmosphere, is the major loss for N(+) ions and maintains these ions in PCSS below 600 km. Solar EUV ionization cannot generate electron densities at the magnitude measured by the RSS experiment and an additional ionization source about 3 x 10 to the 8th ions/sq cm per sec is required. The ionosphere may undergo a transition from PCSS to diffusive control if the N(+) ion production rates were greater than the H2 flux derived from CH4. In this case, the upward flowing H2 is totally converted to H by reaction with N(+) and the remaining N(+) ions recombine radiatively to create an ionosphere under diffusive control above the peak.

Majeed, T.↗

Investigation of the Saturn dust environment from the analysis of energetic charged particle measurements

In order to assist the Cassini project in evaluating risks of collisions with particulate matter in the rotational equatorial plane, available Pioneer 11 and Voyager energetic charged particle data were reinvestigated to constrain the column mass density of absorbing material within several radial ranges. Within the orbit of Mimas, CRAND proton phase space densities maximize near 2.67 R sub s and exhibit secondary maximum at 2.43 R sub s. From the condition that sources must exceed losses near these maxima and using available theoretical models for CRAND proton production rates, upper bounds on the column mass density at these two radial locations are calculated. Detailed fits of radial diffusion models to the observed flux maxima yield somewhat more restrictive upper limits. The upper limits compare to a lower limit on the column mass density, estimated from previous model calculations by Van Allen. Aside from continuous rings, longitudinally limited, low optical depth clouds of particulates may exist in orbit with several of the inner satellites including Mimas and Enceladus. A brief review of Voyager energetic particle microsignatures that suggest the presence of material co-orbiting with these two satellites is presented. Finally, Pioneer 11 and Voyager measurements of low energy electron fluxes exhibit minima near the location of the tenuous E Ring centered on approx. 4 R sub s. Pioneer 11 pitch angle distributions appear to support the possibility that direct absorption by Ring E particulates produced the observed flux decreases.

Hood, L. L.↗

Summer Work Experience: Determining Methane Combustion Mechanisms and Sub-Scale Diffuser Properties for Space Transporation System Engine Testing

To assess engine performance during the testing of Space Shuttle Main Engines (SSMEs), the design of an optimal altitude diffuser is studied for future Space Transportation Systems (STS). For other Space Transportation Systems, rocket propellant using kerosene is also studied. Methane and dodecane have similar reaction schemes as kerosene, and are used to simulate kerosene combustion processes at various temperatures. The equations for the methane combustion mechanism at high temperature are given, and engine combustion is simulated on the General Aerodynamic Simulation Program (GASP). The successful design of an altitude diffuser depends on the study of a sub-scaled diffuser model tested through two-dimensional (2-D) flow-techniques. Subroutines given calculate the static temperature and pressure at each Mach number within the diffuser flow. Implementing these subroutines into program code for the properties of 2-D compressible fluid flow determines all fluid characteristics, and will be used in the development of an optimal diffuser design.

Williams, Powtawche N.↗

Radiation and chemistry in the stratosphere - Sensitivity to O2 absorption cross sections in the Herzberg continuum

It is suggested that the discrepancies between observed and modeled vertical profiles of such halocarbons as CFCl3, as well as the problem of simultaneously fitting N2O, CH4, CF2Cl2 and CFCl3 profiles with a single eddy diffusion model, are due to an overestimation of the molecular oxygen absorption cross sections in the 200-220 nm spectral region. The replacement of current O2 cross sections in this range with values that are in better agreement with results for the compounds cited leads to N2O, CF2Cl2 and CFCl3 concentration reductions of factors 0.70, 0.62 and 0.19, respectively. Profiles of CH4, H2 and CO remain unchanged, and the predicted concentration of HNO3 above 30 km is reduced by about 50% for yet another improved fit with observations. It is noted that the correction proposed produces a 30% ozone increase near the 20-25 km peak.

Froidevaux, L.↗

Evaluation of transport conditions during physical vapor transport growth of opto-electronic crystals

Transport conditions were evaluated during the vapor phase growth of mercurous chloride crystals in a closed tube. Experimentally observed growth rates were much smaller than those calculated by the Hertz-Knudsen (H-K) equation. The Arrhenius behavior of growth rate with the temperature was used to derive the sticking coefficient. A one-dimensional diffusion model was used to calculate the total mass flux and was compared with the condensing flux. It was predicted that growth occurred in the convecto-diffusive range.

Singh, N. B.↗

Multiple Scattering of Laser Pulses in Snow Over Ice: Modeling the Potential Bias in ICESat Altimetry

The primary goal of NASA's current ICESat and future ICESat2 missions is to map the altitude of the Earth's land ice with high accuracy using laser altimetry technology, and to measure sea ice freeboard. Ice however is a highly transparent optical medium with variable scattering and absorption properties. Moreover, it is often covered by a layer of snow with varying depth and optical properties largely dependent on its age. We describe a modeling framework for estimating the potential altimetry bias caused by multiple scattering in the layered medium. We use both a Monte Carlo technique and an analytical diffusion model valid for optically thick media. Our preliminary numerical results are consistent with estimates of the multiple scattering delay from laboratory measurements using snow harvested in Greenland, namely, a few cm. Planned refinements of the models are described.

Davis, A. B.↗

A Study of High Altitude Hypersonic Flow-Field Radiation

This paper presents the results of the works carried out jointly at Stanford University and Ames Research Center under a grant from the Ballistic Missile Defense Organization (BMDO) (formerly the Strategic Defense Initiative Organization) to explain and understand the results of the two flight experiments, Bow Shock Ultra-Violet 1 and Bow Shock Ultra-Violet 2, carried out by the Organization. A portion of the material contained in this paper has been reported elsewhere in open literature. However, this paper provides (1) the details of scientific contents not available in those literature, (2) the links among those and the logical order of the efforts involved, and (3) some materials not contained in any open literature. The first author is responsible for execution of the work; the second author directed the work of the first author. In the two flight experiments mentioned above, the spectra of radiation in the ultraviolet wavelength range incident on the stagnation point of a blunt body were measured at the flight speeds of 3.8 and 5.2 km/sec over a wide range of altitudes. The results were compared first with the calculations made using the original version of the NEQAIR/STRAP codes written earlier by the second author. At low altitudes, the calculations agreed with the measurement. However, at high altitudes, the calculations underestimated the intensity of the radiation by several orders of magnitudes. A shock tube experiment was carried out at CALSPAN and a plasma-torch experiment was carried out at Stanford University to produce experimental data to help explain the discrepancy. In addition, the shock tube experiment at Ames Research Center carried out independently of the BMDO was also found to be relevant to this question. In this paper, several theoretical models are developed and calculations using the models were carried out to explain the results not only of the flight experiments but also of the CALSPAN, Stanford, and Ames experiments. The are: (1) the diffusion model for the rotational mode to explain the slowness of rotational excitation, (2) assignment of different vibrational temperatures and different relaxation rates for different molecules, and (3) the modification of the NEQAIR code to accommodate the new experimental data. This paper shows that the discrepancy between the flight data and calculation is smaller with the present model, but is still substantial.

Moreau, Stephane↗

Numerical modeling of the hybrid plume plasma rocket

This paper reports on ongoing research in the area of numerical simulations of hybrid plume dynamics, and describes the basic mechanisms for energy and particle transport. The solution approach assumes cylindrical geometry and includes a multiplicity of atomic reactions, as well as the presence of a strong magnetic field. The principal reactions are electron impact ionization, and resonant charge exchange between hot ions and cold neutrals. Radial particle and energy transport is mainly by diffusion. The present description includes a modified Bohm diffusion model for plasma in the core of the plume, and classical neutral particle diffusion in the cooler regions of the flow. Neutrals are allowed to 'free stream' in the low density regions, where the collision mean-free-path becomes comparable, or even larger than the characteristic dimensions of the system. A parameteric scan of the various regimes of interest has been completed and the most significant results are presented.

Chang, F. R.↗

Escarpment evolution on high-elevation rifted margins: Insights derived from a surface processes model that combines diffusion, advection and reaction

Experiments with a surface processes model of large-scale (1-1000 km) long-term (1-100 m.y.) erosional denudation are used to establish the controls on the evolution of a model escarpment that is related to the rifting of a continent. The mdoel describes changes in topographic form as a result of sumultaneous short- and long-range mass transport representing hillslope (diffusive) processes and fluvial transport (advection), repsectively. Fluvial entrainment is modeled as a first-order kinetic reaction which reflects the erodibility of the substrate, and therefore the fluvial system is not necessarily carrying at capacity. One dimensional and planform models demonstrate that the principal controls on the evolution of an initially steep model escarpment are (1) antecedent topography/drainage; (2) the timesale (or equivalently a length scale) in the fluvial entrainment reaction; (3) the flexural response of the lithosphere to denudation; and (4) the relative efficiencies of the short- and long-range transport processes. When rainfall and substrate lithology are uniform, a significant amount of discharge draining over the escarpment top causes it to degrade. Only when the top of the model escarpment coincides with a drainage divide can escarpment retreat occur for these conditions. An additional requirement for retreat of a model escarpment without decline is a long reaction time scale for fluvial entrainment. This corresponds to a substrate that is hard to detach by flucial erosion, and therefore to fluvial erosion that is not transport limited. Coninuous backtilting of an escarpment due ot flexural isostatic uplift in response to denudational unloading helps maintain the scarp top as a divide. It is essntial if the escarpment gradient is to be preserved during retreat in a uniform lithology. Low flexural rigidieties propote steep and slowly retreating escarpments. For given rainfall and substrate conditions, the morphology of a retraeating model escarpment is determined by the ratio of hte short-range diffusive and long range advective transport efficiencies. A low ration (which is interpreted to correspond to a more humid, temperate climate) produces a convex upper slope, and concave lower slope morphology and only major escarpments are predicted to preserve a high scarp gradient. Lithological contrasts in hte model produce more complex morphologies and predict the formation of scarps crowned by an erosionaly resistant caprock. However, resistant caprocks are not an essential requirement for model scarps to retreat. We conclude that the inferred controls and model behavior are both consistent with the present-day morphology of rifted continental margins and with modern conceptual models of landscape evolution.

Kooi, Henk↗

Effect of Gravity Level on the Particle Shape and Size During Zeolite Crystal Growth

A microscopic diffusion model is developed to represent solute transport in the boundary layer of a growing zeolite crystal. This model is used to describe the effect of gravity on particle shape and solute distribution. Particle dynamics and crystal growth kinetics serve as the boundary conditions of flow and convection-diffusion equations. A statistical rate theory is used to obtain the rate of solute transport across the growing interface, which is expressed in terms of concentration and velocity of solute species. Microgravity can significantly decrease the solute velocity across the growing interface compared to its earth-based counterpart. The extent of this reduction highly depends on solute diffusion constant in solution. Under gravity, the flow towards the crystal enhances solute transport rate across the growing interface while the flow away from crystals reduces this rate, suggesting a non-uniform growth rate and thus an elliptic final shape. However, microgravity can significantly reduce the influence of flow and obtain a final product with perfect spherical shape. The model predictions compare favorably with the data of space experiment of zeolites grown in space.

Song, Hong-Wei↗

Ion transport and loss in the Earth's quiet ring current. 2: Diffusion and magnetosphere-ionosphere coupling

We have studied the transport and loss of H(+), He(+), and He(++) ions in the Earth's quiet time ring current (1 to 300 keV/e, 3 to 7 R(sub E), Kp less than 2+, absolute value of Dst less than 11, 70 to 110 degs pitchangles, all LT) comparing the standard radial diffusion model developed for the higher-energy radiation belt particles with measurements of the lower energy ring current ions in a previous paper. Large deviations of that model, which fit only 50% of the data to within a factor of 10, suggested that another transport mechanism is operating in the ring current. Here we derive a modified diffusion coefficient corrected for electric field effects on ring current energy ions that fit nearly 80% of the data to within a factor of 2. Thus we infer that electric field fluctuations from the low-latitude to midlatitude ionosphere (ionospheric dynamo) dominated the ring current transport, rather than high-latitude or solar wind fluctuations. Much of the remaining deviation may arise from convective electric field transport of the E less than 30 keV particles. Since convection effects cannot be correctly treated with this azimuthally symmetric model, we defer treatment of the lowest-energy ions to a another paper. We give chi(exp 2) contours for the best fit, showing the dependence of the fit upon the internal/external spectral power of the predicted electric and magnetic field fluctuations.

Sheldon, R. B.↗

Voyager measurements of the isotopic composition of cosmic-ray aluminum and implications for the propagation of cosmic rays

We report a new measurement of the cosmic-ray isotopic composition of aluminum in the low-energy range form 75 to 206 MeV per nucleon.This measurement was made using the high-energy telescope of the CRS experiment on the Voyager 1 and 2 spacecraft during the time period from 1977 to 1993 with an average solar modulation level about 497 MV, roughly the same as at Earth near sunspot minimum. We obtain approximately 430 Al events of which approximately 35 are Al-26 and 395 are Al-27. The Al isotopes were separated with an average mass resolution sigma of 0.35 amu. Our interpretation of the isotopic composition of cosmic-ray aluminum is based on a standard Leaky-Box model for the interstellar propagation of cosmic-ray nuclei using the latest cross sections of the New Mexico-Saclay collaboration as well as a disk-halo diffusion model. From our observed ratio Al-26/Al-27 of 8.3 +/- 2.4 % we deduce an average interstellar density of about 0.52 (+0.26, -0.2) atoms per cu cm. This density is larger than the value of 0.28 (+0.14, -0.11) atoms per cu cm we found from an analysis of the observed abundance of the longer lived Be-10 made using data from the Voyager detectors over almost the same time interval and using essentially the same propagation program.

Lukasiak, A.↗

The Dissolution of an Interfween Miscible Liquids

The disappearance of the surface tension of the interface of a binary mixture, measured using the dynamic surface light scattering technique, is slower for a binary mixture of higher density contrast. A comparison with a naive diffusion model, expected to provide a lower limit for the speed of dissolution in the absence of gravity shows that the interfacial surface tension disappears much slower than even by diffusion with the effect becoming much more pronounced when density contrast between the liquid phases is increased. Thus, the factor most likely to be responsible for this anomalously slow dissolution is gravity. A mechanism could be based on the competition between diffusive relaxation and sedimentation at the dissolving interface.

Vlad, D.H.↗