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 1,801 records · Page 100

A NEW MODEL FOR MIXING BY DOUBLE-DIFFUSIVE CONVECTION (SEMI-CONVECTION). III. THERMAL AND COMPOSITIONAL TRANSPORT THROUGH NON-LAYERED ODDC

Oscillatory double-diffusive convection (ODDC; also known as semi-convection) refers to a type of double-diffusive instability that occurs in regions of planetary and stellar interiors that have a destabilizing thermal stratification and a stabilizing mean molecular weight stratification. In this series of papers, we use an extensive suite of three-dimensional (3D) numerical simulations to quantify the transport of heat and chemical species by ODDC. Rosenblum et al. first showed that ODDC can either spontaneously form layers that significantly enhance the transport of heat and chemical species compared to microscopic transport or remain in a state dominated by large-scale gravity waves, in which there is a more modest enhancement of the turbulent transport rates. Subsequent studies in this series focused on identifying under what conditions layers form and quantifying transport through layered systems. Here we proceed to characterize transport through systems that are unstable to ODDC, but do not undergo spontaneous layer formation. We measure the thermal and compositional fluxes in non-layered ODDC from both two-dimensional (2D) and 3D numerical simulations, and show that 3D simulations are well approximated by similar simulations in a 2D domain. We find that the turbulent mixing rate in this regime is weak and can, to a first-level approximation, be neglected. We conclude by summarizing the findings of papers I through III into a single prescription for transport systems unstable to ODDC.

79 ASTRONOMY AND ASTROPHYSICS↗

COMBINING DIFFUSIVE SHOCK ACCELERATION WITH ACCELERATION BY CONTRACTING AND RECONNECTING SMALL-SCALE FLUX ROPES AT HELIOSPHERIC SHOCKS

Computational and observational evidence is accruing that heliospheric shocks, as emitters of vorticity, can produce downstream magnetic flux ropes and filaments. This led Zank et al. to investigate a new paradigm whereby energetic particle acceleration near shocks is a combination of diffusive shock acceleration (DSA) with downstream acceleration by many small-scale contracting and reconnecting (merging) flux ropes. Using a model where flux-rope acceleration involves a first-order Fermi mechanism due to the mean compression of numerous contracting flux ropes, Zank et al. provide theoretical support for observations that power-law spectra of energetic particles downstream of heliospheric shocks can be harder than predicted by DSA theory and that energetic particle intensities should peak behind shocks instead of at shocks as predicted by DSA theory. In this paper, a more extended formalism of kinetic transport theory developed by le Roux et al. is used to further explore this paradigm. We describe how second-order Fermi acceleration, related to the variance in the electromagnetic fields produced by downstream small-scale flux-rope dynamics, modifies the standard DSA model. The results show that (i) this approach can qualitatively reproduce observations of particle intensities peaking behind the shock, thus providing further support for the new paradigm, and (ii) stochastic acceleration by compressible flux ropes tends to be more efficient than incompressible flux ropes behind shocks in modifying the DSA spectrum of energetic particles.

79 ASTRONOMY AND ASTROPHYSICS↗

Preliminary Evidence for a Virial Shock around the Coma Galaxy Cluster

Galaxy clusters, the largest gravitationally bound objects in the universe, are thought to grow by accreting mass from their surroundings through large-scale virial shocks. Due to electron acceleration in such a shock, it should appear as a γ-ray, hard X-ray, and radio ring, elongated toward the large-scale filaments feeding the cluster, coincident with a cutoff in the thermal Sunyaev–Zel’dovich (SZ) signal. However, no such signature was found until now, and the very existence of cluster virial shocks has remained a theory. We find preliminary evidence for a large γ-ray ring of ∼5 Mpc minor axis around the Coma cluster, elongated toward the large-scale filament connecting Coma and Abell 1367, detected at the nominal 2.7σ confidence level (5.1σ using control signal simulations). The γ-ray ring correlates both with a synchrotron signal and with the SZ cutoff, but not with Galactic tracers. The γ-ray and radio signatures agree with analytic and numerical predictions if the shock deposits ∼1% of the thermal energy in relativistic electrons over a Hubble time and ∼1% in magnetic fields. The implied inverse Compton and synchrotron cumulative emission from similar shocks can contribute significantly to the diffuse extragalactic γ-ray and low-frequency radio backgrounds. Our results, if confirmed, reveal the prolate structure of the hot gas in Coma, the feeding pattern of the cluster, and properties of the surrounding large-scale voids and filaments. The anticipated detection of such shocks around other clusters would provide a powerful new cosmological probe.

