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 451 records · Page 25

Electron pitch-angle diffusion driven by oblique whistler-mode turbulence.

A general description of cyclotron harmonic resonant pitch-angle scattering is presented. Quasi-linear diffusion coefficients are prescribed in terms of the wave normal distribution of plasma wave energy. Numerical computations are performed for the specific case of relativistic electrons interacting with a band of low frequency whistler-mode turbulence. A parametric treatment of the wave energy distribution permits normalized diffusion coefficients to be presented graphically solely as a function of the electron pitch-angle. The diffusion coefficients generally decrease with increasing cyclotron harmonic number. Higher harmonic diffusion is insignificant at very small electron pitch-angles, but becomes increasingly important as the pitch-angle increases. One thus expected the rate of pitch-angle scattering to decrease with increasing electron energy, since the resonant value of the latter varies proportionately with harmonic number. This indicates that, in mirror-type magnet field geometrics, such as the earth's radiation belts, the diffusion losses of high energy electrons are likely to be appreciably slower than those at low energy.

Lyons, L. R.↗

Simulation of the early startup period of high-temperature heat pipes from the frozen state by a rarefied vapor self-diffusion model

The heat pipe startup process is described physically and is divided into five periods for convenience of analysis. The literature survey revealed that none of the previous attempts to simulate the heat pipe startup process numerically were successful, since the rarefied vapor flow in the heat pipe was not considered. Therefore, a rarefied vapor self-diffusion model is proposed, and the early startup periods, in which the rarefied vapor flow is dominant within the heat pipe, are first simulated numerically. The numerical results show that large vapor density gradients existed along the heat pipe length, and the vapor flow reaches supersonic velocities when the density is extremely low. The numerical results are compared with the experimental data of the early startup period with good agreement.

Cao, Y.↗

Reduction of solar cell efficiency by edge defects across the back-surface-field junction - A developed perimeter model

Material imperfections, impurity clusters and fabrication defects across the back-surface-field junction can degrade the performance of high-efficiency solar cells. The degradation from defects appearing on the circumference of a solar cell is analyzed using a two-region developed perimeter device model. The width of the defective perimeter region is characterized by the range or the distance-of-influence of the defective edge and this width is about two diffusion lengths. The defective edge is characterized by a surface recombination velocity. Family of theoretical curves and numerical examples are presented to show that significant reduction of open-circuit voltage can occur in high-efficiency cells which are thin compared with the diffusion length. In one example, the degradation is decreased from 135 mV to 75 mV when the cell size is increased from 10 to 100 times the diffusion length in a thin cell whose thickness is 1% of the diffusion length.

Sah, C. T.↗

Suppression Characteristics of Cup-Burner Flames in Low Gravity

The structure and suppression of laminar methane-air co-flow diffusion flames formed on a cup burner have been studied experimentally and numerically using physically acting fire-extinguishing agents (CO2, N2, He, and Ar) in normal earth (lg) and zero gravity (0g). The computation uses a direct numerical simulation with detailed chemistry and radiative heat-loss models. An initial observation of the flame without agent was also made at the NASA Glenn 2.2-Second Drop Tower. An agent was introduced into a low-speed coflowing oxidizing stream by gradually replacing the air until extinguishment occurred under a fixed minimal fuel velocity. The suppression of cup-burner flames, which resemble real fires, occurred via a blowoff process (in which the flame base drifted downstream) rather than the global extinction phenomenon typical of counterflow diffusion flames. The computation revealed that the peak reactivity spot (the reaction kernel) formed in the flame base was responsible for attachment and blowoff phenomena of the trailing diffusion flame. The thermal and transport properties of the agents affected the flame extinguishment limits.

Takahashi, Fumiaki↗

Two-Flux Green's Function Analysis for Transient Spectral Radiation in a Composite

An analysis is developed for obtaining transient temperatures in a two-layer semitransparent composite with spectrally dependent properties. Each external boundary of the composite is subjected to radiation and convection. The two-flux radiative transfer equations are solved by deriving a Green's function. This yields the local radiative heat source needed to numerically solve the transient energy equation. An advantage of the two-flux method is that isotropic scattering is included without added complexity. The layer refractive indices are larger than one. This produces internal reflections at the boundaries and the internal interface; the reflections are assumed diffuse. Spectral results using the Green's function method are verified by comparing with numerical solutions using the exact radiative transfer equations. Transient temperature distributions are given to illustrate the effect of radiative heating on one side of a composite with external convective cooling. The protection of a material from incident radiation is illustrated by adding scattering to the layer adjacent to the radiative source.

Siegel, Robert↗

Oxidation of Silicon Particles Suspended in Mullite Matrix With Dual-Mode Vacancy Diffusion of Oxygen

