SEARCH · Engineering Papers

Results for “Convection diffusion equation”

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 235 records · Page 13

All sources NASA NTRS DOE OSTI

NASA NTRS1986

Vapor transport processes of Ge-chalcogenides under normal and reduced gravity conditions

Experimental and theoretical studies of the GeSe-GeI4 system revealed the multicomponent, multireaction nature of this system and the existence of a diffusion boundary layer; they also yielded equations for the prediction of diffusion limited mass transport rates. Vapor transport and crystal growth experiments of the GeSe-GeI4 and GeSe-Xe systems on earth and in a microgravity environment demonstrated the effects of convection and of microgravity on mass flux and crystal morphology. The surface and bulk morphology of GeSe crystals obtained for both systems in microgravity is considerably improved relative to ground-based specimens. Unexpected crystal growth phenomena were observed for the GeSe-Xenon system on recent Shuttle flights.

Wiedemeier, H.↗

NASA NTRS1998

Symbolic Computational Approach to the Marangoni Convection Problem With Soret Diffusion

A recently reported solution for stationary stability of a thermosolutal system with Soret diffusion is re-derived and examined using a symbolic computational package. Symbolic computational languages are well suited for such an analysis and facilitate a pragmatic approach that is adaptable to similar problems. Linearization of the equations, normal mode analysis, and extraction of the final solution are performed in a Mathematica notebook format. An exact solution is obtained for stationary stability in the limit of zero gravity. A closed form expression is also obtained for the location of asymptotes in relevant parameter, (Sm(sub c), Mac(sub c)), space. The stationary stability behavior is conveniently examined within the symbolic language environment. An abbreviated version of the Mathematica notebook is given in the Appendix.

Skarda, J. Raymond↗

NASA NTRS1985

Investigation of parabolic computational techniques for internal high-speed viscous flows

A feasibility study was conducted to assess the applicability of an existing parabolic analysis (ADD-Axisymmetric Diffuser Duct), developed previously for subsonic viscous internal flows, to mixed supersonic/subsonic flows with heat addition simulating a SCRAMJET combustor. A study was conducted with the ADD code modified to include additional convection effects in the normal momentum equation when supersonic expansion and compression waves were present. It is concluded from the present study that for the class of problems where strong viscous/inviscid interactions are present a global iteration procedure is required.

Anderson, O. L.↗

NASA NTRS1985

Numerical investigation of internal high-speed viscous flows using a parabolic technique

A feasibility study has been conducted to assess the applicability of an existing parabolic analysis (ADD-Axisymmetric Diffuser Duct), developed previously for subsonic viscous internal flows, to mixed supersonic/subsonic flows with heat addition simulating a SCRAMJET combustor. A study was conducted with the ADD code modified to include additional convection effects in the normal momentum equation when supersonic expansion and compression waves are present. A set of test problems with weak shock and expansion waves have been analyzed with this modified ADD method and stable and accurate solutions were demonstrated provided the streamwise step size was maintained at levels larger than the boundary layer displacement thickness. Calculations made with further reductions in step size encountered departure solutions consistent with strong interaction theory. Calculations were also performed for a flow field with a flame front in which a specific heat release was imposed to simulate a SCRAMJET combustor. In this case the flame front generated relatively thick shear layers which aggravated the departure solution problem. Qualitatively correct results were obtained for these cases using a marching technique with the convective terms in the normal momentum equation suppressed. It is concluded from the present study that for the class of problems where strong viscous/inviscid interactions are present a global iteration procedure is required.

Anderson, O. L.↗

NASA NTRS1991

Diurnal forcing of planetary atmospheres

A free convection parameterization has been introduced into the Mars Planetary Boundary Layer Model (MPBL). Previously, the model would fail to generate turbulence under conditions of zero wind shear, even when statically unstable. This in turn resulted in erroneous results at the equator, for example, when the lack of Coriolis forcing allowed zero wind conditions. The underlying cause of these failures was the level 2 second-order turbulence closure scheme which derived diffusivities as algebraic functions of the Richardson number (the ratio of static stability to wind shear). In the previous formulation, the diffusivities were scaled by the wind shear--a convenient parameter since it is non-negative. This was the drawback that all diffusivities are zero under conditions of zero shear (viz., the free convection case). The new scheme tests for the condition of zero shear in conjunction with static instability and recalculates the diffusivities using a static stability scaling. The results for a simulation of the equatorial boundary layer at autumnal equinox are presented. (Note that after some wind shear is generated, the model reverts to the traditional diffusivity calculation.)