79 ASTRONOMY AND ASTROPHYSICS↗

Diurnal behaviour of water on Mars

A numerical model of the diurnal transport of water across the Martian surface is developed. The atmospheric boundary layer is modeled in terms of local radiative-convective processes, and radiative effects of ice fogs near the surface are included. Diffusion of water in the ground is treated for the cases of adsorption and condensation. The model is applied to the diurnal variation of water vapor in the atmosphere as observed by Barker (1974). The morning rise in the amount of water vapor can be explained in terms of the evaporation of ground fogs. The evening decrease is compatible with the model if adsorption dominates in the soil. The average level of vapor concentration requires that the atmosphere above the boundary layer be relatively dry. The ground fogs persist until midmorning and should be observable. Some consequences of these conclusions are discussed.

Flasar, F. M.↗

Study of the Effects of Impurities on the Properties of Silicon Materials and Performance of Silicon Solar Cell

Numerical solutions were obtained from the exact one dimensional transmission line circuit model to study the following effects on the terrestrial performance of silicon solar cells: interband Auger recombination; surface recombination at the contact interfaces; enhanced metallic impurity solubility; diffusion profiles; and defect-impurity recombination centers. Thermal recombination parameters of titanium impurity in silicon were estimated from recent experimental data. Based on those parameters, computer model calculations showed that titanium concentration must be kept below 6x10 to the 12th power Ti/cu cm in order to achieve 16% AM1 efficiency in a silicon solar cell of 250 micrometers thick and 1.5 ohm-cm resistivity.

Sah, C. T.↗

The interplanetary transport of solar cosmic rays

Numerical solutions are presented for the propagation of solar cosmic rays in interplanetary space, including the effects of pitch-angle scattering and adiabatic focusing. The intensity-time profiles can be well fitted by a simple radial spatial diffusion equation with scattering mean-free path lambda(fit). The radial mean-free path so obtained is significantly larger than the true scattering mean-free path for low-rigidity particles due to both adiabatic focusing and the inapplicability of the diffusive approximation early in the event. The well-known discrepancy between lambda(fit) and the theoretical predictions may be resolved by these calculations.

Gombosi, T. I.↗

Numerical study of solar flare particle propagation in the heliosphere

Numerical solutions are presented for the propagation of solar cosmic rays in interplanetary space, including the effects of pitch-angle scattering and adiabatic focusing. The intensity-time profiles can be well fitted by a simple radial spatial diffusion equation. For low-rigidity particles the radial mean free path so obtained is significantly larger than the mean free path calculated from the scattering coefficient due to the inapplicability of the diffusive approximation early in the event. The well-known discrepancy between the scattering mean free path and the theoretical predictions may be resolved by these calculations.

Gombosi, T. I.↗

Viscous, resistive magnetohydrodynamic stability computed by spectral techniques

Expansions in Chebyshev polynomials are used to study the linear stability of one-dimensional magnetohydrodynamic quasi-equilibria, in the presence of finite resistivity and viscosity. The method is modeled on the one used by Orszag in accurate computation of solutions of the Orr-Sommerfeld equation. Two Reynolds-like numbers involving Alfven speeds, length scales, kinematic viscosity, and magnetic diffusivity govern the stability boundaries, which are determined by the geometric mean of the two Reynolds-like numbers. Marginal stability curves, growth rates versus Reynolds-like numbers, and growth rates versus parallel wave numbers are exhibited. A numerical result that appears general is that instability has been found to be associated with inflection points in the current profile, though no general analytical proof has emerged. It is possible that nonlinear subcritical three-dimensional instabilities may exist, similar to those in Poiseuille and Couette flow.

