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

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 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 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.↗

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↗

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.↗

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.↗

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.↗

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.↗

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.↗

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.↗

Double-diffusive fingering convection in a porous medium

The characteristics of nonlinear two-dimensional horizontally periodic double-diffusive fingering convection in a saturated porous medium is investigated, using the Darcy equation including Brinkman and Forchheimer terms for the momentum equation. To solve the equations and the corresponding initial and boundary conditions, a Galerkin method is applied in the horizontal direction, and a hybrid finite difference method is used in the vertical direction. The developed computer code was used to compute the thermal convection case, and the results were found to be in good agreement with existing results.

Chen, Falin↗

Stellar convection theory. III - Dynamical coupling of the two convection zones in A-type stars by penetrative motions

The thermal convection occurring over many density scale heights in an A-type star outer envelope, encompassing both the hydrogen and helium convectively unstable zones, is examined by means of anelastic modal equations. The single-mode anelastic equations for such compressible convection display strong overshooting of the motions into adjacent radiative zones, which would preclude diffusive separation of elements in the supposedly quiescent region between the two unstable zones. In addition, the anelastic solutions reveal that the two zones of convective instability are dynamically coupled by the overshooting motions. The two solutions that the nonlinear single-mode equations admit for the same horizontal wavelength are distinguished by the sense of the vertical velocity at the center of the three-dimensional cell. It is suggested that strong horizontal shear flows should be present just below the surface of the star, and that the large-scale motions extending into the stable atmosphere would appear mainly as horizontal flows.

Latour, J.↗

Convective instability in superposed fluid and anisotropic porous layers

The onset of thermal convection due to heating from below in a system consisting of a fluid layer overlying a porous layer with anisotropic permeability and thermal diffusivity is studied. The linear perturbation equations are solved numerically. It is found that the effects of anisotropy on the onset of thermal convection are most profound for small values of the depth ratio xi (ratio of fluid layer thickness to porous layer thickness), since in that case, the onset of convection corresponds to significant motion in both layers. For fixed vlues of the vertical permeability in the porous medium, decreasing the value of zeta (ratio of horizontal to vertical permeability) leads to stabilization of the superposed layer configuration because of increased resistance to motion in the porous medium. For larger values of xi, the onset of motion is increasingly confined to the fluid layer, with the transport of heat through the porous layer occurring primarily by conduction. Accordingly, the influence of zeta on the stability characteristics for larger xi is less significant than the effects of an anisotropic thermal conductivity.

Chen, Falin↗

Numerical simulation of a supersonic reacting mixing layer

In order to arrive at physical models that can adequately describe supersonic combustion, and develop accurate and efficient numerical techniques for the solution of such models' governing equations, a computer program has been developed for the study of reacting flows which considers the multicomponent diffusion and convection of important chemical species, as well as their finite state reaction and the interaction of the fluid mechanics and the chemistry that occurred. The code employs a hybrid Chebyshev pseudospectral technique for integration of the models' resulting governing equations; the program is here used to study a spatially developing and reacting mixing layer.

Drummond, J. Philip↗

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.↗

Modeling and calculation of turbulent transport in free-shear flows

In this work the applicability of the combined bulk convection and gradient transport hypotheses for modeling turbulent diffusion is investigated. The resulting model equation, namely the one-equation model, is solved for free-shear flows by an implicit finite-difference method. Results indicate that significant improvements over previous models can be achieved with this new formulation of turbulent diffusion for both heat and momentum transport.

Biringen, S.↗

Cosmic-ray streaming in the Born approximation

The present work invokes the Born approximation to derive a more accurate expression for the streaming of cosmic rays parallel to the mean magnetic field. While all prior results pertaining to the helicity dependence of the diffusion coefficient and convection speed can be recovered as special cases from this streaming equation, it is concluded that a new set of transport parameters presented here is more appropriate for the solar modulation of galactic cosmic rays. In addition, a new parameter related to time variability, which may be a dominant cause of charge sign-dependent transport of solar particles, is introduced.

Bieber, J. W.↗

Forced Convection and Sedimentation Past a Flat Plate

The steady laminar flow of a well-mixed suspension of monodisperse solid spheres, convected steadily past a horizontal flat plate and sedimenting under the action of gravity, is examined. It is shown that, in the limit as Re approaches infinity and epsilon approaches 0, where Re is the bulk Reynolds number and epsilon is the ratio of the particle radius a to the characteristic length scale L, the analysis for determining the particle concentration profile has several aspects in common with that of obtaining the temperature profile in forced-convection heat transfer from a wall to a fluid stream moving at high Reynolds and Prandtl numbers. Specifically, it is found that the particle concentration remains uniform throughout the O(Re(exp -1/2)) thick Blasius boundary layer except for two O(epsilon(exp 2/3)) thin regions on either side of the plate, where the concentration profile becomes non-uniform owing to the presence of shear-induced particle diffusion which balances the particle flux due to convection and sedimentation. The system of equations within this concentration boundary layer admits a similarity solution near the leading edge of the plate, according to which the particle concentration along the top surface of the plate increases from its value in the free stream by an amount proportional to X(exp 5/6), with X measuring the distance along the plate, and decreases in a similar fashion along the underside. But, unlike the case of gravity settling on an inclined plate in the absence of a bulk flow at infinity considered earlier, here the concentration profile remains continuous everywhere. For values of X beyond the region near the leading edge, the particle concentration profile is obtained through the numerical solution of the relevant equations. It is found that, as predicted from the similarity solution, there exists a value of X at which the particle concentration along the top side of the plate attains its maximum value phi(sub m) and that, beyond this point, a stagnant sediment layer will form that grows steadily in time. This critical value of X is computed as a function of phi(sub s), the particle volume fraction in the free stream. In contrast, but again in conformity with the similarity solution, for values of X sufficiently far removed from the leading edge along the underside of the plate, a particle-free region is predicted to form adjacent to the plate. This model, with minor modifications, can be used to describe particle migration in other shear flows, as, for example, in the case of crossflow microfiltration.

Pelekasis, Nikolaos A.↗