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 343 records · Page 19

Comparison of big event with calculations of the air shower development

The incidence of high energy hadrons and electron-photons in air showers at various stages of development is calculated. Numerical calculation is used to solve the diffusion equation for a nuclear cascade and analytical calculation for cascade shower induced gamma rays. From these calculations, one can get the longitudinal development of the high energy hadron and electron-photon components, and the energy spectra of these components at various depths of air shower development. The total number of hadrons (N sub H) and electron-photon components (N sub gamma) are related according to stages of the air shower development and primary energy. The relation of the total energy of hadron and electron-photon component above the threshold energy is given. The energy balance between both components is also a useful parameter to study high energy events accompanying air showers. The relation of N sub H and fractional hadronic energy E (sum E sub H sup gamma/sum E sub H sup gamma + Sum E sub gamma) is calculated. This relation is helpful to understand the stage of air shower development(t) and primary energy (E sub p).

Niwa, M.↗

Monte Carlo study of a model of diffusion-controlled reactions

The Monte Carlo method is used to perform averages over sink configurations in the present study of diffusion-controlled reactions occurring between solute particles and immobile spherical sinks. In order to determine the average steady state solute concentration profile in a locally perturbed solution for sink volume fractions phi of less than 0.3, the diffusion equation in the monopolar plus dipolar approximation of diffusive couplings between the sinks is solved numerically. The Monte Carlo method is shown to be the most accurate and efficient in the phi = 0.001-0.1 region, where a system of only 25 sinks suffices and the monopolar approximation alone is sufficiently accurate.

Beenakker, C. W. J.↗

A multidimensional model for aerosols - Description of computational analogs

The numerical algorithms used to simulate the advection, diffusion, sedimentation, coagulation, and condensational growth of atmospheric aerosols are described. The model can be used in one, two, or three spatial dimensions. The continuity equation in a generalized horizontal and vertical coordinate system is developed, which allows the model to be quickly adapted to a wide variety of dynamical models of global or regional scale. Algorithms are developed to treat the various physical processes, and the results of simulations are presented, which show the strengths and weaknesses of these algorithms. Although the emphasis is on the modeling of aerosols, the work is also applicable to the simulations of the transport of gases.

Toon, O. B.↗

Field-aligned dynamics of chemically induced perturbations to the ionosphere

Ionospheric modification using photochemically reactive vapors is studied with a one-dimensional, multiion, fluid model of plasma flow along magnetic field lines. The magnitudes of ion and electron density changes are determined by considering both chemical processes (i.e., photoionization, ion-molecule reactions, dissociative recombination, electron attachment) and transport processes (i.e., multispecies diffusion, electric currents, ambipolar electric fields). The numerical treatment in the model is not specific to any type of release or any interaction chemistry. It has been used to simulate releases of Ba, CO2, and CF3Br in the ionosphere, but generalization to other species may be easily accomplished. The results of the calculations are found to be in good agreement with experimental observations. The feasibility of modifying the parallel current paths in the auroral F region is examined.

Bernhardt, P. A.↗

Numerical simulation of low Prandtl number turbulent mixing

Numerical simulations of turbulent mixing of strongly diffusive scalar fields were carried out with and without subgrid-scale modeling of the small-scale strain field. For low Reynolds number flows, when the rate of strain field (determined primarily by the small scales) is fully resolved, the scalar microstructure was found to collapse under Batchelor rate-of-strain scaling even for small Prandtl numbers, in agreement with Kerr. For high Reynolds number flows, when small-scale straining is modeled with a subgrid-scale model, the scalar microstructure follows the Batchelor, Howells, and Townsend predictions that the small-scale rate-of-strain is irrelevant.

Gibson, C.↗

High Resolution Radio Imaging of the Merging Galaxies NGC3256 and NGC4194

