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 793 records · Page 44

Modeling the Transport Processes of a Pair of Solar Energetic Particle Events Observed by Parker Solar Probe near Perihelion

We present model calculations of the transport processes of solar energetic particles in the corona and interplanetary medium for two events detected by Parker Solar Probe near its second perihelion on 2019 April 2 and April 4. In the 2019 April 2 event, the <100 keV proton differential intensity measured by the ISOIS EPI-Lo instrument increased by more than a factor of 10 above the pre-event intensity, whereas the ~1 MeV proton differential intensity detected by the EPI-Hi instrument did not show any intensity enhancement. In the 2019 April 4 event, the ~1 MeV proton intensity showed an increase of more than a factor of 100 above the pre-event intensity, but the <100 keV proton intensity enhancement was rather low and gradual. The observational fact that only the 2019 April 4 event displayed a high-energy intensity enhancement indicates that the associated acceleration process in the second event was more energetic than the first event. However, the gradual and low enhancement of the low-energy proton intensity in the second event is at odds with the scenario. In this paper, we investigate the injection and transport processes of protons in the corona and interplanetary magnetic fields with numerical model simulations. Our model calculations suggest that the gradual and low rise of the low-energy protons in the April 4 event was probably due to the different diffusion conditions in these two events and the fact that the spacecraft was moving away from the source region during these two events.

Solar energetic particles↗

Low-Frequency Oscillations and Transport Processes Induced by Multiscale Transverse Structures in the Polar Wind Outflow: A Three-Dimensional Simulation

Multiscale transverse structures in the magnetic-field-aligned flows have been frequently observed in the auroral region by FAST and Freja satellites. A number of multiscale processes, such as broadband low-frequency oscillations and various cross-field transport effects are well correlated with these structures. To study these effects, we have used our three-dimensional multifluid model with multiscale transverse inhomogeneities in the initial velocity profile. Self-consistent-frequency mode driven by local transverse gradients in the generation of the low field-aligned ion flow and associated transport processes were simulated. Effects of particle interaction with the self-consistent time-dependent three-dimensional wave potential have been modeled using a distribution of test particles. For typical polar wind conditions it has been found that even large-scale (approximately 50 - 100 km) transverse inhomogeneities in the flow can generate low-frequency oscillations that lead to significant flow modifications, cross-field particle diffusion, and other transport effects. It has also been shown that even small-amplitude (approximately 10 - 20%) short-scale (approximately 10 km) modulations of the original large-scale flow profile significantly increases low-frequency mode generation and associated cross-field transport, not only at the local spatial scales imposed by the modulations but also on global scales. Note that this wave-induced cross-field transport is not included in any of the global numerical models of the ionosphere, ionosphere-thermosphere, or ionosphere-polar wind. The simulation results indicate that the wave-induced cross-field transport not only affects the ion outflow rates but also leads to a significant broadening of particle phase-space distribution and transverse particle diffusion.

Ganguli, Supriya B.↗

A numerical study of the pitch-angle scattering of cosmic rays

The results are presented of a careful study of finite-difference solutions to the problem of cosmic-ray transport, including pitch-angle scattering. In contrast to some recent studies, the diffusion approximation is confirmed in cases where the scattering mean free path is small compared with other length scales. The reasons for the discrepancy with the conclusions of Gombosi and Owens are discussed.

Kota, J.↗

Numerical method for boundary layers with blowing - The exponential box scheme

The paper describes a new numerical scheme based on exponential difference operator concepts combined with Keller's (1968) box scheme approach to produce a stable second-order accurate finite-difference scheme for convection-diffusion problems arising in boundary layer flows in the presence of massive injection through a porous surface. The technique is demonstrated by application to the self-similar boundary layer equations with massive blowing at the surface.

El-Mistikawy, T. M.↗

Semiconductor photoelectrochemistry

Semiconductor photoelectrochemical reactions are investigated. A model of the charge transport processes in the semiconductor, based on semiconductor device theory, is presented. It incorporates the nonlinear processes characterizing the diffusion and reaction of charge carriers in the semiconductor. The model is used to study conditions limiting useful energy conversion, specifically the saturation of current flow due to high light intensity. Numerical results describing charge distributions in the semiconductor and its effects on the electrolyte are obtained. Experimental results include: an estimate rate at which a semiconductor photoelectrode is capable of converting electromagnetic energy into chemical energy; the effect of cell temperature on the efficiency; a method for determining the point of zero zeta potential for macroscopic semiconductor samples; a technique using platinized titanium dioxide powders and ultraviolet radiation to produce chlorine, bromine, and iodine from solutions containing their respective ions; the photoelectrochemical properties of a class of layered compounds called transition metal thiophosphates; and a technique used to produce high conversion efficiency from laser radiation to chemical energy.

