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.

At least 127 records · Page 7

Direct Coupling Method for Time-Accurate Solution of Incompressible Navier-Stokes Equations

A noniterative finite difference numerical method is presented for the solution of the incompressible Navier-Stokes equations with second order accuracy in time and space. Explicit treatment of convection and diffusion terms and implicit treatment of the pressure gradient give a single pressure Poisson equation when the discretized momentum and continuity equations are combined. A pressure boundary condition is not needed on solid boundaries in the staggered mesh system. The solution of the pressure Poisson equation is obtained directly by Gaussian elimination. This method is tested on flow problems in a driven cavity and a curved duct.

Soh, Woo Y.↗

The structure of premixed particle-cloud flames

The structure of premixed flames propagating in combustible systems, containing uniformly distributed volatile fuel particles, in an oxidizing gas mixture, is analyzed. It is presumed that the fuel particles vaporize first to yield a gaseous fuel of known chemical structure, which is subsequently oxidized in the gas phase. The analysis is performed in the asymptotic limit, where the value of the characteristic Zeldovich number, based on the gas-phase oxidation of the gaseous fuel is large, and for values of phi(u) greater than or equal to 1.0, where phi(u) is the equivalence ratio based on the fuel available in the fuel particles. The structure of the flame is presumed to consist of a preheat vaporization zone where the rate of the gas-phase chemical reaction is small, a reaction zone where convection and the rate of vaporization of the fuel particles are small and a convection zone where diffusive terms in the conservation equations are small. For given values phi(u) the analysis yields results for the burning velocity and phi(g) where phi(g) is the effective equivalence ratio in the reaction zone. The analysis shows that even though phi(u) greater than or equal to 1.0, for certain cases the calculated value of phi(g) is less than unity. This prediction is in agreement with experimental observations.

Seshadri, K.↗

Frozen-Flux Modelling for Epochs 1915 and 1980

The frozen-flux hypothesis for the Earth's liquid core assumes that convective terms dominate diffusive terms in the induction equation governing the behaviour of the magnetic field at the surface of the core. While highly plausible on the basis of estimates of physical parameters, the hypothesis has been questioned. To study this hypothesis, we improve the method which tests the consistency of magnetic observations with the hypothesis by constructing simple, flux-conserving core-field models fitting the data at pairs of epochs. We introduce a new approach that fixes the patch configurations at each of the two epochs before inversion, so that each configuration is consistent with its respective data set but possesses the same patch topology. We expand upon the inversion algorithm, using quadratic programming to maintain the proper flux sign within patches; the modelling calculations are also extended to include data types that depend non-linearly on the model. Every test of a hypothesis depends on the characterization of the observational uncertainties; we undertake a thorough review of this question. For main-field models, the primary source of uncertainty comes from the crustal field. We base our analysis on statistical models of the crustal magnetization, adjusted to bring it into better conformity with our data set. The noise model permits us to take into account the correlations between the measurements and requires that a different weighting be given to horizontal and vertical components. It also indicates that the observations should be fit more closely than has been the practice heretofore. We apply the revised method to Magsat data from 1980 and survey and observatory data from 1915.5, two data sets believed to be particularly difficult to reconcile with the frozen-flux hypothesis. We compute a pair of simple, flux-conserving models that fit the averaged data from each epoch. We therefore conclude that present knowledge of the geomagnetic fields of 1980 and 1915.5 is consistent with the frozen-flux hypothesis.

OBrien, Michael S.↗

Model for premixed turbulent flames

Flame speed and mean degree of reaction profiles are determined for an ensemble mean one-dimensional premixed turbulent flame. Elementary closure models are employed to treat turbulence closure problems arising from the nonlinear convective and chemical reaction rate terms in the time-dependent fuel conservation equation; a constant eddy diffusivity is useed to model correlation terms derived from the convective term. An expression developed for fuel consumption is solved to obtain the turbulent flame speed and the variation of mean fuel mass fraction across the flame. Comparison of the numerical analysis with experimental data indicates the importance of the quenching of the combustion reaction in high-wavenumber eddies.

