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 433 records · Page 24

Seasonal and magnetic storm related changes in the thermosphere induced by eddy mixing

The possibility of explaining the seasonal and magnetic storm related changes in the thermosphere and the ionosphere through variation in eddy diffusion coefficient in the lower thermosphere is investigated theoretically by obtaining simultaneous numerical solutions of the relevant continuity, momentum, and energy balance equations in a self-consistent manner. It is shown that various important features of thermospheric seasonal behavior, including the winter helium bulge phenomenon, can be explained by assuming effective reduction in the winter-time eddy diffusion coefficient by about one order of magnitude. Storm related changes in the thermal structure and composition are seen to arise by assuming a downward shift of turbulence. Physical implications and relative significance of variation in thermospheric turbulence and large scale inter-hemispheric circulation are briefly discussed.

Sinha, A. K.↗

Diffusion effects on the line intensities of He I and He II in the solar transition region

A heuristic treatment of diffusion in the solar chromosphere-corona transition region is developed. It is shown that diffusion becomes increasingly important with steeper temperature gradients, in active and quiet regions relative to coronal holes, and with increasing excitation potential. Numerical calculations are made for the resonance lines of He I and He II and show that diffusion can enhance these lines. Thus the helium lines may appear relatively weak in coronal holes due to a weakening of the enhancement mechanism. Most transition region lines will be less affected by diffusion than He I or He II.

Shine, R.↗

Numerical simulation of electrophoresis separation processes

A new Petrov-Galerkin finite element formulation has been proposed for transient convection-diffusion problems. Most Petrov-Galerkin formulations take into account the spatial discretization, and the weighting functions so developed give satisfactory solutions for steady state problems. Though these schemes can be used for transient problems, there is scope for improvement. The schemes proposed here, which consider temporal as well as spatial discretization, provide improved solutions. Electrophoresis, which involves the motion of charged entities under the influence of an applied electric field, is governed by equations similiar to those encountered in fluid flow problems, i.e., transient convection-diffusion equations. Test problems are solved in electrophoresis and fluid flow. The results obtained are satisfactory. It is also expected that these schemes, suitably adapted, will improve the numerical solutions of the compressible Euler and the Navier-Stokes equations.

Ganjoo, D. K.↗

Relaxation of the loss-cone by quasi-linear diffusion of the electron-cyclotron maser instability in the solar corona

A self-consistent numeric two-dimensional code using kinetic-wave particle equations has been developed and applied to the maser dynamics of the solar corona. Time histories generated by the code reveal details of the evolution of the linear phase, the saturation and subsequent relaxation of the loss cone. Quantitative values are obtained for the saturation time, the amount of the converted free energy, the critical wave energy level for onset of maser diffusion, the velocity and pitch angle range of the resonant particles, and the relative importance of the different magnetoionic modes and harmonics. A wide range of initial conditions is covered by varying the loss-cone distribution and the ambient cold plasma parameters.

Aschwanden, Markus J.↗

Preliminary numerical analysis of improved gas chromatograph model

A mathematical model for the gas chromatograph was developed which incorporates the heretofore neglected transport mechanisms of intraparticle diffusion and rates of adsorption. Because a closed-form analytical solution to the model does not appear realizable, techniques for the numerical solution of the model equations are being investigated. Criteria were developed for using a finite terminal boundary condition in place of an infinite boundary condition used in analytical solution techniques. The class of weighted residual methods known as orthogonal collocation is presently being investigated and appears promising.

Woodrow, P. T.↗

Thermodynamics and kinetics of reactions in protective coating systems

Investigations of fluoride activated packs with Al:Ni ratios greater than 50 a/o prove that the specimen surface is not in equilibrium with the pack at high Al:Ni ratios but that an activity gradient exists between pack and specimen. Therefore, gaseous diffusion and possibly surface reactions play a role in determining the overall rate of Al deposition in such packs. Noticeable differences in coating behavior have been obtained in packs activated with chloride and iodide, and it appears that poorest results are obtained with iodides, better with chlorides, and best with fluorides. A numerical method has been perfected for calculating rates of solid-state diffusion controlled coating formation, allowing for the variation of diffusivity with composition in the NiAl phase. Layer growth rates can now be accurately predicted from a knowledge of the surface and substrate compositions. Furthermore, the correct diffusion profiles are obtained by this method. These differ substantially from the profile obtained when the diffusivity is assumed constant.

Gupta, B.↗

Vortex-scalar element calculations of a diffusion flame stabilized on a plane mixing layer