Houben, Howard C.↗

NASA NTRS1974

The differential rotation of the solar surface

The large-scale flow in the solar convection zone is discussed. The objective is to deduce from observation the principal physical balances in the governing equations. The simplest set of equations that seem potentially realistic are utilized. It is assumed that magnetic fields are negligible, that mixing-length theory gives an accurate representation of the mean structure, and that rigid-body rotation exists at the base. It is deduced that the zonal glow is geostrophic, the meridional flow is controlled by friction, and diffusive heating balances advective cooling due to vertical motion. Next, a detailed calculation of the latitude profile of surface angular velocity is performed, based on approximate equations containing only the principal balances, and agrees well with observation. Finally, it is demonstrated that the amplitude of the flow may be determined by the constraint that the net equatorward angular momentum flux vanish. The conclusions are consistent with the hypothesis that the flow is primarily a single large axisymmetric convection cell in each hemisphere.

Gierasch, P. J.↗

NASA NTRS1993

Numerical modeling of physical vapor transport under microgravity conditions: Effect of thermal creep and stress

One of the most promising applications of microgravity (micro-g) environments is the manufacture of exotic and high-quality crystals in closed cylindrical ampoules using physical vapor transport (PVT) processes. The quality enhancements are believed to be due to the absence of buoyant convection in the weightless environment - resulting in diffusion-limited transport of the vapor. In a typical experiment, solid-phase sample material is initially contained at one end of the ampoule. The sample is made to sublime into the vapor phase and deposit onto the opposite end by maintaining the source at an elevated temperature with respect to the deposit. Identification of the physical factors governing both the rates and uniformity of crystal growth, and the optimization of the micro-g technology, will require an accurate modeling of the vapor transport within the ampoule. Previous micro-g modeling efforts have approached the problem from a 'classical' convective/diffusion formulation, in which convection is driven by the action of buoyancy on thermal and solutal density differences. The general conclusion of these works have been that in low gravity environments the effect of buoyancy on vapor transport is negligible, and vapor transport occurs in a diffusion-limited mode. However, it has been recently recognized than in the non-isothermal (and often low total pressure) conditions encountered in ampoules, the commonly-assumed no-slip boundary condition to the differential equations governing fluid motion can be grossly unrepresentative of the actual situation. Specifically, the temperature gradients can give rise to thermal creep flows at the ampoule side walls. In addition, temperature gradients in the vapor itself can, through the action of thermal stress, lead to bulk fluid convection.

Mackowski, Daniel W.↗

NASA NTRS1984

Reducing numerical diffusion for incompressible flow calculations

A number of approaches for improving the accuracy of incompressible, steady-state flow calculations are examined. Two improved differencing schemes, Quadratic Upstream Interpolation for Convective Kinematics (QUICK) and Skew-Upwind Differencing (SUD), are applied to the convective terms in the Navier-Stokes equations and compared with results obtained using hybrid differencing. In a number of test calculations, it is illustrated that no single scheme exhibits superior performance for all flow situations. However, both SUD and QUICK are shown to be generally more accurate than hybrid differencing.

Claus, R. W.↗

NASA NTRS1977

Theoretical and experimental studies of the deposition of Na2SO4 from seeded combustion gases

A basic point in the hot corrosion of turbine components is the deposition of sodium sulfate from flames containing sodium and sulfur. An experimental study is described which examines a dew point prediction theory based on the local thermochemical equilibrium (LTCE) method, and a means to calculate the deposition rate is suggested. In addition, a convective diffusion theory, based on the assumption of a chemically frozen boundary layer, utilizing the LTCE results, and imposing the additional effects of mass transport, was also successful in predicting dew points for Na2SO4-seeded combustion gases. A multicomponent mass transfer equation was derived to predict NaSO4 deposition rate via vapor transport at temperatures below the deposition onset temperature.