Buoncristiani, A. M.↗

Double-diffusive convection with sidewalls

Stommel et al. (1956) have first described an instability, known as thermosolutal convection, thermohaline convection, or double-diffusive convection. This instability may occur in the case of a fluid in a gravitational field with two diffusing components present. The present study is concerned with the effect of sidewalls on flow in the fingering regime in the absence of applied horizontal gradients. The work was motivated by numerical results obtained on the basis of a simulation of thermosolutal convection occurring during the unidirectional solidification of a binary alloy. In this case, the unperturbed solute field in the liquid ahead of the solidifying planar interface has an exponential vertical profile because of the rejection or preferential incorporation of solute by the solid phase.

Mcfadden, G. B.↗

Effects of impurities on silicon solar-cell performance

Model analyses indicate that sophisticated solar cell designs (back surface fields, optical reflectors, surface passivation, and double layer antireflective coatings) can produce devices with conversion efficiencies above 20%. To realize this potential, the quality of the silicon from which the cells are made must be improved; and these excellent electrical properties must be maintained during device processing. As the cell efficiency rises, the sensitivity to trace contaminants also increases. For example, the threshold Ti impurity concentraion at which cell performance degrades is more than an order of magnitude lower for an 18% cell than for a 16% cell. Similar behavior occurs for numerous other metal species which introduce deep level traps that stimulate the recombination of photogenerated carriers in silicon. Purification via crystal growth in conjunction with gettering steps to preserve the large diffusion length of the as grown material can lead to the production of devices with efficiencies above 18%, as verified experimentally.

Hopkins, R. H.↗

Scaling of the beam plasma discharge for low magnetic fields

A theoretical analysis of the scaling law and the value of the threshold current for beam plasma discharge (BPD) is presented, based on the requirement for an absolute instability near the plasma frequency. It is shown that both the scaling law as well as the numerical values of Ic are consistent with the experimental data, in the low pressure regimes and for weak magnetic field experiments if the dominant particle loss mechanism is due to Bohm diffusion. The implications of the findings to electron injection in space are discussed.

Papadopoulos, K.↗

Impurity effects in silicon for high efficiency solar cells

Model analyses indicate that sophisticated solar cell designs including, e.g., back surface fields, optical reflectors, surface passivation, and double layer antireflective coatings can produce devices with conversion efficiencies above 20 percent (AM1). To realize this potential, the quality of the silicon from which the cells are made must be improved; and these excellent electrical properties must be maintained during device processing. As the cell efficiency rises, the sensitivity to trace contaminants also increases. For example, the threshold Ti impurity concentration at which cell performance degrades is more than an order of magnitude lower for an 18-percent cell. Similar behavior occurs for numerous other metal species which introduce deep level traps that stimulate the recombination of photogenerated carriers in silicon. Purification via crystal growth in conjunction with gettering steps to preserve the large diffusion length of the as-grown material can lead to the production of devices with efficiencies aboved 18 percent, as has been verified experimentally.

Hopkins, R. H.↗

Thermochemical erosion

A quasi-steady multidimensional Stefan problem is treated by means of similarity, asymptotic methods, and numerical analysis. The continuous wearing away of a solid 'erodee' is effected by gradual consumption of a solid 'eroder'. In the process under study, the erodee melts and the (liquid) melt diffuses to the eroder, where an exothermic, indefinitely rapid, diffusion-controlled surface reaction occurs. The products of reaction convey some of the chemically-generated heat to the surface of the erodee, so that further melting may occur, and the process may proceed. The analysis seeks to estimate the quasi-steady rates of consumption of the erodee and eroder, for the (convenient) self-similar case of confocal parabolas defining the boundaries of the melt layer. The analysis is simplified by observing that the flow is of lubrication-theory type: creeping, primarily unidirectional motion owing to the thinness of the melt layer between the erodee and eroder, with diffusion of heat and mass predominantly perpendicular to the principal motion. The analysis is completed by consideration of the force balance on the eroder.

Butler, Gerald↗

Magnetic fields in the accretion disks of cataclysmic variables

The radial transport of magnetic flux in CV accretion disks is considered. Turbulent diffusion in a disk with finite outer radius Rd leads to decay of large-scale magnetic fields, as magnetic flux can leak out at the outer edge of the disk. Numerical computations of the decay rate and magnetic structure are presented. It is shown that the decay time is significantly shorter than the accretion time. Centrifugally driven winds may be possible from the outer parts of the disk, provided there exist efficient dynamo processes which regenerate the field.