Mullite bond coats filled with silicon particles have recently been proposed for the next generation of environmental barrier coatings (EBCs). The sacrificial oxidation of the silicon particles eliminates some of the environmental durability issues associated with traditional EBC designs. To investigate this novel design approach, a mathematical formulation and resulting numerical method have been developed for simulating the oxidation of silicon particles embedded in a mullite matrix. The focus is on the oxidation of the silicon in a dry oxygen environment where the oxidant is presumably oxygen. Oxygen transport from the exterior environment through the mullite is assumed to occur via two lattice defect diffusion mechanisms. As a result, the rate of oxygen transport is dictated by the self-diffusivity of the lattice defects. The numerical method is implemented using the COMSOL Multiphysics® program (COMSOL, Inc.). The method is applied to simulate the oxidation of silicon particles in a mullite pellet in recent oxidation experiments. Numerical results compare favorably with the oxidation behavior observed in the oxidation experiments. The numerical method is also used to investigate the effects of various features of the mullite-silicon particle system on the oxidation behavior.

silicon oxidation↗

Numerical simulation of axisymmetric turbulent flow in combustors and diffusors

A method for predicting turbulent flow in combustors and diffusers is developed. The Navier-Stokes equations, incorporating a turbulence kappa-epsilon model equation, were solved in a nonorthogonal curvilinear coordinate system. The solution applied the finite volume method to discretize the differential equations and utilized the SIMPLE algorithm iteratively to solve the differenced equations. A zonal grid method, wherein the flow field was divided into several subsections, was developed. This approach permitted different computational schemes to be used in the various zones. In addition, grid generation was made a more simple task. However, treatment of the zonal boundaries required special handling. Boundary overlap and interpolating techniques were used and an adjustment of the flow variables was required to assure conservation of mass, momentum and energy fluxes. The numerical accuracy was assessed using different finite differencing methods, i.e., hybrid, quadratic upwind and skew upwind, to represent the convection terms. Flows in different geometries of combustors and diffusers were simulated and results compared with experimental data and good agreement was obtained.

Yung, Chain Nan↗

Passive Rocket Diffuser Testing: Reacting Flow Performance of Advanced Configurations

Upper stage rocket engines are optimized for operation at the edge of space. Ground testing often requires simulation of high-altitude environments to prevent damage from off-design operation. Passive supersonic diffusers are an effective means of establishing the partial vacuum needed around the nozzle exit. The low cost of diffusers relative to active vacuum sources (ejectors, pumps, etc.) provides a perpetual incentive to improve their performance and expand the envelope of passively testable engines. Second-throat diffusers have been the de facto standard for rocket testing since the 1950s, and numerous studies have focused on parametric optimization of the topology. Yet, even the best second-throat designs may require more driving pressure than an upper stage engine can provide. To overcome conventional limits, engineers at NASA’s Stennis Space Center (SSC) adapted an obscure cold-flow diffuser topology to the extremes of immersion in supersonic combustion products. Aerodynamic experimentation was conducted at SSC’s E-3 test stand to quantify the performance of prototypical hot-fire spike diffuser hardware and enable direct comparison to prior states of the art. A LO X /GH 2 thruster was contoured to approximate the regeneratively-cooled nozzle of an upper-stage engine and fired into 28 downstream diffuser configurations: 16 second-throat, 2 centerbody, and 10 spike. Temperature and pressure were recorded along the test article walls for a range of chamber pressures spanning 3.5 - 5.0 MPa. The second-throat configurations with the lowest start and unstart pressure ratios were shown to outperform equivalent systems in literature and taken as baselines for evaluation of advanced diffuser performance. Compared to these exemplars, centerbody diffusers decreased start pressure ratio by 11% but increased unstart pressure ratio by 15%. Spike diffusers reduced start and unstart pressure ratios by 24% and 9%, respectively.

Rocket↗

Numerical solution of compressible viscous flows

The presented review is concerned with the problem of calculating compressible viscous flows. Basic numerical considerations and problems associated with calculating viscous flows are examined and current numerical approaches toward the solution of the Navier-Stokes equations are discussed. It is pointed out that the numerical solution of the full time-dependent equations for turbulent flow is not practical with present computers. Therefore, turbulence effects must be accounted for by modeling. Developments related to turbulence modeling are described. In connection with a discussion of numerical methods for solving viscous flow equations, attention is given to numerical domains of dependence of typical explicit and implicit methods, the diffusion problem, the convection-diffusion problem, and the split-hybrid method.

Maccormack, R. W.↗

The correlation of the cosmic X-ray background with the light of galaxies

Models for the origin of the cosmic X-ray background (CXB) radiation attribute it to numerous weak, unresolved discrete sources (e.g., active galactic nuclei), to some diffuse emission mechanism (e.g., thermal bremmstrahlung from a hot intergalactic plasma), or to an unresolved population of discrete sources which evolve substantially with redshift (e.g., quasars). It is shown that the comparison of the CXB flux variations with flux variations in the light from galaxies can measure the relative contribution of these possible sources of the CXB. Applying the technique to available Uhuru data, it is found that the absence of correlation between the optical flux variations and the X-ray flux variations sets an upper limit of approximately 50% on the fraction of the CXB originating with any classes of X-ray sources substantially represented among bright (apparent photographic magnitude less than or equal to 15.5) galaxies. The data are consistent with nearly all of the CXB being due to diffuse emission or to a class of sources whose density increases rapidly with redshift.

Turner, E. L.↗

The effects of the Asselin time filter on numerical solutions to the linearized shallow-water wave equations

