Discrete and continuum links to a nonlinear coupled transport problem of interacting populations

NASA Astrophysics Data System (ADS)

Duong, M. H.; Muntean, A.; Richardson, O. M.

2017-07-01

We are interested in exploring interacting particle systems that can be seen as microscopic models for a particular structure of coupled transport flux arising when different populations are jointly evolving. The scenarios we have in mind are inspired by the dynamics of pedestrian flows in open spaces and are intimately connected to cross-diffusion and thermo-diffusion problems holding a variational structure. The tools we use include a suitable structure of the relative entropy controlling TV-norms, the construction of Lyapunov functionals and particular closed-form solutions to nonlinear transport equations, a hydrodynamics limiting procedure due to Philipowski, as well as the construction of numerical approximates to both the continuum limit problem in 2D and to the original interacting particle systems.

Ion-Scale Structure in Mercury's Magnetopause Reconnection Diffusion Region

NASA Technical Reports Server (NTRS)

Gershman, Daniel J.; Dorelli, John C.; DiBraccio, Gina A.; Raines, Jim M.; Slavin, James A.; Poh, Gangkai; Zurbuchen, Thomas H.

2016-01-01

The strength and time dependence of the electric field in a magnetopause diffusion region relate to the rate of magnetic reconnection between the solar wind and a planetary magnetic field. Here we use approximately 150 milliseconds measurements of energetic electrons from the Mercury Surface, Space Environment, GEochemistry, and Ranging (MESSENGER) spacecraft observed over Mercury's dayside polar cap boundary (PCB) to infer such small-scale changes in magnetic topology and reconnection rates. We provide the first direct measurement of open magnetic topology in flux transfer events at Mercury, structures thought to account for a significant portion of the open magnetic flux transport throughout the magnetosphere. In addition, variations in PCB latitude likely correspond to intermittent bursts of approximately 0.3 to 3 millivolts per meter reconnection electric fields separated by approximately 5 to10 seconds, resulting in average and peak normalized dayside reconnection rates of approximately 0.02 and approximately 0.2, respectively. These data demonstrate that structure in the magnetopause diffusion region at Mercury occurs at the smallest ion scales relevant to reconnection physics.

New upper limits on the local metagalactic ionizing radiation density

NASA Technical Reports Server (NTRS)

Vogel, Stuart N.; Weymann, Ray; Rauch, Michael; Hamilton, Tom

1995-01-01

We have obtained H-alpha observations with the Maryland-Caltech Fabry-Perot Spectrometer attached to the Cassegrain focus of the 1.5 m telescope at Palomer Observatory in order to set limits on the number of ionizing photons from the local metagalactic radiation field. We have observed the SW component of the Haynes-Giovanelli cloud H I 1225+01, an intergalactic cloud which should be optimum for measuring the metagalactic flux because it is nearly opaque to ionizing photons, it does not appear to be significantly shielded from the metagalactic radiation field, and the limits on embedded or nearby ionizing sources are unusually low. For the area of the cloud with an H I column density greater than 10(exp 19)/sq cm we set a 2 sigma limit of 1.1 x 10(exp -19) ergs/sq cm/s/sq arcsec (20 mR) for the surface brightness of diffuse H-alpha. This implies a 2 sigma upper limit on the incident one-sided ionizing flux of Phi(sub ex) is less than 3 x 10(exp 4)/sq cm/s. For a radiation field of the form J(sub nu) is approximately nu(exp -1.4), this yields a firm 2 sigma upper limit on the local metagalactic photoionization rate of Gamma is less than 2 x 10(exp -13)/s, and an upper limit for the radiation field J(sub nu) at the Lyman limit of J(sub nu0) is less than 8 x 10(exp -23) ergs/sq cm/Hz/sr. We discuss previous efforts to constrain the metagalactic ionizing flux using H-alpha surface brightness observations and also other methods, and conclude that our result places the firmest upper limit on this flux. We also observed the 7 min diameter region centered on 3C 273 in which H-alpha emission at a velocity of approximately 1700 km/s was initially reported by Williams and Schommer. In agreement with T. B. Williams (private communication) we find the initial detection was spurious. We obtain a 2 sigma upper limit of 1.8 x 10(exp -19) ergs/sq cm/s/sq arcsec (32 mR) for the mean surface brightness of diffuse H-alpha, about a factor of 6 below the published value.

TOPICAL REVIEW: Physics of thermoelectric cooling

NASA Astrophysics Data System (ADS)

Gurevich, Yu G.; Logvinov, G. N.

2005-12-01

A new approach is suggested to explain the Peltier effect. It assumes that the Peltier effect is not an isothermal effect. The approach is based on the occurrences of induced thermal fluxes in a structure which consists of two conducting media, through which a dc electric current flows. These induced thermal diffusion fluxes arise to compensate for the change in the thermal flux caused by the electric current (the drift thermal flux) flowing through the junction, in accordance with the general Le Châtelier-Braun principle. The occurrence of these thermal diffusion fluxes leads to temperature heterogeneity in the structure and, as a result, to a cooling or heating of the junction. Within the framework of this concept, the thermoelectric cooling is analysed. It is shown that in the general case the Peltier effect always occurs together with another thermoelectric effect. This thermoelectric effect is predicted for the first time, and we have called it the barrierless thermoelectric effect. Both these effects essentially depend on the junction surface thermal resistance. The Peltier effect disappears in the limiting case of a very large surface thermal resistance, while the barrierless effect disappears in the limiting case of a very small surface thermal resistance. The dependence of thermoelectric cooling on the geometrical dimensions of the structure is noted, and the corresponding interpretation of this fact is discussed. It is shown that the thermoelectric cooling (heating) is a thermodynamically reversible process in the linear approximation of the electric current applied.

Diffuse reflectance of TiO 2 pigmented paints: Spectral dependence of the average pathlength parameter and the forward scattering ratio

NASA Astrophysics Data System (ADS)

Vargas, William E.; Amador, Alvaro; Niklasson, Gunnar A.

2006-05-01

Diffuse reflectance spectra of paint coatings with different pigment concentrations, normally illuminated with unpolarized radiation, have been measured. A four-flux radiative transfer approach is used to model the diffuse reflectance of TiO2 (rutile) pigmented coatings through the solar spectral range. The spectral dependence of the average pathlength parameter and of the forward scattering ratio for diffuse radiation, are explicitly incorporated into this four-flux model from two novel approximations. The size distribution of the pigments has been taken into account to obtain the averages of the four-flux parameters: scattering and absorption cross sections, forward scattering ratios for collimated and isotropic diffuse radiation, and coefficients involved in the expansion of the single particle phase function in terms of Legendre polynomials.

Diffusive growth of a single droplet with three different boundary conditions

NASA Astrophysics Data System (ADS)

Tavassoli, Z.; Rodgers, G. J.

2000-02-01

We study a single, motionless three-dimensional droplet growing by adsorption of diffusing monomers on a 2D substrate. The diffusing monomers are adsorbed at the aggregate perimeter of the droplet with different boundary conditions. Models with both an adsorption boundary condition and a radiation boundary condition, as well as a phenomenological model, are considered and solved in a quasistatic approximation. The latter two models allow particle detachment. In the short time limit, the droplet radius grows as a power of the time with exponents of 1/4, 1/2 and 3/4 for the models with adsorption, radiation and phenomenological boundary conditions, respectively. In the long time limit a universal growth rate as $[t/\\ln(t)]^{1/3}$ is observed for the radius of the droplet for all models independent of the boundary conditions. This asymptotic behaviour was obtained by Krapivsky \\cite{krapquasi} where a similarity variable approach was used to treat the growth of a droplet with an adsorption boundary condition based on a quasistatic approximation. Another boundary condition with a constant flux of monomers at the aggregate perimeter is also examined. The results exhibit a power law growth rate with an exponent of 1/3 for all times.

Non-kinematic Flux-transport Dynamos Including the Effects of Diffusivity Quenching

DOE Office of Scientific and Technical Information (OSTI.GOV)

Ichimura, Chiaki; Yokoyama, Takaaki

2017-04-10

Turbulent magnetic diffusivity is quenched when strong magnetic fields suppress turbulent motion in a phenomenon known as diffusivity quenching. Diffusivity quenching can provide a mechanism for amplifying magnetic field and influencing global velocity fields through Lorentz force feedback. To investigate this effect, we conducted mean field flux-transport dynamo simulations that included the effects of diffusivity quenching in a non-kinematic regime. We found that toroidal magnetic field strength is amplified by up to approximately 1.5 times in the convection zone as a result of diffusivity quenching. This amplification is much weaker than that in kinematic cases as a result of Lorentzmore » force feedback on the system’s differential rotation. While amplified toroidal fields lead to the suppression of equatorward meridional flow locally near the base of the convection zone, large-scale equatorward transport of magnetic flux via meridional flow, which is the essential process of the flux-transport dynamo, is sustainable in our calculations.« less

CRASH: A BLOCK-ADAPTIVE-MESH CODE FOR RADIATIVE SHOCK HYDRODYNAMICS-IMPLEMENTATION AND VERIFICATION

DOE Office of Scientific and Technical Information (OSTI.GOV)

Van der Holst, B.; Toth, G.; Sokolov, I. V.

We describe the Center for Radiative Shock Hydrodynamics (CRASH) code, a block-adaptive-mesh code for multi-material radiation hydrodynamics. The implementation solves the radiation diffusion model with a gray or multi-group method and uses a flux-limited diffusion approximation to recover the free-streaming limit. Electrons and ions are allowed to have different temperatures and we include flux-limited electron heat conduction. The radiation hydrodynamic equations are solved in the Eulerian frame by means of a conservative finite-volume discretization in either one-, two-, or three-dimensional slab geometry or in two-dimensional cylindrical symmetry. An operator-split method is used to solve these equations in three substeps: (1)more » an explicit step of a shock-capturing hydrodynamic solver; (2) a linear advection of the radiation in frequency-logarithm space; and (3) an implicit solution of the stiff radiation diffusion, heat conduction, and energy exchange. We present a suite of verification test problems to demonstrate the accuracy and performance of the algorithms. The applications are for astrophysics and laboratory astrophysics. The CRASH code is an extension of the Block-Adaptive Tree Solarwind Roe Upwind Scheme (BATS-R-US) code with a new radiation transfer and heat conduction library and equation-of-state and multi-group opacity solvers. Both CRASH and BATS-R-US are part of the publicly available Space Weather Modeling Framework.« less

Transfer matrix method for four-flux radiative transfer.

PubMed

Slovick, Brian; Flom, Zachary; Zipp, Lucas; Krishnamurthy, Srini

2017-07-20

We develop a transfer matrix method for four-flux radiative transfer, which is ideally suited for studying transport through multiple scattering layers. The model predicts the specular and diffuse reflection and transmission of multilayer composite films, including interface reflections, for diffuse or collimated incidence. For spherical particles in the diffusion approximation, we derive closed-form expressions for the matrix coefficients and show remarkable agreement with numerical Monte Carlo simulations for a range of absorption values and film thicknesses, and for an example multilayer slab.

Energy and variance budgets of a diffusive staircase with implications for heat flux scaling

NASA Astrophysics Data System (ADS)

Hieronymus, M.; Carpenter, J. R.

2016-02-01

Diffusive convection, the mode of double-diffusive convection that occur when both temperature and salinity increase with increasing depth, is commonplace throughout the high latitude oceans and diffusive staircases constitute an important heat transport process in the Arctic Ocean. Heat and buoyancy fluxes through these staircases are often estimated using flux laws deduced either from laboratory experiments, or from simplified energy or variance budgets. We have done direct numerical simulations of double-diffusive convection at a range of Rayleigh numbers and quantified the energy and variance budgets in detail. This allows us to compare the fluxes in our simulations to those derived using known flux laws and to quantify how well the simplified energy and variance budgets approximate the full budgets. The fluxes are found to agree well with earlier estimates at high Rayleigh numbers, but we find large deviations at low Rayleigh numbers. The close ties between the heat and buoyancy fluxes and the budgets of thermal variance and energy have been utilized to derive heat flux scaling laws in the field of thermal convection. The result is the so called GL-theory, which has been found to give accurate heat flux scaling laws in a very wide parameter range. Diffusive convection has many similarities to thermal convection and an extension of the GL-theory to diffusive convection is also presented and its predictions are compared to the results from our numerical simulations.

PAINeT: An object-oriented software package for simulations of flow-field, transport coefficients and flux terms in non-equilibrium gas mixture flows

NASA Astrophysics Data System (ADS)

Istomin, V. A.

2018-05-01

The software package Planet Atmosphere Investigator of Non-equilibrium Thermodynamics (PAINeT) has been devel-oped for studying the non-equilibrium effects associated with electronic excitation, chemical reactions and ionization. These studies are necessary for modeling process in shock tubes, in high enthalpy flows, in nozzles or jet engines, in combustion and explosion processes, in modern plasma-chemical and laser technologies. The advantages and possibilities of the package implementation are stated. Within the framework of the package implementation, based on kinetic theory approximations (one-temperature and state-to-state approaches), calculations are carried out, and the limits of applicability of a simplified description of shock-heated air flows and any other mixtures chosen by the user are given. Using kinetic theory algorithms, a numerical calculation of the heat fluxes and relaxation terms can be performed, which is necessary for further comparison of engineering simulation with experi-mental data. The influence of state-to-state distributions over electronic energy levels on the coefficients of thermal conductivity, diffusion, heat fluxes and diffusion velocities of the components of various gas mixtures behind shock waves is studied. Using the software package the accuracy of different approximations of the kinetic theory of gases is estimated. As an example state-resolved atomic ionized mixture of N/N+/O/O+/e- is considered. It is shown that state-resolved diffusion coefficients of neutral and ionized species vary from level to level. Comparing results of engineering applications with those given by PAINeT, recommendations for adequate models selection are proposed.

DOE Office of Scientific and Technical Information (OSTI.GOV)

Poutanen, Juri, E-mail: juri.poutanen@utu.fi

Rosseland mean opacity plays an important role in theories of stellar evolution and X-ray burst models. In the high-temperature regime, when most of the gas is completely ionized, the opacity is dominated by Compton scattering. Our aim here is to critically evaluate previous works on this subject and to compute the exact Rosseland mean opacity for Compton scattering over a broad range of temperature and electron degeneracy parameter. We use relativistic kinetic equations for Compton scattering and compute the photon mean free path as a function of photon energy by solving the corresponding integral equation in the diffusion limit. Asmore » a byproduct we also demonstrate the way to compute photon redistribution functions in the case of degenerate electrons. We then compute the Rosseland mean opacity as a function of temperature and electron degeneracy and present useful approximate expressions. We compare our results to previous calculations and find a significant difference in the low-temperature regime and strong degeneracy. We then proceed to compute the flux mean opacity in both free-streaming and diffusion approximations, and show that the latter is nearly identical to the Rosseland mean opacity. We also provide a simple way to account for the true absorption in evaluating the Rosseland and flux mean opacities.« less

Global diffusive fluxes of methane in marine sediments

NASA Astrophysics Data System (ADS)

Egger, Matthias; Riedinger, Natascha; Mogollón, José M.; Jørgensen, Bo Barker

2018-06-01

Anaerobic oxidation of methane provides a globally important, yet poorly constrained barrier for the vast amounts of methane produced in the subseafloor. Here we provide a global map and budget of the methane flux and degradation in diffusion-controlled marine sediments in relation to the depth of the methane oxidation barrier. Our new budget suggests that 45-61 Tg of methane are oxidized with sulfate annually, with approximately 80% of this oxidation occurring in continental shelf sediments (<200 m water depth). Using anaerobic oxidation as a nearly quantitative sink for methane in steady-state diffusive sediments, we calculate that 3-4% of the global organic carbon flux to the seafloor is converted to methane. We further report a global imbalance of diffusive methane and sulfate fluxes into the sulfate-methane transition with no clear trend with respect to the corresponding depth of the methane oxidation barrier. The observed global mean net flux ratio between sulfate and methane of 1.4:1 indicates that, on average, the methane flux to the sulfate-methane transition accounts for only 70% of the sulfate consumption in the sulfate-methane transition zone of marine sediments.

Integral approximations to classical diffusion and smoothed particle hydrodynamics

DOE PAGES

Du, Qiang; Lehoucq, R. B.; Tartakovsky, A. M.

2014-12-31

The contribution of the paper is the approximation of a classical diffusion operator by an integral equation with a volume constraint. A particular focus is on classical diffusion problems associated with Neumann boundary conditions. By exploiting this approximation, we can also approximate other quantities such as the flux out of a domain. Our analysis of the model equation on the continuum level is closely related to the recent work on nonlocal diffusion and peridynamic mechanics. In particular, we elucidate the role of a volumetric constraint as an approximation to a classical Neumann boundary condition in the presence of physical boundary.more » The volume-constrained integral equation then provides the basis for accurate and robust discretization methods. As a result, an immediate application is to the understanding and improvement of the Smoothed Particle Hydrodynamics (SPH) method.« less