Van Ballegooijen, A. A.↗

Effect of ohmic, mass-transfer, and kinetic resistances on linear-sweep voltammetry in a cylindrical-pore electrode

A means is presented to account for the effect of ohmic, mass-transfer, and kinetic resistances on linear-sweep voltammograms by modeling a pore in a porous matrix as a cylindrical-pore electrode, and solving the mass and charge conservation equations in the context of this geometry for the simply redox reaction O + ne(-) yield R where both O and R are soluble species. Both analytical and numerical techniques are used to solve the governing equations. The calculated peak currents and potentials are correlated by empirical formulas to the measurable parameters: sweep rate, concentration of the redox species, diffusion coefficient, conductivity of the electrolyte, and pore dimensions. Using the correlations, a methodology is established for determining if the redox reaction kinetics are irreversible or reversible (Nernstian). If the reaction is irreversible, it is shown how the standard rate constant and the transfer coefficient may be extracted from linear-sweep voltammetry data, or, if the reaction is reversible, how the number of electrons transferred may be deduced.

Weidner, John W.↗

Adjoint operator approach to shape design for internal incompressible flows

The problem of determining the profile of a channel or duct that provides the maximum static pressure rise is solved. Incompressible, laminar flow governed by the steady state Navier-Stokes equations is assumed. Recent advances in computational resources and algorithms have made it possible to solve the direct problem of determining such a flow through a body of known geometry. It is possible to obtain a set of adjoint equations, the solution to which permits the calculation of the direction and relative magnitude of change in the diffuser profile that leads to a higher pressure rise. The solution to the adjoint problem can be shown to represent an artificially constructed flow. This interpretation provides a means to construct numerical solutions to the adjoint equations that do not compromise the fully viscous nature of the problem. The algorithmic and computational aspects of solving the adjoint equations are addressed. The form of these set of equations is similar but not identical to the Navier-Stokes equations. In particular some issues related to boundary conditions and stability are discussed.

Cabuk, H.↗

Effects of a parallel electric field and the geomagnetic field in the topside ionosphere on auroral and photoelectron energy distributions

The consequences of electric field acceleration and an inhomogeneous magnetic field on auroral electron energy distributions in the topside ionosphere are investigated. The one-dimensional, steady state electron transport equation includes elastic and inelastic collisions, an inhomogeneous magnetic field, and a field-aligned electric field. The case of a self-consistent polarization electric field is considered first. The self-consistent field is derived by solving the continuity equation for all ions of importance, including diffusion of O(+) and H(+), and the electron and ion energy equations to derive the electron and ion temperatures. The system of coupled electron transport, continuity, and energy equations is solved numerically. Recognizing observations of parallel electric fields of larger magnitude than the baseline case of the polarization field, the effect of two model fields on the electron distribution function is investigated. In one case the field is increased from the polarization field magnitude at 300 km to a maximum at the upper boundary of 800 km, and in another case a uniform field is added to the polarization field. Substantial perturbations of the low energy portion of the electron flux are produced: an upward directed electric field accelerates the downward directed flux of low-energy secondary electrons and decelerates the upward directed component. Above about 400 km the inhomogeneous magnetic field produces anisotropies in the angular distribution of the electron flux. The effects of the perturbed energy distributions on auroral spectral emission features are noted.

Min, Q.-L.↗

Effects of a Parallel Electric Field and the Geomagnetic Field in the Topside Ionosphere on Auroral and Photoelectron Energy Distributions

The consequences of electric field acceleration and an inhomogencous magnetic field on auroral electron energy distributions in the topside ionosphere are investigated. The one- dimensional, steady state electron transport equation includes elastic and inelastic collisions, an inhomogencous magnetic field, and a field-aligned electric field. The case of a self-consistent polarization electric field is considered first. The self-consistent field is derived by solving the continuity equation for all ions of importance, including diffusion of 0(+) and H(+), and the electron and ion energy equations to derive the electron and ion temperatures. The system of coupled electron transport, continuity, and energy equations is solved numerically. Recognizing observations of parallel electric fields of larger magnitude than the baseline case of the polarization field, the effect of two model fields on the electron distribution function in investigated. In one case the field is increased from the polarization field magnitude at 300 km to a maximum at the upper boundary of 800 km, and in another case a uniform field is added to the polarization field. Substantial perturbations of the low energy portion of the electron flux are produced: an upward directed electric field accelerates the downward directed flux of low-energy secondary electrons and decelerates the upward directed component. Above about 400 km the inhomogencous magnetic field produces anisotropies in the angular distribution of the electron flux. The effects of the perturbed energy distributions on auroral spectral emission features are noted.

