Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Diffusion models”

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,171 records · Page 65

Ion radial diffusion in an electrostatic impulse model for stormtime ring current formation

Guiding-center simulations of stormtime transport of ring-current and radiation-belt ions having first adiabatic invariants mu is approximately greater than 15 MeV/G (E is approximately greater than 165 keV at L is approximately 3) are surprisingly well described (typically within a factor of approximately less than 4) by the quasilinear theory of radial diffusion. This holds even for the case of an individual model storm characterized by substorm-associated impulses in the convection electric field, provided that the actual spectrum of the electric field is incorporated in the quasilinear theory. Correction of the quasilinear diffusion coefficient D(sub LL)(sup ql) for drift-resonance broadening (so as to define D(sub LL)(sup ql)) reduced the typical discrepancy with the diffusion coefficients D(sub LL)(sup sim) deduced from guiding-center simulations of representative-particle trajectories to a factor of approximately 3. The typical discrepancy was reduced to a factor of approximately 1.4 by averaging D(sub LL)(sup sim), D(sub LL)(sup ql), and D(sub LL)(sup rb) over an ensemble of model storms characterized by different (but statistically equivalent) sets of substorm-onset times.

Chen, Margaret W.↗

Interplanetary propagation of flare-associated energetic particles

The basic propagation process of flare-associated energetic particles in interplanetary space is studied on the basis of a model which combines a Gaussian coronal injection profile and interplanetary particle densities found by a theory of focused diffusion. The model is used to describe 30 electron and proton events which originate from the western hemisphere of the sun. A comparison of calculated and observed density profiles shows that the scattering mean free path is 0.1-0.3 AU for 4-80 MeV protons. The value is two or three times smaller for 0.5-1.1 and 3-12 MeV electrons. Thus the scattering mean free path is only slightly rigidity-dependent, contrary to that predicted by the quasi-linear theory of pitch-angle scattering. The rms width is found to be less than an hour for most proton and electron events. This width, which decreases with velocity, is not rigidity-dependent.

Ma Sung, L. S.↗

Cosmic ray albedo gamma rays from the quiet sun

We estimate the flux of gamma-rays that result from collisions of high energy galactic cosmic rays with the solar atmosphere. An important aspect of our model is the propagation of cosmic rays through the magnetic fields of the inner solar systems. We use diffusion to model propagation down to the bottom of the corona. Below the corona we trace particle orbits through the photospheric fields to determine the location of cosmic ray interactions in the solar atmosphere and evolve the resultant cascades. For our nominal choice of parameters, we predict an integrated flux of gamma rays (at 1 AU) of F(E(sub gamma) greater than 100 MeV) approximately = 5 x 10(exp -8)/sq cm sec. This can be an order of magnitude above the galactic background and should be observable by the Energetic Gamma Ray experiment telescope (EGRET).

Seckel, D.↗

The X-ray shadow of the high-latitude molecular cloud MBM 12

ROSAT XRT/PSPC observations show a deep shadow cast by the high-latitude molecular cloud MBM 12 in the 3/4 keV diffuse background. Modeling of the shadow implies that less than 20 percent of the typical high-latitude 3/4 keV diffuse background intensity is emitted in front of the cloud (D = 60-70 pc). A weaker shadow consistent with the lower optical depth at higher energies was observed in the 1.5 keV band. Since little shadowing was seen in the 1/4 keV band, this observation places strong constraints on the amount of 0.5-2 keV emission that is intermixed with the source of the observed 1/4 keV flux.

Snowden, S. L.↗

A soft X-ray map of the Perseus cluster of galaxies

A 0.5-3-keV X-ray map of the Perseus cluster of galaxies is presented. The map shows a region of strong emission centered near NGC 1275 plus a highly elongated emission region which lies along the line of bright galaxies that dominates the core of the cluster. The data are compared with various models that include point and diffuse sources. One model which adequately represents the data is the superposition of a point source at NGC 1275 and an isothermal ellipsoid resulting from the bremsstrahlung emission of cluster gas. The ellipsoid has a major core radius of 20.5 arcmin and a minor core radius of 5.5 arcmin, consistent with the values obtained from galaxy counts. All acceptable models provide evidence for a compact source (less than 3 arcmin FWHM) at NGC 1275 containing about 25% of the total emission. Since the diffuse X-ray and radio components have radically different morphologies, it is unlikely that the emissions arise from a common source, as proposed in inverse-Compton models.