Kohl, F. J.↗

NASA NTRS1988

Continuum theories for fluid-particle flows: Some aspects of lift forces and turbulence

A general framework is outlined for the modeling of fluid particle flows. The momentum exchange between the constituents embodies both lift and drag forces, constitutive equations for which can be made explicit with reference to known single particle analysis. Relevant results for lift are reviewed, and invariant representations are posed. The fluid and particle velocities and the particle volume fraction are then decomposed into mean and fluctuating parts to characterize turbulent motions, and the equations of motion are averaged. In addition to the Reynolds stresses, further correlations between concentration and velocity fluctuations appear. These can be identified with turbulent transport processes such as eddy diffusion of the particles. When the drag force is dominant, the classical convection dispersion model for turbulent transport of particles is recovered. When other interaction forces enter, particle segregation effects can arise. This is illustrated qualitatively by consideration of turbulent channel flow with lift effects included.

Mctigue, David F.↗

NASA NTRS1999

Diffusion, Fluxes, Friction Forces, and Joule Heating in Two-Temperature Multicomponent Magnetohydrodynamics

The relationship between Joule heating, diffusion fluxes, and friction forces has been studied for both total and electron thermal energy equations, using general expressions for multicomponent diffusion in two-temperature plasmas with the velocity dependent Lorentz force acting on charged species in a magnetic field. It is shown that the derivation of Joule heating terms requires both diffusion fluxes and friction between species which represents the resistance experienced by the species moving at different relative velocities. It is also shown that the familiar Joule heating term in the electron thermal energy equation includes artificial effects produced by switching the convective velocity from the species velocity to the mass-weighted velocity, and thus should not be ignored even when there is no net energy dissipation.

Chang, C. H.↗

NASA NTRS2004

Flows in the Solar Convection Zone

Flows within the solar convection zone are the primary drivers of the Sun's magnetic activity cycle. Differential rotation stretches out the magnetic field and converts poloidal fields into toroidal fields. Zones of strong radial shear are found at both the surface and at the base of the convection zone (the tachocline). The poleward meridional flow near the surface transports magnetic flux that is observed to reverse the magnetic poles near the time of cycle maxima. The deeper (and as yet unobserved), equatorward meridional flow should carry magnetic flux toward the equator where it reconnects with oppositely directed fields in the other hemisphere. The non-axisymmetric flows (granules, supergranules, and giant cells) also transport magnetic flux but in a more random, diffusive, manner. Supergranules and giant cells also play significant roles in driving the large-scale, axisymmetric flows themselves. The effect of solar rotation on supergranulation produces the shear layer near the surface and enhances the meridional flow. The effect of solar rotation on giant cells should produce the latitudinal differential rotation, the shear in the tachocline, and the meridional circulation. In this presentation I will describe the observed and theorized characteristics of the flows in the solar convection zone and discuss their connections to the solar activity cycle.

Hathaway, D. H.↗

NASA NTRS2007

Transient Numerical Modeling of Catalytic Channels

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.

Struk, Peter M.↗

NASA NTRS2000

Numerical Simulation of a Chemically Reacting Sorbent Bed for LSS Applications

A detailed numerical model of a chemisorption bed has been developed. The model is based on the constant pressure mass transport equation for gaseous flow through a packed bed, and the equation for diffusion and reaction within a spherical particle. Because there is a wealth of data from the NASA and the Navy bodies of literature, the LiOH-H2O-CO2 system is chosen for application of the model and interpretation of results. Prior models of this system from the life support literature are limited. The current model incorporates many of the features of elaborate models developed for investigation of industrial systems or energy applications (e.g., coal, desulphurization): it distinguishes bulk convection and bed dispersion; mass transport to the particle surface, transport within the particle, and reaction. It uses the nonsteady (not pseudo-steady state) form of the equations. The chemistry is modeled as a multi-step, reversible reaction with evolving solid structure. The resulting system of equations is large. The ODEPACK family of solvers is used to integrate the system. Reaction coefficients are determined by experiment. Typical results of the model are illustrated with mission input parameters. Using the model, an explanation is offered for 1) the varied performance results found after pre-breathing (or after simulated pre-breathe conditions), 2) interrupted use and 3) low temperature use. In addition, options for a reusable canister are explored. The computational resource implications of adding energy equations are discussed briefly, as are applicability to other relevant space and undersea systems.