Min, Q.-L.↗

Numerical Simulation of a Chemically Reacting Sorbent Bed for LSS Applications

A detailed numerical model of a chemisorption bed has been developed. The model is based on the constant pressure mass transport equation for gaseous flow through a packed bed, and the equation for diffusion and reaction within a spherical particle. Because there is a wealth of data from the NASA and the Navy bodies of literature, the LiOH-H2O-CO2 system is chosen for application of the model and interpretation of results. Prior models of this system from the life support literature are limited. The current model incorporates many of the features of elaborate models developed for investigation of industrial systems or energy applications (e.g., coal, desulphurization): it distinguishes bulk convection and bed dispersion; mass transport to the particle surface, transport within the particle, and reaction. It uses the nonsteady (not pseudo-steady state) form of the equations. The chemistry is modeled as a multi-step, reversible reaction with evolving solid structure. The resulting system of equations is large. The ODEPACK family of solvers is used to integrate the system. Reaction coefficients are determined by experiment. Typical results of the model are illustrated with mission input parameters. Using the model, an explanation is offered for 1) the varied performance results found after pre-breathing (or after simulated pre-breathe conditions), 2) interrupted use and 3) low temperature use. In addition, options for a reusable canister are explored. The computational resource implications of adding energy equations are discussed briefly, as are applicability to other relevant space and undersea systems.

Luna, Bernadette↗

Galactic Starburst NGC 3603 from X-Rays to Radio

NGC 3603 is the most massive and luminous visible starburst region in the Galaxy. We present the first Chandra/ACIS-I X-ray image and spectra of this dense, exotic object, accompanied by deep cm-wavelength ATCA radio image at similar or less than 1 inch spatial resolution, and HST/ground-based optical data. At the S/N greater than 3 level, Chandra detects several hundred X-ray point sources (compared to the 3 distinct sources seen by ROSAT). At least 40 of these sources are definitely associated with optically identified cluster O and WR type members, but most are not. A diffuse X-ray component is also seen out to approximately 2 feet (4 pc) form the center, probably arising mainly from the large number of merging/colliding hot stellar winds and/or numerous faint cluster sources. The point-source X-ray fluxes generally increase with increasing bolometric brightnesses of the member O/WR stars, but with very large scatter. Some exceptionally bright stellar X-ray sources may be colliding wind binaries. The radio image shows (1) two resolved sources, one definitely non-thermal, in the cluster core near where the X-ray/optically brightest stars with the strongest stellar winds are located, (2) emission from all three known proplyd-like objects (with thermal and non-thermal components, and (3) many thermal sources in the peripheral regions of triggered star-formation. Overall, NGC 3603 appears to be a somewhat younger and hotter, scaled-down version of typical starbursts found in other galaxies.

Moffat, A. F. J.↗

Combustion Characteristics in a Non-Premixed Cool-Flame Regime of n-Heptane in Microgravity

A series of distinct phenomena have recently been observed in single-fuel-droplet combustion tests performed on the International Space Station (ISS). This study attempts to simulate the observed flame behavior numerically using a gaseous n-heptane fuel source in zero gravity and a time-dependent axisymmetric (2D) code, which includes a detailed reaction mechanism (127 species and 1130 reactions), diffusive transport, and a radiation model (for CH4, CO, CO2, H2O, and soot). The calculated combustion characteristics depend strongly on the air velocity around the fuel source. In a near-quiescent air environment (< or = 2 mm/s), with a sufficiently large fuel injection velocity (1 cm/s), a growing spherical diffusion flame extinguishes at ≈1200 K due to radiative heat losses. This is typically followed by a transition to the low-temperature (cool-flame) regime with a reaction zone (at ≈700 K) in close proximity to the fuel source. The 'cool flame' regime is formed due to the negative temperature coefficient in the low-temperature chemistry. After a relatively long period (≈18 s) of the cool flame regime, a flash re-ignition occurs, associated with flame-edge propagation and subsequent extinction of the re-ignited flame. In a low-speed (≈3 mm/s) airstream (which simulates the slight droplet movement), the diffusion flame is enhanced upstream and experiences a local extinction downstream at ≈1200 K, followed by steady flame pulsations (≈0.4 Hz). At higher air velocities (4-10 mm/s), the locally extinguished flame becomes steady state. The present axisymmetric computational approach helps in revealing the non-premixed 'cool flame' structure and 2D flame-flow interactions observed in recent microgravity droplet combustion experiments.

diffusion flame↗