Cash, W.↗

Far-infrared sources in Cygnus X - An extended emission complex at DR 21 and unresolved sources at S106 and ON 2

A new map of 90 micron emission with 15 arcmin resolution is presented of a complex of sources approximately centered on DR 21. The emission is primarily from diffuse sources. Comparison of the 90 micron map with maps at 20 microns and 2.7 GHz leads to a model for diffuse sources in which there are at least two components of interstellar dust, one at 40 K and one at 165 K, responsible for the 90 micron and 20 micron emission respectively. Both components and the ionized gas must be mixed together; dust in an outer neutral shell is relatively unimportant for these sources. The diffuse sources have small infrared excesses and weak temperature gradients. Unresolved, probably compact, sources were observed at S106 and ON 2 and are briefly discussed.

Campbell, M. F.↗

Reduced-order modeling of parameterized multi-group diffusion k-eigenvalue problems

In this paper, a group-wise Reduced-Order Model, based on a Proper Orthogonal Decomposition (POD) technique, is presented for parameterized multi-group diffusion k-eigenvalue problems. The group-wise approach is shown to be more robust than model-order reduction based on monolithic multi-group POD modes. In addition, fewer POD modes are required. The method of snapshots is employed to generate the offline training set. Two Naïve and a Greedy snapshot-generation strategies are compared. The input parameter space is surveyed using Latin Hypercube Sampling techniques. The effectiveness of the group-wise model-order reduction method is demonstrated on two reactor-physics benchmarks, a two-group 3D Pressurized Water Reactor problem with 26 uncertain input parameters and a seven-group UO 2 -MOX fuel 3D mini-core with 287 uncertain input parameters. In both examples the Reduced-Order Models built using the snapshots generated by the Greedy algorithm proved to be slightly better. Model-order reduction yields results that are within 1 pcm of the full-order model, with speed-up factors of about 1000–3000, depending upon the test case.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Hydrophilic domain structure in polymer exchange membranes: Simulations of NMR spin diffusion experiments to address ability for model discrimination

Here, we detail the development of a flexible simulation program (NMR_DIFFSIM) that solves the nuclear magnetic resonance (NMR) spin diffusion equation for arbitrary polymer architectures. The program was used to explore the proton ( 1 H) NMR spin diffusion behavior predicted for a range of geometrical models describing polymer exchange membranes. These results were also directly compared with the NMR spin diffusion behavior predicted for more complex domain structures obtained from molecular dynamics (MD) simulations. The numerical implementation and capabilities of NMR_DIFFSIM were demonstrated by evaluating the experimental NMR spin diffusion behavior for the hydrophilic domain structure in sulfonated Diels-Alder Poly(Phenylene) (SDAPP) polymer membranes. The impact of morphology variations as a function of sulfonation and hydration level on the resulting NMR spin diffusion behavior were determined. These simulations allowed us to critically address the ability of NMR spin diffusion to discriminate between different structural models, and to highlight the extremely high fidelity experimental data required to accomplish this. A direct comparison of experimental double-quantum-filtered 1 H NMR spin diffusion in SDAPP membranes to the spin diffusion behavior predicted for MD-proposed morphologies revealed excellent agreement, providing experimental support for the MD structures at low to moderate hydration levels.

36 MATERIALS SCIENCE↗

Cosmic Ray Hysteresis as Evidence for Time-dependent Diffusive Processes in the Long Term Solar Modulation

A simple one-dimensional time-dependent diffusion-convection model for the modulation of cosmic rays is presented. This model predicts that the observed intensity at a given time is approximately equal to the intensity given by the time independent diffusion convection solution under interplanetary conditions which existed a time iota in the past, (U(t sub o) = U sub s(t sub o - tau)) where iota is the average time spent by a particle inside the modulating cavity. Delay times in excess of several hundred days are possible with reasonable modulation parameters. Interpretation of phase lags observed during the 1969 to 1970 solar maximum in terms of this model suggests that the modulating region is probably not less than 10 a.u. and maybe as much as 35 a.u. in extent.

