Search NASA⌕ Search

SEARCH · Search NASA

Results for “Numerical diffusion”

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 487 records · Page 27

Bridgman crystal growth in low gravity - A scaling analysis

The results of an order-of-magnitude or scaling analysis are compared with those of numerical simulations of the effects of steady low gravity on compositional nonuniformity in crystals grown by the Bridgman-Stockbarger technique. In particular, the results are examined of numerical simulations of the effect of steady residual acceleration on the transport of solute in a gallium-doped germanium melt during directional solidification under low-gravity conditions. The results are interpreted in terms of the relevant dimensionless groups associated with the process, and scaling techniques are evaluated by comparing their predictions with the numerical results. It is demonstrated that, when convective transport is comparable with diffusive transport, some specific knowledge of the behavior of the system is required before scaling arguments can be used to make reasonable predictions.

Alexander, J. I. D.↗

A flux-limited treatment for the conductive evaporation of spherical interstellar gas clouds

In this work, we present and analyze a new analytic solution for the saturated (flux-limited) thermal evaporation of a spherical cloud. This work is distinguished from earlier analytic studies by allowing the thermal conductivity to change continuously from a diffusive to a saturated form, in a manner usually employed only in numerical calculations. This closed form solution will be of interest as a computational benchmark. Using our calculated temperature profiles and mass-loss rates, we model the thermal evaporation of such a cloud under typical interstellar medium (ISM) conditions, with some restrictions. We examine the ionization structure of the cloud-ISM interface and evaluate column densities of carbon, nitrogen, oxygen, neon, and silicon ions toward the cloud. In accord with other investigations, we find that ionization equilibrium is far from satisfied under the assumed conditions. Since the inclusion of saturation effects in the heat flux narrows the thermal interface relative to its classical structure, we also find that saturation effects tend to lower predicted column densities.

Dalton, William W.↗

Optical observations of NGC 2915: A nearby blue compact dwarf galaxy

This paper presents B and R band Charge Coupled Device (CCD) images and medium resolution spectroscopy of NGC 2915, a relatively isolated BCD (blue compact dwarf) galaxy at a distance of approximately 5 Mpc. NGC 2915 contains two stellar populations: a high surface brightness blue core population and a red diffuse population. The core population contains all of the H II, and numerous embedded objects. It is the locus of current high mass star formation. The brightest embedded objects are likely to be young ionizing clusters, while many of the fainter objects are likely to be individual supergiant stars with masses up to approximately 25 solar mass, or blends of a few such stars. Curious aligned structures on the SE side of the galaxy are seen and their nature discussed. The spectrum of the core is dominated by bright narrow emission lines like that of a high excitation and low metallicity (less than half solar) H II region. The continuum is flat, with Balmer and Ca II features seen in absorption. The velocity of the Ca II features suggest contamination by galactic interstellar absorption. There is a significant velocity gradient in the spectra, probably indicative of rotation. Outside of its core, NGC 2915 resembles a dE (dwarf elliptical) galaxy, in that it has an exponential surface brightness profile, is red ((B-R)(sub 0) = 1.65), and has a low extrapolated central surface brightness (B(0)(sub c) = 22.44). NGC 2915's properties are compared with other BCDs, concentrating on two morphologically similar BCDs that are near enough to resolve into stars: NGC 1705 and NGC 5253. It is noted that the presence of winds in BCDs invalidates closed box chemical evolution models and the remaining constraints on star formation duration are relatively weak. Some BCDs, including NGC 2915, may be able to maintain their present star formation rate for Gyr time scales. This suggests that the overall evolution of these BCDs may be much slower than the approximately 10 Myr burst time scales commonly quoted. However, shortly after the formation of a massive (10(exp 6) solar mass) cluster a BCD will have all the properties of strong starburst galaxy).

Meurer, G. R.↗

Vapor Growth of ZnSe-based Compound Semiconductors in Low Gravity environment on International Space Station