The vortex-scalar element method, a scheme which utilizes vortex elements to discretize the region of high vorticity and scalar elements to represent species or temperature fields, is utilized in the numerical simulations of a two-dimensional reacting mixing layer. Computations are performed for a diffusion flame at high Reynolds and Peclet numbers without resorting to turbulence models. In the nonreacting flow, the mean and fluctuation profiles of a conserved scalar show good agreement with experimental measurements. Results for the reacting flow indicate that for temperature independent kinetics, the chemical reaction begins immediately downstream of the splitter plate where mixing starts. Results for the reacting flow with Arrhenius kinetics show an ignition delay, which depends on reactant temperature, before significant chemical reaction occurs. Harmonic forcing changes the structure of the layer, and concomitantly the rates of mixing and reaction, in accordance with experimental results. Strong stretch within the braids in the nonequilibrium kinetics case causes local flame quenching due to the temperature drop associated with the large convective fluxes.

Ghoniem, Ahmed F.↗

Detailed Studies on the Structure and Dynamics of Reacting Dusty Flows at Normal and Microgravity

Dusty reacting flows are of particular interest for a wide range of applications. Inert particles can alter the flammability and extinction limits of a combustible mixture. Reacting particles can release substantial amount of heat and can be used either for power generation or propulsion. Accumulation of combustible particles in air can result in explosions which, for example, can occur in grain elevators, during lumber milling and in mine galleries. Furthermore, inert particles are used as flow velocity markers in reacting flows, and their velocity is measured by non-intrusive laser diagnostic techniques. Despite their importance, dusty reacting flows have been less studied and understood compared to gas phase as well as sprays. The addition of solid particles in a flowing gas stream can lead to strong couplings between the two phases, which can be of dynamic, thermal, and chemical nature. The dynamic coupling between the two phases is caused by the inertia that causes the phases to move with different velocities. Furthermore, gravitational, thermophoretic, photophoretic, electrophoretic, diffusiophoretic, centrifugal, and magnetic forces can be exerted on the particles. In general, magnetic, electrophoretic, centrifugal, photophoretic, and diffusiophoretic can be neglected. On the other hand, thermophoretic forces, caused by steep temperature gradients, can be important. The gravitational forces are almost always present and can affect the dynamic response of large particles. Understanding and quantifying the chemical coupling between two phases is a challenging task. However, all reacting particles begin this process as inert particles, and they must be heated before they participate in the combustion process. Thus, one must first understand the interactions of inert particles in a combustion environment. The in-detail understanding of the dynamics and structure of dusty flows can be only advanced by considering simple flow geometries such as the opposed-jet, stagnation-type. In such configurations the imposed strain rate is well characterized, and the in-depth understanding of the details of the physico-chemical processes can be systematically obtained. A number of computational and experimental studies on spray and particle flows have been conducted in stagnation-type configurations. Numerically, the need for a hybrid Eulerian-Lagrangian approach has been identified by Continillo and Sirignano, and the use of such approach has allowed for the prediction of the phenomenon of droplet flow reversal. Gomez and Rosner have conducted a detailed study on the particle response in the opposed-jet configuration, and the particle thermophoretic diffusivities were determined experimentally. Sung, Law and co-workers have conducted numerical studies on the effect of strain rate and temperature gradients on the dynamics of inert particles, as a way of understanding potential errors in experimental LDV data that may arise from thermophoretic forces. This investigation is a combined experimental and numerical study on the details of reacting dusty flows. The specific tasks are: (1) Experimental determination of laminar flame speeds, and extinction strain rates of dusty flows at normal- and micro-gravity as functions of the particle type, particle initial diameter, particle initial number density, and gas phase chemical composition; (2) Detailed numerical simulation of the experiments. Results are compared with experiments and the adequacy of theoretical models is assessed; and (3) Provision of enhanced insight into the thermo-chemical coupling between the two phases.

Andac, M. Gurhan↗

Nitric oxide delta band absorption measurements in the lower thermosphere

The first successful measurements of nitric oxide band absorption in the earth's atmosphere have been made by using the high-resolution ultraviolet spectrometers on OSO 8. These measurements were conducted between January and August 1977 and between 30 deg north and south latitudes. The nitric oxide absorption data, at the (1,0) delta band wavelength of 1829.41 A, are numerically inverted to produce nitric oxide density profiles between 85 and 145 km. The major finding of the observations is that for altitudes below 95 km there is more nitric oxide at sunrise than at sunset. A numerical model, incorporating time-dependent ionospheric chemistry and vertical diffusive transport (both eddy and molecular), is used to simulate the OSO 8 nitric oxide profiles. The peak of the model sunset profile lies at 95-km altitude, in agreement with observation, if the vertical component of the eddy diffusion coefficient is 10 million sq cm/s near 100 km.