Ogallagher, J. J.↗

Stochastic simulation of reaction-diffusion systems: A fluctuating-hydrodynamics approach

Here, we develop numerical methods for stochastic reaction-diffusion systems based on approaches used for fluctuating hydrodynamics (FHD). For hydrodynamic systems, the FHD formulation is formally described by stochastic partial differential equations (SPDEs). In the reaction-diffusion systems we consider, our model becomes similar to the reaction-diffusion master equation (RDME) description when our SPDEs are spatially discretized and reactions are modeled as a source term having Poisson fluctuations. However, unlike the RDME, which becomes prohibitively expensive for an increasing number of molecules, our FHD-based description naturally extends from the regime where fluctuations are strong, i.e., each mesoscopic cell has few (reactive) molecules, to regimes with moderate or weak fluctuations, and ultimately to the deterministic limit. By treating diffusion implicitly, we avoid the severe restriction on time step size that limits all methods based on explicit treatments of diffusion and construct numerical methods that are more efficient than RDME methods, without compromising accuracy. Guided by an analysis of the accuracy of the distribution of steady-state fluctuations for the linearized reaction-diffusion model, we construct several two-stage (predictor-corrector) schemes, where diffusion is treated using a stochastic Crank-Nicolson method, and reactions are handled by the stochastic simulation algorithm of Gillespie or a weakly second-order tau leaping method. We find that an implicit midpoint tau leaping scheme attains second-order weak accuracy in the linearized setting and gives an accurate and stable structure factor for a time step size of an order of magnitude larger than the hopping time scale of diffusing molecules. We study the numerical accuracy of our methods for the Schlögl reaction-diffusion model both in and out of thermodynamic equilibrium. We demonstrate and quantify the importance of thermodynamic fluctuations to the formation of a two-dimensional Turing-like pattern and examine the effect of fluctuations on three-dimensional chemical front propagation. Furthermore, by comparing stochastic simulations to deterministic reaction-diffusion simulations, we show that fluctuations accelerate pattern formation in spatially homogeneous systems and lead to a qualitatively different disordered pattern behind a traveling wave.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Calculating Mass Diffusion in High-Pressure Binary Fluids

A comprehensive mathematical model of mass diffusion has been developed for binary fluids at high pressures, including critical and supercritical pressures. Heretofore, diverse expressions, valid for limited parameter ranges, have been used to correlate high-pressure binary mass-diffusion-coefficient data. This model will likely be especially useful in the computational simulation and analysis of combustion phenomena in diesel engines, gas turbines, and liquid rocket engines, wherein mass diffusion at high pressure plays a major role.

Bellan, Josette↗

Plasma diffusion into the wake of Venus

A model of the diffusion of ionosheath plasma into the wake region of Venus is presented. It is shown that particle diffusion, which is assumed to be a consequence of the fluctuating magnetic field observed in the wake of Venus by Pioneer Venus and Veneras 9 and 10, can explain the plasma observations made in the wake by Veneras 9 and 10. The pressure due to these diffusing particles when balanced against the ionospheric pressure yield ionopause heights less than 1000 km for zenith angles less than about 135 deg. The model also predicts significant fluxes of low energy electrons and ions for zenith angles less than 135 deg, which are capable of producing the observed nightside ionosphere.

Gombosi, T.↗

Modeling the formation of methane hydrate-bearing intervals in fine-grained sediments