We present new 6cm and 4cm radio continuum images of the central regions of the merging galaxy systems NGC3256 and NGC4194. NGC3256 is imaged with a resolution of approx. 1 in. or approx. 190pc; NGC4194 is imaged with a resolution of approx. 0.3 in. or approx. 50pc. In both systems, we detect numerous compact radio sources embedded in more diffuse radio emission. We detect 65 compact sources in NGC3256 at 6cm and we detect 46 compact sources in NGC4194, both to a limiting luminosity of approx. 5 x 10(exp 18) W/ Hz or approx. 5 times the luminosity of Cas A. Most of the compact radio sources are loosely associated with active star forming regions but not with specific optical emission sources. Several compact radio sources in NGC3256 are near positions of compact X-ray sources detected by Lira et al.. In both NGC3256 and NGC4194, we are able to measure reliable spectral indices for the stronger sources. We find in NGC3256 approx. 20% have nominally flat radio spectral indices (indicating they are dominated by thermal radio emission from HII regions) while approx. 80% have nominally steep spectral indices (indicating they are dominated by nonthermal emission from supernova remnants). In NGC4194, half the compact radio sources have flat spectral indices and half have steep indices. For the flat-spectrum sources, we estimate the number of young massive stars and the associated ionized gas masses. For the steep-spectrum sources, we estimate supernova rates. We compare these results with those from other well-studied merging galaxy systems. We gratefully acknowledge use of the NRAO Very Large Array (VLA) and the VLA Archive. NRAO is a facility of the National Science Foundation, operated under cooperative agreement by Associated Universities, Inc.

Neff, S. G.↗

Obscured Star-Formation in Merging Galaxies: High Resolution Radio Imaging of a Time-Ordered Sequence