Six low-gravity experiments of crystal growth of ZnSe-based compound semiconductors by physical vapor transport (PVT) are scheduled to be processed inside Low Gradient Furnace (LGF) on Materials Science Research Rack (MSRR) in Material Science Laboratory (MSL) on International Space Station (ISS) starting from February 27th, 2023. The flight samples consist of three unseeded growths of ZnSe, Cr-doped ZnSe and ZnSeTe as well as growth of ZnSe on top of three single crystal ZnSe seeds with orientation of (100), (110) and (111). The main scientific objective of the flight experiments is to establish the effects of gravity-driven fluid flows on Earth to (i) the non-uniform incorporation of impurities and defects and (ii) the deviation from stoichiometry observed in the grown crystals caused by growth interface fluctuations and evolution as the results of buoyancy-driven convection and irregular fluid-flows. The investigation will also evaluate the effects of gravity on the PVT process by examining (i) the growth kinetics on various seed orientations, (ii) the dopant segregation and distribution in the Cr doped ZnSe, and (iii) the compositional segregation and distribution in the ternary ZnSeTe crystal grown by PVT. Additionally, from X-ray diffraction/topography and by measuring the density and distribution of various crystalline defects it will assess the self-induced strain effects developed during processing at elevated temperatures caused by the weight of the crystals. During the ground-based preparations, the following tasks have been accomplished: • Establish the thermodynamic properties of the ZnSe binary and related ternary systems by measuring the partial pressures of the vapor phase coexisting with the condensed phases and perform phase diagram calculations for comparison. • Evaluate the effects of gravity-driven convection in the growth process by performing experiments under various vapor transport orientations relative to gravity vector, i.e., horizontal, vertically stabilized and destabilized configurations. • Establish a quantitative correlation between growth parameters and the characteristics and properties of the grown crystals by characterizing the crystals using a variety of techniques. • Evaluate the fundamentals of the current theories on vapor transport and crystal growth kinetics by performing in-situ and real-time monitoring techniques during growth. • Develop a one-dimensional diffusion limited theory of mass transport and two- and three-dimensional numerical simulations of both mass transport and heat transfer and thus establish the microgravity requirements for the flight experiments and to predict the characteristics of the grown crystals such as the solid-vapor interface shapes… etc. and establish a fundamental understanding of the crystal growth process from the comparison with experimental results.

ZnSe↗

Examination of various turbulence models for application in liquid rocket thrust chambers

There is a large variety of turbulence models available. These models include direct numerical simulation, large eddy simulation, Reynolds stress/flux model, zero equation model, one equation model, two equation k-epsilon model, multiple-scale model, etc. Each turbulence model contains different physical assumptions and requirements. The natures of turbulence are randomness, irregularity, diffusivity and dissipation. The capabilities of the turbulence models, including physical strength, weakness, limitations, as well as numerical and computational considerations, are reviewed. Recommendations are made for the potential application of a turbulence model in thrust chamber and performance prediction programs. The full Reynolds stress model is recommended. In a workshop, specifically called for the assessment of turbulence models for applications in liquid rocket thrust chambers, most of the experts present were also in favor of the recommendation of the Reynolds stress model.

Hung, R. J.↗

Effect of gravity on methane-air combustion

Analytical and numerical techniques dealing with the theoretical description of the influence of zero and reduced gravitational acceleration on diffusion flames, with a view to improving understanding of fires in space vehicles, were developed in support of experimental work performed in this area. This was done in order to confirm qualitative understanding of the process, to determine the quantitative accuracy of numerical predictions, and to establish a mathematical model of the process for subsequent use as a predictive and exploratory tool. The following results were accomplished: (1) derivation of differential equations and boundary conditions describing the system, (2) details of the computations, using a FORTRAN computer program, for calculating the flow and heat and mass transfer in two dimensions (both steady and unsteady). It was shown that the experimental behavior can be reproduced with fair accuracy, provided that the time step is sufficiently short.

Elghobashi, S.↗

3-D Modeling of Double-Diffusive Convection During Directional Solidification of a Non-Dilute Alloy with Application to the HgCdTe Growth Under Microgravity Conditions

A numerical calculation for a non-dilute alloy solidification was performed using the FIDAP finite element code. For low growth velocities plane front solidification occurs. The location and the shape of the interface was determined using melting temperatures from the HgCdTe liquidus curve. The low thermal conductivity of the solid HgCdTe causes thermal short circuit through the ampoule walls, resulting in curved isotherms in the vicinity of the interface. Double-diffusive convection in the melt is caused by radial temperature gradients and by material density inversion with temperature. Cooling from below and the rejection at the solid-melt interface of the heavier HgTe-rich solute each tend to reduce convection. Because of these complicating factors dimensional rather then non-dimensional modeling was performed. Estimates of convection contributions for various gravity conditions was performed parametrically. For gravity levels higher then 1 0 -7 of earth's gravity it was found that the maximum convection velocity is extremely sensitive to gravity vector orientation and can be reduced at least by factor of 50% for precise orientation of the ampoule in the microgravity environment. The predicted interface shape is in agreement with one obtained experimentally by quenching. The results of 3-D modeling are compared with previous 2-D finding. A video film featuring melt convection will be presented.

Bune, Andris V.↗

Analytic Couple Modeling Introducing Device Design Factor, Fin Factor, Thermal Diffusivity Factor, and Inductance Factor

