Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “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 937 records · Page 52

Multigrid method with weighted mean scheme

Multigrid algorithms based on the weighted mean scheme are developed for the solution of the two dimensional incompressible Navier-Stokes equations. They are applied to two typical problems encountered in engineering applications, namely, the convection diffusion problem of the Benard convection cell, and the driven cavity problem. An analysis of the smoothing rates and stability is given. The efficiency of the multigrid method is investigated.

Lustman, L. R.↗

The treatment of convected vortices in compressible potential flow

A method is described for incorporating line vortices into the three dimensional compressible potential flow equation. A modified Biot-Savart law is used to compute a vortical velocity field, which is added to the gradient of the potential to form a total velocity. A rapidly converging approximate factorization (AFZ) scheme is then used to compute a potential such that the modified potential flow equation as well as the appropriate boundary conditions, based on total velocity, are satisfied. As part of a coupled iteration procedure, the positions of the line vortices are computed so that they convect with the total flow. The method is used to compute the field due to a single line vortex convecting past a wing. This represents an approximation of the effect of a canard or other lifting surface ahead of the wing, which sheds a tip vortex. It is seen that the flow field is substantially modified by the passage of the vortex. Unlike Euler equation schemes, which are also used to compute these flows, the solutions exhibit no numerical diffusion: The convected vortices retain their initial upstream width. Euler solutions, on the other hand, involve a vorticity which is numerically convected in an Eulerian frame and, unless extensive adaptive grid refinement is used they result in vortices with spread as they convect. Also, the potential flow method requires approximately two orders of magnitude less computing time and much less computer storage than the Euler methods.

Steinhoff, J.↗

Viscous, resistive magnetohydrodynamic stability computed by spectral techniques

Expansions in Chebyshev polynomials are used to study the linear stability of one-dimensional magnetohydrodynamic quasi-equilibria, in the presence of finite resistivity and viscosity. The method is modeled on the one used by Orszag in accurate computation of solutions of the Orr-Sommerfeld equation. Two Reynolds-like numbers involving Alfven speeds, length scales, kinematic viscosity, and magnetic diffusivity govern the stability boundaries, which are determined by the geometric mean of the two Reynolds-like numbers. Marginal stability curves, growth rates versus Reynolds-like numbers, and growth rates versus parallel wave numbers are exhibited. A numerical result that appears general is that instability has been found to be associated with inflection points in the current profile, though no general analytical proof has emerged. It is possible that nonlinear subcritical three-dimensional instabilities may exist, similar to those in Poiseuille and Couette flow.

Dahlburg, R. B.↗

Numerical Descriptions of Cosmic-Ray Transport

The behavior of energetic particles in the solar system is described by a well known Fokker-Planck equation. Although analytic methods yield insight into the nature of its solutions, especially in the diffusion regime, calculations that go beyond diffusion are very complicated. The reliability of these calculations is of concern, because numerical methods are notorious for their errors and artifacts. The well known Milne problem of classical transport theory was analyzed with the aid of three different numerical methods. These are: (1) The method of eigenfunctions in which the distribution function is approximated by a sum of eigenfunctions of the scattering operator, (2) Numerical solutions of a finite difference aquation; and (3) Direct simulation of the scattering and streaming of individual particles with the aid of Monte Carlo methods.

Earl, J. A.↗

Numerical calculations of turbulent reacting flow in a gas-turbine combustor

A numerical study for confined, axisymmetrical, turbulent diffusion flames is presented. Local mean gas properties are predicted by solving the appropriate conservation equations in the finite difference form with the corresponding boundary conditions. The k-epsilon two-equation turbulence model is employed to describe the turbulent nature of the flow. A two-step kinetic model is assumed to govern the reaction mechanism. The finite reaction rate is the smaller of an Arrhenius type of reaction rate and a modified version of eddy-breakup model. Reasonable agreement is observed between calculations and measurements, but to obtain better agreement, more work is needed on improvements of the above mathematical models. However, the present numerical study offers an improvement in the analysis and design of the gas turbine combustors.

Lin, Chin-Shun↗

A theory for turbulent curved wall jets

