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 145 records · Page 8

All sources NASA NTRS DOE OSTI

NASA NTRS1992

On an origin of numerical diffusion: Violation of invariance under space-time inversion

The invariant properties of the convection equation du/dt + adu/dx = 0 (where d is the partial differential operator) with respect to spatial reflection, time reversal, and space-time inversion are studied. Generally, a finite-difference analog of this equation may possess some or none of these properties. It is shown that, under certain conditions, the von Neumann amplification factor of an analog satisfies a special relation for each invariant property this analog possesses. Particularly, an analog is neutrally stable and thus free of numerical diffusion if it possesses the invariant property related to space-time inversion. It is also explained why generally (1) an upwind scheme possesses neither the invariant property related to spatial reflection nor that related to space-time inversion, and (2) an explicit scheme possesses neither the invariant property related to time reversal nor that related to space-time inversion. Extension to the viscous case and a remarkable connection between the current work and a new numerical framework for solving conservation laws are also discussed.

Chang, Sin-Chung↗

NASA NTRS1986

A spectral element method for the simulation of unsteady incompressible flows with heat transfer

The spectral element method is a high-order finite element technique for solution of the Navier-Stokes and energy equations. In the isoparametric spectral element discretization, the domain is broken up into general brick elements, and the dependent and independent variables represented as high-order tensor-product Lagrangian interpolants through Chebyshev collocation points. The nonlinear and convective terms in the governing equations are treated with explicit collocation, while the pressure and diffusive contributions are handled implicitly using variational projection operators. The method is applied to flow past a cylinder, flow in grooved channels, and natural convection in an enclosure.

Karniadakis, George E.↗

NASA NTRS1996

Relationship Between Solution Shear Viscosity and Density at the Saturation Point

Properties of supersaturated solutions such as the density, viscosity, and solute diffusivity are dependent on the solute concentration. The diffusion-boundary-layer equations are derived and solved for the natural convection case with the viscosity and density dependent on the solute concentration. The solution obtained demonstrates that, at the vicinity of the saturation concentration c(sub s), there is a non-trivial dependence of the solution viscosity eta on its density rho: eta(c(sub s)) = eta(rho(sub s)) varies as rho(sub s)(exp 1/2), where rho(sub s) = rho(c(sub s)). This result has been verified in experiments with aqueous solutions of inorganic and organic salts.

Izmailov, Alexander F.↗

NASA NTRS1989

Hopf bifurcation in the driven cavity

Incompressible two dimensional calculations are reported for the impulsively started lid driven cavity with aspect ratio two. The algorithm is based on the time dependent streamfunction equation, with a Crank-Nicolson differencing scheme for the diffusion terms, and with an Adams-Bashforth scheme for the convection terms. A multigrid method is used to solve the linear implicit equations at each time step. Periodic asymptotic solutions have been found for Re = 10000 and for Re = 5000. The Re = 5000 results are validated by grid refinement calculations. The solutions are shown to be precisely periodic, and care is taken to demonstrate that asymptotic states were reached. A discussion is included about the indicators that are used to show that an asymptotic state was reached, and to show that the asymptotic state is indeed periodic.

Goodrich, John W.↗

NASA NTRS1984

What Would a Dynamo Theorist like to Know About the Dynamics of the Solar Convection Zone?

Observations about the current status of solar dynamo theory are given. The induction equation for magnetic field is solved using assumed velocities and parametric representations of the inductive or diffusive effects of velocities. The equations of motion governing these flow are not solved in parallel. Results from global compressible convection models are discussed. Differential rotation and convection are also investigated.

Gilman, P. A.↗

NASA NTRS2013

Modeling of Particle Acceleration at Multiple Shocks via Diffusive Shock Acceleration: Preliminary Results

Successful forecasting of energetic particle events in space weather models require algorithms for correctly predicting the spectrum of ions accelerated from a background population of charged particles. We present preliminary results from a model that diffusively accelerates particles at multiple shocks. Our basic approach is related to box models in which a distribution of particles is diffusively accelerated inside the box while simultaneously experiencing decompression through adiabatic expansion and losses from the convection and diffusion of particles outside the box. We adiabatically decompress the accelerated particle distribution between each shock by either the method explored in Melrose and Pope (1993) and Pope and Melrose (1994) or by the approach set forth in Zank et al. (2000) where we solve the transport equation by a method analogous to operator splitting. The second method incorporates the additional loss terms of convection and diffusion and allows for the use of a variable time between shocks. We use a maximum injection energy (E(sub max)) appropriate for quasi-parallel and quasi-perpendicular shocks and provide a preliminary application of the diffusive acceleration of particles by multiple shocks with frequencies appropriate for solar maximum (i.e., a non-Markovian process).

Parker, L. Neergaard↗

NASA NTRS1985

Multigrid techniques for the solution of the passive scalar advection-diffusion equation