A set of convenient thermoelectric device solutions have been derived in order to capture a number of factors which are previously only resolved with numerical techniques. The concise conversion efficiency equations derived from governing equations provide intuitive and straight-forward design guidelines. These guidelines allow for better device design without requiring detailed numerical modeling. The analytical modeling accounts for factors such as i) variable temperature boundary conditions, ii) lateral heat transfer, iii) temperature variable material properties, and iv) transient operation. New dimensionless parameters, similar to the figure of merit, are introduced including the device design factor, fin factor, thermal diffusivity factor, and inductance factor. These new device factors allow for the straight-forward description of phenomenon generally only captured with numerical work otherwise. As an example a device design factor of 0.38, which accounts for thermal resistance of the hot and cold shoes, can be used to calculate a conversion efficiency of 2.28 while the ideal conversion efficiency based on figure of merit alone would be 6.15. Likewise an ideal couple with efficiency of 6.15 will be reduced to 5.33 when lateral heat is accounted for with a fin factor of 1.0.

Seebeck effect↗

Waves in space plasmas - The mirror trapping of hot auroral electrons

A brief review is given of the problem of precipitation of auroral electrons by electrostatic Bernstein waves. Since the magnetospheric loss cone is small, only moderately small intense levels of wave turbulence are required to remove any large anisotropy sources of free energy and to maintain a weakly anisotropic electron distribution on strong diffusion precipitation. The electrostatic electron cyclotron harmonic waves are nonconvectively unstable for weak loss cone anisotropies and over a large range of parameters for both the hot and cold distributions. Since the instability is nonconvective, weak wave growth can be maintained independent of the flux level of the hot electrons, i.e., the instability does not have the stably trapped flux limit imposed by convective amplification. Recent plasma numerical simulations show that the nonlinear evolution of this instability involves both the pitch angle diffusion of the hot electrons and the heating of the cold electrons.

Ashour-Abdalla, M.↗

Direct simulation of hypersonic flows over blunt slender bodies

Results of a numerical study of low-density hypersonic flow about cylindrically blunted wedges and spherically blunted cones with body half angles of 0, 5, and 10 deg are presented. Most of the transitional flow regime encountered during entry between the free molecule and continuum regimes is simulated for a reentry velocity of 7.5 km/s by including freestream conditions of 70 to 100 km. The bodies are at zero angle of incidence and have diffuse and finite catalytic surfaces. Translational, thermodynamic, and chemical nonequilibrium effects are considered in the numerical simulation by utilizing the direct simulation Monte Carlo (DSMC) method. The numerical simulations show that noncontinuum effects such as surface temperature jump, and velocity slip are evident for all cases considered. The onset of chemical dissociation occurs at a simulated altitude of 96 km for the two-dimensional configurations. Comparisons between the DSMC and continuum viscous shock-layer calculations highlight the significant difference in flowfield structure predicted by the two methods.

Moss, J. N.↗

International Conference on Numerical Methods in Fluid Dynamics, 11th, Williamsburg, VA, June 27-July 1, 1988, Proceedings

Recent advances in computational fluid dynamics (CFD) are discussed in reviews and reports. Topics addressed include CFD models in plasma dynamics, parallel computation for simulation studies, CFD for hypersonic airbreathing aircraft, multigrid methods for the steady incompressible Navier-Stokes equations, upwind differencing techniques, TV stable schemes for shock-interacting flows, Euler models of hypersonic vortex flows, parallel multilevel adaptive methods, and vortex methods for slightly viscous three-dimensional flows. Consideration is given to the accuracy of node-based solutions on irregular meshes, multigrid calculations for cascades, a finite-volume-element method for planar cavity flow, parallel heterogeneous mesh refinement for advection-diffusion equations, the convergence of the spectral-viscosity method for nonlinear conservation laws, and numerical simulations of Taylor vortices in a spherical gap.

Dwoyer, D. L.↗

NOR-USA Scientific Traverse of East Antarctica: Science and Logistics on a Three-Month Expedition Across Antarctica's Farthest Frontier

Dr. Albert's current research is centered on transfer processes in porous media, including air-snow exchange in the Polar Regions and in soils in temperate areas. Her research includes field measurements, laboratory experiments, and theoretical modeling. Mary conducts field and laboratory measurements of the physical properties of natural terrain surfaces, including permeability, microstructure, and thermal conductivity. Mary uses the measurements to examine the processes of diffusion and advection of heat, mass, and chemical transport through snow and other porous media. She has developed numerical models for investigation of a variety of problems, from interstitial transport to freezing of flowing liquids. These models include a two-dimensional finite element code for air flow with heat, water vapor, and chemical transport in porous media, several multidimensional codes for diffusive transfer, as well as a computational fluid dynamics code for analysis of turbulent water flow in moving-boundary phase change problems.

Albert, Mary R.↗

Implicit Space-Time Conservation Element and Solution Element Schemes