Massie, S. T.↗

Dissolution of multicomponent bubbles

The behavior of an isolated, stationary, multicomponent gas bubble in a glassmelt containing several dissolved gases is considered. The relevant mass-transport equations are formulated and calculations are performed for the case of two diffusing gases using a quasi-stationary model and a numerical solution of the exact mass-transfer equations. The results obtained from these two approaches are compared. The factors which govern the dissolution or growth of a bubble are thermodynamic and kinetic in origin. The tendency of a bubble to grow or shrink at long times is controlled by departure from overall equilibrium, whereas the short-time bubble dynamics may be dominated by kinetic effects. As a result of the existence of these dual influences, maxima and/or minima occur in the functional dependence of the bubble radius on time.

Weinberg, M. C.↗

Generation of auroral Omega bands by shear instability of the neutral winds

Thermospheric neutral wind acceleration via ion drag in the conducting E-region of the ionosphere is greatly increased by electron precipitation associated with auroras. This increased acceleration can lead to the development of significant horizontal wind shears, which were found to be unstable to the Kelvin-Helmholtz shear instability. Numerical simulation of the neutral response to an intense, postmidnight, diffuse aurora shows tne formation of an E-region 'jet stream' within the aurora, with peak winds speeds greather than 700 m/s after one hour. It is proposed that this jet stream produces unstable Kelvin-Helmholtz waves, which can drive waves of discrete aurora along the poleward boundary of the preexisting diffuse aurora. It is suggested that such auroral waves, driven by the neutral winds, form eastward propagating waves (omega bands) occasionally observed along the poleward boundary of postmidnight diffuse auroras. It was found that neutral wind shears that develop in response to discrete auroral arcs are unstable; however, the resulting wind waves are not expected to drive significant auroral waves along discrete arcs.

Lyons, L. R.↗

Transition mixing study

A computer model capable of analyzing the flow field in the transition liner of small gas turbine engines is developed. A FORTRAN code has been assembled from existing codes and physical submodels and used to predict the flow in several test geometries which contain characteristics similar to transition liners, and for which experimental data was available. Comparisons between the predictions and measurements indicate that the code produces qualitative results but that the turbulence models, both K-E and algebraic Reynolds Stress, underestimate the cross-stream diffusion. The code has also been used to perform a numerical experiment to examine the effect of a variety of parameters on the mixing process in transition liners. Comparisons illustrate that geometries with significant curvature show a drift of the jet trajectory toward the convex wall and weaker wake region vortices and decreased penetration for jets located on the convex wall of the liner, when compared to jets located on concave walls. Also shown were the approximate equivalency of angled slots and round holes and a technique by which jet mixing correlations developed for rectangular channels can be used for can geometries.

Reynolds, R.↗

Two dimensional analysis of a high temperature gaseous radiation receiver

The characteristics of the Flowing Gas Radiation Receiver (FGRR), a device that absorbs solar radiation volumetrically in a gas to produce high temperatures for space propulsion and power applications, are analyzed using a two-dimensional axisymmetric numerical model of the flow and radiation fields within a diffusely reflecting channel. The results show that an FGRR system is capable of generating temperatures in excess of 3000 K with collection efficiencies of approximately 75 percent for a channel with a reflectivity of 0.9. For a collinear radiation source, outflow temperatures of 3193 and 3092 K were achieved for axial and radial flow inputs, respectively, with receiver efficiencies of 0.82 and 0.76.

Mcfall, K. A.↗

A numerical study of Li-SF6 wick combustion - Forced and mixed convective burning

A numerical study is conducted to study Li-SF6 wick diffusion flames under mixed convective burning conditions at a pressure of 0.01 MPa. Both planar and cylindrical wicks are considered. The model is based on a conserved scalar approach. The objective of this study is to assess the effects of particular parameters on the burning rate and heat transfer. The flat-plate solution yields a fuel mass burning rate per unit surface area following the x exp -1/2 dependence of the classical similarity solution, where x is the streamwise distance. Cylindrical wick geometries yield enhanced burning rates over planar wicks. For the case of mixed convective burning, the burning rate results approach either the forced or natural convective burning limits as ambient streamwise velocity is changed. Critical Richardson numbers specifying these burning limits are determined for a given condition. Reducing gravity results in a lower burning rate because the influence of natural convection is diminished. Under reduced gravity of 1/1000 of the sea-level value, mixed convective burning nearly resembles forced convection.