Dahlburg, R. B.↗

Third-moment closure of turbulence for predictions of separating and reattaching shear flows: A study of Reynolds-stress closure model

A numerical study of computations in backward-facing steps with flow separation and reattachment, using the Reynolds stress closure is presented. The highlight of this study is the improvement of the Reynold-stress model (RSM) by modifying the diffusive transport of the Reynolds stresses through the formulation, solution and subsequent incorporation of the transport equations of the third moments, bar-u(i)u(j)u(k), into the turbulence model. The diffusive transport of the Reynolds stresses, represented by the gradients of the third moments, attains greater significance in recirculating flows. The third moments evaluated by the development and solution of the complete transport equations are superior to those obtained by existing algebraic correlations. A low-Reynolds number model for the transport equations of the third moments is developed and considerable improvement in the near-wall profiles of the third moments is observed. The values of the empirical constants utilized in the development of the model are recommended. The Reynolds-stress closure is consolidated by incorporating the equations of k and e, containing the modified diffusion coefficients, and the transport equations of the third moments into the Reynolds stress equations. Computational results obtained by the original k-e model, the original RSM and the consolidated and modified RSM are compared with experimental data. Overall improvement in the predictions is seen by consolidation of the RMS and a marked improvement in the profiles of bar-u(i)u(j)u(k) is obtained around the reattachment region.

Amano, R. S.↗

Studies of humid continental haze during SPACE

A concept for a solar radiometer network to provide supporting data during the Satellite Preciptiation and Cloud Experiment (SPACE) was developed. Each of the 9 prime and 10 supplementary SPACE ground sites will be equipped with an upward pointing global solar pyranometer. About half of the sites will also be equipped with upward pointing diffuse (shade ring) solar pyranometers, and a downward pointing global albedo pyranometer. These radiometers will be used to monitor the spatial and temporal variability of solar insolation and haze optical depth. The insolation data will ultimately be input to numerical models of the pre-storm and near-storm boundary layer. The optical depth data will be compared with simultaneous measurements from airborne and satellite-based passive visible radiometers and airborne lidars.

Bowdle, D. A.↗

Joining of ceramics for high temperature applications

Summarized is a literature survey of the methods for joining ceramics to ceramics or ceramics to metals for high temperature applications. Also mechanical properties and potential applications of the joints are considered. The joining of ceramics is usually carried out by brazing or diffusion bonding. Especially the latter has been found useful, increasing the application of bonded ceramics. The possibility of using electron beam and laser beam welding for joining ceramics has also recently been investigated. The bonding of ceramics has found numerous applications typical for high operating temperatures, i.e., sensors and thermocouples.

Vilpas, Martti↗

Magnetohydrodynamic waves and particle acceleration in a coronal loop

A model is presented for the acceleration of electrons in a flaring solar coronal loop. It is suggested that the primary energy release event in the flare establishes a spectrum of MHD waves within the loop which accelerate electrons by wave-particle resonant interactions. The nature of the MHD oscillations are described and numerical and analytic solution of the dispersion relation are examined. It is found that the behavior of a particle depends on whether or not it is in resonance with the wave spectrum. Resonant particles follow chaotic trajectories and can diffuse in velocity space between wave effective potentials, thereby gaining energy.

De La Beaujardiere, Jean-Francois↗

Lateral charge transport from heavy-ion tracks in integrated circuit chips

A 256K DRAM has been used to study the lateral transport of charge (electron-hole pairs) induced by direct ionization from heavy-ion tracks in an IC. The qualitative charge transport has been simulated using a two-dimensional numerical code in cylindrical coordinates. The experimental bit-map data clearly show the manifestation of lateral charge transport in the creation of adjacent multiple-bit errors from a single heavy-ion track. The heavy-ion data further demonstrate the occurrence of multiple-bit errors from single ion tracks with sufficient stopping power. The qualitative numerical simulation results suggest that electric-field-funnel-aided (drift) collection accounts for single error generated by an ion passing through a charge-collecting junction, while multiple errors from a single ion track are due to lateral diffusion of ion-generated charge.