A simple theoretical model is proposed to describe the flow of a turbulent wall jet along a curved surface into a quiescent atmosphere. An integral method is used to solve the momentum equation and identifies three contributions to the spreading rate of the wall jet: (1) turbulent diffusion in the wall jet; (2) wall curvature; and (3) rate of change of wall curvature. Closed from approximate solutions are found for the case of a plane wall, a circular cylinder, and a logarithmic spiral surface. Comparison with experimental data for these three cases is made showing good agreement.

Roberts, Leonard↗

Simulations of high Mach number perpendicular shocks with resistive electrons

A simulation code which models the ions as microparticles and the electrons as a resistive massless fluid is employed to study the structure of high Mach number perpendicular shocks. It is found that stable stationary shock solutions can be obtained for Alfven Mach numbers (M sub A) between 5 and 60 for upstream plasmas where the ratio of the plasma pressure to the magnetic pressure is 1, providing that the upstream resistive diffusion length is much smaller than the ion inertial length. For much larger resistive diffusion lengths, the magnetic field overshoot is damped, and the imbalance in the electron momentum equation results in a periodic fluctuation of the fraction of reflected ions. In the limit of M sub A of less than 10, the magnetic overshoot and the fraction of reflected ions increase with increasing M sub A, while at higher Mach numbers the fraction of reflected ions peaks at about 40 percent and the magnetic field overshoot increases at a much slower rate. Electron inertial effects are also considered.

Quest, K. B.↗

An explicit Runge-Kutta method for turbulent reacting flow calculations

The paper presents a numerical method for the solution of the conservation equations governing steady, reacting, turbulent viscous flow in two-dimensional geometries, in both Cartesian and axisymmetric coordinates. These equations are written in Favre-averaged form and closed with a first order model. A two-equation K-epsilon model, where low Reynolds number and compressibility effects are included, and a modified eddy-break up model are used to simulate fluid mechanics turbulence, chemistry and turbulence-combustion interaction. The solution is obtained by using a pseudo-unsteady method with improved perturbation propagation properties. The equations are discretized in space by using a finite volume formulation. An explicit multi-stage dissipative Runge-Kutta algorithm is then used to advance the flow equations in the pseudo-time. The method is applied to the computation of both diffusion and premixed turbulent reacting flows. The computed temperature distributions compare favorably with experimental data.

Boretti, A. A.↗

Computations of the three-dimensional flow and heat transfer within a coolant passage of a radial turbine blade

A numerical code is developed for computing three-dimensional, turbulent, compressible flow within coolant passages of turbine blades. The code is based on a formulation of the compressible Navier-Stokes equations in a rotating frame of reference in which the velocity dependent variable is specified with respect to the rotating frame instead of the inertial frame. The algorithm employed to obtain solutions to the governing equation is a finite-volume LU algorithm that allows convection, source, as well as diffusion terms to be treated implicitly. In this study, all convection terms are upwind differenced by using flux-vector splitting, and all diffusion terms are centrally differenced. This paper describes the formulation and algorithm employed in the code. Some computed solutions for the flow within a coolant passage of a radial turbine are also presented.

Shih, T. I.-P.↗

Cosmic-ray heating of cooling flows - A critical analysis

It is shown that a combination of MHD wave-mediated cosmic ray heating and thermal conduction could balance cooling in intracluster media and substantially reduce the rate of inflow. The appropriate system of steady state equations is solved, including a new self-consistent formulation for the cosmic-ray diffusivity. Models which can produce substantial positive temperature gradients in static configurations are found when conduction is reduced by a factor of 10 or more. These models have too-flat thermal pressure profiles compared with observations. It is found that cosmic-ray heating is unlikely either to stabilize positive density perturbations against condensation or to contribute appreciably to the powering of the optical filaments.

Loewenstein, Michael↗

A numerical analysis of high-temperature heat pipe startup from the frozen state

Continuum and rarefied vapor flows co-exist along the heat pipe length for most of the startup period. A two-region model is proposed in which the vapor flow in the continuum region is modeled by the compressible Navier-Stokes equations, and the vapor flow in the rarefied region is simulated by a self-diffusion model. The two vapor regions are linked with appropriate boundary conditions, and heat pipe wail, wick, and vapor flow are solved as a conjugate problem. The numerical solutions for the entire heat pipe startup process from the frozen state are compared with the corresponding experimental data with good agreement.