Luna, Bernadette↗

NASA NTRS1975

Cosmic-ray streaming and anisotropies

The paper is concerned with the differential current densities and anisotropies that exist in the interplanetary cosmic-ray gas, and in particular with a correct formulation and simple interpretation of the momentum equation that describes these on a local basis. Two examples of the use of this equation in the interpretation of previous data are given. It is demonstrated that in interplanetary space, the electric-field drifts and convective flow parallel to the magnetic field of cosmic-ray particles combine as a simple convective flow with the solar wind, and that there exist diffusive currents and transverse gradient drift currents. Thus direct reference to the interplanetary electric-field drifts is eliminated, and the study of steady-state and transient cosmic-ray anisotropies is both more systematic and simpler.

Forman, M. A.↗

NASA NTRS1993

A numerical procedure for analysis of finite rate reacting flows

Combustion processes in rocket propulsion systems are characterized by the existence of multiple, vastly differing time and length scales, as well as flow-speeds at wide variation of Mach numbers. The chemical kinetics processes in the highly active reaction zone are characterized by much smaller scales compared to fluid convective and diffusive time scales. An operator splitting procedure for transient finite rate chemistry problems has been developed using a pressure based method, which can be applied to all speed flows without difficulties. The splitting of chemical kinetics terms formed the fluid-mechanical terms of the species equation ameliorated the difficulties associated with the disparate time scales and stiffness in the set of equations which describes highly exothermic combustion. A combined efficient ordinary differential equations (ODE) solver was used to integrate the effective chemical source terms over the residence time at each grid cell. One and two dimensional reacting flow situations were carried out to demonstrate and verify the current procedure. Different chemical kinetics with different degrees of nonlinearity have also been incorporated to test the robustness and generality of the proposed method.

Shang, H. M.↗

NASA NTRS2000

Modeling of a Non-Dilute Alloy Solidification Under Terrestrial and Microgravity Conditions

A general finite element model of non-dilute alloy solidification was used to simulate growth of HgCdTe from the melt. 2-D and 3-D models with and without Oberbeck-Boussinesq (O-B) approximation were used. Verification of the 3-D model was based on experimental data. The location and the shape of the interface were determined using melting temperatures obtained from the HgCdTe liquidus curve. Double-diffusive convection in the melt is caused by radial temperature gradients and by rejection at the solid-melt interface of the heavier HgTe-rich solute. Because of these complicating factors, dimensional rather than non-dimensional modeling was performed. Models based on complete set of Navier-Stokes equations for compressible media are superior to the O-B approximation in microgravity, because they provide more accurate estimation of convection.

Bune, Andris V.↗

NASA NTRS1976

Planetary spin period acceleration of particles in the Jovian magnetosphere

A four-step mechanism is proposed for the acceleration of energetic protons and relativistic electrons in Jupiter's magnetosphere. According to this mechanism, photoelectrons and ions from the Jovian ionosphere are: (1) ejected along magnetic-field lines toward the equator by the centrifugal force of corotation; (2) accelerated by magnetic-field annihiliation in the magnetotail, which process is modulated at Jupiter's rotational frequency; (3) trapped on closed field lines in the reconnection process, convected inward toward Jupiter from the merging region, and subjected to adiabatic compression; and (4) diffused inward by the conventional radial-diffusion process through a violation of the third adiabatic invariant. It is shown that the proposed mechanism produces magnetic moments much larger than those available from inward diffusion of solar-wind particles or motional emf acceleration at the Galilean satellites, provides a natural explanation for the 10-hr periodicity of the energetic particle fluxes observed inside the magnetosphere by the Pioneer spacecraft, and also produces a 10-hr periodicity in the energetic particle flux from the magnetosphere into interplanetary space in such a way that the phase of interplanetary flux variations is locked to the rotational phase of Jupiter

Carbary, J. F.↗