We present new, deep, high resolution 6cm and 4cm radio continuum images of the central regions of a time-ordered sequence of seven large galaxy mergers. The radio observations are able to detect star-forming re- gions that are completely obscured at optical wavelengths. In all systems, we detect numerous compact radio sources embedded in more diffuse ra- dio emission, with limiting luminosities of approx. 1-5 x 10(exp l8) W Hz or approx. 1-5 times the luminosity of Cas A. Many of the compact radio sources are loosely associated with active starforming regions but not with specific optical or W emission sources. Several of the compact radio sources are coincident with Ultra-luminous X-ray objects (ULX's). In most systems, we are able to measure reliable spectral indices for the stronger sources. We find that the fraction of compact radio cources with nominally flat radio spectral indices (indicating they ae dominated by thermal radio emission from HII regions) decreases with merger age, while the fraction of sources with nonimally steep spectral indices (indicating they are dominated by nonthermal emission from supernova remnants) increases. For the flat-spectrum sources, we estimate the numbers of young massive stars, associated ionized gas masses, we estimate supernova rates and required star-formation rates, We compare these results with those from other well-studied merging galaxy systems and from other determinations of star-formation rates. We gratefully acknowledge use of the NRAO Very Large Array (VLA) and the VLA Archive. NRAO is a facility of the National Science Foundation, operated under cooperative agreement by Associated Universities, Inc.

Neff, S. G.↗

Galactic Cosmic-Ray Energy Spectra and Composition during the 2009-2010 Solar Minimum Period

We report new measurements of the elemental energy spectra and composition of galactic cosmic rays during the 2009-2010 solar minimum period using observations from the Cosmic Ray Isotope Spectrometer (CRIS) onboard the Advanced Composition Explorer. This period of time exhibited record-setting cosmic-ray intensities and very low levels of solar activity. Results are given for particles with nuclear charge 5 <= Z <= 28 in the energy range approx. 50-550 MeV / nucleon. Several recent improvements have been made to the earlier CRIS data analysis, and therefore updates of our previous observations for the 1997-1998 solar minimum and 2001-2003 solar maximum are also given here. For most species, the reported intensities changed by less than approx. 7%, and the relative abundances changed by less than approx. 4%. Compared with the 1997-1998 solar minimum relative abundances, the 2009-2010 abundances differ by less than 2sigma, with a trend of fewer secondary species observed in the more recent time period. The new 2009-2010 data are also compared with results of a simple "leaky-box" galactic transport model combined with a spherically symmetric solar modulation model. We demonstrate that this model is able to give reasonable fits to the energy spectra and the secondary-to-primary ratios B/C and (Sc+Ti+V)/Fe. These results are also shown to be comparable to a GALPROP numerical model that includes the effects of diffusive reacceleration in the interstellar medium.

Galactic Cosmic Rays↗

Quasi-linear pitch angle and energy diffusion of pickup ions near Comet Halley

The process of ion pickup in the environment of Halley's comet is studied in order to see if velocity diffusion driven by the observed level of turbulence can explain the observed development of the implanted ion distribution. The theoretical description used is based on a quasi-linear approach and considers the implantation and transport of cometary ions along solar wind flow lines. To make such a study requires some way of extrapolating the measurements on the Giotto trajectory into the upstream region; models for mass loading and turbulence are used. A simplified kinetic equation describing the source, convection, and quasi-linear velocity diffusion of the heavy cometary ions is solved numerically along flow lines parallel to the sun-comet line. Full two-dimensional (pitch angle and velocity) distributions are obtained at positions along the Giotto trajectory, which are compared with measurements. This study finds that quasi-linear theory, with the empirical model for the observed turbulence level, produces the right order of pitch angle diffusion.

Huddleston, D. E.↗

PIV Measurements in Weakly Buoyant Gas Jet Flames

Despite numerous experimental investigations, the characterization of microgravity laminar jet diffusion flames remains incomplete. Measurements to date have included shapes, temperatures, soot properties, radiative emissions and compositions, but full-field quantitative measurements of velocity are lacking. Since the differences between normal-gravity and microgravity diffusion flames are fundamentally influenced by changes in velocities, it is imperative that the associated velocity fields be measured in microgravity flames. Velocity measurements in nonbuoyant flames will be helpful both in validating numerical models and in interpreting past microgravity combustion experiments. Pointwise velocity techniques are inadequate for full-field velocity measurements in microgravity facilities. In contrast, Particle Image Velocimetry (PIV) can capture the entire flow field in less than 1% of the time required with Laser Doppler Velocimetry (LDV). Although PIV is a mature diagnostic for normal-gravity flames , restrictions on size, power and data storage complicate these measurements in microgravity. Results from the application of PIV to gas jet flames in normal gravity are presented here. Ethane flames burning at 13, 25 and 50 kPa are considered. These results are presented in more detail in Wernet et al. (2000). The PIV system developed for these measurements recently has been adapted for on-rig use in the NASA Glenn 2.2-second drop tower.

Sunderland, Peter B.↗

Effects of Lewis Number on Temperatures of Spherical Diffusion Flames

Spherical diffusion flames supported on a porous sphere were studied numerically and experimentally. Experiments were performed in 2.2 s and 5.2 s microgravity facilities. Numerical results were obtained from a Chemkin-based program. The program simulates flow from a porous sphere into a quiescent environment, yields both steady-state and transient results, and accounts for optically thick gas-phase radiation. The low flow velocities and long residence times in these diffusion flames lead to enhanced radiative and diffusive effects. Despite similar adiabatic flame temperatures, the measured and predicted temperatures varied by as much as 700 K. The temperature reduction correlates with flame size but characteristic flow times and, importantly, Lewis number also influence temperature. The numerical results show that the ambient gas Lewis number would have a strong effect on flame temperature if the flames were steady and nonradiating. For example, a 10% decrease in Lewis number would increase the steady-state flame temperature by 200 K. However, for these transient, radiating flames the effect of Lewis number is small. Transient predictions of flame sizes are larger than those observed in microgravity experiments. Close agreement could not be obtained without either increasing the model s thermal and mass diffusion properties by 30% or reducing mass flow rate by 25%.

Santa, K. J.↗

Application of the Space-Time Conservation Element and Solution Element Method to One-Dimensional Advection-Diffusion Problems

Test problems are used to examine the performance of several one-dimensional numerical schemes based on the space-time conservation and solution element (CE/SE) method. Investigated in this paper are the CE/SE schemes constructed previously for solving the linear unsteady advection-diffusion equation and the schemes derived here for solving the nonlinear viscous and inviscid Burgers equations. In comparison with the numerical solutions obtained using several traditional finite-difference schemes with similar accuracy, the CE/SE solutions display much lower numerical dissipation and dispersion errors.

Wang, Xiao-Yen↗

Thermal Transport Model for Heat Sink Design

A document discusses the development of a finite element model for describing thermal transport through microcalorimeter arrays in order to assist in heat-sinking design. A fabricated multi-absorber transition edge sensor (PoST) was designed in order to reduce device wiring density by a factor of four. The finite element model consists of breaking the microcalorimeter array into separate elements, including the transition edge sensor (TES) and the silicon substrate on which the sensor is deposited. Each element is then broken up into subelements, whose surface area subtends 10 10 microns. The heat capacity per unit temperature, thermal conductance, and thermal diffusivity of each subelement are the model inputs, as are the temperatures of each subelement. Numerical integration using the Finite in Time Centered in Space algorithm of the thermal diffusion equation is then performed in order to obtain a temporal evolution of the subelement temperature. Thermal transport across interfaces is modeled using a thermal boundary resistance obtained using the acoustic mismatch model. The document concludes with a discussion of the PoST fabrication. PoSTs are novel because they enable incident x-ray position sensitivity with good energy resolution and low wiring density.

Chervenak, James A.↗

Eddy diffusion coefficients and the variance of the atmosphere 30-60 km

The results of numerical models or of new observational programs are checked by comparing them with past observations. In view of the differing analysis techniques or differing data samples, the eddy diffusivities presented here agree remarkably well with past estimates. However, in the application of K-values to two-dimensional models, the actual magnitude of the diffusivities is no more important than their spatial patterns, i.e., their gradients with height and latitude. It should thus be noted that the present patterns are often much different from those of past results.

Nastrom, G. D.↗

On the origin of steep edges and filaments in vorticity and potential vorticity fields

High-resolution numerical calculations are shown which capture the fundamental process responsible for the intensification of vorticity gradients in an isolated vortex subject to externally imposed distrubances. Imposition of almost any weak strain or shear field in stripping away the relatively weak vorticity at the edge of the vortex and leaves it with gradients four to six orders of magnitude greater than in the initial state. Calculations displaying such enormous gradients have never been reported previously, because of the artificial eddy diffusivities that always limit such gradients in standard numerical models. The present calculations, which have no such limitations, have been made possible by the development of a novel and robust new numerical technique for vortex dynamics called contour surgery.

Dritschel, D. G.↗

Stagnation point heat transfer augmentation due to free stream turbulence

A model has been proposed for the momentum eddy diffusivity induced by free stream turbulence intensity and integral length scale. The eddy diffusivity model is applied to the stagnation point of a cylinder situated in a uniform crossflow in the presence of free stream turbulence. A numerical solution of the governing momentum and energy equations with the proposed eddy diffusivity model yielded results for the skin friction coefficient and the Nusselt number. The numerical predictions of the present work are compared with experimental data and the agreement between the two is seen to be very good.

Gorla, R. S. R.↗

A computer program for the simulation of heat and moisture flow in soils

A computer program that simulates the flow of heat and moisture in soils is described. The space-time dependence of temperature and moisture content is described by a set of diffusion-type partial differential equations. The simulator uses a predictor/corrector to numerically integrate them, giving wetness and temperature profiles as a function of time. The simulator was used to generate solutions to diffusion-type partial differential equations for which analytical solutions are known. These equations include both constant and variable diffusivities, and both flux and constant concentration boundary conditions. In all cases, the simulated and analytic solutions agreed to within the error bounds which were imposed on the integrator. Simulations of heat and moisture flow under actual field conditions were also performed. Ground truth data were used for the boundary conditions and soil transport properties. The qualitative agreement between simulated and measured profiles is an indication that the model equations are reasonably accurate representations of the physical processes involved.

Camillo, P.↗