Engineering Papers⌕ Search

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.

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

277 records · Page 16
All sourcesNASA NTRSDOE OSTI
NASA NTRS1995

Coupled Monte Carlo Probability Density Function/ SPRAY/CFD Code Developed for Modeling Gas-Turbine Combustor Flows

The success of any solution methodology for studying gas-turbine combustor flows depends a great deal on how well it can model various complex, rate-controlling processes associated with turbulent transport, mixing, chemical kinetics, evaporation and spreading rates of the spray, convective and radiative heat transfer, and other phenomena. These phenomena often strongly interact with each other at disparate time and length scales. In particular, turbulence plays an important role in determining the rates of mass and heat transfer, chemical reactions, and evaporation in many practical combustion devices. Turbulence manifests its influence in a diffusion flame in several forms depending on how turbulence interacts with various flame scales. These forms range from the so-called wrinkled, or stretched, flamelets regime, to the distributed combustion regime. Conventional turbulence closure models have difficulty in treating highly nonlinear reaction rates. A solution procedure based on the joint composition probability density function (PDF) approach holds the promise of modeling various important combustion phenomena relevant to practical combustion devices such as extinction, blowoff limits, and emissions predictions because it can handle the nonlinear chemical reaction rates without any approximation. In this approach, mean and turbulence gas-phase velocity fields are determined from a standard turbulence model; the joint composition field of species and enthalpy are determined from the solution of a modeled PDF transport equation; and a Lagrangian-based dilute spray model is used for the liquid-phase representation with appropriate consideration of the exchanges of mass, momentum, and energy between the two phases. The PDF transport equation is solved by a Monte Carlo method, and existing state-of-the-art numerical representations are used to solve the mean gasphase velocity and turbulence fields together with the liquid-phase equations. The joint composition PDF approach was extended in our previous work to the study of compressible reacting flows. The application of this method to several supersonic diffusion flames associated with scramjet combustor flow fields provided favorable comparisons with the available experimental data. A further extension of this approach to spray flames, three-dimensional computations, and parallel computing was reported in a recent paper. The recently developed PDF/SPRAY/computational fluid dynamics (CFD) module combines the novelty of the joint composition PDF approach with the ability to run on parallel architectures. This algorithm was implemented on the NASA Lewis Research Center's Cray T3D, a massively parallel computer with an aggregate of 64 processor elements. The calculation procedure was applied to predict the flow properties of both open and confined swirl-stabilized spray flames.

Source record↗
NASA NTRS1996

Fluid Creep Effects on Near-Wall Solute Transport for Non-Isothermal Ampoules

There is a growing practical and theoretical interest in developing accurate macroscopic modelling for flows arising in chemical or physical vapor transport (VT) crystal growth experiments, including those conducted in reduced gravity environments. Rosner was the first person to point out that previously neglected rarefield gas dynamics phenomena (Stefan and bouyancy-driven flows) become rather important sources of convection. In particular, the combination of rarefaction and strong gradients of temperature (and/or concentration) tangential to the side-walls of the ampoule induces convective flows known as thermal (and concentration) 'creep' respectively. His order-of-magnitude estimates revealed that thermal creep effects can be non-negligible even at normal gravitational levels. On the macroscopic level, the bulk fluid mechanics can be adequately described by the familiar macroscopic equations as long as the boundary conditions are modified to account for the integrated effect of kinetic boundary layers adjacent to solid boundaries. Motivated by the growing importance of these phenomena, we have embarked on a series of computational studies to elucidate these fundamental creep-induced effects for a rarefied gas in simple, two-dimensional confined geometries. However, unlike previous related studies, we resort to a microscopic description of the gas, mathematically expressed by the Boltzmann integro-differential equation. We employ the direct simulation Monte Carlo (DSMC) method of Bird, the theoretical foundations and several practical applications. In the case of thermally induced flows, the no-time counter method of Bird is used, as implemented for a hard-sphere gas. The scheme has been also extended to account for realistic molecular interaction models, an extension necessary if the diffusion physics underlying concentration creep are to be captured.

Papadopoulos, Dimitrios↗
NASA NTRS1993

Influence of test configuration on the combustion characteristics of polymers as ignition sources