The solution of elliptic passive scalar advection-diffusion equations is required in the analysis of many turbulent flow and convective heat transfer problems. The accuracy of the solution may be affected by the presence of regions containing large gradients of the dependent variables. The multigrid concept of local grid refinement is a method for improving the accuracy of the calculations in these problems. In combination with the multilevel acceleration techniques, an accurate and efficient computational procedure is developed. In addition, a robust implementation of the QUICK finite-difference scheme is described. Calculations of a test problem are presented to quantitatively demonstrate the advantages of the multilevel-multigrid method.

Phillips, R. E.↗

NASA NTRS1986

Migration of the separation point on a deforming cylinder

An iterative scheme of solving the Navier-Stokes equations for unsteady two-dimensional flow about a deforming and translating cylinder is given. In the kth iteration, the convected vorticity appears as the source term in an equation of transient diffusion of vorticity. A novel integral transform is used to reduce this transient vorticity equation to a k-dimensional heat equation. The bounded solution of this equation is obtained with a general method of superposition for problems involving a moving boundary. An equation describing the migration of the separation point on a deforming cylinder in unsteady cross flows is derived from the analytically obtained velocity field. The radial expansion of the cylinder surface is shown to hasten the separation time and to increase the separation angle. The results imply that in the steady flow past a body that is moving forward and sinking at constant rates, the locus of points at which the azimuthal component of skin friction changes sign originates on the leeward ray at a point downstream of the front tip.

Lin, S. P.↗

NASA NTRS2015

Badhwar-O'Neill 2014 Galactic Cosmic Ray Flux Model Description

The Badhwar-O'Neill (BON) Galactic Cosmic Ray (GCR) model is based on GCR measurements from particle detectors. The model has mainly been used by NASA to certify microelectronic systems and the analysis of radiation health risks to astronauts in space missions. The BON14 model numerically solves the Fokker-Planck differential equation to account for particle transport in the heliosphere due to diffusion, convection, and adiabatic deceleration under the assumption of a spherically symmetric heliosphere. The model also incorporates an empirical time delay function to account for the lag of the solar activity to reach the boundary of the heliosphere. This technical paper describes the most recent improvements in parameter fits to the BON model (BON14). Using a comprehensive measurement database, it is shown that BON14 is significantly improved over the previous version, BON11.

O'Neill, P. M.↗

NASA NTRS2015

NASA Galactic Cosmic Radiation Environment Model: Badhwar-O'Neill (2014)

The Badhwar‐O'Neill (BON) Galactic Cosmic Ray (GCR) flux model is used by NASA to certify microelectronic systems and in the analysis of radiation health risks for human space flight missions. Of special interest to NASA is the kinetic energy region below 4.0 GeV/n due to the fact that exposure from GCR behind shielding (e.g., inside a space vehicle) is heavily influenced by the GCR particles from this energy domain. The BON model numerically solves the Fokker‐Planck differential equation to account for particle transport in the heliosphere due to diffusion, convection, and adiabatic deceleration under the assumption of a spherically symmetric heliosphere. The model utilizes a GCR measurements database from various particle detectors to determine the boundary conditions. By using an updated GCR database and improved model fit parameters, the new BON model (BON14) is significantly improved over the previous BON models for describing the GCR radiation environment of interest to human space flight.

O'Neill, P. M.↗

NASA NTRS2015

NASA Galactic Cosmic Radiation Environment Model: Badhwar - O'Neill (2014)

The Badhwar-O'Neill (BON) Galactic Cosmic Ray (GCR) flux model has been used by NASA to certify microelectronic systems and in the analysis of radiation health risks for human space flight missions. Of special interest to NASA is the kinetic energy region below 4.0 GeV/n due to the fact that exposure from GCR behind shielding (e.g., inside a space vehicle) is heavily influenced by the GCR particles from this energy domain. The BON model numerically solves the Fokker-Planck differential equation to account for particle transport in the heliosphere due to diffusion, convection, and adiabatic deceleration under the assumption of a spherically symmetric heliosphere. The model utilizes a comprehensive database of GCR measurements from various particle detectors to determine boundary conditions. By using an updated GCR database and improved model fit parameters, the new BON model (BON14) is significantly improved over the previous BON models for describing the GCR radiation environment of interest to human space flight.

Golge, S.↗

NASA NTRS2013

Modeling of Particle Acceleration at Multiple Shocks Via Diffusive Shock Acceleration: Preliminary Results