Damaso, R. C.↗

Soot Volume Fraction Imaging

A new technique is described for the full-field determination of soot volume fractions via laser extinction measurements. This technique differs from previously reported point-wise methods in that a two-dimensional array (i.e., image) of data is acquired simultaneously. In this fashion, the net data rate is increased, allowing the study of time-dependent phenomena and the investigation of spatial and temporal correlations. A telecentric imaging configuration is employed to provide depth-invariant magnification and to permit the specification of the collection angle for scattered light. To improve the threshold measurement sensitivity, a method is employed to suppress undesirable coherent imaging effects. A discussion of the tomographic inversion process is provided, including the results obtained from numerical simulation. Results obtained with this method from an ethylene diffusion flame are shown to be in close agreement with those previously obtained by sequential point-wise interrogation.

Greenberg, Paul S.↗

Detailed Multidimensional Simulations of the Structure and Dynamics of Flames

Numerical simulations in which the various physical and chemical processes can be independently controlled can significantly advance our understanding of the structure, stability, dynamics and extinction of flames. Therefore, our approach has been to use detailed time-dependent, multidimensional, multispecies numerical models to perform carefully designed computational experiments of flames on Earth and in microgravity environments. Some of these computational experiments are complementary to physical experiments performed under the Microgravity Program while others provide a fundamental understanding that cannot be obtained from physical experiments alone. In this report, we provide a brief summary of our recent research highlighting the contributions since the previous microgravity combustion workshop. There are a number of mechanisms that can cause flame instabilities and result in the formation of dynamic multidimensional structures. In the past, we have used numerical simulations to show that it is the thermo-diffusive instability rather than an instability due to preferential diffusion that is the dominant mechanism for the formation of cellular flames in lean hydrogen-air mixtures. Other studies have explored the role of gravity on flame dynamics and extinguishment, multi-step kinetics and radiative losses on flame instabilities in rich hydrogen-air flames, and heat losses on burner-stabilized flames in microgravity. The recent emphasis of our work has been on exploring flame-vortex interactions and further investigating the structure and dynamics of lean hydrogen-air flames in microgravity. These topics are briefly discussed after a brief discussion of our computational approach for solving these problems.

Patnaik, G.↗

Modeling and Validation of Damped Plexiglas Windows for Noise Control

Windows are a significant path for structure-borne and air-borne noise transmission in general aviation aircraft. In this paper, numerical and experimental results are used to evaluate damped plexiglas windows for the reduction of structure-borne and air-borne noise transmitted into the interior of an aircraft. In contrast to conventional homogeneous windows, the damped plexiglas windows were fabricated using two or three layers of plexiglas with transparent viscoelastic damping material sandwiched between the layers. Transmission loss and radiated sound power measurements were used to compare different layups of the damped plexiglas windows with uniform windows of the same nominal thickness. This vibro-acoustic test data was also used for the verification and validation of finite element and boundary element models of the damped plexiglas windows. Numerical models are presented for the prediction of radiated sound power for a point force excitation and transmission loss for diffuse acoustic excitation. Radiated sound power and transmission loss predictions are in good agreement with experimental data. Once validated, the numerical models were used to perform a parametric study to determine the optimum configuration of the damped plexiglas windows for reducing the radiated sound power for a point force excitation.

Buehrle, Ralph D.↗

New Species of Fire Discovered: Fingering Flamelets Form a Dynamic Population

Poets and artists have long used fire as a metaphor for life. At the NASA Glenn Research Center, recent experiments in a subcritical Rayleigh number flow channel demonstrated that this analogy holds up surprisingly well when tools developed to characterize a biological population are applied to a class of fire that occurs in near-extinction, weakly convective environments (such as microgravity) or in vertically confined spaces (such as our apparatus). Under these conditions, the flame breaks into numerous 'flamelets" that form a Turing-type reaction-diffusion fingering pattern as they spread across the fuel. It is standard practice on U.S. spacecraft for the astronaut crew to turn off the ventilation to help extinguish a fire, both to eliminate the fresh oxygen supply and to reduce the distribution of the smoke. When crew members think that the fire is fully extinguished, they reactivate the ventilation system to clear the smoke. However, some flamelets can survive, and our experiments have demonstrated that flamelets quickly grow into a large fire when ventilation increases.

Olson, Sandra L.↗