Diffusion of low energy auroral electrons in the atmosphere.
Diffusion rate calculation of low energy auroral electrons in polar atmosphere by Monte Carlo methods with attention to backscattering coefficient
SEARCH · Engineering Papers
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.
Diffusion rate calculation of low energy auroral electrons in polar atmosphere by Monte Carlo methods with attention to backscattering coefficient
Selected properties of 57Ni-40Ti-3Hf (wt%) were measured to further develop this material for use in aerospace rolling-element bearings. The average grain size and microindentation hardness of the material at room temperature were 90 μm and 549 HV, respectively. From room temperature to 500 °C, electrical resistivity increased from 8.9×10–7 to 11.3×10–7 Ω⋅m, thermal expansion coefficients varied from 1.1×10–5 to 1.2×10–5 mm/mm⋅°C, specific heat increased from 0.44 to 0.53 W⋅sec/g⋅K, thermal diffusivity rose from approximately 0.003 to 0.005 cm2/sec and thermal conductivity increased from approximately 0.09 to 0.19 W/cm⋅K. The behavior of the material can be explained by well-known metallurgical theory. These results are expected to provide a basis for design and modeling of future aerospace components.
A new growth rate formula (NGRF) is developed for the rate of growth of cloud droplets by condensation. The theory used is a modification of the Lees-Shankar theory in which the two-stream Maxwellian distribution function of Lees is used in Maxwell's method of moments to determine the transport of water vapor to and heat away from the droplet. Boundary conditions at the droplet are the usual conditions set in terms of accommodation coefficients, and the solution passes smoothly into diffusion flow in the far region. Comparisons are given between NGRF and the conventional formula showing close agreement (approximately 0.1%) for large radii with significant difference (approximately 5%) for small radii (not greater than 1 micron). Growth times for haze droplets in a Laktionov chamber are computed.
Platform under growth region stabilizes melt for more uniform crystal growth. In new crucible, platform just below growth interface so melt is too shallow to support convection. Critical depth for onset of pertinent instability calculated from heat flux through surface of melt, volume coefficient of thermal expansion, thermal conductivity, thermal diffusivity, and kinematic viscosity.
A diffusion equation (FDT) is presented with a coefficient that reduces to the appropriate limiting form in the streaming and near thermodynamic limits for a moving fluid in which the dominant source of opacity is Thomson scattering. The present results are compared to those obtained with the corresponding equations for an absorptive medium. It is found that FDT for a scattering medium is accurate to better than less than about 17 percent over the range of optical depths of tau in the range of about 0 to 3.
The intensity, time-delay, and frequency content of radar echoes from the Magellan altimetry system are reduced to several parameters that are of great use in addressing many geological issues of the surface of Venus. These parameters include planetary radius, power reflection coefficient (reflectivity, both uncorrected and corrected for diffuse scattering), rms slope, and scattering functions (the behavior of backscatter as a function of incidence angle). Because the surface of Venus often reflects radio energy in unpredictable ways, models of radar scattering and their associated algorithms occasionally fail to accurately solve for the above surface parameters. Methods for identifying possible 'problem' altimetry data footprints, and techniques for resolving some key ambiguities are presented.
This paper presents a transient model of catalytic combustion suitable for isolated channels and monolith reactors. The model is a lumped two-phase (gas and solid) model where the gas phase is quasi-steady relative to the transient solid. Axial diffusion is neglected in the gas phase; lateral diffusion, however, is accounted for using transfer coefficients. The solid phase includes axial heat conduction and external heat loss due to convection and radiation. The combustion process utilizes detailed gas and surface reaction models. The gas-phase model becomes a system of stiff ordinary differential equations while the solid phase reduces, after discretization, into a system of stiff ordinary differential-algebraic equations. The time evolution of the system came from alternating integrations of the quasi-steady gas and transient solid. This work outlines the numerical model and presents some sensitivity studies on important parameters including internal transfer coefficients, catalytic surface site density, and external heat-loss (if applicable). The model is compared to two experiments using CO fuel: (1) steady-state conversion through an isothermal platinum (Pt) tube and (2) transient propagation of a catalytic reaction inside a small Pt tube. The model requires internal mass-transfer resistance to match the experiments at lower residence times. Under mass-transport limited conditions, the model reasonably predicted exit conversion using global mass-transfer coefficients. Near light-off, the model results did not match the experiment precisely even after adjustment of mass-transfer coefficients. Agreement improved for the first case after adjusting the surface kinetics such that the net rate of CO adsorption increased compared to O2. The CO / O2 surface mechanism came from a sub-set of reactions in a popular CH4 / O2 mechanism. For the second case, predictions improved for lean conditions with increased external heat loss or adjustment of the kinetics as in the first case. Finally, the results show that different initial surface-species distribution leads to different steady-states under certain conditions. These results demonstrate the utility of a lumped two-phase model of a transient catalytic combustor with detailed chemistry.
A next-generation docking system is being developed by the National Aeronautics and Space Administration (NASA) to support Constellation Space Exploration Missions to low Earth orbit (LEO), to the Moon, and to Mars. A number of investigations were carried out to quantify the properties of candidate elastomer materials for use in the main interface seal of the Low Impact Docking System (LIDS). This seal forms the gas pressure seal between two mating spacecraft. Three candidate silicone elastomer compounds were examined: Esterline ELA-SA-401, Parker Hannifin S0383-70, and Parker Hannifin S0899-50. All three materials were characterized as low-outgassing compounds, per ASTM E595, so as to minimize the contamination of optical and solar array systems. Important seal properties such as outgas levels, durometer, tensile strength, elongation to failure, glass transition temperature, permeability, compression set, Yeoh strain energy coefficients, coefficients of friction, coefficients of thermal expansion, thermal conductivity and diffusivity were measured and are reported herein.
The fluid-induced forces, both steady and unsteady, acting upon an impeller of a centrifugal pump, and impeller blade-diffuser vane interaction in centrifugal pumps with vaned radial diffusers were evaluated experimentally and theoretically. Knowledge of the steady and unsteady forces, and the associated rotordynamic coefficients are required to effectively model the rotor dynamics of the High Pressure Fuel Turbopump (HPFTP) of the Space Shuttle Main Engine (SSME). These forces and rotordynamic coefficients were investigated using different impellers in combination with volutes and vaned diffusers, and axial inducers. These rotor forces are global. Local forces and pressures are also important in impeller-diffuser interaction, for they may cause cavitation damage and even vane failures. Thus, in a separate investigation, impeller wake, and impeller blade and diffuser vane pressure measurements were made. The nature of the rotordynamic forces is discussed, the experimental facility is described, and the measurements of unsteady forces and pressure are reported together with a brief and incomplete attempt to calculate these flows.
Approximate solutions are sometimes used to determine the heat transfer and temperatures in a semitransparent material in which conduction and thermal radiation are acting. A comparison of the Milne-Eddington two-flux approximation and the diffusion approximation for combined conduction and radiation heat transfer in a ceramic material was preformed to determine the accuracy of the diffusion solution. A plane gray semitransparent layer without a substrate and a non-gray semitransparent plane layer on an opaque substrate were considered. For the plane gray layer the material is semitransparent for all wavelengths and the scattering and absorption coefficients do not vary with wavelength. For the non-gray plane layer the material is semitransparent with constant absorption and scattering coefficients up to a specified wavelength. At higher wavelengths the non-gray plane layer is assumed to be opaque. The layers are heated on one side and cooled on the other by diffuse radiation and convection. The scattering and absorption coefficients were varied. The error in the diffusion approximation compared to the Milne-Eddington two flux approximation was obtained as a function of scattering coefficient and absorption coefficient. The percent difference in interface temperatures and heat flux through the layer obtained using the Milne-Eddington two-flux and diffusion approximations are presented as a function of scattering coefficient and absorption coefficient. The largest errors occur for high scattering and low absorption except for the back surface temperature of the plane gray layer where the error is also larger at low scattering and low absorption. It is shown that the accuracy of the diffusion approximation can be improved for some scattering and absorption conditions if a reflectance obtained from a Kubelka-Munk type two flux theory is used instead of a reflection obtained from the Fresnel equation. The Kubelka-Munk reflectance accounts for surface reflection and radiation scattered back by internal scattering sites while the Fresnel reflection only accounts for surface reflections.
Electron and proton damage coefficients have been measured for several types of silicon cells of recent manufacture using 1 MeV electrons and protons with energies of 8.3, 17.5 and 40 MeV. LPE (AlGa)As-GaAs cells were also irradiated and pseudodamage coefficients derived under the assumption that the irradiation changes the diffusion length in only one of the n- or p-type layers. After irradiation the cells were annealed isochronally up to 450 C. The damage coefficient and annealing data for silicon cells are in substantial agreement with previous work. The GaAs cells have pseudo-damage coefficients which are dependent on fluence, and have about the same energy dependence as Si damage coefficients.
The aerodynamic characteristics of the augmentor wing concept with hypermixing primary nozzles were investigated. A large-scale semispan model in the Ames 40- by 80-Foot Wind Tunnel and Static Test Facility was used. The trailing edge, augmentor flap system occupied 65% of the span and consisted of two fixed pivot flaps. The nozzle system consisted of hypermixing, lobe primary nozzles, and BLC slot nozzles at the forward inlet, both sides and ends of the throat, and at the aft flap. The entire wing leading edge was fitted with a 10% chord slat and a blowing slot. Outboard of the flap was a blown aileron. The model was tested statically and at forward speed. Primary parameters and their ranges included angle of attack from -12 to 32 degrees, flap angles of 20, 30, 45, 60 and 70 degrees, and deflection and diffuser area ratios from 1.16 to 2.22. Thrust coefficients ranged from 0 to 2.73, while nozzle pressure ratios varied from 1.0 to 2.34. Reynolds number per foot varied from 0 to 1.4 million. Analysis of the data indicated a maximum static, gross augmentation of 1.53 at a flap angle of 45 degrees. Analysis also indicated that the configuration was an efficient powered lift device and that the net thrust was comparable with augmentor wings of similar static performance. Performance at forward speed was best at a diffuser area ratio of 1.37.
Angrite LEW 87051 consists of large olivine crystals set in a fine-grained groundmass that clearly represents a crystallized melt. A few olivines contain Ca-poor, Cr-rich cores that crystallized from a very different melt than the outer part of the crystals constituting the majority of olivine in LEW 87051. We evaluate a model in which the cores formed through fractional crystallization of one melt, then were incorporated into a different melt as xenocrysts, whereupon the original zoning patterns were modified by diffusion. Using a similar approach, we calculate zoning patterns for the cores that would result from perfect fractional crystallization, compare them with the observed zoning, and determine whether the differences could result from diffusive modification consistent with known diffusion rates for Ca, Mn, and Cr. Using distribution coefficients from the 1400 C, IW + 1 experiments, we computed CaO, Cr2O3, and MnO abundances in the hypothetical parent melt by inverting the olivine at the centers of the cores. We further assumed that the primary zoning profile for CaO is essentially unmodified, because the diffusion rate of Ca in olivine is slow. We carried out the fractional crystallization calculation until the calculated Ca content was that observed at break in zoning profiles at the outer edge of the core. We then normalized the distance of this calculated profile to the length of the observed profile in the olivine core and calculated profiles for MnO and Cr2O3. The CaO zoning profile agrees well with the observed profile. The observed MnO profile is slightly higher than the calculated profile near the edge. Diffusion calculations indicate that reversal of the general trend of primary zoning through diffusion would require that diffusion of Cr is 5-10x faster than that of Mn.
A theoretical model is presented which accounts for the dopant segregation in liquid-phase electroepitaxy in terms of dopant transport in the liquid phase (by electromigration and diffusion), the growth velocity, and the Peltier effect at the substrate-solution interface. The contribution of dopant electromigration to the magnitude of the effective segregation coefficient is dominant in the absence of convection; the contribution of the Peltier effect becomes significant only in the presence of pronounced convection. Quantitative expressions which relate the segregation coefficient to the growth parameters also permit the determination of the diffusion constant and electromigration mobility of the dopant in the liquid phase. The model was found to be in good agreement with the measured segregation characteristics of Sn in the electroepitaxial growth of GaAs from Ga-As solutions. For Sn in Ga-As solution at 900 C the diffusion constant was found to be 4 x 10 to the -5 sq cm/s and the electromigration velocity (toward the substrate with a positive polarity 2 x 10 to the -5 cm/s current density of 10 A/sq cm.
Wind tunnel diffuser performance is evaluated by comparing experimental data with analytical results predicted by an one-dimensional integration procedure with skin friction coefficient, a two-dimensional interactive boundary layer procedure for analyzing conical diffusers, and a two-dimensional, integral, compressible laminar and turbulent boundary layer code. Pressure, temperature, and velocity data for a 3.25 deg equivalent cone half-angle diffuser (37.3 in., 94.742 cm outlet diameter) was obtained from the one-tenth scale Altitude Wind Tunnel modeling program at the NASA Lewis Research Center. The comparison is performed at Mach numbers of 0.162 (Re = 3.097x19(6)), 0.326 (Re = 6.2737x19(6)), and 0.363 (Re = 7.0129x10(6)). The Reynolds numbers are all based on an inlet diffuser diameter of 32.4 in., 82.296 cm, and reasonable quantitative agreement was obtained between the experimental data and computational codes.
Turbulence measurements were obtained in the nearly separated flow in a 13 deg total angle of divergence conical diffuser coupled to a constant area tailpipe. Air at 207 newtons per square centimeter and 308 K provided an inlet velocity of about 51 meters per second at an inlet unit Reynolds number of 63.7 million per meter. Very high longitudinal turbulence intensities accompanied the diffusion process with peak values approaching 40 percent when normalized by the local centerline velocity. Predictions of the pressure recovery coefficient using a mixing length concept were good in the early stages of diffusion. In the latter stages of diffusion satisfactory predictions of the pressure recovery were obtained with an empirical method.
This work theoretically addressed two subjects: 1) onset of convection, 2) distribution of impurities. Onset of convection was considered analytically and numerically. Crystal growth was characterized by slow surface incorporation kinetics, i.e. growth kinetic coefficient beta (cm/s) small as compared to the typical bulk diffusion rate, D(sub 1)/h, where D(sub 1) is diffusivity of major crystallizing protein and h is the crystal size. Scaling type analysis predicted two laws on how the convection rate, v, essentially the Peclet number, Pe exactly equal to vh/D(sub 1), depends on dimensionless kinetic coefficient a exactly equal to beta h/D(sub 1). Namely: Pe = C(sub 2/5)(aRa(sup 2/5)) and Pe = C(sub 1) aRa. Here, Reynolds number Ra = rho(sub 1)(sup 0)gh(sup 3)(rho(sub p) - rho(sub w))/rho(sup p)rho(sub 1)vD(sub 1), v being solution viscosity. The constants C(sub 2/5), exactly equal to 0.28 and C(sub 1) exactly equal to 10(exp -2) found from the full scale computer simulation for a cylindrical crystal inside big cylindrical vessel. The linear boundary conditions connecting protein and impurity concentration at the interface with the flux to/from the interface was applied. No-slip condition for Navier-Shocker equations was employed. With these conditions, flow and concentration distributions were calculated. Validity of the Pe(Ra) dependencies follows for wide range of parameters for which numerical calculations have been accomplished and presented by various points.
The effect of diffusion on the saturation and absorption properties of molecular gases is considered. The analysis starts with a simple derivation of the multi-level rate equations. Incorporating the diffusion term, calculations are carried out both for the homogeneous and inhomogeneous broadening cases. Both graphical and analytical solutions are obtained for the absorption coefficient and the saturation intensity. It is shown that the effect of diffusion is characterized by a diffusion function, which varies inversely with the radius of the beam at low pressures and small beam sizes. The analysis is done in terms of two parameters k1a and k2a, which are simple functions of the intensity of the beam I, the effective relaxation rate w, and the effective diffusion constant D the a is designated as the radius of the beam.