We present preliminary results from a model that diffusively accelerates particles at multiple shocks. Our basic approach is related to box models (Protheroe and Stanev, 1998; Moraal and Axford, 1983; Ball and Kirk, 1992; Drury et al., 1999) in which a distribution of particles is diffusively accelerated inside the box while simultaneously experiencing decompression through adiabatic expansion and losses from the convection and diffusion of particles outside the box (Melrose and Pope, 1993; Zank et al., 2000). We adiabatically decompress the accelerated particle distribution between each shock by either the method explored in Melrose and Pope (1993) and Pope and Melrose (1994) or by the approach set forth in Zank et al. (2000) where we solve the transport equation by a method analogous to operator splitting. The second method incorporates the additional loss terms of convection and diffusion and allows for the use of a variable time between shocks. We use a maximum injection energy (Emax) appropriate for quasi-parallel and quasi-perpendicular shocks (Zank et al., 2000, 2006; Dosch and Shalchi, 2010) and provide a preliminary application of the diffusive acceleration of particles by multiple shocks with frequencies appropriate for solar maximum (i.e., a non-Markovian process).

Parker, Linda Neergaard↗

NASA NTRS1974

Solar modulation

The important aspects of solar modulation of galactic cosmic rays are outlined. The conditions in the interplanetary medium that affect the modulation are summarized, and equations are developed to describe the behavior of cosmic rays in the solar system, including the effects of diffusion, convection, and energy loss due to adiabatic deceleration in the expanding solar wind. The particle diffusion coefficient is determined from the power spectra of interplanetary magnetic field irregularities. After some solutions for the modulation equations have been derived, the probable forms of the unmodulated cosmic ray spectra are discussed. Other topics include solar cycle variations, energy loss effects, anisotrophies, and gradients.

Fisk, L. A.↗

NASA NTRS1988

On the application of pseudo-spectral FFT technique to non-periodic problems

The reduction-to-periodicity method using the pseudo-spectral Fast Fourier Transform (FFT) technique is applied to the solution of nonperiodic problems including the two-dimensional Navier-Stokes equations. The accuracy of the method is demonstrated by calculating derivatives of given functions, one- and two-dimensional convective-diffusive problems, and by comparing the relative errors due to the FFT method with seocnd order Finite Difference Methods (FDM). Finally, the two-dimensional Navier-Stokes equations are solved by a fractional step procedure using both the FFT and the FDM methods for the driven cavity flow and the backward facing step problems. Comparisons of these solutions provide a realistic assessment of the FFT method indicating its range of applicability.

Biringen, S.↗

NASA NTRS1989

An asymptotic description of transient settling and ultrafiltration of colloidal dispersions

The present model for the sedimentation of colloidal systems includes a diffusion term in the governing equation which, in the regions above the sediment, acts as small perturbation to the Kynch (1952) theory. Within the sediment, diffusion is comparable to convection due to the the high solid volume fraction. An application of the method of matched asymptotic expansions to the conservation equation allows a complete description of the settling process to be formulated, with specific attention to volume-fraction evolution in the sediment.

Davis, K. E.↗

NASA NTRS1971

Transport of solar flare protons: Comparison of a new analytic model with spacecraft measurements

An analytic solution was obtained to the complete Fokker-Planck equation including the effects of convection, interplanetary deceleration and acceleration, corotation, and anisotropic diffusion. The solution yields a time to maximum for the particle flux of about 10 h and an exponential decay time of about 5 h. Several solar flare particle events were observed with the solar and galactic cosmic ray experiment on OGO-6. Comparisons of the calculated time dependence of the fluxes with these observations of 1 to 70 MeV protons show that the model adequately describes both the rise and decay times.

Lupton, J. E.↗

NASA NTRS1972

Galactic cosmic ray modulation from 1965-1970.

Numerical solutions of the cosmic-ray equation of transport within the solar cavity and including the effects of diffusion, convection, and energy losses due to adiabatic deceleration, have been used to reproduce the modulation of galactic electrons, protons, and helium nuclei observed during the period from 1965 to 1970. Kinetic energies between 10 and 10,000 MeV/nucleon are considered. Computed and observed spectra are given for the years 1965, 1968, 1969, and 1970 together with the diffusion coefficients. These diffusion coefficients are assumed to be of separable form in rigidity and radial dependence, and are consistent with the available magnetic-field power spectra. The force-field solutions are given for these diffusion coefficients and galactic spectra and are compared with the numerical solutions. It is shown that the energy losses and convection lead to near-earth nuclei spectra at kinetic energies less than or equal to 100 MeV/nucleon in which the differential intensity is proportional to the kinetic energy with little dependence on the form of the galactic spectrum. This dependence is in agreement with the observed spectra of all species of atomic nuclei and it is argued that this provides strong observational evidence for the presence of energy losses in the propagation process, and for the exclusion of low-energy galactic nuclei from near earth.

Urch, I. H.↗

NASA NTRS1981

Generalized Dufort-Frankel spectral methods

An explicit time-advancing scheme for the spectral solution of parabolic equations is presented. Several two-dimensional examples are considered, including convection-diffusion and nonlinear problems, under various boundary conditions. Numerical evidence demonstrates the efficiency and accuracy of the spectral approach.

Lustman, L.↗