Gouldin, F. C.↗

A Probabilistic Scheme for Semilinear Nonlocal Diffusion Equations with Volume Constraints

This work presents a probabilistic scheme for solving semilinear nonlocal diffusion equations with volume constraints and integrable kernels. The nonlocal model of interest is defined by a time-dependent semilinear partial integro-differential equation (PIDE), in which the integro-differential operator consists of both local convection-diffusion and nonlocal diffusion operators. Here, our numerical scheme is based on the direct approximation of the nonlinear Feynman–Kac formula that establishes a link between nonlinear PIDEs and stochastic differential equations. The exploitation of the Feynman–Kac representation avoids solving dense linear systems arising from nonlocal operators. Compared with existing stochastic approaches, our method can achieve first-order convergence after balancing the temporal and spatial discretization errors, which is a significant improvement of existing probabilistic/stochastic methods for nonlocal diffusion problems. Error analysis of our numerical scheme is established. The effectiveness of our approach is shown in two numerical examples. The first example considers a three-dimensional nonlocal diffusion equation to numerically verify the error analysis results. The second example presents a physics problem motivated by the study of heat transport in magnetically confined fusion plasmas.

97 MATHEMATICS AND COMPUTING↗

Numerical vorticity capturing for vortex-solid body interaction problems

A new numerical method based on adding a term to the Euler/Navier-Stokes equations has been demonstrated it can effectively treat vortex-dominated flows using low-order numerical schemes and coarse grids. The numerical diffusion introduced by convection schemes can be eliminated by modifying the velocity field as a result of solving the modified flow governing equations. The modification of the velocity field, which conserves the total vorticity, essentially convects the vorticity toward its local extreme to offset the numerical diffusion. This method preserves the vortex structure even when vortices travel on coarse grid region. The method is implemented with two distinctive codes: a Navier-Stokes flow solver based on vorticity-velocity formulation and an Euler code based on velocity-pressure formulation. Problems presented in this paper include airfoil dynamic stall, vortex airfoil interaction, and vortex-fuselage interactions.

Wang, Clin M.↗

Badhwar-O'Neill 2011 Galactic Cosmic Ray Model Update and Future Improvements

The Badhwar-O'Neill Galactic Cosmic Ray (GCR) Model based on actual GCR measurements is used by deep space mission planners for the certification of microelectronic systems and the analysis of radiation health risks to astronauts in space missions. The BO GCR Model provides GCR flux in deep space (outside the earth's magnetosphere) for any given time from 1645 to present. The energy spectrum from 50 MeV/n - 20 GeV/n is provided for ions from hydrogen to uranium. This work describes the most recent version of the BO GCR model (BO'11). BO'11 determined the GCR flux at a given time applying an emperical time delay function to past sunspot activity. We describe the GCR measurement data used in the BO'11 update - modern data from BESS, PAMELA, CAPRICE, and ACE emphasized more than the older balloon data used for the previous BO model (BO'10). We look at the GCR flux for the last 24 solar minima and show how much greater the flux was for the cycle 24 minimum in 2010. The BO'11 Model uses the traditional, steady-state Fokker-Planck differential equation to account for particle transport in the heliosphere due to diffusion, convection, and adiabatic deceleration. It assumes a radially symmetrical diffusion coefficient derived from magnetic disturbances caused by sunspots carried outward by a constant solar wind. A more complex differential equation is now being tested to account for particle transport in the heliosphere in the next generation BO model. This new model is time-dependent (no longer a steady state model). In the new model, the dynamics and anti-symmetrical features of the actual heliosphere are accounted for so emperical time delay functions will no longer be required. The new model will be capable of simulating the more subtle features of modulation - such as the Sun's polarity and modulation dependence on gradient and curvature drift. This improvement is expected to significantly improve the fidelity of the BO GCR model. Preliminary results of its performance will be presented.

O'Neill, Pat M.↗

Badhwar-O'Neill 2011 Galactic Cosmic Ray Model Update and Future Improvements