Cao, Y.↗

Boundary layer blockage in expansion tube nozzles

The results of a first order perfect gas correction for the effects of the boundary layer formation within expansion tubes with nozzles are presented. The analytical model developed to describe the boundary layer formation within the expansion tube and an expansion nozzle located at the end of the acceleration tube is based on the Karman integral equations. The results of this analytical model are compared with experimental data from an expansion diffuser. The model provides a useful tool for the preliminary design of nozzles for such facilities.

Sudnitsin, Olga↗

Self-Consistent Simulation of the Brownian Stage of Dust Growth

It is a widely accepted view that in proto-planetary accretion disks the collision and following sticking of dust particles embedded in the gas eventually leads to the formation of planetesimals (coagulation). For the smallest dust grains, Brownian motion is assumed to be the dominant source of their relative velocities leading to collisions between these dust grains. As the dust grains grow they eventually couple to the turbulent motion of the gas which then drives the coagulation much more efficiently. Many numerical coagulation simulations have been carried out to calculate the fractal dimension of the aggregates, which determines the duration of the ineffective Brownian stage of growth. Predominantly on-lattice and off-lattice methods were used. However, both methods require simplification of the astrophysical conditions. The aggregates found by those methods had a fractal dimension of approximately 2 which is equivalent to a constant, mass-independent friction time. If this value were valid for the conditions in an accretion disk, this would mean that the coagulation process would finally 'freeze out' and the growth of a planetesimal would be impossible within the lifetime of an accretion disk. In order to investigate whether this fractal dimension is model independent, we simulate self-consistently the Brownian stage of the coagulation by an N-particle code. This method has the advantage that no further assumptions about homogeneity of the dust have to be made. In our model, the dust grains are considered as aggregates built up of spheres. The equation of motion of the dust grains is based on the probability density for the diffusive transport within the gas atmosphere. Because of the very low number density of the dust grains, only 2-body-collisions have to be considered. As the Brownian stage of growth is very inefficient, the system is to be simulated over long periods of time. In order to find close particle pairs of the system which are most likely to undergo a collision, we use a particle-in-cell (PIC) method for the early stages of the simulation where the system is still very homogeneous and a tree method later when the particles are more clustered.

Kempf, S.↗

Tape-Drop Transient Model for In-Situ Automated Tape Placement of Thermoplastic Composites

Composite parts of nonuniform thickness can be fabricated by in-situ automated tape placement (ATP) if the tape can be started and stopped at interior points of the part instead of always at its edges. This technique is termed start/stop-on-the-part, or, alternatively, tape-add/tape-drop. The resulting thermal transients need to be managed in order to achieve net shape and maintain uniform interlaminar weld strength and crystallinity. Starting-on-the-part has been treated previously. This paper continues the study with a thermal analysis of stopping-on-the-part. The thermal source is switched off when the trailing end of the tape enters the nip region of the laydown/consolidation head. The thermal transient is determined by a Fourier-Laplace transform solution of the two-dimensional, time-dependent thermal transport equation. This solution requires that the Peclet number Pe (the dimensionless ratio of inertial to diffusive heat transport) be independent of time and much greater than 1. Plotted isotherms show that the trailing tape-end cools more rapidly than the downstream portions of tape. This cooling can weaken the bond near the tape end; however the length of the affected region is found to be less than 2 mm. To achieve net shape, the consolidation head must continue to move after cut-off until the temperature on the weld interface decreases to the glass transition temperature. The time and elapsed distance for this condition to occur are computed for the Langley ATP robot applying PEEK/carbon fiber composite tape and for two upgrades in robot performance. The elapsed distance after cut-off ranges from about 1 mm for the present robot to about 1 cm for the second upgrade.

Costen, Robert C.↗

Innermost Van Allen Radiation Belt for High Energy Protons at Saturn