The experimental evaluation of polymers as ignition sources for metals was accomplished at the NASA White Sands Test Facility (WSTF) using a standard promoted combustion test. These tests involve the transient burning of materials in high-pressure oxygen environments. They have provided data from which design decisions can be made; data include video recordings of ignition and non-ignition for specific combinations of metals and polymers. Other tests provide the measured compositions of combustion products for polymers at select burn times and an empirical basis for estimating burn rates. With the current test configuration, the detailed analysis of test results requires modeling a three-dimensional, transient convection process involving fluid motion, thermal conduction and convection, the diffusion of chemical species, and the erosion of sample surface. At the high pressure extremes, it even requires the analysis of turbulent, transient convection where the physics of the problem are not well known and the computation requirements are not practical at this time. An alternative test configuration that can be analyzed with a relatively-simple convection model was developed during the summer period. The principal change constitutes replacing a large-diameter polymer disk at the end of the metal test rod with coaxial polymer cylinders that have a diameter nearer to that of the metal rod. The experimental objective is to assess the importance of test geometries on the promotion of metal ignition by testing with different lengths of the polymer and, with an extended effort, to analyze the surface combustion in the redesigned promoted combustion tests through analytical modeling of the process. The analysis shall use the results of cone-calorimeter tests of the polymer material to model primary chemical reactions and, with proper design of the promoted combustion test, modeling of the convection process could be conveniently limited to a quasi-steady boundary layer analysis where the economical solution of parabolic equations is involved. The products for the summer period are: (1) a conceptual-level redesign of the test apparatus, and (2) the development and use of an approximate integral boundary layer analysis to demonstrate the influence of geometry changes prior to testing. A computer code STAN5, an accurate numerical boundary layer model whose earlier versions were developed for the NASA Lewis Research Center by the Fellow, also was installed and validated on the WSTF and New Mexico State University computer systems as a starting point in the development of a more detailed fluid mechanics and combustion model.

Julien, Howard L.↗
NASA NTRS2014

A Computational and Experimental Study of Coflow Laminar Methane/Air Diffusion Flames: Effects of Fuel Dilution, Inlet Velocity, and Gravity

The influences of fuel dilution, inlet velocity, and gravity on the shape and structure of laminar coflow CH4-air diffusion flames were investigated computationally and experimentally. A series of nitrogen-diluted flames measured in the Structure and Liftoff in Combustion Experiment (SLICE) on board the International Space Station was assessed numerically under microgravity (mu g) and normal gravity (1g) conditions with CH4 mole fraction ranging from 0.4 to 1.0 and average inlet velocity ranging from 23 to 90 cm/s. Computationally, the MC-Smooth vorticity-velocity formulation was employed to describe the reactive gaseous mixture, and soot evolution was modeled by sectional aerosol equations. The governing equations and boundary conditions were discretized on a two-dimensional computational domain by finite differences, and the resulting set of fully coupled, strongly nonlinear equations was solved simultaneously at all points using a damped, modified Newton's method. Experimentally, flame shape and soot temperature were determined by flame emission images recorded by a digital color camera. Very good agreement between computation and measurement was obtained, and the conclusions were as follows. (1) Buoyant and nonbuoyant luminous flame lengths are proportional to the mass flow rate of the fuel mixture; computed and measured nonbuoyant flames are noticeably longer than their 1g counterparts; the effect of fuel dilution on flame shape (i.e., flame length and flame radius) is negligible when the flame shape is normalized by the methane flow rate. (2) Buoyancy-induced reduction of the flame radius through radially inward convection near the flame front is demonstrated. (3) Buoyant and nonbuoyant flame structure is mainly controlled by the fuel mass flow rate, and the effects from fuel dilution and inlet velocity are secondary.

flames↗
NASA NTRS2000

Convective Melting of Particles in Flow Under Microgravity Conditions

Study of melting of dispersed or packed solid particles in a fluid under gravity and microgravity conditions provides benchmark information for many engineering applications such as material processing, environmental assessment and protection and space fire protection. During such processes, packed or dispersed solid particles are interacting with fluid flow at above-melting temperatures. By unmasking the buoyancy effects in coupled flow, phase change and heat transport phenomena, a better understanding of melting rate in non-thermal equilibrium, convective conditions can be studied. A series of flight experiments were conducted onboard the NASA KC-135 Microgravity Research Airplane. The Particle Melting in Plow (PMF) module was designed to allow flow through the initially packed ice particles at controlled temperature and velocity. To achieve this, a close-loop flow system was designed. Video images were taken to record the visualization of the melting process, from which a time variation of packed particle thickness distribution at different times can be obtained by the image analysis method. The fluid temperature distribution within the melting zone is measured by thermocouples. An infrared camera was mounted from the top of the test section to record the ice-water thermal images at a given location. The results from thermal images yield local temperature variation between melting solid and liquid and local Stephan number. Typical results for a number of cases are presented. The mathematical model, describing mass, energy, and momentum balance equations for the liquid and solid phases, is presented. It is found that melting rate is influenced mainly by the ratio of Reynolds number (based on the initial particle diameter) to the Froud number, and Stephan number. At the absence of gravity, Froud number approaches zero, Reynolds number and Stephan number become dominant factors governing the melting rate. The numerically determined results are compared with the experimental ones. It is found that the discrepancy between the predicted and measured melting rate is largely due to the inaccuracy in the constitutive equations for effective thermal physical properties, such as effective thermal conductivity and diffusivity, and transport properties, such as particle interaction coefficient, and local heat transfer coefficient of particles.