The Badhwar-O'Neill Galactic Cosmic Ray (GCR) Model based on actual GR measurements is used by deep space mission planners for the certification of micro-electronic systems and the analysis of radiation health risks to astronauts in space missions. The BO GCR Model provides GCR flux in deep space (outside the earth's magnetosphere) for any given time from 1645 to present. The energy spectrum from 50 MeV/n-20 GeV/n is provided for ions from hydrogen to uranium. This work describes the most recent version of the BO GCR model (BO'11). BO'11 determines the GCR flux at a given time applying an empirical time delay function to past sunspot activity. We describe the GCR measurement data used in the BO'11 update - modern data from BESS, PAMELA, CAPRICE, and ACE emphasized for than the older balloon data used for the previous BO model (BO'10). We look at the GCR flux for the last 24 solar minima and show how much greater the flux was for the cycle 24 minimum in 2010. The BO'11 Model uses the traditional, steady-state Fokker-Planck differential equation to account for particle transport in the heliosphere due to diffusion, convection, and adiabatic deceleration. It assumes a radially symmetrical diffusion coefficient derived from magnetic disturbances caused by sunspots carried onward by a constant solar wind. A more complex differential equation is now being tested to account for particle transport in the heliosphere in the next generation BO model. This new model is time-dependent (no longer a steady state model). In the new model, the dynamics and anti-symmetrical features of the actual heliosphere are accounted for so empirical time delay functions will no longer be required. The new model will be capable of simulating the more subtle features of modulation - such as the Sun's polarity and modulation dependence on the gradient and curvature drift. This improvement is expected to significantly improve the fidelity of the BO GCR model. Preliminary results of its performance will be presented.

O'Neill, Pat M.↗

An improvement of convection fidelity in Euler calculations

A new solution procedure is presented to solve the Euler equations for steady, compressible, rotational, inviscid flows. The approach is aimed at achieving real inviscid solutions in Euler calculations by eliminating numerical diffusion. The variables in the Euler equations are divided into elliptic and convective quantities, using the Clebsch velocity decomposition. The convective quantities are then transported without numerical contamination using an efficient convection operator, while the elliptic quantities are integrated with a relaxation procedure. This approach provides a generalization of the full potential formulation to rotation Euler physics by allowing variations of convective quantities. Results are demonstrated for several transonic flows in two dimensions.

Lee, K. D.↗

Formation and dynamics of large-scale magnetic structures in the ionosphere of Venus

The formation and dynamics of large-scale magnetic structures in the ionosphere of Venus are examined. It is shown that such structures must be the result of steady state convection of interplanetary field lines into the ionosphere by the small amount of solar wind plasma (less than or approximately equal to 1-5 percent) absorbed by the planetary atmosphere below the ionopause, rather than isolated remnants of large fields persisting for long periods without connection to the solar wind induced current and convection pattern. In particular, it is demonstrated that the magnetic diffusion of such structures would result in their dissipation with time scales of 1-10 min, if they were not steady state structures in convective and diffusive equilibriuim. It is shown that the equations governing the diffusion of these magnetic structures are similar to those governing diffusion of a gas out of an enclosed chamber with a porous wall, and a simple analog is illustrated. The application of these results to magnetic fields of astrophysical plasmas is discussed.

Cloutier, P. A.↗

Features of a reattaching turbulent shear layer in divergent channel flow

Experimental data have been obtained in an incompressible turbulent flow over a rearward-facing step in a diverging channel flow. Mean velocities, Reynolds stresses, and triple products that were measured by a laser Doppler velocimeter are presented for two cases of tunnel wall divergence. Eddy viscosities, production, convection, turbulent diffusion, and dissipation (balance of kinetic energy equation) terms are extracted from the data. These data are compared with various eddy-viscosity turbulence models. Numerical calculations incorporating the k-epsilon and algebraic-stress turbulence models are compared with the data. When determining quantities of engineering interest, the modified algebraic-stress model (ASM) is a significant improvement over the unmodified ASM and the unmodified k-epsilon model; however, like the others, it dramatically overpredicts the experimentally determined dissipation rate.

Driver, D. M.↗