Sediment grain size exerts a fundamental control on how methane hydrates are distributed within the pore space. Fine-grained muds are the predominant sediments in continental margins, and hydrates in these sediments have often been observed in semi-vertical veins and fractures. In several instances, these hydrate veins/fractures are found in discrete depth intervals a few tens meters thick within the gas hydrate stability zone (GHSZ) surrounded by hydrate-free sediments above and below. As they are not obviously connected with free gas occurring beneath the base of the GHSZ, these isolated hydrate-bearing intervals have been interpreted as formed by microbial methane generated in situ. To investigate further the formation of these hydrate deposits, we applied a time-dependent advection-diffusion-reaction model that includes the effects of sedimentation, solute diffusion, and microbial methane generation. The microbial methane generation term depends on the amount of metabolizable organic carbon deposited at the seafloor, which is degraded at a prescribed rate resulting in methane formation beneath the sulfate reduction zone. In the model, methane hydrate precipitates once the dissolved methane concentration is greater than solubility, or hydrate dissolves if concentration goes below solubility. If the deposition of organic carbon at the seafloor is kept constant in time, we found that the predicted amounts of hydrate formed in discrete intervals within the GHSZ are much less than those estimated from observations. We then investigated the effect of temporal variations in the deposition of organic carbon. If greater amounts of organic carbon are deposited during some time interval, methane generation is enhanced during burial in the corresponding sediment interval. With variations in organic carbon deposition that are consistent with observations in continental margin sediments, we were able to reproduce the methane hydrate contents estimated in discrete depth intervals. Our results support the suggestion that in situ microbial methane generation is the source for hydrates within fine-grained sediments.

Malinverno, Alberto↗

Diffusion NMR methods applied to xenon gas for materials study

We report initial NMR studies of (i) xenon gas diffusion in model heterogeneous porous media and (ii) continuous flow laser-polarized xenon gas. Both areas utilize the pulsed gradient spin-echo (PGSE) techniques in the gas phase, with the aim of obtaining more sophisticated information than just translational self-diffusion coefficients--a brief overview of this area is provided in the Introduction. The heterogeneous or multiple-length scale model porous media consisted of random packs of mixed glass beads of two different sizes. We focus on observing the approach of the time-dependent gas diffusion coefficient, D(t) (an indicator of mean squared displacement), to the long-time asymptote, with the aim of understanding the long-length scale structural information that may be derived from a heterogeneous porous system. We find that D(t) of imbibed xenon gas at short diffusion times is similar for the mixed bead pack and a pack of the smaller sized beads alone, hence reflecting the pore surface area to volume ratio of the smaller bead sample. The approach of D(t) to the long-time limit follows that of a pack of the larger sized beads alone, although the limiting D(t) for the mixed bead pack is lower, reflecting the lower porosity of the sample compared to that of a pack of mono-sized glass beads. The Pade approximation is used to interpolate D(t) data between the short- and long-time limits. Initial studies of continuous flow laser-polarized xenon gas demonstrate velocity-sensitive imaging of much higher flows than can generally be obtained with liquids (20-200 mm s-1). Gas velocity imaging is, however, found to be limited to a resolution of about 1 mm s-1 owing to the high diffusivity of gases compared with liquids. We also present the first gas-phase NMR scattering, or diffusive-diffraction, data, namely flow-enhanced structural features in the echo attenuation data from laser-polarized xenon flowing through a 2 mm glass bead pack. c2002 John Wiley & Sons, Ltd.

NASA Discipline Life Sciences Technologies↗

Effects of heavy-element settling on solar neutrino fluxes and interior structure

We consider the effects of gravitational settling of both He and heavier elements on the predicted solar neutrino fluxes and interior sound speed and density profiles. We find that while the structural changes that result from the inclusion of both He and heavy-element settling are only slightly larger than the changes resulting from the inclusion of He settling alone, the additional increases in expected neutrino fluxes are of comparable size. Our preferred model with both He and heavy-element settling has neutrino count rates of 9.0 SNU for Cl-37 detectors and 137 SNU for Ga-71 detectors, as compared to 7.1 and 127 SNU for a comparable model without any diffusive separation, or 8.0 and 132 SNU for a model that includes He settling alone. We suggest that the correction factors by which the predicted neutrino fluxes of solar models calculated without including the effects of diffusion should be multiplied are 1.25 +/- 0.08 for Cl detectors, 1.07 +/- 0.02 for Ga detectors, and 1.28 +/- 0.09 for the B-8 flux (1 sigma errors). Comparison of internal sound speed and density profiles strongly suggests that the additional changes in calculated p-mode oscillation frequencies due to the inclusion of heavy-element settling will be small compared to the changes that result from He settling alone, especially for the higher degree modes. All models with diffusive separation give much better agreement with the observed depth of the convection zone than do nondiffusive models. The model that includes both He and heavy-element settling requires an initial He mass fraction Y = 0.280 and has a surface He abundance of Y = 0.251 at the solar age.