Artificial numerical dissipation is in important issue in large Reynolds number computations. In such computations, the artificial dissipation inherent in traditional numerical schemes can overwhelm the physical dissipation and yield inaccurate results on meshes of practical size. In the present work, the space-time conservation element and solution element method is used to construct new and accurate implicit numerical schemes such that artificial numerical dissipation will not overwhelm physical dissipation. Specifically, these schemes have the property that numerical dissipation vanishes when the physical viscosity goes to zero. These new schemes therefore accurately model the physical dissipation even when it is extremely small. The new schemes presented are two highly accurate implicit solvers for a convection-diffusion equation. The two schemes become identical in the pure convection case, and in the pure diffusion case. The implicit schemes are applicable over the whole Reynolds number range, from purely diffusive equations to convection-dominated equations with very small viscosity. The stability and consistency of the schemes are analysed, and some numerical results are presented. It is shown that, in the inviscid case, the new schemes become explicit and their amplification factors are identical to those of the Leapfrog scheme. On the other hand, in the pure diffusion case, their principal amplification factor becomes the amplification factor of the Crank-Nicolson scheme.

Chang, Sin-Chung↗

Transport of absolute angular momentum in quasi-axisymmetric equatorial jet streams

It is well known that prograde equatorial jet stresses cannot occur in an axisymmetric inviscid fluid, owing to the constraints of local angular momentum conservation. For a viscous fluid, the constraints of mass conservation prevent the formation of any local maximum of absolute angular momentum (m) without a means of transferring m against its gradient (delta m) in the meridional plane. The circumstances under which m can be diffused up-gradient by normal molecular viscosity are derived, and illustrated with reference to numerical simulations of axisymmetric flows in a cylindrical annulus. Viscosity is shown to act so as to tend to expel m from the interior outwards from the rotation axis. Such an effect can produce local super-rotation even in a mechanically isolated fluid. The tendency of viscosity to result in the expulsion of m is shown to be analogous in certain respects to a vorticity-mixing hypothesis for the effects of non-axisymmetric eddies of the zonally-averaged flow. It is shown how the advective and diffusive transport of m by non-axisymmetric eddies can be represented by the Transformed Eulerian Mean meridional circulation and the Eliassen-Palm (EP) flux of Andrews and McIntyre respectively, in the zonal mean. Constraints on the form and direction of the EP flux in an advective/diffusive flow for such eddies are derived, by analogy with similar constraints on the diffusive flux of m due to viscosity.

Read, P. L.↗

Pitch angle diffusion in the Jovian magnetodisc

Using a linear approximation to the field variation and a mathematical technique for distorting the orbit integration into the complex gyrophase plane, a formula is derived for the nonadiabatic change in magnetic moment. This change varies as the cosine of the particle gyrophase at the equator, so that if a particle traverses the equator several times with different values of this gyrophase, it suffers randomly correlated values for the nonadiabatic change in magnetic moment. The algebraic formula for this nonadiabatic change is found to be in excellent agreement with previous numerical computations. The phase space density, averaged over several bounce periods, is found to satisfy a diffusion equation in pitch angle. The diffusion coefficient is dependent on particle energy, pitch angle, and the field line along which the particles are moving.

Birmingham, T. J.↗

Thermal analysis of Bridgman-Stockbarger growth

A thermal analysis of a cylindrical HgCdTe sample in a Bridgman-Stockbarger crystal growth configuration was conducted with emphasis on the thermal profile, interface shape and position, and the thermal gradients at the liquid-solid interface. Alloys of HgTe and CdTe with compositions approximating 20 percent CdTe, 80 percent HgTe were used. This composition results in a bandgap suited for the detection of 10.6 micron CO2 radiation. The sensitivity of the sample thermal characteristics to important growth parameters, such as thermal diffusivities, thermal conductivities, furnace temperature profile, ampoule dimensions, and growth velocity was assessed. Numerical techniques and associated computational models necessary to analyze the heat transfer process within the sample and the Bridgman-Stockbarger boundary conditions were developed. This thermal analysis mode was programmed in FORTRAN V, and is currently operational on the MSFC Univac 1100 system.

Knopf, F. W.↗

Numerical solution of the problem of flame propagation by the use of the random element method

A numerical, grid-free algorithm is presented for one-dimensional reaction-diffusion model of laminar flame propagation in premixed gases. It is based on the random element method we developed for the analysis of diffusional processes. The effect of combustion is taken into account by applying the principle of fractional steps to separate the process of diffusion, modeled by the random walk of computational elements, from the exothermic effects of chemical reaction, monitoring their strength. The validity of the algorithm is demonstrated by application to flame propagation problems for which exact solutions exist. The flame speed evaluated by its use oscillates around the exact value at a relatively small amplitude, while the temperature and species concentration profiles are self-correcting in their convergence to the exact solution. A satisfactory resolution is obtained by the use of quite a small number of computational elements which automatically adjust their distribution of fit sharp gradients.

Ghoniem, A. F.↗