Benthic nutrient sources to hypereutrophic upper Klamath Lake, Oregon, USA.

PubMed

Kuwabara, James S; Topping, Brent R; Lynch, Dennis D; Carter, James L; Essaid, Hedeff I

2009-03-01

Three collecting trips were coordinated in April, May, and August 2006 to sample the water column and benthos of hypereutrophic Upper Klamath Lake (OR, USA) through the annual cyanophyte bloom of Aphanizomenon flos-aquae. A pore-water profiler was designed and fabricated to obtain the first high-resolution (centimeter-scale) estimates of the vertical concentration gradients of macro- and micronutrients for diffusive-flux determinations. A consistently positive benthic flux for soluble reactive phosphorus (SRP) was observed with solute release from the sediment, ranging between 0.4 and 6.1 mg/m(2)/d. The mass flux over an approximate 200-km(2) lake area was comparable in magnitude to riverine inputs. An additional concern related to fish toxicity was identified when dissolved ammonium also displayed consistently positive benthic fluxes of 4 to 134 mg/m(2)/d, again comparable to riverine inputs. Although phosphorus was a logical initial choice by water quality managers for the limiting nutrient when nitrogen-fixing cyanophytes dominate, initial trace-element results from the lake and major inflowing tributaries suggested that the role of iron limitation on primary productivity should be investigated. Dissolved iron became depleted in the lake water column during the course of the algal bloom, while dissolved ammonium and SRP increased. Elevated macroinvertebrate densities, at least of the order of 10(4) individuals/m(2), suggested that the diffusive-flux estimates may be significantly enhanced by bioturbation. In addition, heat-flux modeling indicated that groundwater advection of nutrients could also significantly contribute to internal nutrient loading. Accurate environmental assessments of lentic systems and reasonable expectations for point-source management require quantitative consideration of internal solute sources.

Benthic nutrient sources to hypereutrophic Upper Klamath Lake, Oregon, USA

USGS Publications Warehouse

Kuwabara, J.S.; Topping, B.R.; Lynch, D.D.; Carter, J.L.; Essaid, H.I.

2009-01-01

Three collecting trips were coordinated in April, May, and August 2006 to sample the water column and benthos of hypereutrophic Upper Klamath Lake (OR, USA) through the annual cyanophyte bloom of Aphanizomenon flos-aquae. A porewater profiler was designed and fabricated to obtain the first high-resolution (centimeter-scale) estimates of the vertical, concentration gradients of macro- and micronutrients for diffusive-flux determinations. A consistently positive benthic flux for soluble reactive phosphorus (SRP) was observed with solute release from the sediment, ranging between 0.4 and 6.1 mg/m2/d. The mass flux over an approximate 200-km2 lake area was comparable in magnitude to riverine inputs. An additional concern, related to fish toxicity was identified when dissolved ammonium also displayed consistently positive benthic fluxes of 4 to 134 mg/m2/d, again, comparable to riverine inputs. Although phosphorus was a logical initial choice by water quality managers for the limiting nutrient when nitrogen-fixing cyanophytes dominate, initial trace-element results from the lake and major inflowing tributaries suggested that the role of iron limitation on primary productivity should be investigated. Dissolved iron became depleted in the lake water column during the course of the algal bloom, while dissolved ammonium and SRP increased. Elevated macroinvertebrate densities, at least of the order of 104 individuals/m2, suggested, that the diffusive-flux estimates may be significantly enhanced, by bioturbation. In addition, heat-flux modeling indicated that groundwater advection of nutrients could also significantly contribute to internal nutrient loading. Accurate environmental assessments of lentic systems and reasonable expectations for point-source management require quantitative consideration of internal solute sources ?? 2009 SETAC.

Relation between the neutrino flux from Centaurus A and the associated diffuse neutrino flux

DOE Office of Scientific and Technical Information (OSTI.GOV)

Koers, Hylke B. J.; Tinyakov, Peter; Institute for Nuclear Research, 60th October Anniversary Prospect 7a, 117312, Moscow

2008-10-15

Based on recent results obtained by the Pierre Auger Observatory (PAO), it has been hypothesized that Centaurus A (Cen A) is a source of ultrahigh-energy cosmic rays (UHECRs) and associated neutrinos. We point out that the diffuse neutrino flux may be used to constrain the source model if one assumes that the ratio between the UHECR and neutrino fluxes outputted by Cen A is representative for other sources. Under this assumption we investigate the relation between the neutrino flux from Cen A and the diffuse neutrino flux. Assuming furthermore that Cen A is the source of two UHECR events observedmore » by PAO, we estimate the all-sky diffuse neutrino flux to be {approx}200-5000 times larger than the neutrino flux from Cen A. As a result, the diffuse neutrino fluxes associated with some of the recently proposed models of UHECR-related neutrino production in Cen A are above existing limits. Regardless of the underlying source model, our results indicate that the detection of neutrinos from Cen A without the accompanying diffuse flux would mean that Cen A is an exceptionally efficient neutrino source.« less

Quantification of natural vapor fluxes of trichloroethene in the unsaturated zone at Picatinny Arsenal, New Jersey

USGS Publications Warehouse

Smith, James A.; Tisdale, Amy K.; Cho, H. Jean

1996-01-01

The upward flux of trichloroethene (TCE) vapor through the unsaturated zone above a contaminated, water-table aquifer at Picatinny Arsenal, New Jersey, has been studied under natural conditions over a 12-month period. Vertical gas-phase diffusion fluxes were estimated indirectly by measuring the TCE vapor concentration gradient in the unsaturated zone and using Fick's law to calculate the flux. The total gas-phase flux (e.g., the sum of diffusion and advection fluxes) was measured directly with a vertical flux chamber (VFC). In many cases, the upward TCE vapor flux was several orders of magnitude greater than the upward TCE diffusion flux, suggesting that mechanisms other than steady-state vapor diffusion are contributing to the vertical transport of TCE vapors through the unsaturated zone. The measured total flux of TCE vapor from the subsurface to the atmosphere is approximately 50 kg/yr and is comparable in magnitude to the removal rate of TCE from the aquifer by an existing pump-and-treat system and by discharge into a nearby stream. The net upward flux of TCE is reduced significantly during a storm event, presumably due to the mass transfer of TCE from the soil gas to the infiltrating rainwater and its subsequent downward advection. Several potential problems associated with the measurement of total gas-phase fluxes are discussed.

Predicting the drying properties of sludge based on hydrothermal treatment under subcritical conditions.

PubMed

Mäkelä, Mikko; Fraikin, Laurent; Léonard, Angélique; Benavente, Verónica; Fullana, Andrés

2016-03-15

The effects of hydrothermal treatment on the drying properties of sludge were determined. Sludge was hydrothermally treated at 180-260 °C for 0.5-5 h using NaOH and HCl as additives to influence reaction conditions. Untreated sludge and attained hydrochar samples were then dried under identical conditions with a laboratory microdryer and an X-ray microtomograph was used to follow changes in sample dimensions. The effective moisture diffusivities of sludge and hydrochar samples were determined and the effect of process conditions on respective mean diffusivities evaluated using multiple linear regression. Based on the results the drying time of untreated sludge decreased from approximately 80 min to 37-59 min for sludge hydrochar. Drying of untreated sludge was governed by the falling rate period where drying flux decreased continuously as a function of sludge moisture content due to heat and mass transfer limitations and sample shrinkage. Hydrothermal treatment increased the drying flux of sludge hydrochar and decreased the effect of internal heat and mass transfer limitations and sample shrinkage especially at higher treatment temperatures. The determined effective moisture diffusivities of sludge and hydrochar increased as a function of decreasing moisture content and the mean diffusivity of untreated sludge (8.56·10(-9) m(2) s(-1)) and sludge hydrochar (12.7-27.5·10(-9) m(2) s(-1)) were found statistically different. The attained regression model indicated that treatment temperature governed the mean diffusivity of hydrochar, as the effects of NaOH and HCl were statistically insignificant. The attained results enabled prediction of sludge drying properties through mean moisture diffusivity based on hydrothermal treatment conditions. Copyright © 2015 Elsevier Ltd. All rights reserved.

Calculation of Transfer Functions of Multilayer Biotissues in the Problems of Correction of Their Fluorescence Spectra

NASA Astrophysics Data System (ADS)

Lysenko, S. A.

2018-01-01

A method for rapid calculation of a flux of stimulated fluorescence of a multilayer optically dense medium with inhomogeneous distribution of the fluorophore has been developed. The light field in the medium at the excitation wavelength of fluorescence is represented by a superposition of incident collimated, incident diffuse, and reflected diffuse fluxes. A two-stream approximation is used to describe the light field in the medium at the wavelength of emission of the fluorescence. Fluxes in adjacent elementary layers of the medium and on its surface are connected by simple matrix operators that are obtained using a combination of engineering approaches of radiation-transfer theory and single-scattering approximation. The calculations of fluorescence fluxes of a four-layer biotissue that are excited and recorded at 400-800 nm are compared with their Monte Carlo simulation with a discrepancy of 1%. The effect of the propagation medium on the fluorescence spectra of 5-ALA-induced protoporphyrin IX that are recorded from human skin was studied, and a technique for their correction that is based on measurements and quantitative analysis of the diffuse reflectance spectrum of skin was proposed.

COMPARISON OF IMPLICIT SCHEMES TO SOLVE EQUATIONS OF RADIATION HYDRODYNAMICS WITH A FLUX-LIMITED DIFFUSION APPROXIMATION: NEWTON–RAPHSON, OPERATOR SPLITTING, AND LINEARIZATION

DOE Office of Scientific and Technical Information (OSTI.GOV)

Tetsu, Hiroyuki; Nakamoto, Taishi, E-mail: h.tetsu@geo.titech.ac.jp

Radiation is an important process of energy transport, a force, and a basis for synthetic observations, so radiation hydrodynamics (RHD) calculations have occupied an important place in astrophysics. However, although the progress in computational technology is remarkable, their high numerical cost is still a persistent problem. In this work, we compare the following schemes used to solve the nonlinear simultaneous equations of an RHD algorithm with the flux-limited diffusion approximation: the Newton–Raphson (NR) method, operator splitting, and linearization (LIN), from the perspective of the computational cost involved. For operator splitting, in addition to the traditional simple operator splitting (SOS) scheme,more » we examined the scheme developed by Douglas and Rachford (DROS). We solve three test problems (the thermal relaxation mode, the relaxation and the propagation of linear waves, and radiating shock) using these schemes and then compare their dependence on the time step size. As a result, we find the conditions of the time step size necessary for adopting each scheme. The LIN scheme is superior to other schemes if the ratio of radiation pressure to gas pressure is sufficiently low. On the other hand, DROS can be the most efficient scheme if the ratio is high. Although the NR scheme can be adopted independently of the regime, especially in a problem that involves optically thin regions, the convergence tends to be worse. In all cases, SOS is not practical.« less

Preferential flow, diffuse flow, and perching in an interbedded fractured-rock unsaturated zone

NASA Astrophysics Data System (ADS)

Nimmo, John R.; Creasey, Kaitlyn M.; Perkins, Kim S.; Mirus, Benjamin B.

2017-03-01

Layers of strong geologic contrast within the unsaturated zone can control recharge and contaminant transport to underlying aquifers. Slow diffuse flow in certain geologic layers, and rapid preferential flow in others, complicates the prediction of vertical and lateral fluxes. A simple model is presented, designed to use limited geological site information to predict these critical subsurface processes in response to a sustained infiltration source. The model is developed and tested using site-specific information from the Idaho National Laboratory in the Eastern Snake River Plain (ESRP), USA, where there are natural and anthropogenic sources of high-volume infiltration from floods, spills, leaks, wastewater disposal, retention ponds, and hydrologic field experiments. The thick unsaturated zone overlying the ESRP aquifer is a good example of a sharply stratified unsaturated zone. Sedimentary interbeds are interspersed between massive and fractured basalt units. The combination of surficial sediments, basalts, and interbeds determines the water fluxes through the variably saturated subsurface. Interbeds are generally less conductive, sometimes causing perched water to collect above them. The model successfully predicts the volume and extent of perching and approximates vertical travel times during events that generate high fluxes from the land surface. These developments are applicable to sites having a thick, geologically complex unsaturated zone of substantial thickness in which preferential and diffuse flow, and perching of percolated water, are important to contaminant transport or aquifer recharge.

Preferential flow, diffuse flow, and perching in an interbedded fractured-rock unsaturated zone

USGS Publications Warehouse

Nimmo, John R.; Creasey, Kaitlyn M; Perkins, Kimberlie; Mirus, Benjamin B.

2017-01-01

Layers of strong geologic contrast within the unsaturated zone can control recharge and contaminant transport to underlying aquifers. Slow diffuse flow in certain geologic layers, and rapid preferential flow in others, complicates the prediction of vertical and lateral fluxes. A simple model is presented, designed to use limited geological site information to predict these critical subsurface processes in response to a sustained infiltration source. The model is developed and tested using site-specific information from the Idaho National Laboratory in the Eastern Snake River Plain (ESRP), USA, where there are natural and anthropogenic sources of high-volume infiltration from floods, spills, leaks, wastewater disposal, retention ponds, and hydrologic field experiments. The thick unsaturated zone overlying the ESRP aquifer is a good example of a sharply stratified unsaturated zone. Sedimentary interbeds are interspersed between massive and fractured basalt units. The combination of surficial sediments, basalts, and interbeds determines the water fluxes through the variably saturated subsurface. Interbeds are generally less conductive, sometimes causing perched water to collect above them. The model successfully predicts the volume and extent of perching and approximates vertical travel times during events that generate high fluxes from the land surface. These developments are applicable to sites having a thick, geologically complex unsaturated zone of substantial thickness in which preferential and diffuse flow, and perching of percolated water, are important to contaminant transport or aquifer recharge.

Border-Crossing Model for the Diffusive Coarsening of Wet Foams

NASA Astrophysics Data System (ADS)

Durian, Douglas; Schimming, Cody

For dry foams, the transport of gas from small high-pressure bubbles to large low-pressure bubbles is dominated by diffusion across the thin soap films separating neighboring bubbles. For wetter foams, the film areas become smaller as the Plateau borders and vertices inflate with liquid. So-called ``border-blocking'' models can explain some features of wet-foam coarsening based on the presumption that the inflated borders totally block the gas flux; however, this approximation dramatically fails in the wet/unjamming limit where the bubbles become close-packed spheres. Here, we account for the ever-present border-crossing flux by a new length scale defined by the average gradient of gas concentration inside the borders. We argue that it is proportional to the geometric average of film and border thicknesses, and we verify this scaling and the numerical prefactor by numerical solution of the diffusion equation. Then we show how the dA / dt =K0 (n - 6) von Neumann law is modified by the appearance of terms that depend on bubble size and shape as well as the concentration gradient length scale. Finally, we use the modified von Neumann law to compute the growth rate of the average bubble, which is not constant.

Anisotropic diffusion in mesh-free numerical magnetohydrodynamics

NASA Astrophysics Data System (ADS)

Hopkins, Philip F.

2017-04-01

We extend recently developed mesh-free Lagrangian methods for numerical magnetohydrodynamics (MHD) to arbitrary anisotropic diffusion equations, including: passive scalar diffusion, Spitzer-Braginskii conduction and viscosity, cosmic ray diffusion/streaming, anisotropic radiation transport, non-ideal MHD (Ohmic resistivity, ambipolar diffusion, the Hall effect) and turbulent 'eddy diffusion'. We study these as implemented in the code GIZMO for both new meshless finite-volume Godunov schemes (MFM/MFV). We show that the MFM/MFV methods are accurate and stable even with noisy fields and irregular particle arrangements, and recover the correct behaviour even in arbitrarily anisotropic cases. They are competitive with state-of-the-art AMR/moving-mesh methods, and can correctly treat anisotropic diffusion-driven instabilities (e.g. the MTI and HBI, Hall MRI). We also develop a new scheme for stabilizing anisotropic tensor-valued fluxes with high-order gradient estimators and non-linear flux limiters, which is trivially generalized to AMR/moving-mesh codes. We also present applications of some of these improvements for SPH, in the form of a new integral-Godunov SPH formulation that adopts a moving-least squares gradient estimator and introduces a flux-limited Riemann problem between particles.

Inter-cusp Ion and Electron Transport in a Nstar-derivative Ion Thruster

NASA Technical Reports Server (NTRS)

Foster, John E.

2001-01-01