Zoutendyk, J. A.↗

Compact finite volume methods for the diffusion equation

The paper describes an approach to treating initial-boundary-value problems by finite volume methods in which the parallel between differential and difference arguments is closely maintained. By using intrinsic geometrical properties of the volume elements, it is possible to describe discrete versions of the div, curl, and grad operators which lead, using summation-by-parts techniques, to familiar energy equations as well as the div curl = 0 and curl grad = 0 identities. For the diffusion equation, these operators describe compact schemes whose convergence is assured by the energy equations and which yield both the potential and the flux vector with second-order accuracy. A simplified potential form is especially useful for obtaining numerical results by multigrid and ADI methods.

Rose, Milton E.↗

Lidar multiple scattering: Dependence on atmospheric parameters

Backscatter lidars are very useful tools in gaining information about the atmosphere. Inversion of backscatter signals can be used to retrieve the spatial distribution of clouds and aerosols, e.g., vertical extensions of clouds and aerosol layers, detection of 'subvisible' cirrus, etc., very accurately compared to passive radiometers. To retrieve extinction coefficient profiles, in general, the so-called lidar equation is solved. This equation describes a laser pulse which propagates through an atmosphere with transmission exp(- integral of sigma(z')dz'), which is scattered singly under 180 degrees towards a receiver, and which is attenuated on its way back again by exp(- integral of sigma(z')dz'). For inversion, further information on optical properties of the scatterers are required (extinction to backscatter ratio). However, it is evident that, at least for optically thick clouds and/or large receiver's field-of-views, this approach may be erroneous because of multiply scattered photons. Several studies were performed to simulate multiple scattering in case of collimated beams or pulses including Monte Carlo, transport, and diffusion theory. Unfortunately, most of the approaches require an enormous amount of computer time. Thus, it is meaningful to elaborate approximate methods, and one such method is discussed. Our method is based on Bissonnett (1988). We have improved the numerical simulation of limited field-of-view detectors, and compared the results to Monte Carlo calculations and an exact model of Altmann (1989).

Wiegner, Matthias↗

Computational and Experimental Study of Energetic Materials in a Counterflow Microgravity Environment

Counterflow diffusion flames are studied for various fuels flowing against decomposition products from solid ammonium perchlorate (AP) pellets in order to obtain fundamental understanding of composite propellant flame structure and chemistry. We illustrate this approach through a combined experimental and numerical study of a fuel mixture consisting of C2H4 CO + H2, and C2H2 + C2H4 flowing against solid AP. For these particular AP-fuel systems, the resulting flame zone simulates the various flame structures that are ex+ to exist between reaction products from Ap crystals and a hydrocarbon binder. As in all our experimental studies, quantitative species and temperature profiles have been measured between the fuel exit and AP surface. Species measured included CN, NH, NO, OH, N2, CO2, CO, H2, CO, HCl, and H2O. Temperature was measured using a thermocouple at the exit, spontaneous Raman scattering measurements throughout the flame, OH rotational population distributions, and NO vibrational population distributions. The burning rate of AP was also measured as a function of strain rate, given by the separation distance between the AP surface and the gaseous hydrocarbon fuel tube exit plane. This distance was nominally set at 5 mm, although studies have been performed for variations in separation distance. The measured 12 scalars are compared with predictions from a detailed gas-phase kinetics model consisting of 86 species and 531 reactions. Model predictions are found to be in good agreement with experiment and illustrate the type of kinetic features that may be expected to occur in propellants when AP particle size distributions are varied. Furthermore, the results constitute the continued development of a necessary database and validation of a comprehensive model for studying more complex AP-solid fuel systems in microgravity. Exploratory studies have also been performed with liquid and solid fuels at normal gravity. Because of melting (and hence dripping) and deep thermal wave penetration into the liquid, these experiments were found feasible, but not used for obtaining quantitative data. Microgravity experiments are needed to eliminate the dripping and boiling phenomena of these systems at normal gravity. Microgravity tests in the NASA Glenn 2.2 second drop tower were performed (1) to demonstrate the feasibility of performing propellant experiments using the NASA Glenn microgravity facilities, (2) to develop the operational procedures for safe handing of the energetic materials and disposal of their toxic combustion by-products and (3) to obtain initial measurements of the AP burning rate and flame structure under microgravity conditions. Experiments were conducted on the CH4/AP system previously studied at normal gravity using a modified design of the counterflow burner and a NASA Glenn Pig Rig, i.e., one of the existing drop rigs for general-purpose usage. In these experiments, the AP burning rate was measured directly with a linear variable differential transducer (LVDT) and video imaging of the flame structure was recorded ignition was achieved by hot wires stretched across the AP surfaces. Initial drop tower combustion data show that with the same burner separation distance and flow conditions of the normal gravity experiments, the AP burning rate is approximately a factor of two lower. This difference is likely a result of radiation effects, but further tests with longer test times need to be conducted to verify that steady state conditions were achieved under microgravity conditions.

Takahashi, Fumiaki↗

Large-Scale Structure of the Molecular Gas in Taurus Revealed by High Spatial Dynamic Range Spectral Line Mapping

Viewgraph topics include: optical image of Taurus; dust extinction in IR has provided a new tool for probing cloud morphology; observations of the gas can contribute critical information on gas temperature, gas column density and distribution, mass, and kinematics; the Taurus molecular cloud complex; average spectra in each mask region; mas 2 data; dealing with mask 1 data; behavior of mask 1 pixels; distribution of CO column densities; conversion to H2 column density; variable CO/H2 ratio with values much less than 10(exp -4) at low N indicated by UV results; histogram of N(H2) distribution; H2 column density distribution in Taurus; cumulative distribution of mass and area; lower CO fractional abundance in mask 0 and 1 regions greatly increases mass determined in the analysis; masses determined with variable X(CO) and including diffuse regions agrees well with the found from L(CO); distribution of young stars as a function of molecular column density; star formation efficiency; star formation rate and gas depletion; and enlarged images of some of the regions with numerous young stars. Additional slides examine the origin of the Taurus molecular cloud, evolution from HI gas, kinematics as a clue to its origin, and its relationship to star formation.

molecular clouds↗

THE EFFECT OF ROTATION ON OSCILLATORY DOUBLE-DIFFUSIVE CONVECTION (SEMICONVECTION)

Oscillatory double-diffusive convection (ODDC, more traditionally called semiconvection) is a form of linear double-diffusive instability that occurs in fluids that are unstably stratified in temperature (Schwarzschild unstable), but stably stratified in chemical composition (Ledoux stable). This scenario is thought to be quite common in the interiors of stars and giant planets, and understanding the transport of heat and chemical species by ODDC is of great importance to stellar and planetary evolution models. Fluids unstable to ODDC have a tendency to form convective thermocompositional layers that significantly enhance the fluxes of temperature and chemical composition compared with microscopic diffusion. Although a number of recent studies have focused on studying properties of both layered and nonlayered ODDC, few have addressed how additional physical processes such as global rotation affect its dynamics. In this work, we study first how rotation affects the linear stability properties of rotating ODDC. Using direct numerical simulations, we then analyze the effect of rotation on properties of layered and nonlayered ODDC, and we study how the angle of the rotation axis with respect to the direction of gravity affects layering. We find that rotating systems can be broadly grouped into two categories based on the strength of rotation. The qualitative behavior in the more weakly rotating group is similar to nonrotating ODDC, but strongly rotating systems become dominated by vortices that are invariant in the direction of the rotation vector and strongly influence transport. We find that whenever layers form, rotation always acts to reduce thermal and compositional transport.

79 ASTRONOMY AND ASTROPHYSICS↗

Refine your search to explore more results.