Experimental study of a three-dimensional, shear-driven, turbulent boundary layer using a three-dimensional laser Doppler velocimeter

Experimental data have been obtained in a three-dimensional, shear-driven, turbulent boundary layer subjected to sudden transverse strain. Measurements made with a newly developed, three-dimensional laser Doppler velocimeter include all components of mean flow, turbulent Reynolds stresses, and triple-product correlations. Eddy viscosities, production, convection, turbulent diffusion, and dissipation (balance of kinetic energy equation) are extracted from the data. It was shown experimentally that eddy viscosity is strongly direction-dependent (not a scalar), indicating a need for full Reynolds stress modeling. Turbulence models for eddy viscosity, stress diffusion, pressure rate-of-strain, and anisotropy of Reynolds stresses are tested using the data. Calculations using a Reynolds stress equation model are compared with the data.

Driver, D. M.↗

The telegraph equation in charged particle transport

We present a new derivation of the telegraph equation which modifies its coefficients. First, an infinite order partial differential equation is obtained for the velocity space solid angle-averaged phase-space distribution of particles which underwent at least a few collisions. It is shown that, in the lowest order asymptotic expansion, this equation simplifies to the well-known diffusion equation. The second-order asymptotic expansion for isotropic small-angle scattering results in a modified telegraph equation with a signal propagation speed of v(5/11) exp 1/2 instead of the usual v/3 exp 1/2. Our derivation of a modified telegraph equation follows from an expansion of the Boltzmann equation in the relevant smallness parameters and not from a truncation of an eigenfunction expansion. This equation is consistent with causality. It is shown that, under steady state conditions in a convecting plasma, the telegraph equation may be regarded as a diffusion equation with a modified transport coefficient, which describes a combination of diffusion and cosmic-ray inertia.

Gombosi, T. I.↗

Quasi-linear pitch angle and energy diffusion of pickup ions near Comet Halley

The process of ion pickup in the environment of Halley's comet is studied in order to see if velocity diffusion driven by the observed level of turbulence can explain the observed development of the implanted ion distribution. The theoretical description used is based on a quasi-linear approach and considers the implantation and transport of cometary ions along solar wind flow lines. To make such a study requires some way of extrapolating the measurements on the Giotto trajectory into the upstream region; models for mass loading and turbulence are used. A simplified kinetic equation describing the source, convection, and quasi-linear velocity diffusion of the heavy cometary ions is solved numerically along flow lines parallel to the sun-comet line. Full two-dimensional (pitch angle and velocity) distributions are obtained at positions along the Giotto trajectory, which are compared with measurements. This study finds that quasi-linear theory, with the empirical model for the observed turbulence level, produces the right order of pitch angle diffusion.

Huddleston, D. E.↗

A method for calculating unsteady turbulent boundary layers in two- and three-dimensional flows.

The governing equations for an unsteady turbulent boundary layer on a swept infinite cylinder, composed of a continuity equation, a pair of momentum equations and a pair of rate equations for the shear stress, based on the turbulent kinetic-energy equation, are solved numerically. These rate equations for the shear-stress express a balance between the convection, production, dissipation and diffusion of the turbulent shear stress components. Due to the physical model of turbulence employed, this system of equations is hyperbolic. An explicit, second-order accurate, conditionally stable difference scheme for this system, together with appropriate boundary and initial conditions, is formulated and developed into a computer program. Calculations for oscillating free stream flows with no pressure gradient show significant unsteady effects on the turbulent boundary layer.

Singleton, R. E.↗

Development of a New Membrane Casting Apparatus for Studying Macrovoid Defects in Low-G