Diffusion of electrons and ions to anode surfaces between the magnetic cusps of a NASA Solar Electric Propulsion Technology Application Readiness ion thruster has been characterized. Ion flux measurements were made at the anode and at the screen grid electrode. The measurements indicated that the average ion current density at the anode and at the screen grid were approximately equal. Additionally, it was found that the electron flux to the anode between cusps is best described by the classical cross-field diffusion coefficient.

Ion and Electron Transport in an Nstar-derivative Ion Thruster. Revised

NASA Technical Reports Server (NTRS)

Foster, John E.

2001-01-01

Diffusion of electrons and ions to anode surfaces between the magnetic cusps of a NASA Solar Electric Propulsion Technology Application Readiness ion thruster has been characterized. Ion flux measurements were made at the anode and at the screen grid electrode. The measurements indicated that the average ion current density at the anode and at the screen grid were approximately equal. Additionally, it was found that the electron flux to the anode between cusps is best described by the classical cross-field diffusion coefficient.

Nonlinear restrictions on dynamo action. [in magnetic fields of astrophysical objects

NASA Technical Reports Server (NTRS)

Vainshtein, Samuel I.; Cattaneo, Fausto

1992-01-01

Astrophysical dynamos operate in the limit of small magnetic diffusivity. In order for magnetic reconnection to occur, very small magnetic structures must form so that diffusion becomes effective. The formation of small-scale fields is accompanied by the stretching of the field lines and therefore by an amplification of the magnetic field strength. The back reaction of the magnetic field on the motions leads to the eventual saturation of the dynamo process, thus posing a constraint on the amount of magnetic flux that can be generated by dynamo action, It is argued that in the limit of small diffusivity only a small amount of flux, many orders of magnitude less than the observed fluxes, can be created by dynamo processes.

The diffuse extreme-ultraviolet background - Constraints on hot coronal plasma

NASA Technical Reports Server (NTRS)

Paresce, F.; Stern, R.

1981-01-01

The Apollo-Soyuz data and data reported by Cash et al. (1976) have been reanalyzed in terms of both isothermal models and temperature distribution models. In the latter case, a power-law form is assumed for the relation between emission measure and temperature. A new upper limit on diffuse flux in the 20-73 eV band derived from Apollo-Soyuz observations made in the earth's shadow has been incorporated in the calculation. In the considered investigation the results of the new analysis are presented and the implications for the physical properties of the hot component of the interstellar medium are discussed. The analysis of the Berkeley extreme ultraviolet (EUV) diffuse background measurements using either isothermal or power law temperature distribution models for the emitting plasma indicates excellent qualitative agreement with hard X-ray data that suggest the sun to be immersed in a hot plasma that pervades most of space out to approximately 100 pc.

Monte Carlo charged-particle tracking and energy deposition on a Lagrangian mesh.

PubMed

Yuan, J; Moses, G A; McKenty, P W

2005-10-01

A Monte Carlo algorithm for alpha particle tracking and energy deposition on a cylindrical computational mesh in a Lagrangian hydrodynamics code used for inertial confinement fusion (ICF) simulations is presented. The straight line approximation is used to follow propagation of "Monte Carlo particles" which represent collections of alpha particles generated from thermonuclear deuterium-tritium (DT) reactions. Energy deposition in the plasma is modeled by the continuous slowing down approximation. The scheme addresses various aspects arising in the coupling of Monte Carlo tracking with Lagrangian hydrodynamics; such as non-orthogonal severely distorted mesh cells, particle relocation on the moving mesh and particle relocation after rezoning. A comparison with the flux-limited multi-group diffusion transport method is presented for a polar direct drive target design for the National Ignition Facility. Simulations show the Monte Carlo transport method predicts about earlier ignition than predicted by the diffusion method, and generates higher hot spot temperature. Nearly linear speed-up is achieved for multi-processor parallel simulations.

Spatial Fluctuations in the Diffuse Cosmic X-Ray Background. Ph.D. Thesis

NASA Technical Reports Server (NTRS)

Shafer, R. A.

1983-01-01

The bright, essentially isotropic, X-ray sky flux above 2 keV yields information on the universe at large distances. However, a definitive understanding of the origin of the flux is lacking. Some fraction of the total flux is contributed by active galactic nuclei and clusters of galaxies, but less than one percent of the total is contributed by the or approximately 3 keV band resolved sources, which is the band where the sky flux is directly observed. Parametric models of AGN (quasar) luminosity function evolution are examined. Most constraints are by the total sky flux. The acceptability of particular models hinges on assumptions currently not directly testable. The comparison with the Einstein Observatory 1 to keV low flux source counts is hampered by spectral uncertainties. A tentative measurement of a large scale dipole anisotropy is consistent with the velocity and direction derived from the dipole in the microwave background. The impact of the X-ray anisotropy limits for other scales on studies of large-scale structure in the universe is sketched. Models of the origins of the X-ray sky flux are reviewed, and future observational programs outlined.

Local bounds preserving stabilization for continuous Galerkin discretization of hyperbolic systems

NASA Astrophysics Data System (ADS)

Mabuza, Sibusiso; Shadid, John N.; Kuzmin, Dmitri

2018-05-01

The objective of this paper is to present a local bounds preserving stabilized finite element scheme for hyperbolic systems on unstructured meshes based on continuous Galerkin (CG) discretization in space. A CG semi-discrete scheme with low order artificial dissipation that satisfies the local extremum diminishing (LED) condition for systems is used to discretize a system of conservation equations in space. The low order artificial diffusion is based on approximate Riemann solvers for hyperbolic conservation laws. In this case we consider both Rusanov and Roe artificial diffusion operators. In the Rusanov case, two designs are considered, a nodal based diffusion operator and a local projection stabilization operator. The result is a discretization that is LED and has first order convergence behavior. To achieve high resolution, limited antidiffusion is added back to the semi-discrete form where the limiter is constructed from a linearity preserving local projection stabilization operator. The procedure follows the algebraic flux correction procedure usually used in flux corrected transport algorithms. To further deal with phase errors (or terracing) common in FCT type methods, high order background dissipation is added to the antidiffusive correction. The resulting stabilized semi-discrete scheme can be discretized in time using a wide variety of time integrators. Numerical examples involving nonlinear scalar Burgers equation, and several shock hydrodynamics simulations for the Euler system are considered to demonstrate the performance of the method. For time discretization, Crank-Nicolson scheme and backward Euler scheme are utilized.

Multiyear search for a diffuse flux of muon neutrinos with AMANDA-II

DOE Office of Scientific and Technical Information (OSTI.GOV)

Achterberg, A.; Duvoort, M. R.; Heise, J.

2007-08-15

A search for TeV-PeV muon neutrinos from unresolved sources was performed on AMANDA-II data collected between 2000 and 2003 with an equivalent live time of 807 days. This diffuse analysis sought to find an extraterrestrial neutrino flux from sources with nonthermal components. The signal is expected to have a harder spectrum than the atmospheric muon and neutrino backgrounds. Since no excess of events was seen in the data over the expected background, an upper limit of E{sup 2}{phi}{sub 90percentC.L.}<7.4x10{sup -8} GeV cm{sup -2} s{sup -1} sr{sup -1} is placed on the diffuse flux of muon neutrinos with a {phi}{proportional_to}E{sup -2}more » spectrum in the energy range 16 TeV to 2.5 PeV. This is currently the most sensitive {phi}{proportional_to}E{sup -2} diffuse astrophysical neutrino limit. We also set upper limits for astrophysical and prompt neutrino models, all of which have spectra different from {phi}{proportional_to}E{sup -2}.« less

Martian CH(4): sources, flux, and detection.

PubMed

Onstott, T C; McGown, D; Kessler, J; Lollar, B Sherwood; Lehmann, K K; Clifford, S M

2006-04-01

Recent observations have detected trace amounts of CH(4) heterogeneously distributed in the martian atmosphere, which indicated a subsurface CH(4) flux of ~2 x 10(5) to 2 x 10(9) cm(2) s(1). Four different origins for this CH(4) were considered: (1) volcanogenic; (2) sublimation of hydrate- rich ice; (3) diffusive transport through hydrate-saturated cryosphere; and (4) microbial CH(4) generation above the cryosphere. A diffusive flux model of the martian crust for He, H(2), and CH(4) was developed based upon measurements of deep fracture water samples from South Africa. This model distinguishes between abiogenic and microbial CH(4) sources based upon their isotopic composition, and couples microbial CH(4) production to H(2) generation by H(2)O radiolysis. For a He flux of approximately 10(5) cm(2) s(1) this model yields an abiogenic CH(4) flux and a microbial CH(4) flux of approximately 10(6) and approximately 10(9) cm(2) s(1), respectively. This flux will only reach the martian surface if CH(4) hydrate is saturated in the cryosphere; otherwise it will be captured within the cryosphere. The sublimation of a hydrate-rich cryosphere could generate the observed CH(4) flux, whereas microbial CH(4) production in a hypersaline environment above the hydrate stability zone only seems capable of supplying approximately 10(5) cm(2) s(1) of CH(4). The model predicts that He/H(2)/CH(4)/C(2)H(6) abundances and the C and H isotopic values of CH(4) and the C isotopic composition of C(2)H(6) could reveal the different sources. Cavity ring-down spectrometers represent the instrument type that would be most capable of performing the C and H measurements of CH(4) on near future rover missions and pinpointing the cause and source of the CH(4) emissions.

Bounded fractional diffusion in geological media: Definition and Lagrangian approximation

NASA Astrophysics Data System (ADS)

Zhang, Yong; Green, Christopher T.; LaBolle, Eric M.; Neupauer, Roseanna M.; Sun, HongGuang

2016-11-01

Spatiotemporal fractional-derivative models (FDMs) have been increasingly used to simulate non-Fickian diffusion, but methods have not been available to define boundary conditions for FDMs in bounded domains. This study defines boundary conditions and then develops a Lagrangian solver to approximate bounded, one-dimensional fractional diffusion. Both the zero-value and nonzero-value Dirichlet, Neumann, and mixed Robin boundary conditions are defined, where the sign of Riemann-Liouville fractional derivative (capturing nonzero-value spatial-nonlocal boundary conditions with directional superdiffusion) remains consistent with the sign of the fractional-diffusive flux term in the FDMs. New Lagrangian schemes are then proposed to track solute particles moving in bounded domains, where the solutions are checked against analytical or Eulerian solutions available for simplified FDMs. Numerical experiments show that the particle-tracking algorithm for non-Fickian diffusion differs from Fickian diffusion in relocating the particle position around the reflective boundary, likely due to the nonlocal and nonsymmetric fractional diffusion. For a nonzero-value Neumann or Robin boundary, a source cell with a reflective face can be applied to define the release rate of random-walking particles at the specified flux boundary. Mathematical definitions of physically meaningful nonlocal boundaries combined with bounded Lagrangian solvers in this study may provide the only viable techniques at present to quantify the impact of boundaries on anomalous diffusion, expanding the applicability of FDMs from infinite domains to those with any size and boundary conditions.

3D Monte-Carlo study of toroidally discontinuous limiter SOL configurations of Aditya tokamak

NASA Astrophysics Data System (ADS)

Sahoo, Bibhu Prasad; Sharma, Devendra; Jha, Ratneshwar; Feng, Yühe

2017-08-01

The plasma-neutral transport in the scrape-off layer (SOL) region formed by toroidally discontinuous limiters deviates from usual uniform SOL approximations when 3D effects caused by limiter discreteness begin to dominate. In an upgrade version of the Aditya tokamak, originally having a toroidally localized poloidal ring-like limiter, the newer outboard block and inboard belt limiters are expected to have smaller connection lengths and a multiple fold toroidal periodicity. The characteristics of plasma discharges may accordingly vary from the original observations of large diffusivity, and a net improvement and the stability of the discharges are desired. The estimations related to 3D effects in the ring limiter plasma transport are also expected to be modified and are updated by predictive simulations of transport in the new block limiter configuration. A comparison between the ring limiter results and those from new simulations with block limiter SOL shows that for the grids produced using same core plasma equilibrium, the modified SOL plasma flows and flux components have enhanced poloidal periodicity in the block limiter case. These SOL modifications result in a reduced net recycling for the equivalent edge density values. Predictions are also made about the relative level of the diffusive transport and its impact on the factors limiting the operational regime.

Calculation of the neutron diffusion equation by using Homotopy Perturbation Method

DOE Office of Scientific and Technical Information (OSTI.GOV)

Koklu, H., E-mail: koklu@gantep.edu.tr; Ozer, O.; Ersoy, A.

The distribution of the neutrons in a nuclear fuel element in the nuclear reactor core can be calculated by the neutron diffusion theory. It is the basic and the simplest approximation for the neutron flux function in the reactor core. In this study, the neutron flux function is obtained by the Homotopy Perturbation Method (HPM) that is a new and convenient method in recent years. One-group time-independent neutron diffusion equation is examined for the most solved geometrical reactor core of spherical, cubic and cylindrical shapes, in the frame of the HPM. It is observed that the HPM produces excellent resultsmore » consistent with the existing literature.« less

Contribution of sediment fluxes and transformations to the summer nitrogen budget of an Upper Mississippi River backwater system

USGS Publications Warehouse

James, W.F.; Richardson, W.B.; Soballe, D.M.

2008-01-01

Routing nitrate through backwaters of regulated floodplain rivers to increase retention could decrease loading to nitrogen (N)-sensitive coastal regions. Sediment core determinations of N flux were combined with inflow-outflow fluxes to develop mass balance approximations of N uptake and transformations in a flow-controlled backwater of the Upper Mississippi River (USA). Inflow was the dominant nitrate source (>95%) versus nitrification and varied as a function of source water concentration since flow was constant. Nitrate uptake length increased linearly, while uptake velocity decreased linearly, with increasing inflow concentration to 2 mg l-1, indicating limitation of N uptake by loading. N saturation at higher inflow concentration coincided with maximum uptake capacity, 40% uptake efficiency, and an uptake length 2 times greater than the length of the backwater. Nitrate diffusion and denitrification in sediment accounted for 27% of the backwater nitrate retention, indicating that assimilation by other biota or denitrification on other substrates were the dominant uptake mechanisms. Ammonium export from the backwater was driven by diffusive efflux from the sediment. Ammonium increased from near zero at the inflow to a maximum mid-lake, then declined slightly toward the outflow due to uptake during transport. Ammonium export was small compared to nitrate retention. ?? 2007 Springer Science+Business Media B.V.

Gradual Diffusion and Punctuated Phase Space Density Enhancements of Highly Relativistic Electrons: Van Allen Probes Observations

NASA Technical Reports Server (NTRS)

Baker, D. N.; Jaynes, A. N.; Li, X.; Henderson, M. G.; Kanekal, S. G.; Reeves, G. D.; Spence, H. E.; Claudepierre, S. G.; Fennell, J. F.; Hudson, M. K.

2014-01-01

The dual-spacecraft Van Allen Probes mission has provided a new window into mega electron volt (MeV) particle dynamics in the Earth's radiation belts. Observations (up to E (is) approximately 10MeV) show clearly the behavior of the outer electron radiation belt at different timescales: months-long periods of gradual inward radial diffusive transport and weak loss being punctuated by dramatic flux changes driven by strong solar wind transient events. We present analysis of multi-MeV electron flux and phase space density (PSD) changes during March 2013 in the context of the first year of Van Allen Probes operation. This March period demonstrates the classic signatures both of inward radial diffusive energization and abrupt localized acceleration deep within the outer Van Allen zone (L (is) approximately 4.0 +/- 0.5). This reveals graphically that both 'competing' mechanisms of multi-MeV electron energization are at play in the radiation belts, often acting almost concurrently or at least in rapid succession.

Bounded fractional diffusion in geological media: Definition and Lagrangian approximation

USGS Publications Warehouse

Zhang, Yong; Green, Christopher T.; LaBolle, Eric M.; Neupauer, Roseanna M.; Sun, HongGuang

2016-01-01

Spatiotemporal Fractional-Derivative Models (FDMs) have been increasingly used to simulate non-Fickian diffusion, but methods have not been available to define boundary conditions for FDMs in bounded domains. This study defines boundary conditions and then develops a Lagrangian solver to approximate bounded, one-dimensional fractional diffusion. Both the zero-value and non-zero-value Dirichlet, Neumann, and mixed Robin boundary conditions are defined, where the sign of Riemann-Liouville fractional derivative (capturing non-zero-value spatial-nonlocal boundary conditions with directional super-diffusion) remains consistent with the sign of the fractional-diffusive flux term in the FDMs. New Lagrangian schemes are then proposed to track solute particles moving in bounded domains, where the solutions are checked against analytical or Eularian solutions available for simplified FDMs. Numerical experiments show that the particle-tracking algorithm for non-Fickian diffusion differs from Fickian diffusion in relocating the particle position around the reflective boundary, likely due to the non-local and non-symmetric fractional diffusion. For a non-zero-value Neumann or Robin boundary, a source cell with a reflective face can be applied to define the release rate of random-walking particles at the specified flux boundary. Mathematical definitions of physically meaningful nonlocal boundaries combined with bounded Lagrangian solvers in this study may provide the only viable techniques at present to quantify the impact of boundaries on anomalous diffusion, expanding the applicability of FDMs from infinite do mains to those with any size and boundary conditions.

Survey of upper band chorus and ECH waves: Implications for the diffuse aurora

NASA Astrophysics Data System (ADS)

Meredith, Nigel; Horne, Richard; Thorne, Richard; Anderson, Roger

2010-05-01

The origin of the diffuse aurora has been a source of controversy for many years. More recently the question has taken a new significance in view of the associated changes in atmospheric chemistry which may affect the middle atmosphere. Here we use CRRES data to assess the importance of upper band chorus and electron cyclotron harmonic (ECH) waves in the production of the diffuse aurora. Both wave modes increase with increasing geomagnetic activity, suggesting they are related to periods of enhanced convection and/or substorm activity. They are confined to the near-equatorial region which excludes the pre-noon sector from the wave survey. During active conditions intense ECH waves and upper band chorus, with amplitudes exceeding 1 mVm-1, are observed in the region 4 < L < 7 from 2100 to 0600 MLT approximately 20% and 6% of the time respectively. This suggests that both wave modes can put electrons on strong diffusion, but only during active conditions and not at all local times. Scattering rates fall below the strong diffusion limit at other times when the wave amplitudes are weaker. Fluxes of low energy electrons (100 eV < E < 30 keV) also increase with increasing geomagnetic activity in approximately the same region of geospace as the waves, suggesting that these electrons are responsible for the generation of the waves. The patterns of the upper band chorus, ECH waves and low energy electrons are similar to the global morphology of the diffuse aurora, suggesting that both wave modes play significant roles in the production of the diffuse aurora.

The arbitrary order mixed mimetic finite difference method for the diffusion equation

DOE PAGES

Gyrya, Vitaliy; Lipnikov, Konstantin; Manzini, Gianmarco

2016-05-01

Here, we propose an arbitrary-order accurate mimetic finite difference (MFD) method for the approximation of diffusion problems in mixed form on unstructured polygonal and polyhedral meshes. As usual in the mimetic numerical technology, the method satisfies local consistency and stability conditions, which determines the accuracy and the well-posedness of the resulting approximation. The method also requires the definition of a high-order discrete divergence operator that is the discrete analog of the divergence operator and is acting on the degrees of freedom. The new family of mimetic methods is proved theoretically to be convergent and optimal error estimates for flux andmore » scalar variable are derived from the convergence analysis. A numerical experiment confirms the high-order accuracy of the method in solving diffusion problems with variable diffusion tensor. It is worth mentioning that the approximation of the scalar variable presents a superconvergence effect.« less

Improved limit to the diffuse flux of ultrahigh energy neutrinos from the Pierre Auger Observatory

NASA Astrophysics Data System (ADS)

Aab, A.; Abreu, P.; Aglietta, M.; Ahn, E. J.; Al Samarai, I.; Albuquerque, I. F. M.; Allekotte, I.; Allison, P.; Almela, A.; Alvarez Castillo, J.; Alvarez-Muñiz, J.; Alves Batista, R.; Ambrosio, M.; Aminaei, A.; Anchordoqui, L.; Andringa, S.; Aramo, C.; Aranda, V. M.; Arqueros, F.; Arsene, N.; Asorey, H.; Assis, P.; Aublin, J.; Ave, M.; Avenier, M.; Avila, G.; Awal, N.; Badescu, A. M.; Barber, K. B.; Bäuml, J.; Baus, C.; Beatty, J. J.; Becker, K. H.; Bellido, J. A.; Berat, C.; Bertaina, M. E.; Bertou, X.; Biermann, P. L.; Billoir, P.; Blaess, S. G.; Blanco, A.; Blanco, M.; Bleve, C.; Blümer, H.; Boháčová, M.; Boncioli, D.; Bonifazi, C.; Borodai, N.; Brack, J.; Brancus, I.; Bridgeman, A.; Brogueira, P.; Brown, W. C.; Buchholz, P.; Bueno, A.; Buitink, S.; Buscemi, M.; Caballero-Mora, K. S.; Caccianiga, B.; Caccianiga, L.; Candusso, M.; Caramete, L.; Caruso, R.; Castellina, A.; Cataldi, G.; Cazon, L.; Cester, R.; Chavez, A. G.; Chiavassa, A.; Chinellato, J. A.; Chudoba, J.; Cilmo, M.; Clay, R. W.; Cocciolo, G.; Colalillo, R.; Coleman, A.; Collica, L.; Coluccia, M. R.; Conceição, R.; Contreras, F.; Cooper, M. J.; Cordier, A.; Coutu, S.; Covault, C. E.; Cronin, J.; Dallier, R.; Daniel, B.; Dasso, S.; Daumiller, K.; Dawson, B. R.; de Almeida, R. M.; de Jong, S. J.; De Mauro, G.; de Mello Neto, J. R. T.; De Mitri, I.; de Oliveira, J.; de Souza, V.; del Peral, L.; Deligny, O.; Dembinski, H.; Dhital, N.; Di Giulio, C.; Di Matteo, A.; Diaz, J. C.; Díaz Castro, M. L.; Diogo, F.; Dobrigkeit, C.; Docters, W.; D'Olivo, J. C.; Dorofeev, A.; Dorosti Hasankiadeh, Q.; Dova, M. T.; Ebr, J.; Engel, R.; Erdmann, M.; Erfani, M.; Escobar, C. O.; Espadanal, J.; Etchegoyen, A.; Falcke, H.; Fang, K.; Farrar, G.; Fauth, A. C.; Fazzini, N.; Ferguson, A. P.; Fernandes, M.; Fick, B.; Figueira, J. M.; Filevich, A.; Filipčič, A.; Fox, B. D.; Fratu, O.; Freire, M. M.; Fuchs, B.; Fujii, T.; García, B.; Garcia-Pinto, D.; Gate, F.; Gemmeke, H.; Gherghel-Lascu, A.; Ghia, P. L.; Giaccari, U.; Giammarchi, M.; Giller, M.; Głas, D.; Glaser, C.; Glass, H.; Golup, G.; Gómez Berisso, M.; Gómez Vitale, P. F.; González, N.; Gookin, B.; Gordon, J.; Gorgi, A.; Gorham, P.; Gouffon, P.; Griffith, N.; Grillo, A. F.; Grubb, T. D.; Guardincerri, Y.; Guarino, F.; Guedes, G. P.; Hampel, M. R.; Hansen, P.; Harari, D.; Harrison, T. A.; Hartmann, S.; Harton, J. L.; Haungs, A.; Hebbeker, T.; Heck, D.; Heimann, P.; Herve, A. E.; Hill, G. C.; Hojvat, C.; Hollon, N.; Holt, E.; Homola, P.; Hörandel, J. R.; Horvath, P.; Hrabovský, M.; Huber, D.; Huege, T.; Insolia, A.; Isar, P. G.; Jandt, I.; Jansen, S.; Jarne, C.; Johnsen, J. A.; Josebachuili, M.; Kääpä, A.; Kambeitz, O.; Kampert, K. H.; Kasper, P.; Katkov, I.; Kégl, B.; Keilhauer, B.; Keivani, A.; Kemp, E.; Kieckhafer, R. M.; Klages, H. O.; Kleifges, M.; Kleinfeller, J.; Krause, R.; Krohm, N.; Krömer, O.; Kuempel, D.; Kunka, N.; LaHurd, D.; Latronico, L.; Lauer, R.; Lauscher, M.; Lautridou, P.; Le Coz, S.; Lebrun, D.; Lebrun, P.; Leigui de Oliveira, M. A.; Letessier-Selvon, A.; Lhenry-Yvon, I.; Link, K.; Lopes, L.; López, R.; López Casado, A.; Louedec, K.; Lu, L.; Lucero, A.; Malacari, M.; Maldera, S.; Mallamaci, M.; Maller, J.; Mandat, D.; Mantsch, P.; Mariazzi, A. G.; Marin, V.; Mariş, I. C.; Marsella, G.; Martello, D.; Martin, L.; Martinez, H.; Martínez Bravo, O.; Martraire, D.; Masías Meza, J. J.; Mathes, H. J.; Mathys, S.; Matthews, J.; Matthews, J. A. J.; Matthiae, G.; Maurel, D.; Maurizio, D.; Mayotte, E.; Mazur, P. O.; Medina, C.; Medina-Tanco, G.; Meissner, R.; Mello, V. B. B.; Melo, D.; Menshikov, A.; Messina, S.; Meyhandan, R.; Micheletti, M. I.; Middendorf, L.; Minaya, I. A.; Miramonti, L.; Mitrica, B.; Molina-Bueno, L.; Mollerach, S.; Montanet, F.; Morello, C.; Mostafá, M.; Moura, C. A.; Muller, M. A.; Müller, G.; Müller, S.; Mussa, R.; Navarra, G.; Navarro, J. L.; Navas, S.; Necesal, P.; Nellen, L.; Nelles, A.; Neuser, J.; Nguyen, P. H.; Niculescu-Oglinzanu, M.; Niechciol, M.; Niemietz, L.; Niggemann, T.; Nitz, D.; Nosek, D.; Novotny, V.; Nožka, L.; Ochilo, L.; Oikonomou, F.; Olinto, A.; Pacheco, N.; Pakk Selmi-Dei, D.; Palatka, M.; Pallotta, J.; Papenbreer, P.; Parente, G.; Parra, A.; Paul, T.; Pech, M.; PÈ©kala, J.; Pelayo, R.; Pepe, I. M.; Perrone, L.; Petermann, E.; Peters, C.; Petrera, S.; Petrov, Y.; Phuntsok, J.; Piegaia, R.; Pierog, T.; Pieroni, P.; Pimenta, M.; Pirronello, V.; Platino, M.; Plum, M.; Porcelli, A.; Porowski, C.; Prado, R. R.; Privitera, P.; Prouza, M.; Purrello, V.; Quel, E. J.; Querchfeld, S.; Quinn, S.; Rautenberg, J.; Ravel, O.; Ravignani, D.; Revenu, B.; Ridky, J.; Riggi, S.; Risse, M.; Ristori, P.; Rizi, V.; Rodrigues de Carvalho, W.; Rodriguez Fernandez, G.; Rodriguez Rojo, J.; Rodríguez-Frías, M. D.; Rogozin, D.; Rosado, J.; Roth, M.; Roulet, E.; Rovero, A. C.; Saffi, S. J.; Saftoiu, A.; Salamida, F.; Salazar, H.; Saleh, A.; Salesa Greus, F.; Salina, G.; Sánchez, F.; Sanchez-Lucas, P.; Santos, E.; Santos, E. M.; Sarazin, F.; Sarkar, B.; Sarmento, R.; Sato, R.; Scarso, C.; Schauer, M.; Scherini, V.; Schieler, H.; Schiffer, P.; Schmidt, D.; Scholten, O.; Schoorlemmer, H.; Schovánek, P.; Schröder, F. G.; Schulz, A.; Schulz, J.; Schumacher, J.; Sciutto, S. J.; Segreto, A.; Settimo, M.; Shadkam, A.; Shellard, R. C.; Sidelnik, I.; Sigl, G.; Sima, O.; Śmiałkowski, A.; Šmída, R.; Snow, G. R.; Sommers, P.; Sorokin, J.; Squartini, R.; Srivastava, Y. N.; Stanca, D.; Stanič, S.; Stapleton, J.; Stasielak, J.; Stephan, M.; Stutz, A.; Suarez, F.; Suomijärvi, T.; Supanitsky, A. D.; Sutherland, M. S.; Swain, J.; Szadkowski, Z.; Taborda, O. A.; Tapia, A.; Tepe, A.; Theodoro, V. M.; Tiffenberg, J.; Timmermans, C.; Todero Peixoto, C. J.; Toma, G.; Tomankova, L.; Tomé, B.; Tonachini, A.; Torralba Elipe, G.; Torres Machado, D.; Travnicek, P.; Ulrich, R.; Unger, M.; Urban, M.; Valdés Galicia, J. F.; Valiño, I.; Valore, L.; van Aar, G.; van Bodegom, P.; van den Berg, A. M.; van Velzen, S.; van Vliet, A.; Varela, E.; Vargas Cárdenas, B.; Varner, G.; Vasquez, R.; Vázquez, J. R.; Vázquez, R. A.; Veberič, D.; Verzi, V.; Vicha, J.; Videla, M.; Villaseñor, L.; Vlcek, B.; Vorobiov, S.; Wahlberg, H.; Wainberg, O.; Walz, D.; Watson, A. A.; Weber, M.; Weidenhaupt, K.; Weindl, A.; Werner, F.; Widom, A.; Wiencke, L.; Wilczyński, H.; Winchen, T.; Wittkowski, D.; Wundheiler, B.; Wykes, S.; Yang, L.; Yapici, T.; Yushkov, A.; Zas, E.; Zavrtanik, D.; Zavrtanik, M.; Zepeda, A.; Zhu, Y.; Zimmermann, B.; Ziolkowski, M.; Zuccarello, F.; Pierre Auger Collaboration

2015-05-01

Neutrinos in the cosmic ray flux with energies near 1 EeV and above are detectable with the Surface Detector array (SD) of the Pierre Auger Observatory. We report here on searches through Auger data from 1 January 2004 until 20 June 2013. No neutrino candidates were found, yielding a limit to the diffuse flux of ultrahigh energy neutrinos that challenges the Waxman-Bahcall bound predictions. Neutrino identification is attempted using the broad time structure of the signals expected in the SD stations, and is efficiently done for neutrinos of all flavors interacting in the atmosphere at large zenith angles, as well as for "Earth-skimming" neutrino interactions in the case of tau neutrinos. In this paper the searches for downward-going neutrinos in the zenith angle bins 60°-75° and 75°-90° as well as for upward-going neutrinos, are combined to give a single limit. The 90% C.L. single-flavor limit to the diffuse flux of ultrahigh energy neutrinos with an E-2 spectrum in the energy range 1.0 ×1 017 eV - 2.5 ×1 019 eV is Eν2d Nν/d Eν<6.4 ×10-9 GeV cm-2 s-1 sr-1 .

Bioturbation, advection, and diffusion of a conserved tracer in a laboratory flume

NASA Astrophysics Data System (ADS)

Work, P. A.; Moore, P. R.; Reible, D. D.

2002-06-01

Laboratory experiments indicating the relative influences of advection, diffusion, and bioturbation on transport of NaCl tracer between a stream and streambed are described. Data were collected in a recirculating flume housing a box filled with test sediments. Peclet numbers ranged from 0 to 1.5. Sediment components included a medium sand (d50 = 0.31 mm), kaolinite, and topsoil. Lumbriculus variegatus were introduced as bioturbators. Conductivity probes were employed to document the flux of the tracer solution out of the bed. Measurements are compared to one-dimensional effective diffusion models assuming one or two horizontal sediment layers. These simple models provide a good indication of tracer half-life in the bed if a suitable effective diffusion coefficient is chosen but underpredict initial flux and overpredict flux at long times. Organism activity was limited to the upper reaches of the sediment test box but eventually exerts a secondary influence on flux from deeper regions.

Volatile Emissions from Hot Spring Basin, Yellowstone National Park, USA

NASA Astrophysics Data System (ADS)

Werner, C.; Hurwitz, S.; Bergfeld, D.; Evans, W. C.; Lowenstern, J. B.; Jaworowski, C.; Heasler, H.

2007-12-01

The flux and composition of magmatic volatiles were characterized for Hot Spring Basin (HSB), Yellowstone National Park, in August 2006. Diffuse fluxes of CO2 (228 sites) from thermal soil were elevated, with a population distribution similar to that of other acid-sulfate areas in Yellowstone. Thus the estimated diffuse emission rate at HSB is proportionately larger than other areas due to its large area, and could be as high as 1000 td-1 CO2. The diffuse flux of H2S was only above detection limits at 20 of the 31 sites measured. The estimated diffuse H2S emission rate was ~ 4 td-1. Good correlation exists between the log of CO2 flux and shallow soil temperatures, indicating linked steam and gas upflow in the subsurface. The correlation between CO2 and H2S fluxes is weak, and the CO2 / H2S diffuse flux ratio was higher than in fumarolic ratios of CO2 to H2S. This suggests that various reactions, e.g., native sulfur deposition, act to remove H2S from the original gas stream in the diffuse low- temperature environment. Dissolved sulfate flux through Shallow Creek, which drains part of HSB, was ~ 4 td-1. Comparing dissolved sulfate flux to estimates of primary emission of H2S based on fumarolic gas geochemistry gives first order estimates of the sulfur consumed in surficial or subsurface mineral deposition. Total C and S outputs from HSB are comparable to other active volcanic systems.

A third-order computational method for numerical fluxes to guarantee nonnegative difference coefficients for advection-diffusion equations in a semi-conservative form

NASA Astrophysics Data System (ADS)

Sakai, K.; Watabe, D.; Minamidani, T.; Zhang, G. S.

2012-10-01

According to Godunov theorem for numerical calculations of advection equations, there exist no higher-order schemes with constant positive difference coefficients in a family of polynomial schemes with an accuracy exceeding the first-order. We propose a third-order computational scheme for numerical fluxes to guarantee the non-negative difference coefficients of resulting finite difference equations for advection-diffusion equations in a semi-conservative form, in which there exist two kinds of numerical fluxes at a cell surface and these two fluxes are not always coincident in non-uniform velocity fields. The present scheme is optimized so as to minimize truncation errors for the numerical fluxes while fulfilling the positivity condition of the difference coefficients which are variable depending on the local Courant number and diffusion number. The feature of the present optimized scheme consists in keeping the third-order accuracy anywhere without any numerical flux limiter. We extend the present method into multi-dimensional equations. Numerical experiments for advection-diffusion equations showed nonoscillatory solutions.

A model of recovering the parameters of fast nonlocal heat transport in magnetic fusion plasmas

NASA Astrophysics Data System (ADS)

Kukushkin, A. B.; Kulichenko, A. A.; Sdvizhenskii, P. A.; Sokolov, A. V.; Voloshinov, V. V.

2017-12-01

A model is elaborated for interpreting the initial stage of the fast nonlocal transport events, which exhibit immediate response, in the diffusion time scale, of the spatial profile of electron temperature to its local perturbation, while the net heat flux is directed opposite to ordinary diffusion (i.e. along the temperature gradient). We solve the inverse problem of recovering the kernel of the integral equation, which describes nonlocal (superdiffusive) transport of energy due to emission and absorption of electromagnetic (EM) waves with long free path and strong reflection from the vacuum vessel’s wall. To allow for the errors of experimental data, we use the method based on the regularized (in the framework of an ill-posed problem, using the parametric models) approximation of available experimental data. The model is applied to interpreting the data from stellarator LHD and tokamak TFTR. The EM wave transport is considered here in the single-group approximation, however the limitations of the physics model enable us to identify the spectral range of the EM waves which might be responsible for the observed phenomenon.

Revisiting the Fundamental Analytical Solutions of Heat and Mass Transfer: The Kernel of Multirate and Multidimensional Diffusion

NASA Astrophysics Data System (ADS)

Zhou, Quanlin; Oldenburg, Curtis M.; Rutqvist, Jonny; Birkholzer, Jens T.

2017-11-01

There are two types of analytical solutions of temperature/concentration in and heat/mass transfer through boundaries of regularly shaped 1-D, 2-D, and 3-D blocks. These infinite-series solutions with either error functions or exponentials exhibit highly irregular but complementary convergence at different dimensionless times, td. In this paper, approximate solutions were developed by combining the error-function-series solutions for early times and the exponential-series solutions for late times and by using time partitioning at the switchover time, td0. The combined solutions contain either the leading term of both series for normal-accuracy approximations (with less than 0.003 relative error) or the first two terms for high-accuracy approximations (with less than 10-7 relative error) for 1-D isotropic (spheres, cylinders, slabs) and 2-D/3-D rectangular blocks (squares, cubes, rectangles, and rectangular parallelepipeds). This rapid and uniform convergence for rectangular blocks was achieved by employing the same time partitioning with individual dimensionless times for different directions and the product of their combined 1-D slab solutions. The switchover dimensionless time was determined to minimize the maximum approximation errors. Furthermore, the analytical solutions of first-order heat/mass flux for 2-D/3-D rectangular blocks were derived for normal-accuracy approximations. These flux equations contain the early-time solution with a three-term polynomial in √td and the late-time solution with the limited-term exponentials for rectangular blocks. The heat/mass flux equations and the combined temperature/concentration solutions form the ultimate kernel for fast simulations of multirate and multidimensional heat/mass transfer in porous/fractured media with millions of low-permeability blocks of varying shapes and sizes.

Revisiting the Fundamental Analytical Solutions of Heat and Mass Transfer: The Kernel of Multirate and Multidimensional Diffusion

DOE Office of Scientific and Technical Information (OSTI.GOV)

Zhou, Quanlin; Oldenburg, Curtis M.; Rutqvist, Jonny

There are two types of analytical solutions of temperature/concentration in and heat/mass transfer through boundaries of regularly shaped 1D, 2D, and 3D blocks. These infinite-series solutions with either error functions or exponentials exhibit highly irregular but complementary convergence at different dimensionless times, t d0. In this paper, approximate solutions were developed by combining the error-function-series solutions for early times and the exponential-series solutions for late times and by using time partitioning at the switchover time, t d0. The combined solutions contain either the leading term of both series for normal-accuracy approximations (with less than 0.003 relative error) or the firstmore » two terms for high-accuracy approximations (with less than 10-7 relative error) for 1D isotropic (spheres, cylinders, slabs) and 2D/3D rectangular blocks (squares, cubes, rectangles, and rectangular parallelepipeds). This rapid and uniform convergence for rectangular blocks was achieved by employing the same time partitioning with individual dimensionless times for different directions and the product of their combined 1D slab solutions. The switchover dimensionless time was determined to minimize the maximum approximation errors. Furthermore, the analytical solutions of first-order heat/mass flux for 2D/3D rectangular blocks were derived for normal-accuracy approximations. These flux equations contain the early-time solution with a three-term polynomial in √td and the late-time solution with the limited-term exponentials for rectangular blocks. The heat/mass flux equations and the combined temperature/concentration solutions form the ultimate kernel for fast simulations of multirate and multidimensional heat/mass transfer in porous/fractured media with millions of low-permeability blocks of varying shapes and sizes.« less

Revisiting the Fundamental Analytical Solutions of Heat and Mass Transfer: The Kernel of Multirate and Multidimensional Diffusion

DOE PAGES

Zhou, Quanlin; Oldenburg, Curtis M.; Rutqvist, Jonny; ...

2017-10-24

There are two types of analytical solutions of temperature/concentration in and heat/mass transfer through boundaries of regularly shaped 1D, 2D, and 3D blocks. These infinite-series solutions with either error functions or exponentials exhibit highly irregular but complementary convergence at different dimensionless times, t d0. In this paper, approximate solutions were developed by combining the error-function-series solutions for early times and the exponential-series solutions for late times and by using time partitioning at the switchover time, t d0. The combined solutions contain either the leading term of both series for normal-accuracy approximations (with less than 0.003 relative error) or the firstmore » two terms for high-accuracy approximations (with less than 10-7 relative error) for 1D isotropic (spheres, cylinders, slabs) and 2D/3D rectangular blocks (squares, cubes, rectangles, and rectangular parallelepipeds). This rapid and uniform convergence for rectangular blocks was achieved by employing the same time partitioning with individual dimensionless times for different directions and the product of their combined 1D slab solutions. The switchover dimensionless time was determined to minimize the maximum approximation errors. Furthermore, the analytical solutions of first-order heat/mass flux for 2D/3D rectangular blocks were derived for normal-accuracy approximations. These flux equations contain the early-time solution with a three-term polynomial in √td and the late-time solution with the limited-term exponentials for rectangular blocks. The heat/mass flux equations and the combined temperature/concentration solutions form the ultimate kernel for fast simulations of multirate and multidimensional heat/mass transfer in porous/fractured media with millions of low-permeability blocks of varying shapes and sizes.« less

A screening tool for delineating subregions of steady recharge within groundwater models

USGS Publications Warehouse

Dickinson, Jesse; Ferré, T.P.A.; Bakker, Mark; Crompton, Becky

2014-01-01

We have developed a screening method for simplifying groundwater models by delineating areas within the domain that can be represented using steady-state groundwater recharge. The screening method is based on an analytical solution for the damping of sinusoidal infiltration variations in homogeneous soils in the vadose zone. The damping depth is defined as the depth at which the flux variation damps to 5% of the variation at the land surface. Groundwater recharge may be considered steady where the damping depth is above the depth of the water table. The analytical solution approximates the vadose zone diffusivity as constant, and we evaluated when this approximation is reasonable. We evaluated the analytical solution through comparison of the damping depth computed by the analytic solution with the damping depth simulated by a numerical model that allows variable diffusivity. This comparison showed that the screening method conservatively identifies areas of steady recharge and is more accurate when water content and diffusivity are nearly constant. Nomograms of the damping factor (the ratio of the flux amplitude at any depth to the amplitude at the land surface) and the damping depth were constructed for clay and sand for periodic variations between 1 and 365 d and flux means and amplitudes from nearly 0 to 1 × 10−3 m d−1. We applied the screening tool to Central Valley, California, to identify areas of steady recharge. A MATLAB script was developed to compute the damping factor for any soil and any sinusoidal flux variation.

Search for diffuse neutrino flux from astrophysical sources with MACRO

NASA Astrophysics Data System (ADS)

MACRO Collaboration; Ambrosio, M.; Antolini, R.; Auriemma, G.; Bakari, D.; Baldini, A.; Barbarino, G. C.; Barish, B. C.; Battistoni, G.; Becherini, Y.; Bellotti, R.; Bemporad, C.; Bernardini, P.; Bilokon, H.; Bloise, C.; Bower, C.; Brigida, M.; Bussino, S.; Cafagna, F.; Calicchio, M.; Campana, D.; Carboni, M.; Caruso, R.; Cecchini, S.; Cei, F.; Chiarella, V.; Choudhary, B. C.; Coutu, S.; Cozzi, M.; de Cataldo, G.; Dekhissi, H.; de Marzo, C.; de Mitri, I.; Derkaoui, J.; de Vincenzi, M.; di Credico, A.; Erriquez, O.; Favuzzi, C.; Forti, C.; Fusco, P.; Giacomelli, G.; Giannini, G.; Giglietto, N.; Giorgini, M.; Grassi, M.; Grillo, A.; Guarino, F.; Gustavino, C.; Habig, A.; Hanson, K.; Heinz, R.; Iarocci, E.; Katsavounidis, E.; Katsavounidis, I.; Kearns, E.; Kim, H.; Kyriazopoulou, S.; Lamanna, E.; Lane, C.; Levin, D. S.; Lipari, P.; Longley, N. P.; Longo, M. J.; Loparco, F.; Maaroufi, F.; Mancarella, G.; Mandrioli, G.; Margiotta, A.; Marini, A.; Martello, D.; Marzari-Chiesa, A.; Mazziotta, M. N.; Michael, D. G.; Monacelli, P.; Montaruli, T.; Monteno, M.; Mufson, S.; Musser, J.; Nicolò, D.; Nolty, R.; Orth, C.; Osteria, G.; Palamara, O.; Patera, V.; Patrizii, L.; Pazzi, R.; Peck, C. W.; Perrone, L.; Petrera, S.; Pistilli, P.; Popa, V.; Rainò, A.; Reynoldson, J.; Ronga, F.; Rrhioua, A.; Satriano, C.; Scapparone, E.; Scholberg, K.; Sciubba, A.; Serra, P.; Sioli, M.; Sirri, G.; Sitta, M.; Spinelli, P.; Spinetti, M.; Spurio, M.; Steinberg, R.; Stone, J. L.; Sulak, L. R.; Surdo, A.; Tarlè, G.; Togo, V.; Vakili, M.; Walter, C. W.; Webb, R.

2003-04-01

Many galactic and extragalactic astrophysical sources are currently considered promising candidates as high-energy neutrino emitters. Astrophysical neutrinos can be detected as upward-going muons produced in charged-current interactions with the medium surrounding the detector. The expected neutrino fluxes from various models start to dominate on the atmospheric neutrino background at neutrino energies above some tens of TeV. We present the results of a search for an excess of high-energy upward-going muons among the sample of data collected by MACRO during ~5.8 years of effective running time. No significant evidence for this signal was found. As a consequence, an upper limit on the flux of upward-going muons from high-energy neutrinos was set at the level of 1.7×10-14 cm-2s-1sr-1. The corresponding upper limit for the diffuse neutrino flux was evaluated assuming a neutrino power law spectrum. Our result was compared with theoretical predictions and upper limits from other experiments.

Problems with heterogeneous and non-isotropic media or distorted grids

DOE Office of Scientific and Technical Information (OSTI.GOV)

Hyman, J.; Shashkov, M.; Steinberg, S.

1996-08-01

This paper defines discretizations of the divergence and flux operators that produce symmetric, positive-definite, and accurate approximations to steady-state diffusion problems. Because discontinuous material properties and highly distorted grids are allowed, the flux operator, rather than the gradient, is used as a fundamental operator to be discretized. Resulting finite-difference scheme is similar to those obtained from the mixed finite-element method.

Improved limit to the diffuse flux of ultrahigh energy neutrinos from the Pierre Auger Observatory

DOE PAGES

Aab, Alexander

2015-05-26

Neutrinos in the cosmic ray flux with energies near 1 EeV and above are detectable with the Surface Detector array (SD) of the Pierre Auger Observatory. We report here on searches through Auger data from 1 January 2004 until 20 June 2013. No neutrino candidates were found, yielding a limit to the diffuse flux of ultrahigh energy neutrinos that challenges the Waxman-Bahcall bound predictions. Neutrino identification is attempted using the broad time structure of the signals expected in the SD stations, and is efficiently done for neutrinos of all flavors interacting in the atmosphere at large zenith angles, as wellmore » as for “Earth-skimming” neutrino interactions in the case of tau neutrinos. In this paper the searches for downward-going neutrinos in the zenith angle bins 60°–75° and 75°–90° as well as for upward-going neutrinos, are combined to give a single limit. In addition, the 90% C.L. single-flavor limit to the diffuse flux of ultrahigh energy neutrinos with an E –2 spectrum in the energy range 1.0 × 10 17 eV – 2.5 × 10 19 eV is E 2 νdN ν/dE ν < 6.4 × 10 –9 GeV cm –2 s –1 sr –1.« less

Quasi-analytical treatment of spatially averaged radiation transfer in complex terrain

NASA Astrophysics Data System (ADS)

LöWe, H.; Helbig, N.

2012-10-01

We provide a new quasi-analytical method to compute the subgrid topographic influences on the shortwave radiation fluxes and the effective albedo in complex terrain as required for large-scale meteorological, land surface, or climate models. We investigate radiative transfer in complex terrain via the radiosity equation on isotropic Gaussian random fields. Under controlled approximations we derive expressions for domain-averaged fluxes of direct, diffuse, and terrain radiation and the sky view factor. Domain-averaged quantities can be related to a type of level-crossing probability of the random field, which is approximated by long-standing results developed for acoustic scattering at ocean boundaries. This allows us to express all nonlocal horizon effects in terms of a local terrain parameter, namely, the mean-square slope. Emerging integrals are computed numerically, and fit formulas are given for practical purposes. As an implication of our approach, we provide an expression for the effective albedo of complex terrain in terms of the Sun elevation angle, mean-square slope, the area-averaged surface albedo, and the ratio of atmospheric direct beam to diffuse radiation. For demonstration we compute the decrease of the effective albedo relative to the area-averaged albedo in Switzerland for idealized snow-covered and clear-sky conditions at noon in winter. We find an average decrease of 5.8% and spatial patterns which originate from characteristics of the underlying relief. Limitations and possible generalizations of the method are discussed.

Resonant thickening of self-gravitating discs: imposed or self-induced orbital diffusion in the tightly wound limit

NASA Astrophysics Data System (ADS)

Fouvry, Jean-Baptiste; Pichon, Christophe; Chavanis, Pierre-Henri; Monk, Laura

2017-11-01

The secular thickening of a self-gravitating stellar galactic disc is investigated using the dressed collisionless Fokker-Planck equation and the inhomogeneous multicomponent Balescu-Lenard equation. The thick WKB limits for the diffusion fluxes are found using the epicyclic approximation, while assuming that only radially tightly wound transient spirals are sustained by the disc. This yields simple quadratures for the drift and diffusion coefficients, providing a clear understanding of the positions of maximum vertical orbital diffusion within the disc, induced by fluctuations either external or due to the finite number of particles. These thick limits also offer a consistent derivation of a thick disc Toomre parameter, which is shown to be exponentially boosted by the ratio of the vertical to radial scaleheights. Dressed potential fluctuations within the disc statistically induce a vertical bending of a subset of resonant orbits, triggering the corresponding increase in vertical velocity dispersion. When applied to a tepid stable tapered disc perturbed by shot noise, these two frameworks reproduce qualitatively the formation of ridges of resonant orbits towards larger vertical actions, as found in direct numerical simulations, but overestimates the time-scale involved in their appearance. Swing amplification is likely needed to resolve this discrepancy, as demonstrated in the case of razor-thin discs. Other sources of thickening are also investigated, such as fading sequences of slowing bars, or the joint evolution of a population of giant molecular clouds within the disc.

Determination of the reaction rate coefficient of sulphide mine tailings deposited under water.

PubMed

Awoh, Akué Sylvette; Mbonimpa, Mamert; Bussière, Bruno

2013-10-15

The efficiency of a water cover to limit dissolved oxygen (DO) availability to underlying acid-generating mine tailings can be assessed by calculating the DO flux at the tailings-water interface. Fick's equations, which are generally used to calculate this flux, require knowing the effective DO diffusion coefficient (Dw) and the reaction (consumption) rate coefficient (Kr) of the tailings, or the DO concentration profile. Whereas Dw can be accurately estimated, few studies have measured the parameter Kr for submerged sulphide tailings. The objective of this study was to determine Kr for underwater sulphide tailings in a laboratory experiment. Samples of sulphide mine tailings (an approximately 6 cm layer) were placed in a cell under a water cover (approximately 2 cm) maintained at constant DO concentration. Two tailings were studied: TA1 with high sulphide content (83% pyrite) and TA2 with low sulphide content (2.8% pyrite). DO concentration was measured with a microelectrode at various depths above and below the tailings-water interface at 1 mm intervals. Results indicate that steady-state condition was rapidly attained. As expected, a diffusive boundary layer (DBL) was observed in all cases. An iterative back-calculation process using the numerical code POLLUTEv6 and taking the DBL into account provided the Kr values used to match calculated and experimental concentration profiles. Kr obtained for tailings TA1 and TA2 was about 80 d(-1) and 6.5 d(-1), respectively. For comparison purposes, Kr obtained from cell tests on tailings TA1 was lower than Kr calculated from the sulphate production rate obtained from shake-flask tests. Steady-state DO flux at the water-tailings interface was then calculated with POLLUTEv6 using tailings characteristics Dw and Kr. For the tested conditions, DO flux ranged from 608 to 758 mg O2/m(2)/d for tailings TA1 and from 177 to 221 mg O2/m(2)/d for tailings TA2. The impact of placing a protective layer of inert material over the tailings was also investigated for tailings TA1 (with high sulphide content). A protective layer of only 5 cm reduced the DO flux into the tailings at about 5 mg/m(2)/d, compared to 608 mg O2/m(2)/d without a protective layer, or an approximately 99% reduction in flux. Copyright © 2013 Elsevier Ltd. All rights reserved.

Border-crossing model for the diffusive coarsening of two-dimensional and quasi-two-dimensional wet foams

NASA Astrophysics Data System (ADS)

Schimming, C. D.; Durian, D. J.

2017-09-01

For dry foams, the transport of gas from small high-pressure bubbles to large low-pressure bubbles is dominated by diffusion across the thin soap films separating neighboring bubbles. For wetter foams, the film areas become smaller as the Plateau borders and vertices inflate with liquid. So-called "border-blocking" models can explain some features of wet-foam coarsening based on the presumption that the inflated borders totally block the gas flux; however, this approximation dramatically fails in the wet or unjamming limit where the bubbles become close-packed spheres and coarsening proceeds even though there are no films. Here, we account for the ever-present border-crossing flux by a new length scale defined by the average gradient of gas concentration inside the borders. We compute that it is proportional to the geometric average of film and border thicknesses, and we verify this scaling by numerical solution of the diffusion equation. We similarly consider transport across inflated vertices and surface Plateau borders in quasi-two-dimensional foams. And we show how the d A /d t =K0(n -6 ) von Neumann law is modified by the appearance of terms that depend on bubble size and shape as well as the concentration gradient length scales. Finally, we use the modified von Neumann law to compute the growth rate of the average bubble area, which is not constant.

Coupling compositional liquid gas Darcy and free gas flows at porous and free-flow domains interface

DOE Office of Scientific and Technical Information (OSTI.GOV)

Masson, R., E-mail: roland.masson@unice.fr; Team COFFEE INRIA Sophia Antipolis Méditerranée; Trenty, L., E-mail: laurent.trenty@andra.fr

This paper proposes an efficient splitting algorithm to solve coupled liquid gas Darcy and free gas flows at the interface between a porous medium and a free-flow domain. This model is compared to the reduced model introduced in [6] using a 1D approximation of the gas free flow. For that purpose, the gas molar fraction diffusive flux at the interface in the free-flow domain is approximated by a two point flux approximation based on a low-frequency diagonal approximation of a Steklov–Poincaré type operator. The splitting algorithm and the reduced model are applied in particular to the modelling of the massmore » exchanges at the interface between the storage and the ventilation galleries in radioactive waste deposits.« less

Limits on diffuse fluxes of high energy extraterrestrial neutrinos with the AMANDA-B10 detector

DOE Office of Scientific and Technical Information (OSTI.GOV)

Ahrens, J.; Bai, X.; Barwick, S.W.

2003-03-11

Data from the AMANDA-B10 detector taken during the austral winter of 1997 have been searched for a diffuse flux of high energy extraterrestrial muon-neutrinos, as predicted from, e.g., the sum of all active galaxies in the universe. This search yielded no excess events above those expected from the background atmospheric neutrinos, leading to upper limits on the extraterrestrial neutrino flux. For an assumed E{sup -2} spectrum, a 90 percent classical confidence level upper limit has been placed at a level E{sup 2} Phi(E) = 8.4 x 10{sup -7} GeV cm{sup -2} s{sup -1}1 sr{sup -1} (for a predominant neutrino energymore » range 6-1000 TeV) which is the most restrictive bound placed by any neutrino detector. When specific predicted spectral forms are considered, it is found that some are excluded.« less

Effects of stomata clustering on leaf gas exchange.

PubMed

Lehmann, Peter; Or, Dani

2015-09-01

A general theoretical framework for quantifying the stomatal clustering effects on leaf gaseous diffusive conductance was developed and tested. The theory accounts for stomatal spacing and interactions among 'gaseous concentration shells'. The theory was tested using the unique measurements of Dow et al. (2014) that have shown lower leaf diffusive conductance for a genotype of Arabidopsis thaliana with clustered stomata relative to uniformly distributed stomata of similar size and density. The model accounts for gaseous diffusion: through stomatal pores; via concentration shells forming at pore apertures that vary with stomata spacing and are thus altered by clustering; and across the adjacent air boundary layer. Analytical approximations were derived and validated using a numerical model for 3D diffusion equation. Stomata clustering increases the interactions among concentration shells resulting in larger diffusive resistance that may reduce fluxes by 5-15%. A similar reduction in conductance was found for clusters formed by networks of veins. The study resolves ambiguities found in the literature concerning stomata end-corrections and stomatal shape, and provides a new stomata density threshold for diffusive interactions of overlapping vapor shells. The predicted reduction in gaseous exchange due to clustering, suggests that guard cell function is impaired, limiting stomatal aperture opening. © 2015 The Authors. New Phytologist © 2015 New Phytologist Trust.

You Don't Need Richards'... A New General 1-D Vadose Zone Solution Method that is Reliable

NASA Astrophysics Data System (ADS)

Ogden, F. L.; Lai, W.; Zhu, J.; Steinke, R. C.; Talbot, C. A.

2015-12-01

Hydrologic modelers and mathematicians have strived to improve 1-D Richards' equation (RE) solution reliability for predicting vadose zone fluxes. Despite advances in computing power and the numerical solution of partial differential equations since Richards first published the RE in 1931, the solution remains unreliable. That is to say that there is no guarantee that for a particular set of soil constitutive relations, moisture profile conditions, or forcing input that a numerical RE solver will converge to an answer. This risk of non-convergence renders prohibitive the use of RE solvers in hydrological models that need perhaps millions of infiltration solutions. In lieu of using unreliable numerical RE solutions, researchers have developed a wide array of approximate solutions that more-or-less mimic the behavior of the RE, with some notable deficiencies such as parameter insensitivity or divergence over time. The improved Talbot-Ogden (T-O) finite water-content scheme was shown by Ogden et al. (2015) to be an extremely good approximation of the 1-D RE solution, with a difference in cumulative infiltration of only 0.2 percent over an 8 month simulation comparing the improved T-O scheme with a RE numerical solver. The reason is that the newly-derived fundamental flow equation that underpins the improved T-O method is equivalent to the RE minus a term that is equal to the diffusive flux divided by the slope of the wetting front. Because the diffusive flux has zero mean, this term is not important in calculating the mean flux. The wetting front slope is near infinite (sharp) in coarser soils that produce more significant hydrological interactions between surface and ground waters, which also makes this missing term 1) disappear in the limit, and, 2) create stability challenges for the numerical solution of RE. The improved T-O method is a replacement for the 1-D RE in soils that can be simulated as homogeneous layers, where the user is willing to neglect the effects of soil water diffusivity. This presentation emphasizes the transformative nature of the improved T-O finite water-content solution, and highlights the benefits of the methods' reliability in high-resolution large watershed simulations in the high performance computing environment, and discusses coupling of the soil matrix and non-Darcian macropores.

Mass transport at rotating disk electrodes: effects of synthetic particles and nerve endings.

PubMed

Chiu, Veronica M; Lukus, Peter A; Doyle, Jamie L; Schenk, James O

2011-11-01

An unstirred layer (USL) exists at the interface of solids with solutions. Thus, the particles in brain tissue preparations possess a USL as well as at the surface of a rotating disk electrode (RDE) used to measure chemical fluxes. Time constraints for observing biological kinetics based on estimated thicknesses of USLs at the membrane surface in real samples of nerve endings were estimated. Liposomes, silica, and Sephadex were used separately to model the tissue preparation particles. Within a solution stirred by the RDE, both diffusion and hydrodynamic boundary layers are formed. It was observed that the number and size of particles decreased the following: the apparent diffusion coefficient excluding Sephadex, boundary layer thicknesses excluding silica, sensitivity excluding diluted liposomes (in agreement with results from other laboratories), limiting current potentially due to an increase in the path distance, and mixing time. They have no effect on the detection limit (6 ± 2 nM). The RDE kinetically resolves transmembrane transport with a timing of approximately 30 ms. Copyright © 2011 Elsevier Inc. All rights reserved.

Understanding of flux-limited behaviors of heat transport in nonlinear regime

NASA Astrophysics Data System (ADS)

Guo, Yangyu; Jou, David; Wang, Moran

2016-01-01

The classical Fourier's law of heat transport breaks down in highly nonequilibrium situations as in nanoscale heat transport, where nonlinear effects become important. The present work is aimed at exploring the flux-limited behaviors based on a categorization of existing nonlinear heat transport models in terms of their theoretical foundations. Different saturation heat fluxes are obtained, whereas the same qualitative variation trend of heat flux versus exerted temperature gradient is got in diverse nonlinear models. The phonon hydrodynamic model is proposed to act as a standard to evaluate other heat flux limiters because of its more rigorous physical foundation. A deeper knowledge is thus achieved about the phenomenological generalized heat transport models. The present work provides deeper understanding and accurate modeling of nonlocal and nonlinear heat transport beyond the diffusive limit.

Determination of pollutant diffusion coefficients in naturally formed biofilms using a single tube extractive membrane bioreactor

DOE Office of Scientific and Technical Information (OSTI.GOV)

Zhang, S.F.; Splendiani, A.; Freitas dos Santos, L.M.

A novel technique has been used to determine the effective diffusion coefficients for 1,1,2-trichloroethane (TCE), a nonreacting tracer, in biofilms growing on the external surface of a silicone rubber membrane tube during degradation of 1,2-dichloroethane (DCE) by Xanthobacter autotrophicus GJ10 and monochlorobenzene (MCB) by Pseudomonas JS150. Experiments were carried out in a single tube extractive membrane bioreactor (STEMB), whose configuration makes it possible to measure the transmembrane flux of substrates. A video imaging technique (VIT) was employed for in situ biofilm thickness measurement and recording. Diffusion coefficients of TCE in the biofilms and TCE mass transfer coefficients in the liquidmore » films adjacent to the biofilms were determined simultaneously using a resistances-in-series diffusion model. It was found that the flux and overall mass transfer coefficient of TCE decrease with increasing biofilm thickness, showing the importance of biofilm diffusion on the mass transfer process. Similar fluxes were observed for the nonreacting tracer (TCE) and the reactive substrates (MCB or DCE), suggesting that membrane-attached biofilm systems can be rate controlled primarily by substrate diffusion. The TCE diffusion coefficient in the JS150 biofilm appeared to be dependent on biofilm thickness, decreasing markedly for biofilm thicknesses of >1 mm. The values of the TCE diffusion coefficients in the JS150 biofilms <1-mm thick are approximately twice those in water and fall to around 30% of the water value for biofilms >1-mm thick.« less

New constraints on all flavor Galactic diffuse neutrino emission with the ANTARES telescope

NASA Astrophysics Data System (ADS)

Albert, A.; André, M.; Anghinolfi, M.; Anton, G.; Ardid, M.; Aubert, J.-J.; Avgitas, T.; Baret, B.; Barrios-Martí, J.; Basa, S.; Belhorma, B.; Bertin, V.; Biagi, S.; Bormuth, R.; Bourret, S.; Bouwhuis, M. C.; Bruijn, R.; Brunner, J.; Busto, J.; Capone, A.; Caramete, L.; Carr, J.; Celli, S.; Cherkaoui El Moursli, R.; Chiarusi, T.; Circella, M.; Coelho, J. A. B.; Coleiro, A.; Coniglione, R.; Costantini, H.; Coyle, P.; Creusot, A.; Díaz, A. F.; Deschamps, A.; de Bonis, G.; Distefano, C.; di Palma, I.; Domi, A.; Donzaud, C.; Dornic, D.; Drouhin, D.; Eberl, T.; El Bojaddaini, I.; El Khayati, N.; Elsässer, D.; Enzenhöfer, A.; Ettahiri, A.; Fassi, F.; Felis, I.; Fusco, L. A.; Galatà, S.; Gay, P.; Giordano, V.; Glotin, H.; Grégoire, T.; Gracia Ruiz, R.; Graf, K.; Hallmann, S.; van Haren, H.; Heijboer, A. J.; Hello, Y.; Hernández-Rey, J. J.; Hößl, J.; Hofestädt, J.; Hugon, C.; Illuminati, G.; James, C. W.; de Jong, M.; Jongen, M.; Kadler, M.; Kalekin, O.; Katz, U.; Kießling, D.; Kouchner, A.; Kreter, M.; Kreykenbohm, I.; Kulikovskiy, V.; Lachaud, C.; Lahmann, R.; Lefèvre, D.; Leonora, E.; Lotze, M.; Loucatos, S.; Marcelin, M.; Margiotta, A.; Marinelli, A.; Martínez-Mora, J. A.; Mele, R.; Melis, K.; Michael, T.; Migliozzi, P.; Moussa, A.; Navas, S.; Nezri, E.; Organokov, M.; Pǎvǎlaş, G. E.; Pellegrino, C.; Perrina, C.; Piattelli, P.; Popa, V.; Pradier, T.; Quinn, L.; Racca, C.; Riccobene, G.; Sánchez-Losa, A.; Saldaña, M.; Salvadori, I.; Samtleben, D. F. E.; Sanguineti, M.; Sapienza, P.; Schüssler, F.; Sieger, C.; Spurio, M.; Stolarczyk, Th.; Taiuti, M.; Tayalati, Y.; Trovato, A.; Turpin, D.; Tönnis, C.; Vallage, B.; van Elewyck, V.; Versari, F.; Vivolo, D.; Vizzoca, A.; Wilms, J.; Zornoza, J. D.; Zúñiga, J.; Gaggero, D.; Grasso, D.; ANTARES Collaboration

2017-09-01

The flux of very high-energy neutrinos produced in our Galaxy by the interaction of accelerated cosmic rays with the interstellar medium is not yet determined. The characterization of this flux will shed light on Galactic accelerator features, gas distribution morphology and Galactic cosmic ray transport. The central Galactic plane can be the site of an enhanced neutrino production, thus leading to anisotropies in the extraterrestrial neutrino signal as measured by the IceCube Collaboration. The ANTARES neutrino telescope, located in the Mediterranean Sea, offers a favorable view of this part of the sky, thereby allowing for a contribution to the determination of this flux. The expected diffuse Galactic neutrino emission can be obtained, linking a model of generation and propagation of cosmic rays with the morphology of the gas distribution in the Milky Way. In this paper, the so-called "gamma model" introduced recently to explain the high-energy gamma-ray diffuse Galactic emission is assumed as reference. The neutrino flux predicted by the "gamma model" depends on the assumed primary cosmic ray spectrum cutoff. Considering a radially dependent diffusion coefficient, this proposed scenario is able to account for the local cosmic ray measurements, as well as for the Galactic gamma-ray observations. Nine years of ANTARES data are used in this work to search for a possible Galactic contribution according to this scenario. All flavor neutrino interactions are considered. No excess of events is observed, and an upper limit is set on the neutrino flux of 1.1 (1.2) times the prediction of the "gamma model," assuming the primary cosmic ray spectrum cutoff at 5 (50) PeV. This limit excludes the diffuse Galactic neutrino emission as the major cause of the "spectral anomaly" between the two hemispheres measured by IceCube.

The narrow pulse approximation and long length scale determination in xenon gas diffusion NMR studies of model porous media

NASA Technical Reports Server (NTRS)

Mair, R. W.; Sen, P. N.; Hurlimann, M. D.; Patz, S.; Cory, D. G.; Walsworth, R. L.

2002-01-01

We report a systematic study of xenon gas diffusion NMR in simple model porous media, random packs of mono-sized glass beads, and focus on three specific areas peculiar to gas-phase diffusion. These topics are: (i) diffusion of spins on the order of the pore dimensions during the application of the diffusion encoding gradient pulses in a PGSE experiment (breakdown of the narrow pulse approximation and imperfect background gradient cancellation), (ii) the ability to derive long length scale structural information, and (iii) effects of finite sample size. We find that the time-dependent diffusion coefficient, D(t), of the imbibed xenon gas at short diffusion times in small beads is significantly affected by the gas pressure. In particular, as expected, we find smaller deviations between measured D(t) and theoretical predictions as the gas pressure is increased, resulting from reduced diffusion during the application of the gradient pulse. The deviations are then completely removed when water D(t) is observed in the same samples. The use of gas also allows us to probe D(t) over a wide range of length scales and observe the long time asymptotic limit which is proportional to the inverse tortuosity of the sample, as well as the diffusion distance where this limit takes effect (approximately 1-1.5 bead diameters). The Pade approximation can be used as a reference for expected xenon D(t) data between the short and the long time limits, allowing us to explore deviations from the expected behavior at intermediate times as a result of finite sample size effects. Finally, the application of the Pade interpolation between the long and the short time asymptotic limits yields a fitted length scale (the Pade length), which is found to be approximately 0.13b for all bead packs, where b is the bead diameter. c. 2002 Elsevier Sciences (USA).

The narrow pulse approximation and long length scale determination in xenon gas diffusion NMR studies of model porous media.

PubMed

Mair, R W; Sen, P N; Hürlimann, M D; Patz, S; Cory, D G; Walsworth, R L

2002-06-01

We report a systematic study of xenon gas diffusion NMR in simple model porous media, random packs of mono-sized glass beads, and focus on three specific areas peculiar to gas-phase diffusion. These topics are: (i) diffusion of spins on the order of the pore dimensions during the application of the diffusion encoding gradient pulses in a PGSE experiment (breakdown of the narrow pulse approximation and imperfect background gradient cancellation), (ii) the ability to derive long length scale structural information, and (iii) effects of finite sample size. We find that the time-dependent diffusion coefficient, D(t), of the imbibed xenon gas at short diffusion times in small beads is significantly affected by the gas pressure. In particular, as expected, we find smaller deviations between measured D(t) and theoretical predictions as the gas pressure is increased, resulting from reduced diffusion during the application of the gradient pulse. The deviations are then completely removed when water D(t) is observed in the same samples. The use of gas also allows us to probe D(t) over a wide range of length scales and observe the long time asymptotic limit which is proportional to the inverse tortuosity of the sample, as well as the diffusion distance where this limit takes effect (approximately 1-1.5 bead diameters). The Padé approximation can be used as a reference for expected xenon D(t) data between the short and the long time limits, allowing us to explore deviations from the expected behavior at intermediate times as a result of finite sample size effects. Finally, the application of the Padé interpolation between the long and the short time asymptotic limits yields a fitted length scale (the Padé length), which is found to be approximately 0.13b for all bead packs, where b is the bead diameter. c. 2002 Elsevier Sciences (USA).

Hybrid simplified spherical harmonics with diffusion equation for light propagation in tissues.

PubMed

Chen, Xueli; Sun, Fangfang; Yang, Defu; Ren, Shenghan; Zhang, Qian; Liang, Jimin

2015-08-21

Aiming at the limitations of the simplified spherical harmonics approximation (SPN) and diffusion equation (DE) in describing the light propagation in tissues, a hybrid simplified spherical harmonics with diffusion equation (HSDE) based diffuse light transport model is proposed. In the HSDE model, the living body is first segmented into several major organs, and then the organs are divided into high scattering tissues and other tissues. DE and SPN are employed to describe the light propagation in these two kinds of tissues respectively, which are finally coupled using the established boundary coupling condition. The HSDE model makes full use of the advantages of SPN and DE, and abandons their disadvantages, so that it can provide a perfect balance between accuracy and computation time. Using the finite element method, the HSDE is solved for light flux density map on body surface. The accuracy and efficiency of the HSDE are validated with both regular geometries and digital mouse model based simulations. Corresponding results reveal that a comparable accuracy and much less computation time are achieved compared with the SPN model as well as a much better accuracy compared with the DE one.

Vaporization of a solid surface in an ambient gas

NASA Astrophysics Data System (ADS)

Benilov, M. S.; Jacobsson, S.; Kaddani, A.; Zahrai, S.

2001-07-01

The net flux of vapour from a solid surface in an ambient gas is analysed with the aim to estimate the effect of vaporization cooling on the energy balance of an arc cathode under conditions typical for a high-power current breaker. If the ratio of the equilibrium vapour pressure pv to the ambient pressure p∞ is smaller than unity, the removal of vapour from the surface is due to diffusion into the bulk of the gas. As a consequence, the net flux of the vapour from the surface is much smaller than the emitted flux. An estimate of the diffusion rate under conditions typical for a high-power current breaker indicates that vaporization cooling plays a minor role in the energy balance of the cathode in this case. If ratio pv/p∞ is above unity, the flow of the vapour from the surface appears and the net flux is comparable to the emitted flux. A simple analytical solution has been obtained for this case, which is in a good agreement with results of the Monte Carlo modelling of preceding authors. If pv/p∞ exceeds approximately 4.5, vaporization occurs as into vacuum and the net flux is about 0.82 of the emitted flux.

Wyman's equation and oxygen flux through the red cell.

PubMed

McCabe, Michael; Maguire, David J

2007-01-01

Wyman's equation of 1966 describes the facilitation of flux of a reversibly bound substrate such as oxygen, consequent on the translational diffusion of the binding protein (the carrier). While Wyman's equation, or some modification of it such as that by Murray 2, may provide a realistic description of the flux of oxygen through a dilute solution of haemoglobin (see also Wittenburg), it is unlikely to be the complete explanation, nor even the basis, for oxygen transport through the intact red cell. The mature erythrocyte contains approximately 350 g/l haemoglobin, and while this suggests that only 35% of the available water volume is actually occupied by the protein, the remaining 65% is unavailable for protein translational diffusion due to the mutual exclusion of the haemoglobin molecules. For this reason we have examined other possible mechanisms whereby haemoglobin may facilitate the translational diffusion of oxygen within the erythrocyte. Possible alternatives include rotational diffusion by the haemoglobins, intracellular shuffling of haemoglobins due to shape changes by the erythrocyte, and haemoglobin rotations and oxygen exchange consequent on the charge change which accompanies substration and desubstration of the haemoglobin molecule. Finally the dipole interactions are shown to generate significant intermolecular attractions between adjacent haemoglobins.

Multi-year search for a diffuse flxu of muon neutrinos with AMANDA-II

DOE Office of Scientific and Technical Information (OSTI.GOV)

IceCube Collaboration; Klein, Spencer; Achterberg, A.

2008-04-13

A search for TeV-PeV muon neutrinos from unresolved sources was performed on AMANDA-II data collected between 2000 and 2003 with an equivalent livetime of 807 days. This diffuse analysis sought to find an extraterrestrial neutrino flux from sources with non-thermal components. The signal is expected to have a harder spectrum than the atmospheric muon and neutrino backgrounds. Since no excess of events was seen in the data over the expected background, an upper limit of E{sup 2}{Phi}{sub 90%C.L.} < 7.4 x 10{sup -8} GeV cm{sup -2} s{sup -1} sr{sup -1} is placed on the diffuse flux of muon neutrinos withmore » a {Phi} {proportional_to} E{sup -2} spectrum in the energy range 16 TeV to 2.5 PeV. This is currently the most sensitive {Phi} {proportional_to} E{sup -2} diffuse astrophysical neutrino limit. We also set upper limits for astrophysical and prompt neutrino models, all of which have spectra different than {Phi} {proportional_to} E{sup -2}.« less

A Flux-Corrected Transport Based Hydrodynamic Model for the Plasmasphere Refilling Problem following Geomagnetic Storms

NASA Astrophysics Data System (ADS)

Chatterjee, K.; Schunk, R. W.

2017-12-01

The refilling of the plasmasphere following a geomagnetic storm remains one of the longstanding problems in the area of ionosphere-magnetosphere coupling. Both diffusion and hydrodynamic approximations have been adopted for the modeling and solution of this problem. The diffusion approximation neglects the nonlinear inertial term in the momentum equation and so this approximation is not rigorously valid immediately after the storm. Over the last few years, we have developed a hydrodynamic refilling model using the flux-corrected transport method, a numerical method that is extremely well suited to handling nonlinear problems with shocks and discontinuities. The plasma transport equations are solved along 1D closed magnetic field lines that connect conjugate ionospheres and the model currently includes three ion (H+, O+, He+) and two neutral (O, H) species. In this work, each ion species under consideration has been modeled as two separate streams emanating from the conjugate hemispheres and the model correctly predicts supersonic ion speeds and the presence of high levels of Helium during the early hours of refilling. The ultimate objective of this research is the development of a 3D model for the plasmasphere refilling problem and with additional development, the same methodology can potentially be applied to the study of other complex space plasma coupling problems in closed flux tube geometries. Index Terms: 2447 Modeling and forecasting [IONOSPHERE] 2753 Numerical modeling [MAGNETOSPHERIC PHYSICS] 7959 Models [SPACE WEATHER

The OPGT/MJS plasma wave science team

NASA Technical Reports Server (NTRS)

Scarf, F. L.

1972-01-01

Some properties of a model magnetosphere for Saturn were studied in order to determine the bounds that can be set on surface field strength and trapped particle population. The primary observational constraint was that nonthermal radiation similar to the Jovian radio emissions must be undetectable from Earth. It is argued that for a Saturn surface field of approximately one gauss, particles that are energized as they diffuse in from the magnetopause with conservation of magnetic moment will produce synchrotron radiation levels that are undetectable at a range of 9.5 AU. The plasma instabilities that heat the oncoming wind particles at the bow shock and others that can limit the stably-trapped flux levels are also discussed.

Comparison of fluid neutral models for one-dimensional plasma edge modeling with a finite volume solution of the Boltzmann equation

DOE Office of Scientific and Technical Information (OSTI.GOV)

Horsten, N., E-mail: niels.horsten@kuleuven.be; Baelmans, M.; Dekeyser, W.

2016-01-15

We derive fluid neutral approximations for a simplified 1D edge plasma model, suitable to study the neutral behavior close to the target of a nuclear fusion divertor, and compare its solutions to the solution of the corresponding kinetic Boltzmann equation. The plasma is considered as a fixed background extracted from a detached 2D simulation. We show that the Maxwellian equilibrium distribution is already obtained very close to the target, justifying the use of a fluid approximation. We compare three fluid neutral models: (i) a diffusion model; (ii) a pressure-diffusion model (i.e., a combination of a continuity and momentum equation) assumingmore » equal neutral and ion temperatures; and (iii) the pressure-diffusion model coupled to a neutral energy equation taking into account temperature differences between neutrals and ions. Partial reflection of neutrals reaching the boundaries is included in both the kinetic and fluid models. We propose two methods to obtain an incident neutral flux boundary condition for the fluid models: one based on a diffusion approximation and the other assuming a truncated Chapman-Enskog distribution. The pressure-diffusion model predicts the plasma sources very well. The diffusion boundary condition gives slightly better results overall. Although including an energy equation still improves the results, the assumption of equal ion and neutral temperature already gives a very good approximation.« less

Comparison of transport properties models for numerical simulations of Mars entry vehicles

NASA Astrophysics Data System (ADS)

Hao, Jiaao; Wang, Jingying; Gao, Zhenxun; Jiang, Chongwen; Lee, Chunhian

2017-01-01

Effects of two different models for transport properties, including the approximate model and the collision integral model, on hypersonic flow simulations of Mars entry vehicles are numerically investigated. A least square fitting is firstly performed using the best-available data of collision integrals for Martian atmosphere species within the temperature range of 300-20,000 K. Then, the performance of these two transport properties models are compared for an equilibrium Martian atmosphere gas mixture at 10 kPa and temperatures ranging from 1000 to 10,000 K. Finally, four flight conditions chosen from the trajectory of the Mars Pathfinder entry vehicle are numerically simulated. It is indicated that the approximate model is capable of accurately providing the distributions of species mass fractions and temperatures in the flowfield. Both models give similar translational-rotational and vibrational heat fluxes. However, the chemical diffusion heat fluxes predicted by the approximate model are significantly larger than the results computed by the collision integral model, particularly in the vicinity of the forebody stagnation point, whose maximum relative error of 15% for the super-catalytic case. The diffusion model employed in the approximate model is responsible to the discrepancy. In addition, the wake structure is largely unaffected by the transport properties models.

Search for neutrino-induced cascades with five years of AMANDA data

NASA Astrophysics Data System (ADS)

Abbasi, R.; Abdou, Y.; Abu-Zayyad, T.; Actis, O.; Adams, J.; Aguilar, J. A.; Ahlers, M.; Andeen, K.; Auffenberg, J.; Bai, X.; Baker, M.; Barwick, S. W.; Bay, R.; Bazo Alba, J. L.; Beattie, K.; Beatty, J. J.; Bechet, S.; Becker, J. K.; Becker, K.-H.; Benabderrahmane, M. L.; Berdermann, J.; Berghaus, P.; Berley, D.; Bernardini, E.; Bertrand, D.; Besson, D. Z.; Bissok, M.; Blaufuss, E.; Boersma, D. J.; Bohm, C.; Böser, S.; Botner, O.; Bradley, L.; Braun, J.; Buitink, S.; Carson, M.; Chirkin, D.; Christy, B.; Clem, J.; Clevermann, F.; Cohen, S.; Colnard, C.; Cowen, D. F.; D'Agostino, M. V.; Danninger, M.; Davis, J. C.; de Clercq, C.; Demirörs, L.; Depaepe, O.; Descamps, F.; Desiati, P.; de Vries-Uiterweerd, G.; Deyoung, T.; Díaz-Vélez, J. C.; Dreyer, J.; Dumm, J. P.; Duvoort, M. R.; Ehrlich, R.; Eisch, J.; Ellsworth, R. W.; Engdegård, O.; Euler, S.; Evenson, P. A.; Fadiran, O.; Fazely, A. R.; Feusels, T.; Filimonov, K.; Finley, C.; Foerster, M. M.; Fox, B. D.; Franckowiak, A.; Franke, R.; Gaisser, T. K.; Gallagher, J.; Ganugapati, R.; Geisler, M.; Gerhardt, L.; Gladstone, L.; Glüsenkamp, T.; Goldschmidt, A.; Goodman, J. A.; Grant, D.; Griesel, T.; Groß, A.; Grullon, S.; Gunasingha, R. M.; Gurtner, M.; Ha, C.; Hallgren, A.; Halzen, F.; Han, K.; Hanson, K.; Helbing, K.; Herquet, P.; Hickford, S.; Hill, G. C.; Hoffman, K. D.; Homeier, A.; Hoshina, K.; Hubert, D.; Huelsnitz, W.; Hülß, J.-P.; Hulth, P. O.; Hultqvist, K.; Hussain, S.; Imlay, R. L.; Ishihara, A.; Jacobsen, J.; Japaridze, G. S.; Johansson, H.; Joseph, J. M.; Kampert, K.-H.; Kappes, A.; Karg, T.; Karle, A.; Kelley, J. L.; Kemming, N.; Kenny, P.; Kiryluk, J.; Kislat, F.; Klein, S. R.; Knops, S.; Köhne, J.-H.; Kohnen, G.; Kolanoski, H.; Köpke, L.; Koskinen, D. J.; Kowalski, M.; Kowarik, T.; Krasberg, M.; Krings, T.; Kroll, G.; Kuehn, K.; Kuwabara, T.; Labare, M.; Lafebre, S.; Laihem, K.; Landsman, H.; Lauer, R.; Lehmann, R.; Lennarz, D.; Lünemann, J.; Madsen, J.; Majumdar, P.; Maruyama, R.; Mase, K.; Matis, H. S.; Matusik, M.; Meagher, K.; Merck, M.; Mészáros, P.; Meures, T.; Middell, E.; Milke, N.; Miller, J.; Montaruli, T.; Morse, R.; Movit, S. M.; Nahnhauer, R.; Nam, J. W.; Naumann, U.; Nießen, P.; Nygren, D. R.; Odrowski, S.; Olivas, A.; Olivo, M.; Ono, M.; Panknin, S.; Paul, L.; Pérez de Los Heros, C.; Petrovic, J.; Piegsa, A.; Pieloth, D.; Porrata, R.; Posselt, J.; Price, P. B.; Prikockis, M.; Przybylski, G. T.; Rawlins, K.; Redl, P.; Resconi, E.; Rhode, W.; Ribordy, M.; Rizzo, A.; Rodrigues, J. P.; Roth, P.; Rothmaier, F.; Rott, C.; Roucelle, C.; Ruhe, T.; Rutledge, D.; Ruzybayev, B.; Ryckbosch, D.; Sander, H.-G.; Sarkar, S.; Schatto, K.; Schlenstedt, S.; Schmidt, T.; Schneider, D.; Schukraft, A.; Schultes, A.; Schulz, O.; Schunck, M.; Seckel, D.; Semburg, B.; Seo, S. H.; Sestayo, Y.; Seunarine, S.; Silvestri, A.; Slipak, A.; Spiczak, G. M.; Spiering, C.; Stamatikos, M.; Stanev, T.; Stephens, G.; Stezelberger, T.; Stokstad, R. G.; Stoyanov, S.; Strahler, E. A.; Straszheim, T.; Sullivan, G. W.; Swillens, Q.; Taboada, I.; Tamburro, A.; Tepe, A.; Ter-Antonyan, S.; Tilav, S.; Toale, P. A.; Tosi, D.; Turčan, D.; van Eijndhoven, N.; Vandenbroucke, J.; van Overloop, A.; van Santen, J.; Voigt, B.; Walck, C.; Waldenmaier, T.; Wallraff, M.; Walter, M.; Wendt, C.; Westerhoff, S.; Whitehorn, N.; Wiebe, K.; Wiebusch, C. H.; Wikström, G.; Williams, D. R.; Wischnewski, R.; Wissing, H.; Woschnagg, K.; Xu, C.; Xu, X. W.; Yodh, G.; Yoshida, S.; Zarzhitsky, P.

2011-01-01

We report on the search for electromagnetic and hadronic showers (“cascades”) produced by a diffuse flux of extraterrestrial neutrinos in the AMANDA neutrino telescope. Data for this analysis were recorded during 1001 days of detector livetime in the years 2000-2004. The observed event rates are consistent with the background expectation from atmospheric neutrinos and muons. An upper limit is derived for the diffuse flux of neutrinos of all flavors assuming a flavor ratio of νe:νμ:ντ = 1:1:1 at the detection site. The all-flavor flux of neutrinos with an energy spectrum Φ ∝ E-2 is less than 5.0 × 10-7 GeV s-1 sr-1 cm-2 at a 90% C.L. Here, 90% of the simulated signal would fall within the energy range 40 TeV to 9 PeV. We discuss flux limits in the context of several specific models of extraterrestrial and prompt atmospheric neutrino production.

Studies of the extreme ultraviolet/soft x-ray background

DOE Office of Scientific and Technical Information (OSTI.GOV)

Stern, R.A.

1978-01-01

The results of an extensive sky survey of the extreme ultraviolet (EUV)/soft x-ray background are reported. The data were obtained with a focusing telescope designed and calibrated at U.C. Berkeley which observed EUV sources and the diffuse background as part of the Apollo-Soyuz mission in July, 1975. With a primary field-of-view of 2.3 + 0.1/sup 0/ FWHM and four EUV bandpass filters (16 to 25, 20 to 73, 80 to 108, and 80 to 250 eV) the EUV telescope obtained background data included in the final observational sample for 21 discrete sky locations and 11 large angular scans, as wellmore » as for a number of shorter observations. Analysis of the data reveals as intense flux above 80 eV energy, with upper limits to the background intensity given for the lower energy filters Ca 2 x 10/sup 4/ and 6 x 10/sup 2/ ph cm/sup -2/ sec/sup -1/ ster/sup -1/ eV/sup -1/ at 21 and 45 eV respectively). The 80 to 108 eV flux agrees within statistical errors with the earlier results of Cash, Malina and Stern (1976): the Apollo-Soyuz average reported intensity is 4.0 +- 1.3 ph cm/sup -2/ sec/sup -1/ ster/sup -1/ eV/sup -1/ at Ca 100 eV, or roughly a factor of ten higher than the corresponding 250 eV intensity. The uniformity of the background flux is uncertain due to limitations in the statistical accuracy of the data; upper limits to the point-to-point standard deviation of the background intensity are (..delta..I/I approximately less than 0.8 +- 0.4 (80 to 108 eV) and approximately less than 0.4 +- 0.2 (80 to 250 eV). No evidence is found for a correlation between the telescope count rate and earth-based parameters (zenith angle, sun angle, etc.) for E approximately greater than 80 eV (the lower energy bandpasses are significantly affected by scattered solar radiation. Unlike some previous claims for the soft x-ray background, no simple dependence upon galactic latitude is seen.« less

Radiation pressure driving of a dusty atmosphere

NASA Astrophysics Data System (ADS)

Tsang, Benny T.-H.; Milosavljević, Miloš

2015-10-01

Radiation pressure can be dynamically important in star-forming environments such as ultra-luminous infrared and submillimetre galaxies. Whether and how radiation drives turbulence and bulk outflows in star formation sites is still unclear. The uncertainty in part reflects the limitations of direct numerical schemes that are currently used to simulate radiation transfer and radiation-gas coupling. An idealized setup in which radiation is introduced at the base of a dusty atmosphere in a gravitational field has recently become the standard test for radiation-hydrodynamics methods in the context of star formation. To a series of treatments featuring the flux-limited diffusion approximation as well as a short-characteristics tracing and M1 closure for the variable Eddington tensor approximation, we here add another treatment that is based on the implicit Monte Carlo radiation transfer scheme. Consistent with all previous treatments, the atmosphere undergoes Rayleigh-Taylor instability and readjusts to a near-Eddington-limited state. We detect late-time net acceleration in which the turbulent velocity dispersion matches that reported previously with the short-characteristics-based radiation transport closure, the most accurate of the three preceding treatments. Our technical result demonstrates the importance of accurate radiation transfer in simulations of radiative feedback.

Cuticular gas exchange by Antarctic sea spiders.

PubMed

Lane, Steven J; Moran, Amy L; Shishido, Caitlin M; Tobalske, Bret W; Woods, H Arthur

2018-04-25

Many marine organisms and life stages lack specialized respiratory structures, like gills, and rely instead on cutaneous respiration, which they facilitate by having thin integuments. This respiratory mode may limit body size, especially if the integument also functions in support or locomotion. Pycnogonids, or sea spiders, are marine arthropods that lack gills and rely on cutaneous respiration but still grow to large sizes. Their cuticle contains pores, which may play a role in gas exchange. Here, we examined alternative paths of gas exchange in sea spiders: (1) oxygen diffuses across pores in the cuticle, a common mechanism in terrestrial eggshells, (2) oxygen diffuses directly across the cuticle, a common mechanism in small aquatic insects, or (3) oxygen diffuses across both pores and cuticle. We examined these possibilities by modeling diffusive oxygen fluxes across all pores in the body of sea spiders and asking whether those fluxes differed from measured metabolic rates. We estimated fluxes across pores using Fick's law parameterized with measurements of pore morphology and oxygen gradients. Modeled oxygen fluxes through pores closely matched oxygen consumption across a range of body sizes, which means the pores facilitate oxygen diffusion. Furthermore, pore volume scaled hypermetrically with body size, which helps larger species facilitate greater diffusive oxygen fluxes across their cuticle. This likely presents a functional trade-off between gas exchange and structural support, in which the cuticle must be thick enough to prevent buckling due to external forces but porous enough to allow sufficient gas exchange. © 2018. Published by The Company of Biologists Ltd.

Bright X-ray arcs and the emergence of solar magnetic flux

NASA Technical Reports Server (NTRS)

Chapman, G. A.; Broussard, R. M.

1977-01-01

The Skylab S-056 and S-082A experiments and ground-based magnetograms have been used to study the role of bright X-ray arcs and the emergence of solar magnetic flux in the McMath region 12476. The S-056 X-ray images show a system of one or sometimes two bright arcs within a diffuse emitting region. The arcs seem to directly connect regions of opposite magnetic polarity in the photosphere. Magnetograms suggest the possible emergence of a magnetic flux. The width of the main arc is approximately 6 arcsec when most clearly defined, and the length is approximately 30-50 arcsec. Although the arc system is observed to vary in brightness over a period exceeding 24 hours, it remains fixed in orientation. The temperature of the main arc is approximately 3 x 10 to the 6th K. It is suggested that merging magnetic fields may provide the primary energy source, perhaps accompanied by resistive heating from a force-free current.

Molecules in motion: influences of diffusion on metabolic structure and function in skeletal muscle.

PubMed

Kinsey, Stephen T; Locke, Bruce R; Dillaman, Richard M

2011-01-15

Metabolic processes are often represented as a group of metabolites that interact through enzymatic reactions, thus forming a network of linked biochemical pathways. Implicit in this view is that diffusion of metabolites to and from enzymes is very fast compared with reaction rates, and metabolic fluxes are therefore almost exclusively dictated by catalytic properties. However, diffusion may exert greater control over the rates of reactions through: (1) an increase in reaction rates; (2) an increase in diffusion distances; or (3) a decrease in the relevant diffusion coefficients. It is therefore not surprising that skeletal muscle fibers have long been the focus of reaction-diffusion analyses because they have high and variable rates of ATP turnover, long diffusion distances, and hindered metabolite diffusion due to an abundance of intracellular barriers. Examination of the diversity of skeletal muscle fiber designs found in animals provides insights into the role that diffusion plays in governing both rates of metabolic fluxes and cellular organization. Experimental measurements of metabolic fluxes, diffusion distances and diffusion coefficients, coupled with reaction-diffusion mathematical models in a range of muscle types has started to reveal some general principles guiding muscle structure and metabolic function. Foremost among these is that metabolic processes in muscles do, in fact, appear to be largely reaction controlled and are not greatly limited by diffusion. However, the influence of diffusion is apparent in patterns of fiber growth and metabolic organization that appear to result from selective pressure to maintain reaction control of metabolism in muscle.

Computational and analytical comparison of flux discretizations for the semiconductor device equations beyond Boltzmann statistics

NASA Astrophysics Data System (ADS)

Farrell, Patricio; Koprucki, Thomas; Fuhrmann, Jürgen

2017-10-01

We compare three thermodynamically consistent numerical fluxes known in the literature, appearing in a Voronoï finite volume discretization of the van Roosbroeck system with general charge carrier statistics. Our discussion includes an extension of the Scharfetter-Gummel scheme to non-Boltzmann (e.g. Fermi-Dirac) statistics. It is based on the analytical solution of a two-point boundary value problem obtained by projecting the continuous differential equation onto the interval between neighboring collocation points. Hence, it serves as a reference flux. The exact solution of the boundary value problem can be approximated by computationally cheaper fluxes which modify certain physical quantities. One alternative scheme averages the nonlinear diffusion (caused by the non-Boltzmann nature of the problem), another one modifies the effective density of states. To study the differences between these three schemes, we analyze the Taylor expansions, derive an error estimate, visualize the flux error and show how the schemes perform for a carefully designed p-i-n benchmark simulation. We present strong evidence that the flux discretization based on averaging the nonlinear diffusion has an edge over the scheme based on modifying the effective density of states.

The equilibrium-diffusion limit for radiation hydrodynamics

DOE PAGES

Ferguson, J. M.; Morel, J. E.; Lowrie, R.

2017-07-27

The equilibrium-diffusion approximation (EDA) is used to describe certain radiation-hydrodynamic (RH) environments. When this is done the RH equations reduce to a simplified set of equations. The EDA can be derived by asymptotically analyzing the full set of RH equations in the equilibrium-diffusion limit. Here, we derive the EDA this way and show that it and the associated set of simplified equations are both first-order accurate with transport corrections occurring at second order. Having established the EDA’s first-order accuracy we then analyze the grey nonequilibrium-diffusion approximation and the grey Eddington approximation and show that they both preserve this first-order accuracy.more » Further, these approximations preserve the EDA’s first-order accuracy when made in either the comoving-frame (CMF) or the lab-frame (LF). And while analyzing the Eddington approximation, we found that the CMF and LF radiation-source equations are equivalent when neglecting O(β 2) terms and compared in the LF. Of course, the radiation pressures are not equivalent. It is expected that simplified physical models and numerical discretizations of the RH equations that do not preserve this first-order accuracy will not retain the correct equilibrium-diffusion solutions. As a practical example, we show that nonequilibrium-diffusion radiative-shock solutions devolve to equilibrium-diffusion solutions when the asymptotic parameter is small.« less

Multi-Component Diffusion with Application To Computational Aerothermodynamics

NASA Technical Reports Server (NTRS)

Sutton, Kenneth; Gnoffo, Peter A.

1998-01-01

The accuracy and complexity of solving multicomponent gaseous diffusion using the detailed multicomponent equations, the Stefan-Maxwell equations, and two commonly used approximate equations have been examined in a two part study. Part I examined the equations in a basic study with specified inputs in which the results are applicable for many applications. Part II addressed the application of the equations in the Langley Aerothermodynamic Upwind Relaxation Algorithm (LAURA) computational code for high-speed entries in Earth's atmosphere. The results showed that the presented iterative scheme for solving the Stefan-Maxwell equations is an accurate and effective method as compared with solutions of the detailed equations. In general, good accuracy with the approximate equations cannot be guaranteed for a species or all species in a multi-component mixture. 'Corrected' forms of the approximate equations that ensured the diffusion mass fluxes sum to zero, as required, were more accurate than the uncorrected forms. Good accuracy, as compared with the Stefan- Maxwell results, were obtained with the 'corrected' approximate equations in defining the heating rates for the three Earth entries considered in Part II.

A Novel Diffuse Fraction-Based Two-Leaf Light Use Efficiency Model: An Application Quantifying Photosynthetic Seasonality across 20 AmeriFlux Flux Tower Sites

NASA Astrophysics Data System (ADS)

Yan, Hao; Wang, Shao-Qiang; Yu, Kai-Liang; Wang, Bin; Yu, Qin; Bohrer, Gil; Billesbach, Dave; Bracho, Rosvel; Rahman, Faiz; Shugart, Herman H.

2017-10-01

Diffuse radiation can increase canopy light use efficiency (LUE). This creates the need to differentiate the effects of direct and diffuse radiation when simulating terrestrial gross primary production (GPP). Here, we present a novel GPP model, the diffuse-fraction-based two-leaf model (DTEC), which includes the leaf response to direct and diffuse radiation, and treats maximum LUE for shaded leaves (ɛmsh defined as a power function of the diffuse fraction (Df)) and sunlit leaves (ɛmsu defined as a constant) separately. An Amazonian rainforest site (KM67) was used to calibrate the model by simulating the linear relationship between monthly canopy LUE and Df. This showed a positive response of forest GPP to atmospheric diffuse radiation, and suggested that diffuse radiation was more limiting than global radiation and water availability for Amazon rainforest GPP on a monthly scale. Further evaluation at 20 independent AmeriFlux sites showed that the DTEC model, when driven by monthly meteorological data and MODIS leaf area index (LAI) products, explained 70% of the variability observed in monthly flux tower GPP. This exceeded the 51% accounted for by the MODIS 17A2 big-leaf GPP product. The DTEC model's explicit accounting for the impacts of diffuse radiation and soil water stress along with its parameterization for C4 and C3 plants was responsible for this difference. The evaluation of DTEC at Amazon rainforest sites demonstrated its potential to capture the unique seasonality of higher GPP during the diffuse radiation-dominated wet season. Our results highlight the importance of diffuse radiation in seasonal GPP simulation.