In the present investigation, a one-dimensional linearized analysis is used to determine the effect of Asselin's (1972) time filter on both the computational stability and phase error of numerical solutions for the shallow water wave equations, in cases with diffusion but without rotation. An attempt has been made to establish the approximate optimal values of the filtering parameter nu for each of the 'lagged', Dufort-Frankel, and Crank-Nicholson diffusion schemes, suppressing the computational wave mode without materially altering the physical wave mode. It is determined that in the presence of diffusion, the optimum filter length depends on whether waves are undergoing significant propagation. When moderate propagation is present, with or without diffusion, the Asselin filter has little effect on the spatial phase lag of the physical mode for the leapfrog advection scheme of the three diffusion schemes considered.

Schlesinger, R. E.↗

Effects of a wavy neutral sheet on cosmic ray anisotropies

The first results of a three-dimensional numerical code calculating cosmic ray anisotropies is presented. The code includes diffusion, convection, adiabatic cooling, and drift in an interplanetary magnetic field model containing a wavy neutral sheet. The 3-D model can reproduce all the principal observations for a reasonable set of parameters.

Kota, J.↗

Convergence of infinite dimensional sampled LQR problems - Theory and numerical results

A theory is developed for the convergence of the closed-loop solution to infinite-dimensional discrete-time linear-quadratic regulator (LQR) problems on the infinite time interval to the solution of a corresponding continuous-time LQR problem as the length of the sampling interval tends toward zero. Convergence of solutions to the operator algebraic Riccati equation and corresponding optimal feedback control gains is guaranteed under appropriate uniform stabilizability and detectability conditions and consistent sampling. Also presented are numerical results involving the optimal LQ control of a heat or diffusion equation, a hereditary or delay differential equation, and a hybrid system of ordinary and partial differential equations describing the transverse vibration of a cantilevered Voigt-Kelvin viscoelastic beam with tip mass.

Rosen, I. G.↗

Renormalization group analysis of anisotropic diffusion in turbulent shear flows

The renormalization group is applied to compute anisotropic corrections to the scalar eddy diffusivity representation of turbulent diffusion of a passive scalar. The corrections are linear in the mean velocity gradients. All model constants are computed theoretically. A form of the theory valid at arbitrary Reynolds number is derived. The theory applies only when convection of the velocity-scalar correlation can be neglected. A ratio of diffusivity components, found experimentally to have a nearly constant value in a variety of shear flows, is computed theoretically for flows in a certain state of equilibrium. The theoretical value is well within the fairly narrow range of experimentally observed values. Theoretical predictions of this diffusivity ratio are also compared with data from experiments and direct numerical simulations of homogeneous shear flows with constant velocity and scalar gradients.

Rubinstein, Robert↗

Composition Effects on Phase Formation and Stability

In this report, results from experimental studies of the composition dependence of nucleation are presented. A model for nucleation that takes simultaneous account of the interfacial attachment processes at the growing cluster interface and diffusion into the region surrounding the developing cluster is presented and numerical results are discussed.

Kelton, Kenneth F.↗

Large Eddy Simulation Of Gravitational Effects In Transitional And Turbulent Gas-Jet Diffusion Flames

The influence of gravity on the spatial and the compositional structures of transitional and turbulent hydrocarbon diffusion flames are studies via large eddy simulation (LES) and direct numerical simulation (DNS) of round and planar jets. The subgrid-scale (SGS) closures in LES are based on the filtered mass density function (FMDF) methodology. The FMDF represents the joint probability density function (PDF) of the SGS scalars, and is obtained by solving its transport equation. The fundamental advantage of LES/FMDF is that it accounts for the effects of chemical reaction and buoyancy exactly. The methodology is employed for capturing some of the fundamental influences of gravity in equilibrium flames via realistic chemical kinetic schemes. Some preliminary investigation of the gravity effects in non-equilibrium flames is also conducted, but with idealized chemical kinetics models.

Jaberi, Farhad A.↗

AGN Obscuration Through Dusty Infrared Dominated Flows: Radiation-Hydrodynamics Solution for the Wind - 1

We construct a radiation-hydrodynamics model for the obscuring toroidal structure in active galactic nuclei. In this model the obscuration is produced at parsec scale by a dense, dusty wind which is supported by infrared radiation pressure on dust grains. To find the distribution of radiation pressure, we numerically solve the 2D radiation transfer problem in a flux limited diffusion approximation. We iteratively couple the solution with calculations of stationary 1D models for the wind, and obtain the z-component of the velocity. Our results demonstrate that for AGN luminosities greater than 0.1 L(sub edd) external illumination can support a geometrically thick obscuration via outflows driven by infrared radiation pressure. The terminal velocity of marginally Compton-thin models (0.2 < tau(sub T) < 0.6), is comparable to or greater than the escape velocity. In Compton thick models the maximum value of the vertical component of the velocity is lower than the escape velocity, suggesting that a significant part of our torus is in the form of failed wind. The results demonstrate that obscuration via normal or failed infrared-driven winds is a viable option for the AGN torus problem and AGN unification models. Such winds can also provide an important channel for AGN feedback.

Dorodnitsyn, A.↗