The high energy proton radiation belts of Saturn are energetically dominated by the source from cosmic ray albedo neutron decay (CRAND), trapping of protons from beta decay of neutrons emitted from galactic cosmic ray nuclear interactions with the main rings. These belts were originally discovered in wide gaps between the A-ring, Janus/Epimetheus, Mimas, and Enceladus. The narrow F and G rings significant affected the CRAND protons but did not produce total depletion. Voyager 2 measurements subsequently revealed an outermost CRAND proton belt beyond Enceladus. Although the source rate is small, the trapping times limited by radial magnetospheric diffusion are very long, about ten years at peak measured flux inwards of the G ring, so large fluxes can accumulate unless otherwise limited in the trapping region by neutral gas, dust, and ring body interactions. One proposed final extension of the Cassini Orbiter mission would place perikrone in a 3000-km gap between the inner D ring and the upper atmosphere of Saturn. Experience with CRAND in the Earth's inner Van Allen proton belt suggests that a similar innermost belt might be found in this comparably wide region at Saturn. Radial dependence of magnetospheric diffusion, proximity to the ring neutron source, and northward magnetic offset of Saturn's magnetic equator from the ring plane could potentially produce peak fluxes several orders of magnitude higher than previously measured outside the main rings. Even brief passes through such an intense environment of highly penetrating protons would be a significant concern for spacecraft operations and science observations. Actual fluxes are limited by losses in Saturn's exospheric gas and in a dust environment likely comparable to that of the known CRAND proton belts. The first numerical model of this unexplored radiation belt is presented to determine limits on peak magnitude and radial profile of the proton flux distribution.

Cooper, John F.↗

The role of gravity wave generated advection and diffusion in transport of tracers in the mesosphere

The roles of advection and turbulent diffusion in controlling tracer transport in the meridional plane of the mesosphere are compared using a version of the global primitive equation model for the middle atmosphere developed by Holton and Wehrbein (1980). The results of model sensitivity tests demonstrated that transport in the mesosphere for long-lived species with large-scale heights, such as H2O, CO, and NO(x), is primarily an advective process; on a global scale, small-scale turbulent diffusion produced by gravity wave breaking plays a secondary role. It is concluded that, for such tracers, the eddy diffusion coefficient used in one-dimensional chemical models of the mesosphere should be regarded as a parameterization of the transport by the meridional circulation, not as a parameterization of mixing by small-scale diffusion.

Holton, James R.↗

Thermofluid analysis of the SSME preburner using a gas-gas diffusion model for oxygen and hydrogen combustion at supercritical pressures

The paper discusses the thermofluid analysis of the Space Shuttle Main Engine (SSME) fuelside preburner. The governing equations have been solved numerically to predict flow, heat transfer, mixing, and combustion. A two-fluid approach is adopted in which oxygen is regarded as one fluid and hydrogen is regarded as the other fluid. The chemical kinetics is assumed to be very fast so that combustion is primarily controlled by the rate of mixing between oxygen and hydrogen. The preburner pressure is much greater than the critical pressures of oxygen and hydrogen; hence, a gas-gas diffusion model (rather than an evaporation model) has been developed to compute the rate of interphase mixing. Empirical correlations have been incorporated to account for the effect of slip on the interphase exchange. A sensitivity study has been performed with various model parameters. It is observed that the model can predict possibility of incomplete combustion and local regions of high temperatures under steady operating conditions. Some of these anomalies have been observed in actual tests, and the numerical model is useful for understanding possible causes and remedies. At least some measurements are needed for quantitative verification of the model.

Prakash, C.↗

A two-dimensional numerical simulation of a supersonic, chemically reacting mixing layer

Research has been undertaken to achieve an improved understanding of physical phenomena present when a supersonic flow undergoes chemical reaction. A detailed understanding of supersonic reacting flows is necessary to successfully develop advanced propulsion systems now planned for use late in this century and beyond. In order to explore such flows, a study was begun to create appropriate physical models for describing supersonic combustion, and to develop accurate and efficient numerical techniques for solving the governing equations that result from these models. From this work, two computer programs were written to study reacting flows. Both programs were constructed to consider the multicomponent diffusion and convection of important chemical species, the finite rate reaction of these species, and the resulting interaction of the fluid mechanics and the chemistry. The first program employed a finite difference scheme for integrating the governing equations, whereas the second used a hybrid Chebyshev pseudospectral technique for improved accuracy.

Drummond, J. Philip↗