A new membrane-casting apparatus is developed for studying macrovoid defects in polymeric membranes made by the wet- and dry-casting process in low-gravity. Macrovoids are large (10-50 micron), open cavities interspersed among the smaller pores in the substructure under the gelled skin surface layer of the cast membrane. Although their occurrence is considered endemic to the wet- and dry-casting process since they can lead to compaction or skin rupture in the membrane process, recent studies suggest several useful applications such as transdermal and osmotic drug delivery systems, miniature bioreactors, etc. However, lack of knowledge about the macrovoid formation mechanism is an obstacle to further development of applications using them. An on-going debate is the role of the surface-tension-driven solutocapillary convection during macrovoid formation. The rapid growth of macrovoids within 1-5 seconds and the high polymer concentration in and near macrovoids make it difficult to explain the mechanism of macrovoid growth by diffusion alone, which is the widely accepted hypothesis proposed by Reuvers et al. The hypothesis advanced by our research group can explain this rapid growth via a mechanism that involves diffusion from the casting solution in the meta-stable region to the macrovoid enhanced by solutocapillary convection induced by the steep nonsolvent concentration gradient in the vicinity of the macrovoid. Since macrovoid growth is hypothesized to be the interplay of a solutocapillary-induced driving force counteracted by viscous drag and buoyancy, eliminate the latter provides a means for testing this hypothesis. Moreover, free convection mass transfer in the nonsolvent immersion bath used to cause phase-separation in membrane casting complicates developing a model for both the wet-casting process and macrovoid growth. The low-g environment minimizes gravitationally induced free convection thereby permitting a tractable solution to the ternary diffusion equations that characterize membrane formation. NASA's Parabolic Flight Research Aircraft provides a small window of low-g (approximately 25 s) that can be used to study macrovoid development in both wet- and dry-cast membranes if an appropriate casting apparatus is used. This casting apparatus should be able to cast the membrane in both low- and high-g in a manner so that essential one-dimensional mass transfer conditions are achieved to insure lateral uniformity in the membrane. The apparatus used in previous research on membrane casting in low-gravity was operated with the plunger driven mechanism. The spring-loaded plunger pushes the bottom block containing the polymer casting solution well directly under the absorbent chamber located in the upper stationary block. However, membranes made via this casting apparatus often displayed lateral nonuniformities that precluded obtaining quantitative information on the macrovoid growth process. Thus, it was necessary to determine the reason for these structural irregularities observed in the low-g casting apparatus. Both experimental as well as computer simulation studies of the low-g casting apparatus established that the impulsive action of the plunger caused the undesired structural nonuniformities. The simulation results showed that the width-to-depth aspect ratio of the shallow well that contains the casting solution in this apparatus was not an important factor in minimizing this problem. Even for a 40:1 (width : depth) aspect ratio, any convection induced by the horizontal motion of the interface of the casting solution will be damped out within 6.25x10(exp 4) seconds. However, the experimental studies revealed that the impulsive motion of the plunger caused a 'sloshing' of the casting solution that had to be eliminated. Therefore, the plungerdriven mechanism was changed to a cam-driven mechanism that did not cause any impulsive motion of the casting solution. Other refinements to this new membrane-casting apparatus include provision for removing the membranes from the casting wells in a less destructive manner. This was accomplished by using a slit geometry for the casting well that permitted disassembly for removal of the cast membrane. The materials used in the construction of this casting apparatus were chosen to insure wetting at the side walls and to maintain precise control of the thickness of the polymer solution in the casting well. An additional provision in this new casting apparatus is the ability to carry out both wet- as well as dry-casting. As such, this apparatus permitted the first studies of the wet-casting of polymeric membranes in low-g. Both wet- and dry-casting experiments on NASA's KC-135 research aircraft employing this new membrane-casting apparatus are scheduled in July 2002. The morphology of the resulting membranes will be characterized using an environmental scanning electron microscope (ESEM). The results of these low-g studies will be reported later.

Lee, Hanyong↗

Eulerian Mapping Closure Approach for Probability Density Function of Concentration in Shear Flows

The Eulerian mapping closure approach is developed for uncertainty propagation in computational fluid mechanics. The approach is used to study the Probability Density Function (PDF) for the concentration of species advected by a random shear flow. An analytical argument shows that fluctuation of the concentration field at one point in space is non-Gaussian and exhibits stretched exponential form. An Eulerian mapping approach provides an appropriate approximation to both convection and diffusion terms and leads to a closed mapping equation. The results obtained describe the evolution of the initial Gaussian field, which is in agreement with direct numerical simulations.

He, Guowei↗