Proffitt, Charles R.↗

Factors Controlling Short-Range Methane Migration of Gas Hydrate Accumulations in Thin Coarse-Grained Layers

Natural gas hydrate is often found in marine sediment in heterogeneous distributions in different sediment types. Diffusion may be a dominant mechanism for methane migration and affect hydrate distribution. We use a 1–D advection–diffusion–reaction model to understand hydrate distribution in and surrounding thin coarse–grained layers to examine the sensitivity of four controlling factors in a diffusion–dominant gas hydrate system. These factors are the particulate organic carbon content at seafloor, the microbial reaction rate constant, the sediment grading pattern, and the cementation factor of the coarse–grained layer. We use available data at Walker Ridge 313 in the northern Gulf of Mexico where two ~3–m–thick hydrate–bearing coarse–grained layers were observed at different depths. The results show that the hydrate volume and the total amount of methane within thin, coarse–grained layers are most sensitive to the particulate organic carbon of fine–grained sediments when deposited at the seafloor. The thickness of fine–grained hydrate free zones surrounding the coarse–grained layers is most sensitive to the microbial reaction rate constant. Furthermore, it may be possible to estimate microbial reaction rate constants at other locations by studying the thickness of the hydrate free zones using the Damköhler number. In addition, we note that sediment grading patterns have a strong influence on gas hydrate occurrence within coarse–grained layers.

54 ENVIRONMENTAL SCIENCES↗

The spectral energy distribution of the scattered light from dark clouds

A dark cloud is exposed to the ambient radiation field of integrated starlight in the Galaxy. Scattering of starlight by the dust particles gives rise to a diffuse surface brightness of the dark nebula. The intensity and the spectrum of this diffuse radiation can be used to investigate, e.g., the scattering parameters of the dust, the optical thickness of the cloud, and as a probe of the ambient radiation field at the location of the cloud. An understanding of the scattering process is also a prerequisite for the isolation of broad spectral features due to fluorescence or to any other non-scattering origin of the diffuse light. Model calculations are presented for multiple scattering in a spherical cloud. These calculations show that the different spectral shapes of the observed diffuse light can be reproduced with standard dust parameters. The possibility to use the observed spectrum as a diagnostic tool for analyzing the thickness of the cloud and the dust particle is discussed.

Mattila, Kalevi↗

Fermi Large Area Telescope Second Source Catalog

We present the second catalog of high-energy gamma-ray sources detected by the Large Area Telescope (LAT), the primary science instrument on the Fermi Gamma-ray Space Telescope (Fermi), derived from data taken during the first 24 months of the science phase of the mission, which began on 2008 August 4. Source detection is based on the average flux over the 24-month period. The Second Fermi-LAT catalog (2FGL) includes source location regions, defined in terms of elliptical fits to the 95% confidence regions and spectral fits in terms of power-law, exponentially cutoff power-law, or log-normal forms. Also included are flux measurements in 5 energy bands and light curves on monthly intervals for each source. Twelve sources in the catalog are modeled as spatially extended. We provide a detailed comparison of the results from this catalog with those from the first Fermi-LAT catalog (1FGL). Although the diffuse Galactic and isotropic models used in the 2FGL analysis are improved compared to the 1FGL catalog, we attach caution flags to 162 of the sources to indicate possible confusion with residual imperfections in the diffuse model. The 2FGL catalog contains 1873 sources detected and characterized in the 100 11eV to 100 GeV range of which we consider 127 as being firmly identified and 1171 as being reliably associated with counterparts of known or likely gamma-ray-producing source classes.

Nolan, P. L.↗