Stewart, C.↗
NASA NTRS1999

Control of Meridional Flow in Circular Cylinders by a Travelling Axial Magnetic Field

Convective flow in a Bridgman or float zone configuration significantly affects the interface shape and segregation phenomena. While the primary causative factor for this flow is buoyancy induced convection in an enclosed Bridgman melt, the presence of a free surface gives rise to surface tension driven flows in the floating zone processing of melts. It is of interest to curtail these flows in order to realize near quiescent growth conditions that have shown to result in crystals with good longitudinal and radial homogeneity and thereby of better overall quality. While buoyancy effects can be reduced by careful processing in a low gravity (space) environment, the reduction of Marangoni flows due to surface tension variations is not that straight forward. Attempts have been made with some limited success with the use of external fields to affect the melt thermo-fluid behavior. The use of a static magnetic field that reduces convective contamination through the effects of a non-intrusively induced, dissipative Lorentz force in an electrically conducting melt is one such approach. Experiments have shown that axial fields of the order of 5 Tesla can significantly eliminate convection and yield close to diffusion limited crystal growth conditions. The generation and use of such high magnetic fields require substantial hardware and incur significant costs for its operation. Lately, the use of rotating magnetic fields has been tested in semiconductor crystal growth. The method is fairly well known and commonly used in metal processing but its adaptation to crystal growth of semiconductors is fairly recent. The elegance of the technique rests in its low power requirement (typically 10-20 milli-Tesla at 50-400 Hz) and its efficacy in curtailing deleterious temperature fluctuations in the melt. A rotating magnetic field imposes a rotational force and thereby induces a circulation within the melt that tends to dominate other sporadic convective effects. Thus a known low level of convective flow is introduced into the system. A new novel variation of the Lorentz force mechanism is proposed and investigated in this study. Since one of the desired process conditions in melt crystal growth is the minimization of convective effects, this investigation examines the use of an external field of magnetic origin to counteract existing convective flow within the melt. This is accomplished by utilizing a running or traveling axial magnetic wave in the system. The concept is similar to the use of vibrational means in order to induce streaming flows that oppose buoyant or surface tension driven convection in the system. The rotation direction as well as the magnitude (strength) of this circulation can be easily controlled by external inputs thus affording a direct means of controlling the developing shape of the crystallizing front (interface). The theoretical model of this technique is fully developed and presented in this paper. Results from the solution of the developed governing equations and boundary conditions are also presented. An experimental demonstration of the concept is presented through the suppression of natural convective flow in a mercury column. Implications to crystal growth systems will be fully explored in the final manuscript.

Mazuruk, K.↗
NASA NTRS2001

Experimental and Computational Studies of the Control of Convection of Non-Conducting Liquids During solidification by Means of a Magnetic Field Gradient

The elimination of convection is essential in experimental investigations of diffusive transport (of heat and matter) during solidification. One classical approach to damping convection in a conducting liquid is the application of a magnetic field. The damping phenomenon is the induction, by the motion of a conductor in a magnetic field, of currents which interact with the field to produce Lorentz forces that oppose the flow. However, there are many liquids which are not sufficiently conducting to exploit this phenomenon; examples include the transparent liquids (such as succinonitrile-acetone) that are used as "model alloys" in fundamental solidification studies. There have been several investigations of the solidification of these liquids that have been carried out in orbiting laboratories to eliminate natural convection. The paper describes an investigation of an alternative approach whereby a magnetic field gradient is applied to the liquid. A magnetic body force then arises which is dependent on the susceptibility of the liquid and thereby on the temperature and or concentration. With the field gradient aligned vertically and of correct magnitude, the variation of gravitational body force due to temperature/concentration dependent density can be counterbalanced by a variation in magnetic body force. Experiments have been carried out in a super-conducting magnet at Marshall Space Flight Center to measure velocities in an aqueous manganese chloride solution. The solution was contained in a chamber with temperature controlled end walls and glass side walls. Velocities were measured by particle image velocimetry. Starting from zero current in the magnet (zero field gradient) flow driven by the temperature difference between the end walls was measured. At a critical current the flow was halted. At higher currents the normal convection was reversed. The experiments included ones where the solution was solidified and were accompanied by solution of the flow/transport equations using the software package FLUENT.

Seybert, C. D.↗