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.

Heat transfer in suspensions of rigid particles

NASA Astrophysics Data System (ADS)

Brandt, Luca; Niazi Ardekani, Mehdi; Abouali, Omid

2016-11-01

We study the heat transfer in laminar Couette flow of suspensions of rigid neutrally buoyant particles by means of numerical simulations. An Immersed Boundary Method is coupled with a VOF approach to simulate the heat transfer in the fluid and solid phase, enabling us to fully resolve the heat diffusion. First, we consider spherical particles and show that the proposed algorithm is able to reproduce the correlations between heat flux across the channel, the particle volume fraction and the heat diffusivity obtained in laboratory experiments and recently proposed in the literature, results valid in the limit of vanishing inertia. We then investigate the role of inertia on the heat transfer and show an increase of the suspension diffusivity at finite particle Reynolds numbers. Finally, we vary the relativity diffusivity of the fluid and solid phase and investigate its effect on the effective heat flux across the channel. The data are analyzed by considering the ensemble averaged energy equation and decomposing the heat flux in 4 different contributions, related to diffusion in the solid and fluid phase, and the correlations between wall-normal velocity and temperature fluctuations. Results for non-spherical particles will be examined before the meeting. Supported by the European Research Council Grant No. ERC-2013- CoG-616186, TRITOS. The authors acknowledge computer time provided by SNIC (Swedish National Infrastructure for Computing).

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.

Horizontal convection with mechanical stirring

NASA Astrophysics Data System (ADS)

Griffiths, Ross; Stewart, Kial; Hughes, Graham

2012-11-01

The effects of turbulent mixing on convective circulation forced by a horizontal gradient of buoyancy at the surface is examined using laboratory experiments in which a salt flux is introduced at the surface, at one end of a box, and a freshwater buoyancy condition is applied over the rest of the surface. Horizontal rods are oscillated and yo-yoed continuously through the water column, providing a diffusivity that can be calibrated. The convection reaches a stationary state having zero net salt flux. We find that for small stirring rates the small but finite volume flux from the dense source is significant and a virtual source correction is required to take this into account. The density stratification and overturning volume transport are consistent with a theoretical model for high Rayleigh numbers: the transport ψ increases with diffusivity κ (ψg ~ gκ 1 / 4) . The results show that vertical mixing in the boundary layer is important, particularly in setting the density of the interior and the overturning rate. However, interior mixing is unimportant, which raises an interesting question over whether abyssal mixing rates in the ocean play any significant role in setting the abyssal ocean density or the transport in the Meridional Overturning Circulation.

Plasma density injection and flow during coaxial helicity injection in a tokamak

NASA Astrophysics Data System (ADS)

Hooper, E. B.

2018-02-01

Whole device, resistive MHD simulations of spheromaks and tokamaks have used a large diffusion coefficient that maintains a nearly constant density throughout the device. In the present work, helicity and plasma are coinjected into a low-density plasma in a tokamak with a small diffusion coefficient. As in previous simulations [Hooper et al., Phys. Plasmas 20, 092510 (2013)], a flux bubble is formed, which expands to fill the tokamak volume. The injected plasma is non-uniform inside the bubble. The flow pattern is analyzed; when the simulation is not axisymmetric, an n = 1 mode on the surface of the bubble generates leakage of plasma into the low-density volume. Closed flux is generated following injection, as in experiments and previous simulations. The result provides a more detailed physics analysis of the injection, including density non-uniformities in the plasma that may affect its use as a startup plasma [Raman et al., Phys. Rev. Lett. 97, 175002 (2006)].

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.

DUMAND Summer Workshop, University of California, La Jolla, Calif., July 24-September 2, 1978, Proceedings. Volume 2 - UHE interactions, neutrino astronomy

NASA Technical Reports Server (NTRS)

Roberts, A.

1979-01-01

The volume covers categories on inelastic neutrino scattering and the W-boson, and other ultra-high-energy processes, on pulsars, quasars and galactic nuclei, as well as other point sources and constants from gamma ray astronomy. Individual subjects include weak intermediate vector bosons and DUMAND, the Monte Carlo simulation of inelastic neutrino scattering in DUMAND, and Higgs boson production by very high-energy neutrinos. The observability of the neutrino flux from the inner region of the galactic disk, the diffuse fluxes of high-energy neutrinos, as well as the significance of gamma ray observations for neutrino astronomy are also among the topics covered.

SToRM: A Model for Unsteady Surface Hydraulics Over Complex Terrain

USGS Publications Warehouse

Simoes, Francisco J.

2014-01-01

A two-dimensional (depth-averaged) finite volume Godunov-type shallow water model developed for flow over complex topography is presented. The model is based on an unstructured cellcentered finite volume formulation and a nonlinear strong stability preserving Runge-Kutta time stepping scheme. The numerical discretization is founded on the classical and well established shallow water equations in hyperbolic conservative form, but the convective fluxes are calculated using auto-switching Riemann and diffusive numerical fluxes. The model’s implementation within a graphical user interface is discussed. Field application of the model is illustrated by utilizing it to estimate peak flow discharges in a flooding event of historic significance in Colorado, U.S.A., in 2013.

Estimation of nocturnal CO2 and N2O soil emissions from changes in surface boundary layer mass storage

NASA Astrophysics Data System (ADS)

Grant, Richard H.; Omonode, Rex A.

2018-04-01

Annual budgets of greenhouse and other trace gases require knowledge of the emissions throughout the year. Unfortunately, emissions into the surface boundary layer during stable, calm nocturnal periods are not measurable using most micrometeorological methods due to non-stationarity and uncoupled flow. However, during nocturnal periods with very light winds, carbon dioxide (CO2) and nitrous oxide (N2O) frequently accumulate near the surface and this mass accumulation can be used to determine emissions. Gas concentrations were measured at four heights (one within and three above canopy) and turbulence was measured at three heights above a mature 2.5 m maize canopy from 23 July to 10 September 2015. Nocturnal CO2 and N2O fluxes from the canopy were determined using the accumulation of mass within a 6.3 m control volume and out the top of the control volume within the nocturnal surface boundary layer. Diffusive fluxes were estimated by flux gradient method. The total accumulative and diffusive fluxes during near-calm nights (friction velocities < 0.05 ms-1) averaged 1.16 µmol m-2 s-1 CO2 and 0.53 nmol m-2 s-1 N2O. Fluxes were also measured using chambers. Daily mean CO2 fluxes determined by the accumulation method were 90 to 130 % of those determined using soil chambers. Daily mean N2O fluxes determined by the accumulation method were 60 to 80 % of that determined using soil chambers. The better signal-to-noise ratios of the chamber method for CO2 over N2O, non-stationary flow, assumed Schmidt numbers, and anemometer tilt were likely contributing reasons for the differences in chambers versus accumulated nocturnal mass flux estimates. Near-surface N2O accumulative flux measurements in more homogeneous regions and with greater depth are needed to confirm the conclusion that mass accumulation can be effectively used to estimate soil emissions during nearly calm nights.

Geometry of the diffusive propagation region in the August 14, 1982 solar electron event

NASA Technical Reports Server (NTRS)

Evenson, P. A.

1985-01-01

On August 14, 1982, relativistic electrons arrived promptly after an impulsive gamma ray flare, indicating that very little scattering was taking place in interplanetary space. By ignoring anisotropy data the time profile of the event is well described by interplanetary diffusion except for the derived particle injection time. This discrepancy provides independent evidence that the particles are diffusing in a volume close to the Sun rather than in interplanetary space. The flux at maximum method of determining the number of particles produced is still a good approximation when appropriately applied.

Cell-centered high-order hyperbolic finite volume method for diffusion equation on unstructured grids

NASA Astrophysics Data System (ADS)

Lee, Euntaek; Ahn, Hyung Taek; Luo, Hong

2018-02-01

We apply a hyperbolic cell-centered finite volume method to solve a steady diffusion equation on unstructured meshes. This method, originally proposed by Nishikawa using a node-centered finite volume method, reformulates the elliptic nature of viscous fluxes into a set of augmented equations that makes the entire system hyperbolic. We introduce an efficient and accurate solution strategy for the cell-centered finite volume method. To obtain high-order accuracy for both solution and gradient variables, we use a successive order solution reconstruction: constant, linear, and quadratic (k-exact) reconstruction with an efficient reconstruction stencil, a so-called wrapping stencil. By the virtue of the cell-centered scheme, the source term evaluation was greatly simplified regardless of the solution order. For uniform schemes, we obtain the same order of accuracy, i.e., first, second, and third orders, for both the solution and its gradient variables. For hybrid schemes, recycling the gradient variable information for solution variable reconstruction makes one order of additional accuracy, i.e., second, third, and fourth orders, possible for the solution variable with less computational work than needed for uniform schemes. In general, the hyperbolic method can be an effective solution technique for diffusion problems, but instability is also observed for the discontinuous diffusion coefficient cases, which brings necessity for further investigation about the monotonicity preserving hyperbolic diffusion method.

New approach to effective diffusion coefficient evaluation in the nanostructured two-phase media

NASA Astrophysics Data System (ADS)

Lyashenko, Yu. O.; Liashenko, O. Y.; Morozovich, V. V.

2018-03-01

Most widely used basic and combined models for evaluation of the effective diffusion parameters of inhomogeneous two-phase zone are reviewed. A new combined model of effective medium is analyzed for the description of diffusion processes in the two-phase zones. In this model the effective diffusivity depends on the growth kinetic coefficients of each phase, the volume fractions of phases and on the additional parameter that generally characterizes the structure type of the two-phase zone. Our combined model describes two-phase zone evolution in the binary systems based on consideration of the diffusion fluxes through both phases. The Lattice Monte Carlo method was used to test the validity of different phenomenological models for evaluation of the effective diffusivity in nanostructured two-phase zones with different structural morphology.

Water transport, free volume, and polymer dynamics in crosslinked polymer networks

NASA Astrophysics Data System (ADS)

Frieberg, Bradley; Soles, Christopher

Many technologies rely on amorphous polymer membranes that selectively transport small molecules or ions, which has led to a significant scientific interest in elucidating the mechanisms of transport. A recurring theme among several different materials systems is that free volume and polymer chain dynamics facilitate transport. In order to understand the interplay between free volume, transport and polymer dynamics we quantify these properties for a model epoxy network. The epoxy chemistry allows for systematically varying both the structural rigidity of the network as well as the cross-link density. We performed positron annihilation lifetime spectroscopy measurements to characterize the unoccupied volume and correlated the unoccupied volume to the equilibrium moisture uptake and effective diffusion coefficient. We have recently extended this work to include polymer dynamics measured by quasi-elastic neutron scattering on the NIST High Flux Backscatter Spectrometer. These measurements reveal a strong correlation between the MSD and the transport kinetics, which was even stronger than the correlation previously observed between free volume and water diffusion. These observations challenge previous theories that suggest free volume governs transport.

Quantitative Microplate-Based Respirometry with Correction for Oxygen Diffusion

PubMed Central

2009-01-01

Respirometry using modified cell culture microplates offers an increase in throughput and a decrease in biological material required for each assay. Plate based respirometers are susceptible to a range of diffusion phenomena; as O2 is consumed by the specimen, atmospheric O2 leaks into the measurement volume. Oxygen also dissolves in and diffuses passively through the polystyrene commonly used as a microplate material. Consequently the walls of such respirometer chambers are not just permeable to O2 but also store substantial amounts of gas. O2 flux between the walls and the measurement volume biases the measured oxygen consumption rate depending on the actual [O2] gradient. We describe a compartment model-based correction algorithm to deconvolute the biological oxygen consumption rate from the measured [O2]. We optimize the algorithm to work with the Seahorse XF24 extracellular flux analyzer. The correction algorithm is biologically validated using mouse cortical synaptosomes and liver mitochondria attached to XF24 V7 cell culture microplates, and by comparison to classical Clark electrode oxygraph measurements. The algorithm increases the useful range of oxygen consumption rates, the temporal resolution, and durations of measurements. The algorithm is presented in a general format and is therefore applicable to other respirometer systems. PMID:19555051

Volume Diffusion Growth Kinetics and Step Geometry in Crystal Growth

NASA Technical Reports Server (NTRS)

Mazuruk, Konstantin; Ramachandran, Narayanan

1998-01-01

The role of step geometry in two-dimensional stationary volume diff4sion process used in crystal growth kinetics models is investigated. Three different interface shapes: a) a planar interface, b) an equidistant hemispherical bumps train tAx interface, and c) a train of right angled steps, are used in this comparative study. The ratio of the super-saturation to the diffusive flux at the step position is used as a control parameter. The value of this parameter can vary as much as 50% for different geometries. An approximate analytical formula is derived for the right angled steps geometry. In addition to the kinetic models, this formula can be utilized in macrostep growth models. Finally, numerical modeling of the diffusive and convective transport for equidistant steps is conducted. In particular, the role of fluid flow resulting from the advancement of steps and its contribution to the transport of species to the steps is investigated.

Asynchronous discrete event schemes for PDEs

NASA Astrophysics Data System (ADS)

Stone, D.; Geiger, S.; Lord, G. J.

2017-08-01

A new class of asynchronous discrete-event simulation schemes for advection-diffusion-reaction equations is introduced, based on the principle of allowing quanta of mass to pass through faces of a (regular, structured) Cartesian finite volume grid. The timescales of these events are linked to the flux on the face. The resulting schemes are self-adaptive, and local in both time and space. Experiments are performed on realistic physical systems related to porous media flow applications, including a large 3D advection diffusion equation and advection diffusion reaction systems. The results are compared to highly accurate reference solutions where the temporal evolution is computed with exponential integrator schemes using the same finite volume discretisation. This allows a reliable estimation of the solution error. Our results indicate a first order convergence of the error as a control parameter is decreased, and we outline a framework for analysis.

The path to COVIS: A review of acoustic imaging of hydrothermal flow regimes

NASA Astrophysics Data System (ADS)

Bemis, Karen G.; Silver, Deborah; Xu, Guangyu; Light, Russ; Jackson, Darrell; Jones, Christopher; Ozer, Sedat; Liu, Li

2015-11-01

Acoustic imaging of hydrothermal flow regimes started with the incidental recognition of a plume on a routine sonar scan for obstacles in the path of the human-occupied submersible ALVIN. Developments in sonar engineering, acoustic data processing and scientific visualization have been combined to develop technology which can effectively capture the behavior of focused and diffuse hydrothermal discharge. This paper traces the development of these acoustic imaging techniques for hydrothermal flow regimes from their conception through to the development of the Cabled Observatory Vent Imaging Sonar (COVIS). COVIS has monitored such flow eight times a day for several years. Successful acoustic techniques for estimating plume entrainment, bending, vertical rise, volume flux, and heat flux are presented as is the state-of-the-art in diffuse flow detection.

Concentration dependency in nicotine skin penetration flux from aqueous solutions reflects vehicle induced changes in nicotine stratum corneum retention.

PubMed

Kuswahyuning, Rina; Roberts, Michael S

2014-06-01

This study sought to understand the mechanism by which the steady state flux of nicotine across the human skin from aqueous solutions is markedly decreased at higher nicotine concentrations. Nicotine's steady state flux through human epidermis and its amount in the stratum corneum for a range of aqueous nicotine solutions was determined using Franz diffusion cells, with the nicotine analysed by high performance liquid chromatography (HPLC). Nicotine's thermodynamic activity in the various solutions was estimated from its partial vapour pressure and stratum corneum hydration was determined using a corneometer. The amount of nicotine retained in the stratum corneum was estimated from the nicotine amount found in individual stratum corneum tape strips and a D-Squame determined weight for each strip. The observed steady state flux of nicotine across human epidermis was found to show a parabolic dependence on nicotine concentration, with the flux proportional to its thermodynamic activity up to a concentration of 48% w/w. The nicotine retention in the stratum corneum showed a similar dependency on concentration whereas the diffusivity of nicotine in the stratum corneum appeared to be concentration independent. This retention, in turn, could be estimated from the extent of stratum corneum hydration and the nicotine concentration in the applied solution and volume of water in the skin. Nonlinear dependency of nicotine skin flux on its concentration results from a dehydration induced decrease in its stratum corneum retention at higher concentration and not dehydration induced changes nicotine diffusivity in the stratum corneum.

The mathematical model of radon-222 accumulation in underground mines

NASA Astrophysics Data System (ADS)

Klimshin, A.

2012-04-01

Necessity to control underground mine air radon level arises during building and operating mines as well as auto and railway tunnels including those for metros. Calculation of underground mine air radon level can be fulfilled for estimation of potential radon danger of area for underground structure building. In this work the new mathematical model of radon accumulation in underground mines has been suggested. It takes into consideration underground mine dimensions, air exchange factor and soils ability to emanate radon. The following assumptions have been taken for model development. It is assumed that underground mine is a cylinder of length L and of base area S. Due to ventilation atmosphere air of volume activity Catm, is coming in through one cylinder base and is going out of volume activity Cind from underground mine. Diffusion radon flux is coming in through side surfaces of underground mine. The sources of this flux are radium-226 atoms distributed evenly in rock. For simplification of the task it considered possible to disregard radon emanation by loosened rock and underground waters. As a result of solution of the radon diffusion equation the following expression for calculation of radon volume activity in underground space air has been got: 2·r0 ·λv ·Catm-·l·K0(r0/l)-+D-·K1(r0/l)·C0- Cind = 2·(λ+ λv)·r0 ·l·K0 (r0/l)+ D ·K1(r0/l) . The following designations are used in this expression: Kν(r) - the second genus modified Bessel's function, C0 - equilibrium radon volume activity in soil air, l - diffusion radon length in soil, D - radon diffusion factor, r0 - radius of underground tunnel, λv - factor of air exchange. Expression found may be used for calculation of the minimum factor of necessary air exchange for ensuring safe radon levels in underground spaces. With this worked out model expected levels of radon volume activity were calculated for air in the second metro line underground spaces in the city of Yekaterinburg, Russia.

A finite volume method for trace element diffusion and partitioning during crystal growth

NASA Astrophysics Data System (ADS)

Hesse, Marc A.

2012-09-01

A finite volume method on a uniform grid is presented to compute the polythermal diffusion and partitioning of a trace element during the growth of a porphyroblast crystal in a uniform matrix and in linear, cylindrical and spherical geometry. The motion of the crystal-matrix interface and the thermal evolution are prescribed functions of time. The motion of the interface is discretized and it advances from one cell boundary to next as the prescribed interface position passes the cell center. The appropriate conditions for the flux across the crystal-matrix interface are derived from discrete mass conservation. Numerical results are benchmarked against steady and transient analytic solutions for isothermal diffusion with partitioning and growth. Two applications illustrate the ability of the model to reproduce observed rare-earth element patterns in garnets (Skora et al., 2006) and water concentration profiles around spherulites in obsidian (Watkins et al., 2009). Simulations with diffusion inside the growing crystal show complex concentration evolutions for trace elements with high diffusion coefficients, such as argon or hydrogen, but demonstrate that rare-earth element concentrations in typical metamorphic garnets are not affected by intracrystalline diffusion.

Transport equations of electrodiffusion processes in the laboratory reference frame.

PubMed

Garrido, Javier

2006-02-23

The transport equations of electrodiffusion processes use three reference frames for defining the fluxes: Fick's reference in diffusion, solvent-fixed reference in transference numbers, and laboratory fluxes in electric conductivity. The convenience of using only one reference frame is analyzed here from the point of view of the thermodynamics of irreversible processes. A relation between the fluxes of ions and solvent and the electric current density is deduced first from a mass and volume balance. This is then used to show that (i) the laboratory and Fick's diffusion coefficients are identical and (ii) the transference numbers of both the solvent and the ion in the laboratory reference frame are related. Finally, four experimental methods for the measurement of ion transference numbers are analyzed critically. New expressions for evaluating transference numbers for the moving boundary method and the chronopotentiometry technique are deduced. It is concluded that the ion transport equation in the laboratory reference frame plays a key role in the description of electrodiffusion processes.

A pyramid scheme for three-dimensional diffusion equations on polyhedral meshes

NASA Astrophysics Data System (ADS)

Wang, Shuai; Hang, Xudeng; Yuan, Guangwei

2017-12-01

In this paper, a new cell-centered finite volume scheme is proposed for three-dimensional diffusion equations on polyhedral meshes, which is called as pyramid scheme (P-scheme). The scheme is designed for polyhedral cells with nonplanar cell-faces. The normal flux on a nonplanar cell-face is discretized on a planar face, which is determined by a simple optimization procedure. The resulted discrete form of the normal flux involves only cell-centered and cell-vertex unknowns, and is free from face-centered unknowns. In the case of hexahedral meshes with skewed nonplanar cell-faces, a quite simple expression is obtained for the discrete normal flux. Compared with the second order accurate O-scheme [31], the P-scheme is more robust and the discretization cost is reduced remarkably. Numerical results are presented to show the performance of the P-scheme on various kinds of distorted meshes. In particular, the P-scheme is shown to be second order accurate.

Predicted sensitivity of the KM3NeT/ARCA detector to a diffuse flux of cosmic neutrinos

NASA Astrophysics Data System (ADS)

Coniglione, R.; Fusco, L. A.; Stransky, D.

2016-04-01

The KM3NeT Collaboration has started the construction of a research infrastructure hosting a network of underwater neutrino detectors in the Mediterranean Sea. Two instruments based on the same technology are being built: KM3NeT/ORCA to measure the neutrino mass hierarchy and to study atmospheric neutrino oscillations and KM3NeT/ARCA to detect high-energy cosmic neutrinos both in diffuse and point source mode. The excellent angular resolution of the ARCA detector, with an instrumented volume of about one Gton, will allow for an unprecedented exploration of the neutrino sky searching for neutrinos coming from defined sources of sky regions, like the Galactic Plane and the Fermi Bubbles. It will also look for diffuse high energy neutrino fluxes following the indication provided by the IceCube signal. This contribution will report on the sensitivity of the KM3NeT/ARCA telescope with particular attention to the region of the Galactic Plane. Comparisons with theoretical expectations are also discussed.

Simulations of Plasmasheet Electrons in a Model Magnetosphere with AMIE Potentials: Implications for Diffuse Aurora

NASA Astrophysics Data System (ADS)

Chen, M. W.; Schulz, M.; Lu, G.

2001-12-01

We obtain distributions of precipitating electrons by tracing drift shells of plasmasheet electrons in the limit of strong pitch angle diffusion in Dungey's model magnetosphere, which consists of a dipolar magnetic field plus a uniform southward field. Under strong pitch-angle diffusion particles drift so as to conserve an adiabatic invariant Λ equal to the enclosed phase-space volume (i.e., the cube of the particle momentum p times the occupied flux-tube volume per unit magnetic flux). In the past we applied a quiescent Stern-Volland electric-field model with a cross-tail potential drop of 25 kV and added to it a storm-associated Brice-Nishida cross-magnetospheric electric field with impulses to represent substorm effects. For the present study we use the more realistic Assimilative Model of Ionospheric Electrodynamics (AMIE). We use an analytical expansion to express the AMIE ionospheric potential as a function of latitude and magnetic local time. We map this AMIE potential to latitudes >= 50^o to magnetospheric field lines with (L \\ge 2.5) in Dungey's magnetic field model. We trace the bounce-averaged drift motion of representative plasmasheet electrons for values of \\Lambda corresponding to energies of 0.25-64 keV on field lines of equatorial radial distance r = 6 R_E (L = 5.7), which maps to \\approx 65^o$ latitude in the ionosphere. We use the simulation results to map stormtime phase space distributions taking into account loss due to precipitation. We consider 2 models of electron scattering: (1) the limit of strong scattering everywhere, and (2) an MLT-dependent scattering that is less than everywhere strong in the plasma sheet. From the phase space distributions we calculate the total precipitating electron energy flux into the ionosphere. For this study we focus on the October 19, 1998, storm. We compare qualitatively the simulated energy flux with X-ray intensity from Polar/PIXIE images during this storm.

Gas Flux and Density Surrounding a Cylindrical Aperture in the Free Molecular Flow Regime

NASA Technical Reports Server (NTRS)

Soulas, George C.

2011-01-01

The equations for rigorously calculating the particle flux and density surrounding a cylindrical aperture in the free molecular flow regime are developed and presented. The fundamental equations for particle flux and density from a reservoir and a diffusely reflecting surface will initially be developed. Assumptions will include a Maxwell-Boltzmann speed distribution, equal particle and wall temperatures, and a linear flux distribution along the cylindrical aperture walls. With this information, the equations for axial flux and density surrounding a cylindrical aperture will be developed. The cylindrical aperture will be divided into multiple volumes and regions to rigorously determine the surrounding axial flux and density, and appropriate limits of integration will be determined. The results of these equations will then be evaluated. The linear wall flux distribution assumption will be assessed. The axial flux and density surrounding a cylindrical aperture with a thickness-to-radius ratio of 1.25 will be presented. Finally, the equations determined in this study will be verified using multiple methods.

[Factors affecting benzene diffusion from contaminated soils to the atmosphere and flux characteristics].

PubMed

Du, Ping; Wang, Shi-Jie; Zhao, Huan-Huan; Wu, Bin; Han, Chun-Mei; Fang, Ji-Dun; Li, Hui-Ying; Hosomi, Masaaki; Li, Fa-Sheng

2013-12-01

The influencing factors of benzene diffusion fluxes from sand and black soil to atmosphere were investigated using a flux chamber (30.0 cm x 17.5 cm x 29.0 cm). In this study, the benzene diffusion fluxes were estimated by measuring the benzene concentrations both in the headspace of the chamber and in the soils of different layers. The results indicated that the soil water content played an important role in benzene diffusion fluxes. The diffusion flux showed positive correlation with the initial benzene concentration and the benzene dissolution concentration for both soil types. The changes of air flow rate from 300 to 900 mL x min(-1) and temperature from 20 degrees C to 40 degrees C resulted in increases of the benzene diffusion flux. Our study of benzene diffusion fluxes from contaminated soils will be beneficial for the predicting model, and emergency management and precautions.

Finite element modelling of creep cavity filling by solute diffusion

NASA Astrophysics Data System (ADS)

Versteylen, C. D.; Szymański, N. K.; Sluiter, M. H. F.; van Dijk, N. H.

2018-04-01

In recently discovered self healing creep steels, open-volume creep cavities are filled by the precipitation of supersaturated solute. These creep cavities form on the grain boundaries oriented perpendicular to the applied stress. The presence of a free surface triggers a flux of solute from the matrix, over the grain boundaries towards the creep cavities. We studied the creep cavity filling by finite element modelling and found that the filling time critically depends on (i) the ratio of diffusivities in the grain boundary and the bulk, and (ii) on the ratio of the intercavity distance and the cavity size. For a relatively large intercavity spacing 3D transport is observed when the grain boundary and volume diffusivities are of a similar order of magnitude, while a 2D behaviour is observed when the grain boundary diffusivity is dominant. Instead when the intercavity distance is small, the transport behaviour tends to a 1D behaviour in all cases, as the amount of solute available in the grain boundary is insufficient. A phase diagram with the transition lines is constructed.

Predictive model to describe water migration in cellular solid foods during storage.

PubMed

Voogt, Juliën A; Hirte, Anita; Meinders, Marcel B J

2011-11-01

Water migration in cellular solid foods during storage causes loss of crispness. To improve crispness retention, physical understanding of this process is needed. Mathematical models are suitable tools to gain this physical knowledge. Water migration in cellular solid foods involves migration through both the air cells and the solid matrix. For systems in which the water migration distance is large compared with the cell wall thickness of the solid matrix, the overall water flux through the system is dominated by the flux through the air. For these systems, water migration can be approximated well by a Fickian diffusion model. The effective diffusion coefficient can be expressed in terms of the material properties of the solid matrix (i.e. the density, sorption isotherm and diffusion coefficient of water in the solid matrix) and the morphological properties of the cellular structure (i.e. water vapour permeability and volume fraction of the solid matrix). The water vapour permeability is estimated from finite element method modelling using a simplified model for the cellular structure. It is shown that experimentally observed dynamical water profiles of bread rolls that differ in crust permeability are predicted well by the Fickian diffusion model. Copyright © 2011 Society of Chemical Industry.

Numerical simulation of conservation laws

NASA Technical Reports Server (NTRS)

Chang, Sin-Chung; To, Wai-Ming

1992-01-01

A new numerical framework for solving conservation laws is being developed. This new approach differs substantially from the well established methods, i.e., finite difference, finite volume, finite element and spectral methods, in both concept and methodology. The key features of the current scheme include: (1) direct discretization of the integral forms of conservation laws, (2) treating space and time on the same footing, (3) flux conservation in space and time, and (4) unified treatment of the convection and diffusion fluxes. The model equation considered in the initial study is the standard one dimensional unsteady constant-coefficient convection-diffusion equation. In a stability study, it is shown that the principal and spurious amplification factors of the current scheme, respectively, are structurally similar to those of the leapfrog/DuFort-Frankel scheme. As a result, the current scheme has no numerical diffusion in the special case of pure convection and is unconditionally stable in the special case of pure diffusion. Assuming smooth initial data, it will be shown theoretically and numerically that, by using an easily determined optimal time step, the accuracy of the current scheme may reach a level which is several orders of magnitude higher than that of the MacCormack scheme, with virtually identical operation count.

All-flavor Search for a Diffuse Flux of Cosmic Neutrinos with Nine Years of ANTARES Data

NASA Astrophysics Data System (ADS)

Albert, A.; André, M.; Anghinolfi, M.; Anton, G.; Ardid, M.; Aubert, J.-J.; Aublin, J.; Avgitas, T.; Baret, B.; Barrios-Martí, J.; Basa, S.; Belhorma, B.; Bertin, V.; Biagi, S.; Bormuth, R.; Bourret, S.; Bouwhuis, M. C.; Brânzaş, H.; 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.; Gay, P.; Giordano, V.; Glotin, H.; Grégoire, T.; Ruiz, R. Gracia; Graf, K.; Hallmann, S.; van Haren, H.; Heijboer, A. J.; Hello, Y.; Hernández-Rey, J. J.; Hößl, J.; Hofestädt, J.; 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.; (The ANTARES Collaboration

2018-01-01

The ANTARES detector is at present the most sensitive neutrino telescope in the northern hemisphere. The highly significant cosmic neutrino excess observed by the Antarctic IceCube detector can be studied with ANTARES, exploiting its complementing field of view, exposure, and lower energy threshold. Searches for an all-flavor diffuse neutrino signal, covering nine years of ANTARES data taking, are presented in this Letter. Upward-going events are used to reduce the atmospheric muon background. This work includes for the first time in ANTARES both track-like (mainly {ν }μ ) and shower-like (mainly {ν }e) events in this kind of analysis. Track-like events allow for an increase of the effective volume of the detector thanks to the long path traveled by muons in rock and/or sea water. Shower-like events are well reconstructed only when the neutrino interaction vertex is close to, or inside, the instrumented volume. A mild excess of high-energy events over the expected background is observed in nine years of ANTARES data in both samples. The best fit for a single power-law cosmic neutrino spectrum, in terms of per-flavor flux at 100 TeV, is {{{Φ }}}01f(100 {{TeV}})=(1.7+/- 1.0)× 10‑18 GeV‑1 cm‑2 s‑1 sr‑1 with spectral index {{Γ }}={2.4}-0.4+0.5. The null cosmic flux assumption is rejected with a significance of 1.6σ.

On turbulent diffusion of magnetic fields and the loss of magnetic flux from stars

NASA Technical Reports Server (NTRS)

Vainshtein, Samuel I.; Rosner, Robert

1991-01-01

The turbulent diffusion of magnetic fields in astrophysical objects, and the processes leading to magnetic field flux loss from such objects are discussed with attention to the suppression of turbulent diffusion by back-reaction of magnetic fields on small spatial scales, and on the constraint imposed on magnetic flux loss by flux-freezing within stars. Turbulent magnetic diffusion can be suppressed even for very weak large-scale magnetic fields, so that 'standard' turbulent diffusion is incapable of significant magnetic flux destruction within a star. Finally, magnetic flux loss via winds is shown to be generally ineffective, no matter what the value of the effective magnetic Reynolds number is.

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

Chloride Transport in Porous Lipid Bilayer Membranes

PubMed Central

Andreoli, Thomas E.; Watkins, Mary L.

1973-01-01

This paper describes dissipative Cl- transport in "porous" lipid bilayer membranes, i.e., cholesterol-containing membranes exposed to 1–3 x 10-7 M amphotericin B. P DCl (cm·s-1), the diffusional permeability coefficient for Cl-, estimated from unidirectional 36Cl- fluxes at zero volume flow, varied linearly with the membrane conductance (Gm, Ω-1·cm-2) when the contributions of unstirred layers to the resistance to tracer diffusion were relatively small with respect to the membranes; in 0.05 M NaCl, P DCl was 1.36 x 10-4 cm·s-1 when Gm was 0.02 Ω-1·cm-2. Net chloride fluxes were measured either in the presence of imposed concentration gradients or electrical potential differences. Under both sets of conditions: the values of P DCl computed from zero volume flow experiments described net chloride fluxes; the net chloride fluxes accounted for ∼90–95% of the membrane current density; and, the chloride flux ratio conformed to the Ussing independence relationship. Thus, it is likely that Cl- traversed aqueous pores in these anion-permselective membranes via a simple diffusion process. The zero current membrane potentials measured when the aqueous phases contained asymmetrical NaCl solutions could be expressed in terms of the Goldman-Hodgkin-Katz constant field equation, assuming that the P DNa/P DCl ratio was 0.05. In symmetrical salt solutions, the current-voltage properties of these membranes were linear; in asymmetrical NaCl solutions, the membranes exhibited electrical rectification consistent with constant-field theory. It seems likely that the space charge density in these porous membranes is sufficiently low that the potential gradient within the membranes is approximately linear; and, that the pores are not electrically neutral, presumably because the Debye length within the membrane phase approximates the membrane thickness. PMID:4708408

Fundamental Flux Equations for Fracture-Matrix Interactions with Linear Diffusion

NASA Astrophysics Data System (ADS)

Oldenburg, C. M.; Zhou, Q.; Rutqvist, J.; Birkholzer, J. T.

2017-12-01

The conventional dual-continuum models are only applicable for late-time behavior of pressure propagation in fractured rock, while discrete-fracture-network models may explicitly deal with matrix blocks at high computational expense. To address these issues, we developed a unified-form diffusive flux equation for 1D isotropic (spheres, cylinders, slabs) and 2D/3D rectangular matrix blocks (squares, cubes, rectangles, and rectangular parallelepipeds) by partitioning the entire dimensionless-time domain (Zhou et al., 2017a, b). For each matrix block, this flux equation consists of the early-time solution up until a switch-over time after which the late-time solution is applied to create continuity from early to late time. The early-time solutions are based on three-term polynomial functions in terms of square root of dimensionless time, with the coefficients dependent on dimensionless area-to-volume ratio and aspect ratios for rectangular blocks. For the late-time solutions, one exponential term is needed for isotropic blocks, while a few additional exponential terms are needed for highly anisotropic blocks. The time-partitioning method was also used for calculating pressure/concentration/temperature distribution within a matrix block. The approximate solution contains an error-function solution for early times and an exponential solution for late times, with relative errors less than 0.003. These solutions form the kernel of multirate and multidimensional hydraulic, solute and thermal diffusion in fractured reservoirs.

SToRM: A Model for 2D environmental hydraulics

USGS Publications Warehouse

Simões, Francisco J. M.

2017-01-01

A two-dimensional (depth-averaged) finite volume Godunov-type shallow water model developed for flow over complex topography is presented. The model, SToRM, is based on an unstructured cell-centered finite volume formulation and on nonlinear strong stability preserving Runge-Kutta time stepping schemes. The numerical discretization is founded on the classical and well established shallow water equations in hyperbolic conservative form, but the convective fluxes are calculated using auto-switching Riemann and diffusive numerical fluxes. Computational efficiency is achieved through a parallel implementation based on the OpenMP standard and the Fortran programming language. SToRM’s implementation within a graphical user interface is discussed. Field application of SToRM is illustrated by utilizing it to estimate peak flow discharges in a flooding event of the St. Vrain Creek in Colorado, U.S.A., in 2013, which reached 850 m3/s (~30,000 f3 /s) at the location of this study.

Fem Simulation of Triple Diffusive Natural Convection Along Inclined Plate in Porous Medium: Prescribed Surface Heat, Solute and Nanoparticles Flux

NASA Astrophysics Data System (ADS)

Goyal, M.; Goyal, R.; Bhargava, R.

2017-12-01

In this paper, triple diffusive natural convection under Darcy flow over an inclined plate embedded in a porous medium saturated with a binary base fluid containing nanoparticles and two salts is studied. The model used for the nanofluid is the one which incorporates the effects of Brownian motion and thermophoresis. In addition, the thermal energy equations include regular diffusion and cross-diffusion terms. The vertical surface has the heat, mass and nanoparticle fluxes each prescribed as a power law function of the distance along the wall. The boundary layer equations are transformed into a set of ordinary differential equations with the help of group theory transformations. A wide range of parameter values are chosen to bring out the effect of buoyancy ratio, regular Lewis number and modified Dufour parameters of both salts and nanofluid parameters with varying angle of inclinations. The effects of parameters on the velocity, temperature, solutal and nanoparticles volume fraction profiles, as well as on the important parameters of heat and mass transfer, i.e., the reduced Nusselt, regular and nanofluid Sherwood numbers, are discussed. Such problems find application in extrusion of metals, polymers and ceramics, production of plastic films, insulation of wires and liquid packaging.

Fluctuation of Ultrafiltration Coefficient of Hemodialysis Membrane During Reuse

NASA Astrophysics Data System (ADS)

Arif, Idam; Christin

2010-12-01

Hemodialysis treatment for patient with kidney failure is to regulate body fluid and to excrete waste products of metabolism. The patient blood and the dialyzing solution (dialysate) are flowed counter currently in a dialyzer to allow volume flux of fluid and diffusion of solutes from the blood to the dialysate through a semipermiable membrane. The volume flux of fluid depends on the hydrostatic and the osmotic pressure difference between the blood and the dialysate. It also depends on the membrane parameter that represents how the membrane allows the fluid and the solutes to move across as a result of the pressure difference, known as the ultrafiltration coefficient Kuf. The coefficient depends on the number and the radius of membrane pores for the movement of the fluids and the solutes across the membrane. The measured membrane ultrafiltration coefficient of reused dialyzer shows fluctuation between one uses to another without any significant trend of change. This indicates that the cleaning process carried out before reuse does not cause perfect removal of clots that happen in the previous use. Therefore the unblocked pores are forced to work hardly to obtain targeted volume flux in a certain time of treatment. This may increase the unblocked pore radius. Reuse is stopped when there is indication of blood leakage during the hemodialysis treatment.

An efficient fully-implicit multislope MUSCL method for multiphase flow with gravity in discrete fractured media

NASA Astrophysics Data System (ADS)

Jiang, Jiamin; Younis, Rami M.

2017-06-01

The first-order methods commonly employed in reservoir simulation for computing the convective fluxes introduce excessive numerical diffusion leading to severe smoothing of displacement fronts. We present a fully-implicit cell-centered finite-volume (CCFV) framework that can achieve second-order spatial accuracy on smooth solutions, while at the same time maintain robustness and nonlinear convergence performance. A novel multislope MUSCL method is proposed to construct the required values at edge centroids in a straightforward and effective way by taking advantage of the triangular mesh geometry. In contrast to the monoslope methods in which a unique limited gradient is used, the multislope concept constructs specific scalar slopes for the interpolations on each edge of a given element. Through the edge centroids, the numerical diffusion caused by mesh skewness is reduced, and optimal second order accuracy can be achieved. Moreover, an improved smooth flux-limiter is introduced to ensure monotonicity on non-uniform meshes. The flux-limiter provides high accuracy without degrading nonlinear convergence performance. The CCFV framework is adapted to accommodate a lower-dimensional discrete fracture-matrix (DFM) model. Several numerical tests with discrete fractured system are carried out to demonstrate the efficiency and robustness of the numerical model.

Hybrid Multiscale Finite Volume Method for Advection-Diffusion Equations Subject to Heterogeneous Reactive Boundary Conditions

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

Barajas-Solano, David A.; Tartakovsky, A. M.

2016-10-13

We present a hybrid scheme for the coupling of macro and microscale continuum models for reactive contaminant transport in fractured and porous media. The transport model considered is the advection-dispersion equation, subject to linear heterogeneous reactive boundary conditions. The Multiscale Finite Volume method (MsFV) is employed to define an approximation to the microscale concentration field defined in terms of macroscopic or \\emph{global} degrees of freedom, together with local interpolator and corrector functions capturing microscopic spatial variability. The macroscopic mass balance relations for the MsFV global degrees of freedom are coupled with the macroscopic model, resulting in a global problem for the simultaneous time-stepping of all macroscopic degrees of freedom throughout the domain. In order to perform the hybrid coupling, the micro and macroscale models are applied over overlapping subdomains of the simulation domain, with the overlap denoted as the handshake subdomainmore » $$\\Omega^{hs}$$, over which continuity of concentration and transport fluxes between models is enforced. Continuity of concentration is enforced by posing a restriction relation between models over $$\\Omega^{hs}$$. Continuity of fluxes is enforced by prolongating the macroscopic model fluxes across the boundary of $$\\Omega^{hs}$$ to microscopic resolution. The microscopic interpolator and corrector functions are solutions to local microscopic advection-diffusion problems decoupled from the global degrees of freedom and from each other by virtue of the MsFV decoupling ansatz. The error introduced by the decoupling ansatz is reduced iteratively by the preconditioned GMRES algorithm, with the hybrid MsFV operator serving as the preconditioner.« less

Design and laboratory testing of a chamber device to measure total flux of volatile organic compounds from the unsaturated zone under natural conditions.

PubMed

Tillman, Fred D; Smith, James A

2004-11-01

To determine if an aquifer contaminated with volatile organic compounds (VOCs) has potential for natural remediation, all natural processes affecting the fate and transport of VOCs in the subsurface must be identified and quantified. This research addresses the quantification of air-phase volatile organic compounds (VOCs) leaving the unsaturated zone soil gas and entering the atmosphere-including the additional flux provided by advective soil-gas movement induced by barometric pumping. A simple and easy-to-use device for measuring VOC flux under natural conditions is presented. The vertical flux chamber (VFC) was designed using numerical simulations and evaluated in the laboratory. Mass-balance numerical simulations based on continuously stirred tank reactor equations (CSTR) provided information on flux measurement performance of several sampling configurations with the final chamber configuration measuring greater than 96% of model-simulated fluxes. A laboratory device was constructed to evaluate the flux chamber under both diffusion-only and advection-plus-diffusion transport conditions. The flux chamber measured an average of 82% of 15 diffusion-only fluxes and an average of 95% of 15 additional advection-plus-diffusion flux experiments. The vertical flux chamber has the capability of providing reliable measurement of VOC flux from the unsaturated zone under both diffusion and advection transport conditions.

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.

Optimizing laboratory-based radon flux measurements for sediments.

PubMed

Chanyotha, Supitcha; Kranrod, Chutima; Kritsananuwat, Rawiwan; Lane-Smith, Derek; Burnett, William C

2016-07-01

Radon flux via diffusion from sediments and other materials may be determined in the laboratory by circulating air through the sample and a radon detector in a closed loop. However, this approach is complicated by the necessity of having to determine the total air volume in the system and accounting for any small air leaks that can arise if using extended measurement periods. We designed a simple open-loop configuration that includes a measured mass of wet sediment and water inside a gas-tight reaction flask connected to a drying system and a radon-in-air analyzer. Ambient air flows through two charcoal columns before entering the reaction vessel to eliminate incoming radon. After traveling through the reaction flask, the air passes the drier and the radon analyzer and is then vented. After some time, the radon activity will reach a steady state depending upon the airflow rate. With this approach, the radon flux via diffusion is simply the product of the steady-state radon activity (Bq/m(3)) multiplied by the airflow rate (mL/min). We demonstrated that this setup could produce good results for materials that produce relatively high radon fluxes. We also show that a modified closed system approach, including radon removal of the incoming air by charcoal filtration in a bypass, can produce very good results including samples with very low emission rates. Copyright © 2016 Elsevier Ltd. All rights reserved.

Gyrokinetic modelling of the quasilinear particle flux for plasmas with neutral-beam fuelling

NASA Astrophysics Data System (ADS)

Narita, E.; Honda, M.; Nakata, M.; Yoshida, M.; Takenaga, H.; Hayashi, N.

2018-02-01

A quasilinear particle flux is modelled based on gyrokinetic calculations. The particle flux is estimated by determining factors, namely, coefficients of off-diagonal terms and a particle diffusivity. In this paper, the methodology to estimate the factors is presented using a subset of JT-60U plasmas. First, the coefficients of off-diagonal terms are estimated by linear gyrokinetic calculations. Next, to obtain the particle diffusivity, a semi-empirical approach is taken. Most experimental analyses for particle transport have assumed that turbulent particle fluxes are zero in the core region. On the other hand, even in the stationary state, the plasmas in question have a finite turbulent particle flux due to neutral-beam fuelling. By combining estimates of the experimental turbulent particle flux and the coefficients of off-diagonal terms calculated earlier, the particle diffusivity is obtained. The particle diffusivity should reflect a saturation amplitude of instabilities. The particle diffusivity is investigated in terms of the effects of the linear instability and linear zonal flow response, and it is found that a formula including these effects roughly reproduces the particle diffusivity. The developed framework for prediction of the particle flux is flexible to add terms neglected in the current model. The methodology to estimate the quasilinear particle flux requires so low computational cost that a database consisting of the resultant coefficients of off-diagonal terms and particle diffusivity can be constructed to train a neural network. The development of the methodology is the first step towards a neural-network-based particle transport model for fast prediction of the particle flux.

Edge-based nonlinear diffusion for finite element approximations of convection-diffusion equations and its relation to algebraic flux-correction schemes.

PubMed

Barrenechea, Gabriel R; Burman, Erik; Karakatsani, Fotini

2017-01-01

For the case of approximation of convection-diffusion equations using piecewise affine continuous finite elements a new edge-based nonlinear diffusion operator is proposed that makes the scheme satisfy a discrete maximum principle. The diffusion operator is shown to be Lipschitz continuous and linearity preserving. Using these properties we provide a full stability and error analysis, which, in the diffusion dominated regime, shows existence, uniqueness and optimal convergence. Then the algebraic flux correction method is recalled and we show that the present method can be interpreted as an algebraic flux correction method for a particular definition of the flux limiters. The performance of the method is illustrated on some numerical test cases in two space dimensions.

Enhanced Diffusion of Chlorinated Organic Compounds into Aquitards due to Cracking

NASA Astrophysics Data System (ADS)

Ayral, D.; Otero, M.; Chung, S.; Goltz, M. N.; Huang, J.; Demond, A. H.

2012-12-01

Despite great efforts, remediation of sites contaminated with dense non-aqueous phase liquids (DNAPLs) is very challenging because, even at residual saturations, DNAPLs can act as a long-term source for a dissolved phase contaminant plume. Current models consider the possibility of diffusion and storage of these compounds in unfractured low permeability layers. However, there is a need to consider the impact of cracks, whether naturally occurring or induced by the interaction between low permeable layers and DNAPLs. To evaluate the impact on diffusive fluxes, diffusion coefficients were measured in low permeability materials representative of aquitards at steady-state using the time-lag method. The experimental setup comprised silty soil, packed into a retaining ring, sandwiched in between two reservoirs. The analytical solution for the time-lag method requires constant conditions in the upper and lower reservoirs. The lower reservoir contained pure trichloroethylene (TCE), while the upper reservoir was maintained at a concentration of zero by bubbling air through it, sweeping TCE into toluene trap. In order to predict the flux, the experimental effective diffusion coefficients were used to calculate the flux through uncracked matrix whereas bulk diffusion coefficient was used to calculate flux through the cracks. By using the experimentally-obtained diffusion coefficients and experimentally-measured crack intensity factors (the ratio of the area of cracks to the uncracked area), the total flux was estimated over extended time periods. These calculations, based on experimental data, were used to evaluate if diffusive-based fluxes in the presence of cracks were significantly greater than in the case of diffusion into an uncracked matrix. The enhanced diffusive fluxes were evaluated to determine whether there is the potential for significantly greater storage in the low permeable layers in the case of cracks, or whether the possibility of advective fluxes into the cracks needs to be considered as well.

Biophysical mechanisms of trichloroethene uptake and loss in baldcypress growing in shallow contaminated groundwater

USGS Publications Warehouse

Nietch, C.T.; Morris, J.T.; Vroblesky, D.A.

1999-01-01

Wetland vegetation may be useful in the remediation of shallow contaminated aquifers. Mesocosm experiments were conducted to describe the regulatory mechanisms affecting trichloroethene (TCE) removal rates from groundwater by flood-adapted wetland trees at a contaminated site. TCE flux through baldcypress [Taxodium distichum (L) Rich] seedlings grown in glass- carboys decreased from day to night and from August to December. The diel fluctuation coincided with changes in leaf-level physiology, as the daytime flux was significantly correlated with net photosynthesis but not with respiration at night. A decrease in seedling water use from summer to winter explained the large seasonal difference in TCE flux. A simple model that simulates gas-phase diffusion through aerenchyma tested the importance of diffusion of TCE vapor from roots to the stem. The modeled diffusive flux was within 64% of the observed value during the winter but could only explain 8% of the summer flux. Seedling water use was a good estimator of flux during the summer. Hence, evapotranspiration (ET) in the summer may serve as a good predictor for the potential of TCE removal by baldcypress trees, while diffusive flux may better approximate potential contaminant loss in the winter.Wetland vegetation may be useful in the remediation of shallow contaminated aquifers. Mesocosm experiments were conducted to describe the regulatory mechanisms affecting trichloroethene (TCE) removal rates from groundwater by flood-adapted wetland trees at a contaminated site. TCE flux through baldcypress [Taxodium distichum (L) Rich] seedlings grown in glass-carboys decreased from day to night and from August to December. The diel fluctuation coincided with changes in leaf-level physiology, as the daytime flux was significantly correlated with net photosynthesis but not with respiration at night. A decrease in seedling water use from summer to winter explained the large seasonal difference in TCE flux. A simple model that simulates gas-phase diffusion through aerenchyma tested the importance of diffusion of TCE vapor from roots to the stem. The modeled diffusive flux was within 64% of the observed value during the winter but could only explain 8% of the summer flux. Seedling water use was a good estimator of flux during the summer. Hence, evapotranspiration (ET) in the summer may serve as a good predictor for the potential of TCE removal by baldcypress trees, while diffusive flux may better approximate potential contaminant loss in the winter.

Discontinuous Galerkin method for multicomponent chemically reacting flows and combustion

NASA Astrophysics Data System (ADS)

Lv, Yu; Ihme, Matthias

2014-08-01

This paper presents the development of a discontinuous Galerkin (DG) method for application to chemically reacting flows in subsonic and supersonic regimes under the consideration of variable thermo-viscous-diffusive transport properties, detailed and stiff reaction chemistry, and shock capturing. A hybrid-flux formulation is developed for treatment of the convective fluxes, combining a conservative Riemann-solver and an extended double-flux scheme. A computationally efficient splitting scheme is proposed, in which advection and diffusion operators are solved in the weak form, and the chemically stiff substep is advanced in the strong form using a time-implicit scheme. The discretization of the viscous-diffusive transport terms follows the second form of Bassi and Rebay, and the WENO-based limiter due to Zhong and Shu is extended to multicomponent systems. Boundary conditions are developed for subsonic and supersonic flow conditions, and the algorithm is coupled to thermochemical libraries to account for detailed reaction chemistry and complex transport. The resulting DG method is applied to a series of test cases of increasing physico-chemical complexity. Beginning with one- and two-dimensional multispecies advection and shock-fluid interaction problems, computational efficiency, convergence, and conservation properties are demonstrated. This study is followed by considering a series of detonation and supersonic combustion problems to investigate the convergence-rate and the shock-capturing capability in the presence of one- and multistep reaction chemistry. The DG algorithm is then applied to diffusion-controlled deflagration problems. By examining convergence properties for polynomial order and spatial resolution, and comparing these with second-order finite-volume solutions, it is shown that optimal convergence is achieved and that polynomial refinement provides advantages in better resolving the localized flame structure and complex flow-field features associated with multidimensional and hydrodynamic/thermo-diffusive instabilities in deflagration and detonation systems. Comparisons with standard third- and fifth-order WENO schemes are presented to illustrate the benefit of the DG scheme for application to detonation and multispecies flow/shock-interaction problems.

Contrast-enhanced dynamic and diffusion-weighted magnetic resonance imaging at 3.0 T to assess early-stage nasopharyngeal carcinoma.

PubMed

Ni, Liangping; Liu, Ying

2018-04-01

The present study aimed to assess early-stage nasopharyngeal carcinoma (NPC) with dynamic contrast-enhanced magnetic resonance imaging (DCE-MRI) and diffusion-weighted imaging (DWI) at 3.0 T. A total of 44 patients newly diagnosed with NPC were included in the present study. All patients underwent MR examination at 3.0 T using DCE-MRI and DWI. The volume transfer constant ( K trans ), flux rate constant between extravascular extracellular space and plasma ( K ep ), the volume of extravascular extracellular space per unit volume of tissue ( V e ) and the apparent diffusion coefficient (ADC) of tumours were investigated. Furthermore, the correlation between clinical stages and ADC value and K trans were analysed. The diagnostic accuracy of K trans and ADC were estimated using receiver operating characteristic curves. NPC stage correlated positively with K trans and negatively with ADC values. Additionally, tumour K trans negatively correlated with ADC value. The sensitivity and accuracy of combined K trans and ADC in distinguishing between stage II and stage III and stage III and IV were higher than the values of either measurement used separately. The present study suggested that K trans and ADC derived from DCE-MRI and DWI may be useful to detect stage early NPC accurately. K trans and ADC in combination were superior than either alone.

Will open ocean oxygen stress intensify under climate change?

NASA Astrophysics Data System (ADS)

Gnanadesikan, A.; Dunne, J. P.; John, J.

2011-07-01

Global warming is expected to reduce oxygen solubility and vertical exchange in the ocean, changes which would be expected to result in an increase in the volume of hypoxic waters. A simulation made with a full earth system model with dynamical atmosphere, ocean, sea ice and biogeochemical cycling shows that this holds true if the condition for hypoxia is set relatively high. However, the volume of the most hypoxic waters does not increase under global warming, as these waters actually become more oxygenated. We show that the rise in oxygen is associated with a drop in ventilation time. A term-by-term analysis within the least oxygenated waters shows an increased supply of oxygen due to lateral diffusion. compensating an increase in remineralization within these highly hypoxic waters. This lateral diffusive flux is the result of an increase of ventilation along the Chilean coast, as a drying of the region under global warming opens up a region of wintertime convection in our model.

A Case for Radio Galaxies as the Sources of IceCube's Astrophysical Neutrino Flux

DOE PAGES

Hooper, Dan

2016-09-01

Here, we present an argument that radio galaxies (active galaxies with mis-aligned jets) are likely to be the primary sources of the high-energy astrophysical neutrinos observed by IceCube. In particular, if the gamma-ray emission observed from radio galaxies is generated through the interactions of cosmic-ray protons with gas, these interactions can also produce a population of neutrinos with a flux and spectral shape similar to that measured by IceCube. We present a simple physical model in which high-energy cosmic rays are confined within the volumes of radio galaxies, where they interact with gas to generate the observed diffuse fluxes ofmore » neutrinos and gamma rays. In addition to simultaneously accounting for the observations of Fermi and IceCube, radio galaxies in this model also represent an attractive class of sources for the highest energy cosmic rays.« less

Reactive solid surface morphology variation via ionic diffusion.

PubMed

Sun, Zhenchao; Zhou, Qiang; Fan, Liang-Shih

2012-08-14

In gas-solid reactions, one of the most important factors that determine the overall reaction rate is the solid morphology, which can be characterized by a combination of smooth, convex and concave structures. Generally, the solid surface structure varies in the course of reactions, which is classically noted as being attributed to one or more of the following three mechanisms: mechanical interaction, molar volume change, and sintering. Here we show that if a gas-solid reaction involves the outward ionic diffusion of a solid-phase reactant then this outward ionic diffusion could eventually smooth the surface with an initial concave and/or convex structure. Specifically, the concave surface is filled via a larger outward diffusing surface pointing to the concave valley, whereas the height of the convex surface decreases via a lower outward diffusion flux in the vertical direction. A quantitative 2-D continuum diffusion model is established to analyze these two morphological variation processes, which shows consistent results with the experiments. This surface morphology variation by solid-phase ionic diffusion serves to provide a fourth mechanism that supplements the traditionally acknowledged solid morphology variation or, in general, porosity variation mechanisms in gas-solid reactions.

Finite Volume Methods: Foundation and Analysis

NASA Technical Reports Server (NTRS)

Barth, Timothy; Ohlberger, Mario

2003-01-01

Finite volume methods are a class of discretization schemes that have proven highly successful in approximating the solution of a wide variety of conservation law systems. They are extensively used in fluid mechanics, porous media flow, meteorology, electromagnetics, models of biological processes, semi-conductor device simulation and many other engineering areas governed by conservative systems that can be written in integral control volume form. This article reviews elements of the foundation and analysis of modern finite volume methods. The primary advantages of these methods are numerical robustness through the obtention of discrete maximum (minimum) principles, applicability on very general unstructured meshes, and the intrinsic local conservation properties of the resulting schemes. Throughout this article, specific attention is given to scalar nonlinear hyperbolic conservation laws and the development of high order accurate schemes for discretizing them. A key tool in the design and analysis of finite volume schemes suitable for non-oscillatory discontinuity capturing is discrete maximum principle analysis. A number of building blocks used in the development of numerical schemes possessing local discrete maximum principles are reviewed in one and several space dimensions, e.g. monotone fluxes, E-fluxes, TVD discretization, non-oscillatory reconstruction, slope limiters, positive coefficient schemes, etc. When available, theoretical results concerning a priori and a posteriori error estimates are given. Further advanced topics are then considered such as high order time integration, discretization of diffusion terms and the extension to systems of nonlinear conservation laws.

Hydrodynamic theory of diffusion in two-temperature multicomponent plasmas

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

Ramshaw, J.D.; Chang, C.H.

Detailed numerical simulations of multicomponent plasmas require tractable expressions for species diffusion fluxes, which must be consistent with the given plasma current density J{sub q} to preserve local charge neutrality. The common situation in which J{sub q} = 0 is referred to as ambipolar diffusion. The use of formal kinetic theory in this context leads to results of formidable complexity. We derive simple tractable approximations for the diffusion fluxes in two-temperature multicomponent plasmas by means of a generalization of the hydrodynamical approach used by Maxwell, Stefan, Furry, and Williams. The resulting diffusion fluxes obey generalized Stefan-Maxwell equations that contain drivingmore » forces corresponding to ordinary, forced, pressure, and thermal diffusion. The ordinary diffusion fluxes are driven by gradients in pressure fractions rather than mole fractions. Simplifications due to the small electron mass are systematically exploited and lead to a general expression for the ambipolar electric field in the limit of infinite electrical conductivity. We present a self-consistent effective binary diffusion approximation for the diffusion fluxes. This approximation is well suited to numerical implementation and is currently in use in our LAVA computer code for simulating multicomponent thermal plasmas. Applications to date include a successful simulation of demixing effects in an argon-helium plasma jet, for which selected computational results are presented. Generalizations of the diffusion theory to finite electrical conductivity and nonzero magnetic field are currently in progress.« less

Protein gradients in single cells induced by their coupling to "morphogen"-like diffusion

NASA Astrophysics Data System (ADS)

Nandi, Saroj Kumar; Safran, Sam A.

2018-05-01

One of the many ways cells transmit information within their volume is through steady spatial gradients of different proteins. However, the mechanism through which proteins without any sources or sinks form such single-cell gradients is not yet fully understood. One of the models for such gradient formation, based on differential diffusion, is limited to proteins with large ratios of their diffusion constants or to specific protein-large molecule interactions. We introduce a novel mechanism for gradient formation via the coupling of the proteins within a single cell with a molecule, that we call a "pronogen," whose action is similar to that of morphogens in multi-cell assemblies; the pronogen is produced with a fixed flux at one side of the cell. This coupling results in an effectively non-linear diffusion degradation model for the pronogen dynamics within the cell, which leads to a steady-state gradient of the protein concentration. We use stability analysis to show that these gradients are linearly stable with respect to perturbations.

Constraining Gas Diffusivity-Soil Water Content Relationships in Forest Soils Using Surface Chamber Fluxes and Depth Profiles of Multiple Trace Gases

NASA Astrophysics Data System (ADS)

Dore, J. E.; Kaiser, K.; Seybold, E. C.; McGlynn, B. L.

2012-12-01

Forest soils are sources of carbon dioxide (CO2) to the atmosphere and can act as either sources or sinks of methane (CH4) and nitrous oxide (N2O), depending on redox conditions and other factors. Soil moisture is an important control on microbial activity, redox conditions and gas diffusivity. Direct chamber measurements of soil-air CO2 fluxes are facilitated by the availability of sensitive, portable infrared sensors; however, corresponding CH4 and N2O fluxes typically require the collection of time-course physical samples from the chamber with subsequent analyses by gas chromatography (GC). Vertical profiles of soil gas concentrations may also be used to derive CH4 and N2O fluxes by the gradient method; this method requires much less time and many fewer GC samples than the direct chamber method, but requires that effective soil gas diffusivities are known. In practice, soil gas diffusivity is often difficult to accurately estimate using a modeling approach. In our study, we apply both the chamber and gradient methods to estimate soil trace gas fluxes across a complex Rocky Mountain forested watershed in central Montana. We combine chamber flux measurements of CO2 (by infrared sensor) and CH4 and N2O (by GC) with co-located soil gas profiles to determine effective diffusivity in soil for each gas simultaneously, over-determining the diffusion equations and providing constraints on both the chamber and gradient methodologies. We then relate these soil gas diffusivities to soil type and volumetric water content in an effort to arrive at empirical parameterizations that may be used to estimate gas diffusivities across the watershed, thereby facilitating more accurate, frequent and widespread gradient-based measurements of trace gas fluxes across our study system. Our empirical approach to constraining soil gas diffusivity is well suited for trace gas flux studies over complex landscapes in general.

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

Han, Longtao; Krstic, Predrag; Kaganovich, Igor

Here we find that negative charges on an armchair single-walled carbon nanotube (SWCNT) can significantly enhance the migration of a carbon adatom on the external surfaces of SWCNTs, along the direction of the tube axis. Nanotube charging results in stronger binding of adatoms to SWCNTs and consequent longer lifetimes of adatoms before desorption, which in turn increases their migration distance several orders of magnitude. These results support the hypothesis of diffusion enhanced SWCNT growth in the volume of arc plasma. This process could enhance effective carbon flux to the metal catalyst.

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.

Improved modeling of turbulent forced convection heat transfer in straight ducts

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

Rokni, M.; Sunden, B.

1999-08-01

This investigation concerns numerical calculation of turbulent forced convective heat transfer and fluid flow in their fully developed state at low Reynolds number. The authors have developed a low Reynolds number version of the nonlinear {kappa}-{epsilon} model combined with the heat flux models of simple eddy diffusivity (SED), low Reynolds number version of generalized gradient diffusion hypothesis (GGDH), and wealth {proportional_to} earning {times} time (WET) in general three-dimensional geometries. The numerical approach is based on the finite volume technique with a nonstaggered grid arrangement and the SIMPLEC algorithm. Results have been obtained with the nonlinear {kappa}-{epsilon} model, combined with themore » Lam-Bremhorst and the Abe-Kondoh-Nagano damping functions for low Reynolds numbers.« less

An experimental and theoretical study of radiative extinction of diffusion flames

NASA Technical Reports Server (NTRS)

Wichman, Indrek S.; Atreya, A.

1994-01-01

Our work was primarily theoretical and numerical. We investigated the simplified modeling of heat losses in diffusion flames, then we 'ramped up' the level of complexity in each successive study until the final chapter discussed the general problem of soot/flame interaction. With regard to the specific objective of studying radiative extinction, we conclude that in the steady case a self-extinguishing zero-g flame is unlikely to occur. The soot volume fractions are too small. On the other hand, our work does provide rational means for assessing the mixture of chemical energy release and radiative heat release. It also provides clues for suitable 'tailoring' this balance. Thus heat fluxes to surrounding surfaces can be substantially increased by exploiting and modifying its sooting capability.

Effects of cloudiness on global and diffuse UV irradiance in a high-mountain area

NASA Astrophysics Data System (ADS)

Blumthaler, M.; Ambach, W.; Salzgeber, M.

1994-03-01

At the high-mountain station Jungfraujoch (3576 m a.s.l., Switzerland), measurements of the radiation fluxes were made during 16 periods of six to eight weeks by means of a Robertson—Berger sunburn meter (UVB data), an Eppley UVA radiometer and an Eppley pyranometer. Cloudiness, opacity and altitude of clouds were recorded at 30-minute intervals. A second set of instruments was employed for separate measurement of the diffuse radiation fluxes using shadow bands. The global and diffuse UVA- and UVB radiation fluxes change less with cloudiness than the corresponding total radiation fluxes. When the sun is covered by clouds, the global UVA- and UVB radiation fluxes are also affected less than the global total radiation flux. The roughly equal influence of cloudiness on the UVA- and UVB radiation fluxes suggests that the reduction is influenced more by scattering than by ozone. Also, the share of diffuse irradiance in global irradiance is considerably higher for UVA- and UVB irradiance than for total irradiance. At 50° solar elevation and 0/10 cloudiness, the share is 39% for UVB irradiance, 34% for UVA irradiance and 11% for total irradiance. The increased aerosol turbidity after the eruptions of El Chichon and Pinatubo has caused a significant increase in diffuse total irradiance but has not produced any significant changes in diffuse UVA- and UVB irradiances.

Energetics of multicomponent diffusion in molten CaO-Al 2O 3-SiO 2

NASA Astrophysics Data System (ADS)

Liang, Yan; Davis, Andrew M.

2002-02-01

The energetics of multicomponent diffusion in molten CaO-Al2O3-SiO2 (CAS) were examined experimentally at 1440 to 1650°C and 0.5 to 2 GPa. Two melt compositions were investigated: a haplodacitic melt (25 wt.% CaO, 15% Al2O3, and 60% SiO2) and a haplobasaltic melt (35% CaO, 20% Al2O3, and 45% SiO2). Diffusion matrices were measured in a mass-fixed frame of reference with simple oxides as end-member components and Al2O3 as a dependent variable. Chemical diffusion in molten CAS shows clear evidence of diffusive coupling among the components. The diffusive flux of SiO2 is significantly enhanced whenever there is a large CaO gradient that is oriented in a direction opposite to the SiO2 gradient. This coupling effect is more pronounced in the haplodacitic melt and is likely to be significant in natural magmas of rhyolitic to andesitic compositions. The relative magnitude of coupled chemical diffusion is not very sensitive to changes in temperature and pressure. To a good approximation, the measured diffusion matrices follow well-defined Arrhenius relationships with pressure and reciprocal temperature. Typically, a change in temperature of 100°C results in a relative change in the elements of diffusion matrix of 50 to 100%, whereas a change in pressure of 1 GPa introduces a relative change in elements of diffusion matrix of 4 to 6% for the haplobasalt, and less than 5% for the haplodacite. At a pressure of 1 GPa, the ratios between the major and minor eigenvalues of the diffusion matrix λ1/λ2 are not very sensitive to temperature variations, with an average of 5.5 ± 0.2 for the haplobasalt and 3.7 ± 0.6 for the haplodacite. The activation energies for the major and minor eigenvalues of the diffusion matrix are 215 ± 12 and 240 ± 21 kJ mol-1, respectively, for the haplodacite and 192 ± 8 and 217 ± 14 kJ mol-1 for the haplobasalt. These values are comparable to the activation energies for self-diffusion of calcium and silicon at the same melt compositions and pressure. At a fixed temperature of 1500°C, the ratios λ1/λ2 increase with the increase of pressure, with λ1/λ2 varying from 2.5 to 4.1 (0.5 to 1.3 GPa) for the haplodacite and 4 to 6.5 (0.5 to 2.0 GPa) for the haplobasalt. The activation volumes for the major and minor eigenvalues of the diffusion matrix are 0.31 ± 0.44 and 2.3 ± 0.8 cm3 mol-1, respectively, for the haplodacite and -1.48 ± 0.18 and -0.42 ± 0.24 cm3 mol-1 for the haplobasalt. These values are quite different from the activation volumes for self-diffusion of calcium and silicon at the same melt compositions and temperature. These differences in activation volumes between the two melts likely result from a difference in the structure and thermodynamic properties of the melt between the two compositions (e.g., partial molar volume). Applications of the measured diffusion matrices to quartz crystal dissolution in molten CAS reveal that the activation energy and activation volume for quartz dissolution are almost identical to the activation energy and activation volume for diffusion of the minor or slower eigencomponent of the diffusion matrix. This suggests that the diffusion rate of slow eigencomponent is the rate-limiting factor in isothermal crystal dissolution, a conclusion that is likely to be valid for crystal growth and dissolution in natural magmas when diffusion in liquid is the rate-limiting factor.

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.

Viscosity and viscoelasticity of two-phase systems having diffuse interfaces

NASA Technical Reports Server (NTRS)

Hopper, R. W.

1976-01-01

The equilibrium stability criterion for diffuse interfaces in a two-component solution with a miscibility gap requires that the interdiffusion flux vanish. If the system is continuously deformed, convective fluxes disrupt the equilibrium in the interface regions and induce a counter diffusive flux, which is dissipative and contributes to the apparent viscosity of the mixture. Chemical free energy is recoverably stored, causing viscoelastic phenomena. Both effects are significant.

Communication: A method to compute the transport coefficient of pure fluids diffusing through planar interfaces from equilibrium molecular dynamics simulations.

PubMed

Vermorel, Romain; Oulebsir, Fouad; Galliero, Guillaume

2017-09-14

The computation of diffusion coefficients in molecular systems ranks among the most useful applications of equilibrium molecular dynamics simulations. However, when dealing with the problem of fluid diffusion through vanishingly thin interfaces, classical techniques are not applicable. This is because the volume of space in which molecules diffuse is ill-defined. In such conditions, non-equilibrium techniques allow for the computation of transport coefficients per unit interface width, but their weak point lies in their inability to isolate the contribution of the different physical mechanisms prone to impact the flux of permeating molecules. In this work, we propose a simple and accurate method to compute the diffusional transport coefficient of a pure fluid through a planar interface from equilibrium molecular dynamics simulations, in the form of a diffusion coefficient per unit interface width. In order to demonstrate its validity and accuracy, we apply our method to the case study of a dilute gas diffusing through a smoothly repulsive single-layer porous solid. We believe this complementary technique can benefit to the interpretation of the results obtained on single-layer membranes by means of complex non-equilibrium methods.

Effect of particle- and specimen-level transport on product state in compacted-powder combustion synthesis and thermal debinding of polymers from molded powders

NASA Astrophysics Data System (ADS)

Oliveira, Amir Antonio Martins

The existence of large gradients within particles and fast temporal variations in the temperature and species concentration prevents the use of asymptotic approximations for the closure of the volume-averaged, specimen-level formulations. In this case a solution of the particle-level transport problem is needed to complement the specimen-level volume-averaged equations. Here, the use of combined specimen-level and particle-level models for transport in reactive porous media is demonstrated with two examples. For the gasless compacted-powder combustion synthesis, a three-scale model is developed. The specimen-level model is based on the volume-averaged equations for species and temperature. Local thermal equilibrium is assumed and the macroscopic mass diffusion and convection fluxes are neglected. The particle-level model accounts for the interparticle diffusion (i.e., the liquid migration from liquid-rich to liquid-lean regions) and the intraparticle diffusion (i.e., the species mass diffusion within the product layer formed at the surface of the high melting temperature component). It is found that the interparticle diffusion controls the extent of conversion to the final product, the maximum temperature, and to a smaller degree the propagation velocity. The intraparticle diffusion controls the propagation velocity and to a smaller degree the maximum temperature. The initial stages of thermal degradation of EVA from molded specimens is modeled using volume-averaged equations for the species and empirical models for the kinetics of the thermal degradation, the vapor-liquid equilibrium, and the diffusion coefficient of acetic acid in the molten polymer. It is assumed that a bubble forms when the partial pressure of acetic acid exceeds the external ambient pressure. It is found that the removal of acetic acid is characterized by two regimes, a pre-charge dominated regime and a generation dominated regime. For the development of an optimum debinding schedule, the heating rate is modulated to avoid bubbling, while the concentration and temperature follow the bubble-point line for the mixture. The results show a strong dependence on the presence of a pre-charge. It is shown that isolation of the pre-charge effect by using temporary lower heating rates results in an optimum schedule for which the process time is reduced by over 70% when compared to a constant heating rate schedule.

A 1-D Cryothermal Model of Ceres’ Megaregolith: Predictions for Surface Vapor Flux, Subsurface Temperatures and Pore Ice Distribution

NASA Astrophysics Data System (ADS)

Reynolds, Dylan; Wood, Stephen E.; Bapst, Jonathan; Mehlhaff, Joshua; Griffiths, Stephen G.

2014-11-01

We have applied a self-consistent 1-D model for heat diffusion, vapor diffusion, and ice condensation/sublimation, and surface energy balance to investigate our hypothesis for the source of the recently observed water vapor around Ceres [1]. As described in a companion presentation [2], we find that the estimated global flux of 6 kg/s can be produced by steady-state sublimation of subsurface ice driven by the “geothermal” temperature gradient for a heat flux of 1 mW/m2 - the value estimated for a chondritic abundance of heat-producing elements [3,4]. We will present a detailed description of our Ceres cryothermal diffusion model and comparisons with previous models. One key difference is the use of a new physics-based analytic model (‘MaxRTCM’) for calculating the thermal conductivity (Kth) of planetary regolith [5] that has been validated by comparisons to a wide range of laboratory data [6]. MaxRTCM predicts much lower Kth values in the upper regolith than those in previous work [3]. It also accounts for a process first modeled in a study of unstable equatorial ground ice on Mars [7,8], where vapor diffusing up from a receding ice table toward the surface can recondense at shallower depths - eventually forming a steady-state profile of pore ice volume fraction that increases with depth and maintains a constant flux of vapor at all depths [7]. Using MaxRTCM we calculate the corresponding Kth(z) profiles and will present predictions and implications of the resulting temperature profile in the upper few kilometers of Ceres’ megaregolith.References: [1] Küppers et al. (2014), Nature, 505(7484), 525-527. [2] Wood et al., 2014, this meeting. [3] Fanale & Salvail (1989) Icarus 82, 97-110. [4] McCord and Sotin (2005) JGR 110, E05009. [5] Wood (2013) LPSC Abs. 44, 3077. [6] Wood (2014), Icarus, in revision. [7] Mellon et al. (1997), JGR, 102, 19357-69. [8] Clifford (1993), JGR, 98, 10973-11016.

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

The development of a peak-time criterion for designing controlled-release devices.

PubMed

Simon, Laurent; Ospina, Juan

2016-08-25

This work consists of estimating dynamic characteristics for topically-applied drugs when the magnitude of the flux increases to a maximum value, called peak flux, before declining to zero. This situation is typical of controlled-released systems with a finite donor or vehicle volume. Laplace transforms were applied to the governing equations and resulted in an expression for the flux in terms of the physical characteristics of the system. After approximating this function by a second-order model, three parameters of this reduced structure captured the essential features of the original process. Closed-form relationships were then developed for the peak flux and time-to-peak based on the empirical representation. Three case studies that involve mechanisms, such as diffusion, partitioning, dissolution and elimination, were selected to illustrate the procedure. The technique performed successfully as shown by the ability of the second-order flux to match the prediction of the original transport equations. A main advantage of the proposed method is that it does not require a solution of the original partial differential equations. Less accurate results were noted for longer lag times. Copyright © 2016 Elsevier B.V. All rights reserved.

Inventory of File gdas1.t06z.sfluxgrbf06.grib2

Science.gov Websites

hour ave Visible Diffuse Downward Solar Flux [W/m^2] 036 surface NBDSF 0-6 hour ave Near IR Beam Downward Solar Flux [W/m^2] 037 surface NDDSF 0-6 hour ave Near IR Diffuse Downward Solar Flux [W/m^2] 038

Inventory of File gfs.t06z.sfluxgrbf06.grib2

Science.gov Websites

hour ave Visible Diffuse Downward Solar Flux [W/m^2] 036 surface NBDSF 0-6 hour ave Near IR Beam Downward Solar Flux [W/m^2] 037 surface NDDSF 0-6 hour ave Near IR Diffuse Downward Solar Flux [W/m^2] 038

Influence of heat and particle fluxes nonlocality on spatial distribution of plasma density in two-chamber inductively coupled plasma sources

NASA Astrophysics Data System (ADS)

Kudryavtsev, A. A.; Serditov, K. Yu.

2012-07-01

This study presents 2D simulations of the two-chamber inductively coupled plasma source where power is supplied in the small discharge chamber and extends by electron thermal conductivity mechanism to the big diffusion chamber. Depending on pressure, two main scenarios of plasma density and its spatial distribution behavior were identified. One case is characterized by the localization of plasma in the small driver chamber where power is deposed. Another case describes when the diffusion chamber becomes the main source of plasma with maximum of the electron density. The differences in spatial distribution are caused by local or non-local behavior of electron energy transport in the discharge volume due to different characteristic scale of heat transfer with electronic conductivity.

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

Dynamics and mass transport of solutal convection in a closed porous media system

NASA Astrophysics Data System (ADS)

Wen, Baole; Akhbari, Daria; Hesse, Marc

2016-11-01

Most of the recent studies of CO2 sequestration are performed in open systems where the constant partial pressure of CO2 in the vapor phase results in a time-invariant saturated concentration of CO2 in the brine (Cs). However, in some closed natural CO2 reservoirs, e.g., Bravo Dome in New Mexico, the continuous dissolution of CO2 leads to a pressure drop in the gas that is accompanied by a reduction of Cs and thereby affects the dynamics and mass transport of convection in the brine. In this talk, I discuss the characteristics of convective CO2 dissolution in a closed system. The gas is assumed to be ideal and its solubility given by Henry's law. An analytical solution shows that the diffusive base state is no longer self-similar and that diffusive mass transfer declines rapidly. Scaling analysis reveals that the volume ratio of brine and gas η determines the behavior of the system. DNS show that no constant flux regime exists for η > 0 nevertheless, the quantity F /Cs2 remains constant, where F is the dissolution flux. The onset time is only affected by η when the Rayleigh number Ra is small. In this case, the drop in Cs during the initial diffusive regime significantly reduces the effective Ra and therefore delays the onset.

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.

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.

Mass Transfer in a Nanoscale Material Enhanced by an Opposing Flux

NASA Astrophysics Data System (ADS)

Chmelik, Christian; Bux, Helge; Caro, Jürgen; Heinke, Lars; Hibbe, Florian; Titze, Tobias; Kärger, Jörg

2010-02-01

Diffusion is known to be quantified by measuring the rate of molecular fluxes in the direction of falling concentration. In contrast with intuition, considering methanol diffusion in a novel type of nanoporous material (MOF ZIF-8), this rate has now been found to be enhanced rather than slowed down by an opposing flux of labeled molecules. In terms of the key quantities of random particle movement, this result means that the self-diffusivity exceeds the transport diffusivity. It is rationalized by considering the strong intermolecular interaction and the dominating role of intercage hopping in mass transfer in the systems under study.

Roles of Vacancy/Interstitial Diffusion and Segregation in the Microchemistry at Grain Boundaries of Irradiated Fe-Cr-Ni alloys

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

Yang, Ying; Field, Kevin G.; Allen, Todd R.

2016-02-23

A detailed analysis of the diffusion fluxes near and at grain boundaries of irradiated Fe–Cr–Ni alloys, induced by preferential atom-vacancy and atom-interstitial coupling, is presented. The diffusion flux equations were based on the Perks model formulated through the linear theory of the thermodynamics of irreversible processes. The preferential atom-vacancy coupling was described by the mobility model, whereas the preferential atom-interstitial coupling was described by the interstitial binding model. The composition dependence of the thermodynamic factor was modeled using the CALPHAD approach. The calculated fluxes up to 10 dpa suggested the dominant diffusion mechanism for chromium and iron is via vacancy,more » while that for nickel can swing from the vacancy to the interstitial dominant mechanism. The diffusion flux in the vicinity of a grain boundary was found to be greatly modified by the segregation induced by irradiation, leading to the oscillatory behavior of alloy compositions in this region.« less

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 flux in a model of stochastically gated oxygen transport in insect respiration

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

Berezhkovskii, Alexander M.; Shvartsman, Stanislav Y.

Oxygen delivery to insect tissues is controlled by transport through a branched tubular network that is connected to the atmosphere by valve-like gates, known as spiracles. In certain physiological regimes, the spiracles appear to be randomly switching between open and closed states. Quantitative analysis of this regime leads a reaction-diffusion problem with stochastically switching boundary condition. We derive an expression for the diffusive flux at long times in this problem. Our approach starts with the derivation of the passage probability for a single particle that diffuses between a stochastically gated boundary, which models the opening and closing spiracle, and themore » perfectly absorbing boundary, which models oxygen absorption by the tissue. This passage probability is then used to derive an expression giving the diffusive flux as a function of the geometric parameters of the tube and characteristic time scales of diffusion and gate dynamics.« less

Influence of crystal orientation and ion bombardment on the nitrogen diffusivity in single-crystalline austenitic stainless steel

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

Martinavicius, A.; Abrasonis, G.; Moeller, W.

2011-10-01

The nitrogen diffusivity in single-crystalline AISI 316L austenitic stainless steel (ASS) during ion nitriding has been investigated at different crystal orientations ((001), (110), (111)) under variations of ion flux (0.3-0.7 mA cm{sup -2}), ion energy (0.5-1.2 keV), and temperature (370-430 deg. C). The nitrogen depth profiles obtained from nuclear reaction analysis are in excellent agreement with fits using the model of diffusion under the influence of traps, from which diffusion coefficients were extracted. At fixed ion energy and flux, the diffusivity varies by a factor up to 2.5 at different crystal orientations. At (100) orientation, it increases linearly with increasingmore » ion flux or energy. The findings are discussed on the basis of atomistic mechanisms of interstitial diffusion, potential lattice distortions, local decomposition, and ion-induced lattice vibrational excitations.« less

Diffusive flux in a model of stochastically gated oxygen transport in insect respiration.

PubMed

Berezhkovskii, Alexander M; Shvartsman, Stanislav Y

2016-05-28

Oxygen delivery to insect tissues is controlled by transport through a branched tubular network that is connected to the atmosphere by valve-like gates, known as spiracles. In certain physiological regimes, the spiracles appear to be randomly switching between open and closed states. Quantitative analysis of this regime leads a reaction-diffusion problem with stochastically switching boundary condition. We derive an expression for the diffusive flux at long times in this problem. Our approach starts with the derivation of the passage probability for a single particle that diffuses between a stochastically gated boundary, which models the opening and closing spiracle, and the perfectly absorbing boundary, which models oxygen absorption by the tissue. This passage probability is then used to derive an expression giving the diffusive flux as a function of the geometric parameters of the tube and characteristic time scales of diffusion and gate dynamics.

a Cell Vertex Algorithm for the Incompressible Navier-Stokes Equations on Non-Orthogonal Grids

NASA Astrophysics Data System (ADS)

Jessee, J. P.; Fiveland, W. A.

1996-08-01

The steady, incompressible Navier-Stokes (N-S) equations are discretized using a cell vertex, finite volume method. Quadrilateral and hexahedral meshes are used to represent two- and three-dimensional geometries respectively. The dependent variables include the Cartesian components of velocity and pressure. Advective fluxes are calculated using bounded, high-resolution schemes with a deferred correction procedure to maintain a compact stencil. This treatment insures bounded, non-oscillatory solutions while maintaining low numerical diffusion. The mass and momentum equations are solved with the projection method on a non-staggered grid. The coupling of the pressure and velocity fields is achieved using the Rhie and Chow interpolation scheme modified to provide solutions independent of time steps or relaxation factors. An algebraic multigrid solver is used for the solution of the implicit, linearized equations.A number of test cases are anlaysed and presented. The standard benchmark cases include a lid-driven cavity, flow through a gradual expansion and laminar flow in a three-dimensional curved duct. Predictions are compared with data, results of other workers and with predictions from a structured, cell-centred, control volume algorithm whenever applicable. Sensitivity of results to the advection differencing scheme is investigated by applying a number of higher-order flux limiters: the MINMOD, MUSCL, OSHER, CLAM and SMART schemes. As expected, studies indicate that higher-order schemes largely mitigate the diffusion effects of first-order schemes but also shown no clear preference among the higher-order schemes themselves with respect to accuracy. The effect of the deferred correction procedure on global convergence is discussed.

Magnetic flux concentration and zonal flows in magnetorotational instability turbulence

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

Bai, Xue-Ning; Stone, James M., E-mail: xbai@cfa.harvard.edu

2014-11-20

Accretion disks are likely threaded by external vertical magnetic flux, which enhances the level of turbulence via the magnetorotational instability (MRI). Using shearing-box simulations, we find that such external magnetic flux also strongly enhances the amplitude of banded radial density variations known as zonal flows. Moreover, we report that vertical magnetic flux is strongly concentrated toward low-density regions of the zonal flow. Mean vertical magnetic field can be more than doubled in low-density regions, and reduced to nearly zero in high-density regions in some cases. In ideal MHD, the scale on which magnetic flux concentrates can reach a few diskmore » scale heights. In the non-ideal MHD regime with strong ambipolar diffusion, magnetic flux is concentrated into thin axisymmetric shells at some enhanced level, whose size is typically less than half a scale height. We show that magnetic flux concentration is closely related to the fact that the turbulent diffusivity of the MRI turbulence is anisotropic. In addition to a conventional Ohmic-like turbulent resistivity, we find that there is a correlation between the vertical velocity and horizontal magnetic field fluctuations that produces a mean electric field that acts to anti-diffuse the vertical magnetic flux. The anisotropic turbulent diffusivity has analogies to the Hall effect, and may have important implications for magnetic flux transport in accretion disks. The physical origin of magnetic flux concentration may be related to the development of channel flows followed by magnetic reconnection, which acts to decrease the mass-to-flux ratio in localized regions. The association of enhanced zonal flows with magnetic flux concentration may lead to global pressure bumps in protoplanetary disks that helps trap dust particles and facilitates planet formation.« less

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.

Numerical flux formulas for the Euler and Navier-Stokes equations. 2: Progress in flux-vector splitting

NASA Technical Reports Server (NTRS)

Coirier, William J.; Vanleer, Bram

1991-01-01

The accuracy of various numerical flux functions for the inviscid fluxes when used for Navier-Stokes computations is studied. The flux functions are benchmarked for solutions of the viscous, hypersonic flow past a 10 degree cone at zero angle of attack using first order, upwind spatial differencing. The Harten-Lax/Roe flux is found to give a good boundary layer representation, although its robustness is an issue. Some hybrid flux formulas, where the concepts of flux-vector and flux-difference splitting are combined, are shown to give unsatisfactory pressure distributions; there is still room for improvement. Investigations of low diffusion, pure flux-vector splittings indicate that a pure flux-vector splitting can be developed that eliminates spurious diffusion across the boundary layer. The resulting first-order scheme is marginally stable and not monotone.

Photospheric Magnetic Diffusion by Measuring Moments of Active Regions

NASA Astrophysics Data System (ADS)

Engell, Alexander; Longcope, D.

2013-07-01

Photospheric magnetic surface diffusion is an important constraint for the solar dynamo. The HMI Active Region Patches (HARPs) program automatically identify all magnetic regions above a certain flux. In our study we measure the moments of ARs that are no longer actively emerging and can thereby give us good statistical constraints on photospheric diffusion. We also present the diffusion properties as a function of latitude, flux density, and single polarity (leading or following) within each HARP.

Asymptotic, multigroup flux reconstruction and consistent discontinuity factors

DOE PAGES

Trahan, Travis J.; Larsen, Edward W.

2015-05-12

Recent theoretical work has led to an asymptotically derived expression for reconstructing the neutron flux from lattice functions and multigroup diffusion solutions. The leading-order asymptotic term is the standard expression for flux reconstruction, i.e., it is the product of a shape function, obtained through a lattice calculation, and the multigroup diffusion solution. The first-order asymptotic correction term is significant only where the gradient of the diffusion solution is not small. Inclusion of this first-order correction term can significantly improve the accuracy of the reconstructed flux. One may define discontinuity factors (DFs) to make certain angular moments of the reconstructed fluxmore » continuous across interfaces between assemblies in 1-D. Indeed, the standard assembly discontinuity factors make the zeroth moment (scalar flux) of the reconstructed flux continuous. The inclusion of the correction term in the flux reconstruction provides an additional degree of freedom that can be used to make two angular moments of the reconstructed flux continuous across interfaces by using current DFs in addition to flux DFs. Thus, numerical results demonstrate that using flux and current DFs together can be more accurate than using only flux DFs, and that making the second angular moment continuous can be more accurate than making the zeroth moment continuous.« less

Hydrogen mitigation in submerged arc welding

NASA Astrophysics Data System (ADS)

Klimowicz, Steven

With the role of hydrogen in weld metal well understood in its relation to cold cracking, there has been a push to produce welds with lower and lower diffusible hydrogen contents. The push for lower diffusible hydrogen contents has placed pressure on consumables manufactures to create consumables that can achieve the requirements for lower diffusible hydrogen content. Currently EM12K flux is produced so that it can achieve below 4 ml of diffusible hydrogen for every 100g of weld metal deposited (ml/100g) for submerged arc welding (SAW). The recent trend for industry is to preferentially achieve diffusible hydrogen contents below 3 ml/100g. Making it necessary to find a way to modify the flux to achieve a lower diffusible hydrogen content for the welds it produces. To achieve this goal a two phase plan was developed. The first phase was to characterize the entire welding system for hydrogen. Since the goal of the project is hydrogen mitigation, any amount of hydrogen that could be reduced is helpful and therefore must first be discovered. Sources of hydrogen may be found by analyzing the welding wire and base metal, as well as breaking the flux down into its components and production steps. The wire was analyzed for total hydrogen content as was the base metal. The flux and its components were analyzed using differential thermal analysis-simultaneous thermal analysis (DTA-STA) and later vacuum degassing for moisture content. The analysis of the wire showed that the copper coating on the wire was the largest contributor of hydrogen. There was lubricant present on the wire surface as well, but it did not contribute as much as the copper coating. It was found that a simple low temperature baking of the wire was enough to remove the lubricant and coating moisture. The base metal was found to have a similar total hydrogen content to that of the wire. The breakdown of the flux and production process for moisture content analysis revealed that the production process removes the moisture that is added by the water based binder. The second phase of the project was to modify the flux with fluoride additions to remove hydrogen from the arc while welding. The introduction of fluorine into the arc would lower the amount of hydrogen that may be absorbed as diffusible hydrogen by the weld metal. To select the fluorides a series of thermodynamic calculations were performed as well as simple tests to determine the fluorides behavior in a welding arc and flux. From these tests the following fluorides were selected to be used to be added to EM12K flux as oneweight percent additions: SrF 2, K2TiF6, K2SiF6, and LiF. Welds were then run with the experimental fluxes according to AWS A4.3 standard for diffusible hydrogen testing. From these tests it was found that none experimental fluxes were able to achieve a diffusible hydrogen content lower than the original EM12K flux. It was also found that fluoride reduction in a simple flux is a better predictor of fluoride effectiveness than decomposition temperature.

Diffusion-driven fluid dynamics in ideal gases and plasmas

NASA Astrophysics Data System (ADS)

Vold, E. L.; Yin, L.; Taitano, W.; Molvig, K.; Albright, B. J.

2018-06-01

The classical transport theory based on Chapman-Enskog methods provides self-consistent approximations for the kinetic flux of mass, heat, and momentum in a fluid limit characterized with a small Knudsen number. The species mass fluxes relative to the center of mass, or "diffusive fluxes," are expressed as functions of known gradient quantities with kinetic coefficients evaluated using similar analyses for mixtures of gases or plasma components. The sum over species of the diffusive mass fluxes is constrained to be zero in the Lagrange frame, and thus results in a non-zero molar flux leading to a pressure perturbation. At an interface between two species initially in pressure equilibrium, the pressure perturbation driven by the diffusive molar flux induces a center of mass velocity directed from the species of greater atomic mass towards the lighter atomic mass species. As the ratio of the species particle masses increases, this center of mass velocity carries an increasingly greater portion of the mass across the interface and for a particle mass ratio greater than about two, the center of mass velocity carries more mass than the gradient driven diffusion flux. Early time transients across an interface between two species in a 1D plasma regime and initially in equilibrium are compared using three methods; a fluid code with closure in a classical transport approximation, a particle in cell simulation, and an implicit Fokker-Planck solver for the particle distribution functions. The early time transient phenomenology is shown to be similar in each of the computational simulation methods, including a pressure perturbation associated with the stationary "induced" component of the center of mass velocity which decays to pressure equilibrium during diffusion. At early times, the diffusive process generates pressure and velocity waves which propagate outward from the interface and are required to maintain momentum conservation. The energy in the outgoing waves dissipates as heat in viscous regions, and it is hypothesized that these diffusion driven waves may sustain fluctuations in less viscid finite domains after reflections from the boundaries. These fluid dynamic phenomena are similar in gases or plasmas and occur in flow transients with a moderate Knudsen number. The analysis and simulation results show how the kinetic flux, represented in the fluid transport closure, directly modifies the mass averaged flow described with the Euler equations.

Evaporation Flux Distribution of Drops on a Hydrophilic or Hydrophobic Flat Surface by Molecular Simulations.

PubMed

Xie, Chiyu; Liu, Guangzhi; Wang, Moran

2016-08-16

The evaporation flux distribution of sessile drops is investigated by molecular dynamic simulations. Three evaporating modes are classified, including the diffusion dominant mode, the substrate heating mode, and the environment heating mode. Both hydrophilic and hydrophobic drop-substrate interactions are considered. To count the evaporation flux distribution, which is position dependent, we proposed an azimuthal-angle-based division method under the assumption of spherical crown shape of drops. The modeling results show that the edge evaporation, i.e., near the contact line, is enhanced for hydrophilic drops in all the three modes. The surface diffusion of liquid molecular absorbed on solid substrate for hydrophilic cases plays an important role as well as the space diffusion on the enhanced evaporation rate at the edge. For hydrophobic drops, the edge evaporation flux is higher for the substrate heating mode, but lower than elsewhere of the drop for the diffusion dominant mode; however, a nearly uniform distribution is found for the environment heating mode. The evidence shows that the temperature distribution inside drops plays a key role in the position-dependent evaporation flux.

Conservative and bounded volume-of-fluid advection on unstructured grids

NASA Astrophysics Data System (ADS)

Ivey, Christopher B.; Moin, Parviz

2017-12-01

This paper presents a novel Eulerian-Lagrangian piecewise-linear interface calculation (PLIC) volume-of-fluid (VOF) advection method, which is three-dimensional, unsplit, and discretely conservative and bounded. The approach is developed with reference to a collocated node-based finite-volume two-phase flow solver that utilizes the median-dual mesh constructed from non-convex polyhedra. The proposed advection algorithm satisfies conservation and boundedness of the liquid volume fraction irrespective of the underlying flux polyhedron geometry, which differs from contemporary unsplit VOF schemes that prescribe topologically complicated flux polyhedron geometries in efforts to satisfy conservation. Instead of prescribing complicated flux-polyhedron geometries, which are prone to topological failures, our VOF advection scheme, the non-intersecting flux polyhedron advection (NIFPA) method, builds the flux polyhedron iteratively such that its intersection with neighboring flux polyhedra, and any other unavailable volume, is empty and its total volume matches the calculated flux volume. During each iteration, a candidate nominal flux polyhedron is extruded using an iteration dependent scalar. The candidate is subsequently intersected with the volume guaranteed available to it at the time of the flux calculation to generate the candidate flux polyhedron. The difference in the volume of the candidate flux polyhedron and the actual flux volume is used to calculate extrusion during the next iteration. The choice in nominal flux polyhedron impacts the cost and accuracy of the scheme; however, it does not impact the methods underlying conservation and boundedness. As such, various robust nominal flux polyhedron are proposed and tested using canonical periodic kinematic test cases: Zalesak's disk and two- and three-dimensional deformation. The tests are conducted on the median duals of a quadrilateral and triangular primal mesh, in two-dimensions, and on the median duals of a hexahedral, wedge and tetrahedral primal mesh, in three-dimensions. Comparisons are made with the adaptation of a conventional unsplit VOF advection scheme to our collocated node-based flow solver. Depending on the choice in the nominal flux polyhedron, the NIFPA scheme presented accuracies ranging from zeroth to second order and calculation times that differed by orders of magnitude. For the nominal flux polyhedra which demonstrate second-order accuracy on all tests and meshes, the NIFPA method's cost was comparable to the traditional topologically complex second-order accurate VOF advection scheme.

The Role of Nonlocal Heat Flow in Hohlraums

NASA Astrophysics Data System (ADS)

Town, R. P. J.; Short, R. W.; Verdon, C. P.; Afeyan, B. B.; Glenzer, S. H.; Suter, L. J.

1997-11-01

Glenzer,(Submitted to Physical Review Letters.)* using the Thomson scattering technique, has measured the time evolution of the electron temperature in scale-1 hohlraums. The measured peak electron temperature was 5 keV. Lasnex simulations, using a flux-limited Spitzer heat diffusion model with the standard sharp-cutoff flux limiter of 0.05, gave a peak electron temperature of only 3 keV. Good agreement between simulation and experiment was found when Lasnex simulations employed a time-varying flux limiter, which had a value of 0.01 when the main drive came on. The need to severly inhibit heat transport over the entire volume of hot plasma at late time suggests that nonlocal heat flow could be important in explaining these experimental observations. In this presentation we will report on Fokker--Planck calculations of idealized hohlraums and compare them to standard hydrodynamic calculations using flux-limited Spitzer heat flow. This work was supported by the U.S. Department of Energy Office of Inertial Confinement Fusion under Cooperative Agreement No. DE-FC03-92SF19460. Also, work performed under the auspices of the U.S. Department of Energy by Lawrence Livermore National Laboratory under Contract W-7405-ENG-48.

The pump and leak steady-state concept with a variety of regulated leak pathways.

PubMed

Hoffmann, E K

2001-12-01

The paper will reflect on how Ussing has affected my own scientific work and how he created much of the framework within which I have been working. I have used five examples: (i) The first description of a 1:1 exchange diffusion was introduced by Ussing in 1947 and has been found to be of great physiological significance in most cells. We found that Cl- transport in Ehrlich ascites tumor cells (EATC) was completely dominated by an exchange diffusion process, as defined by Ussing, and, thus, the Cl- conductance was much lower than previously estimated from measurements of unidirectional tracer fluxes. This had a major influence on my later description of a swelling-activated Cl- conductance. (ii) The pump-leak steady-state concept for cell volume control was introduced by Krogh in 1946, but it was developed in detail by Leaf and Ussing in 1959. This concept was the basis for me and others, when we later found that the passive ion leaks play an active role in cell volume control. (iii) The use of isotopes and Ussing's famous flux ratio equation provided an ingenious instrument for distinguishing the various transport routes. We used this to identify the Na,K,2Cl cotransport system as accounting for maintaining a [Cl-]i in the EATC far above thermodynamic equilibrium, as well as accounting for the ion uptake during a regulatory volume increase (RVI) in EATC, similar to what Ussing had found in frog skin. (iv) Short-circuit current setup in the Ussing chamber is still used in laboratories around the world to study ion transport across epithelia. A few results on Cl- transport across the operculum epithelium of the small eurohaline fish Fundulus heteroclitus mounted in small Ussing chambers are presented. (v) Shrinkage-activated Na+ conductance and its possible role in isotonic secretion in frog skin glands is finally discussed.

Transport of Internetwork Magnetic Flux Elements in the Solar Photosphere

NASA Astrophysics Data System (ADS)

Agrawal, Piyush; Rast, Mark P.; Gošić, Milan; Bellot Rubio, Luis R.; Rempel, Matthias

2018-02-01

The motions of small-scale magnetic flux elements in the solar photosphere can provide some measure of the Lagrangian properties of the convective flow. Measurements of these motions have been critical in estimating the turbulent diffusion coefficient in flux-transport dynamo models and in determining the Alfvén wave excitation spectrum for coronal heating models. We examine the motions of internetwork flux elements in Hinode/Narrowband Filter Imager magnetograms and study the scaling of their mean squared displacement and the shape of their displacement probability distribution as a function of time. We find that the mean squared displacement scales super-diffusively with a slope of about 1.48. Super-diffusive scaling has been observed in other studies for temporal increments as small as 5 s, increments over which ballistic scaling would be expected. Using high-cadence MURaM simulations, we show that the observed super-diffusive scaling at short increments is a consequence of random changes in barycenter positions due to flux evolution. We also find that for long temporal increments, beyond granular lifetimes, the observed displacement distribution deviates from that expected for a diffusive process, evolving from Rayleigh to Gaussian. This change in distribution can be modeled analytically by accounting for supergranular advection along with granular motions. These results complicate the interpretation of magnetic element motions as strictly advective or diffusive on short and long timescales and suggest that measurements of magnetic element motions must be used with caution in turbulent diffusion or wave excitation models. We propose that passive tracer motions in measured photospheric flows may yield more robust transport statistics.

Evidence for Secondary Flux Rope Generated by the Electron Kelvin-Helmholtz Instability in a Magnetic Reconnection Diffusion Region

NASA Astrophysics Data System (ADS)

Zhong, Z. H.; Tang, R. X.; Zhou, M.; Deng, X. H.; Pang, Y.; Paterson, W. R.; Giles, B. L.; Burch, J. L.; Tobert, R. B.; Ergun, R. E.; Khotyaintsev, Y. V.; Lindquist, P.-A.

2018-02-01

Secondary flux ropes are suggested to play important roles in energy dissipation and particle acceleration during magnetic reconnection. However, their generation mechanism is not fully understood. In this Letter, we present the first direct evidence that a secondary flux rope was generated due to the evolution of an electron vortex, which was driven by the electron Kelvin-Helmholtz instability in an ion diffusion region as observed by the Magnetospheric Multiscale mission. The subion scale (less than the ion inertial length) flux rope was embedded within the electron vortex, which contained a secondary electron diffusion region at the trailing edge of the flux rope. We propose that intense electron shear flow produced by reconnection generated the electron Kelvin-Helmholtz vortex, which induced a secondary reconnection in the exhaust of the primary X line and then led to the formation of the flux rope. This result strongly suggests that secondary electron Kelvin-Helmholtz instability is important for reconnection dynamics.

Evidence for Secondary Flux Rope Generated by the Electron Kelvin-Helmholtz Instability in a Magnetic Reconnection Diffusion Region.

PubMed

Zhong, Z H; Tang, R X; Zhou, M; Deng, X H; Pang, Y; Paterson, W R; Giles, B L; Burch, J L; Tobert, R B; Ergun, R E; Khotyaintsev, Y V; Lindquist, P-A

2018-02-16

Secondary flux ropes are suggested to play important roles in energy dissipation and particle acceleration during magnetic reconnection. However, their generation mechanism is not fully understood. In this Letter, we present the first direct evidence that a secondary flux rope was generated due to the evolution of an electron vortex, which was driven by the electron Kelvin-Helmholtz instability in an ion diffusion region as observed by the Magnetospheric Multiscale mission. The subion scale (less than the ion inertial length) flux rope was embedded within the electron vortex, which contained a secondary electron diffusion region at the trailing edge of the flux rope. We propose that intense electron shear flow produced by reconnection generated the electron Kelvin-Helmholtz vortex, which induced a secondary reconnection in the exhaust of the primary X line and then led to the formation of the flux rope. This result strongly suggests that secondary electron Kelvin-Helmholtz instability is important for reconnection dynamics.

Size effects in non-linear heat conduction with flux-limited behaviors

NASA Astrophysics Data System (ADS)

Li, Shu-Nan; Cao, Bing-Yang

2017-11-01

Size effects are discussed for several non-linear heat conduction models with flux-limited behaviors, including the phonon hydrodynamic, Lagrange multiplier, hierarchy moment, nonlinear phonon hydrodynamic, tempered diffusion, thermon gas and generalized nonlinear models. For the phonon hydrodynamic, Lagrange multiplier and tempered diffusion models, heat flux will not exist in problems with sufficiently small scale. The existence of heat flux needs the sizes of heat conduction larger than their corresponding critical sizes, which are determined by the physical properties and boundary temperatures. The critical sizes can be regarded as the theoretical limits of the applicable ranges for these non-linear heat conduction models with flux-limited behaviors. For sufficiently small scale heat conduction, the phonon hydrodynamic and Lagrange multiplier models can also predict the theoretical possibility of violating the second law and multiplicity. Comparisons are also made between these non-Fourier models and non-linear Fourier heat conduction in the type of fast diffusion, which can also predict flux-limited behaviors.

Migration of a carbon adatom on a charged single-walled carbon nanotube

DOE PAGES

Han, Longtao; Krstic, Predrag; Kaganovich, Igor; ...

2017-02-02

Here we find that negative charges on an armchair single-walled carbon nanotube (SWCNT) can significantly enhance the migration of a carbon adatom on the external surfaces of SWCNTs, along the direction of the tube axis. Nanotube charging results in stronger binding of adatoms to SWCNTs and consequent longer lifetimes of adatoms before desorption, which in turn increases their migration distance several orders of magnitude. These results support the hypothesis of diffusion enhanced SWCNT growth in the volume of arc plasma. This process could enhance effective carbon flux to the metal catalyst.

Surface flux transport simulations: Effect of inflows toward active regions and random velocities on the evolution of the Sun's large-scale magnetic field

NASA Astrophysics Data System (ADS)

Martin-Belda, D.; Cameron, R. H.

2016-02-01

Aims: We aim to determine the effect of converging flows on the evolution of a bipolar magnetic region (BMR), and to investigate the role of these inflows in the generation of poloidal flux. We also discuss whether the flux dispersal due to turbulent flows can be described as a diffusion process. Methods: We developed a simple surface flux transport model based on point-like magnetic concentrations. We tracked the tilt angle, the magnetic flux and the axial dipole moment of a BMR in simulations with and without inflows and compared the results. To test the diffusion approximation, simulations of random walk dispersal of magnetic features were compared against the predictions of the diffusion treatment. Results: We confirm the validity of the diffusion approximation to describe flux dispersal on large scales. We find that the inflows enhance flux cancellation, but at the same time affect the latitudinal separation of the polarities of the bipolar region. In most cases the latitudinal separation is limited by the inflows, resulting in a reduction of the axial dipole moment of the BMR. However, when the initial tilt angle of the BMR is small, the inflows produce an increase in latitudinal separation that leads to an increase in the axial dipole moment in spite of the enhanced flux destruction. This can give rise to a tilt of the BMR even when the BMR was originally aligned parallel to the equator.

Magnetic Turbulence, Fast Magnetic Field line Diffusion and Small Magnetic Structures in the Solar Wind

NASA Astrophysics Data System (ADS)

Zimbardo, G.; Pommois, P.; Veltri, P.

2003-09-01

The influence of magnetic turbulence on magnetic field line diffusion has been known since the early days of space and plasma physics. However, the importance of ``stochastic diffusion'' for energetic particles has been challenged on the basis of the fact that sharp gradients of either energetic particles or ion composition are often observed in the solar wind. Here we show that fast transverse field line and particle diffusion can coexist with small magnetic structures, sharp gradients, and with long lived magnetic flux tubes. We show, by means of a numerical realization of three dimensional magnetic turbulence and by use of the concepts of deterministic chaos and turbulent transport, that turbulent diffusion is different from Gaussian diffusion, and that transport can be inhomogeneous even if turbulence homogeneously fills the heliosphere. Several diagnostics of field line transport and flux tube evolution are shown, and the size of small magnetic structures in the solar wind, like gradient scales and flux tube thickness, are estimated and compared to the observations.

Heat Diffusion in Gases, Including Effects of Chemical Reaction

NASA Technical Reports Server (NTRS)

Hansen, C. Frederick

1960-01-01

The diffusion of heat through gases is treated where the coefficients of thermal conductivity and diffusivity are functions of temperature. The diffusivity is taken proportional to the integral of thermal conductivity, where the gas is ideal, and is considered constant over the temperature interval in which a chemical reaction occurs. The heat diffusion equation is then solved numerically for a semi-infinite gas medium with constant initial and boundary conditions. These solutions are in a dimensionless form applicable to gases in general, and they are used, along with measured shock velocity and heat flux through a shock reflecting surface, to evaluate the integral of thermal conductivity for air up to 5000 degrees Kelvin. This integral has the properties of a heat flux potential and replaces temperature as the dependent variable for problems of heat diffusion in media with variable coefficients. Examples are given in which the heat flux at the stagnation region of blunt hypersonic bodies is expressed in terms of this potential.

Variability in radial sap flux density patterns and sapwood area among seven co-occurring temperate broad-leaved tree species.

PubMed

Gebauer, Tobias; Horna, Viviana; Leuschner, Christoph

2008-12-01

Forest transpiration estimates are frequently based on xylem sap flux measurements in the outer sections of the hydro-active stem sapwood. We used Granier's constant-heating technique with heating probes at various xylem depths to analyze radial patterns of sap flux density in the sapwood of seven broad-leaved tree species differing in wood density and xylem structure. Study aims were to (1) compare radial sap flux density profiles between diffuse- and ring-porous trees and (2) analyze the relationship between hydro-active sapwood area and stem diameter. In all investigated species except the diffuse-porous beech (Fagus sylvatica L.) and ring-porous ash (Fraxinus excelsior L.), sap flux density peaked at a depth of 1 to 4 cm beneath the cambium, revealing a hump-shaped curve with species-specific slopes. Beech and ash reached maximum sap flux densities immediately beneath the cambium in the youngest annual growth rings. Experiments with dyes showed that the hydro-active sapwood occupied 70 to 90% of the stem cross-sectional area in mature trees of diffuse-porous species, whereas it occupied only about 21% in ring-porous ash. Dendrochronological analyses indicated that vessels in the older sapwood may remain functional for 100 years or more in diffuse-porous species and for up to 27 years in ring-porous ash. We conclude that radial sap flux density patterns are largely dependent on tree species, which may introduce serious bias in sap-flux-derived forest transpiration estimates, if non-specific sap flux profiles are assumed.

An adjoint view on flux consistency and strong wall boundary conditions to the Navier–Stokes equations

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

Stück, Arthur, E-mail: arthur.stueck@dlr.de

2015-11-15

Inconsistent discrete expressions in the boundary treatment of Navier–Stokes solvers and in the definition of force objective functionals can lead to discrete-adjoint boundary treatments that are not a valid representation of the boundary conditions to the corresponding adjoint partial differential equations. The underlying problem is studied for an elementary 1D advection–diffusion problem first using a node-centred finite-volume discretisation. The defect of the boundary operators in the inconsistently defined discrete-adjoint problem leads to oscillations and becomes evident with the additional insight of the continuous-adjoint approach. A homogenisation of the discretisations for the primal boundary treatment and the force objective functional yieldsmore » second-order functional accuracy and eliminates the defect in the discrete-adjoint boundary treatment. Subsequently, the issue is studied for aerodynamic Reynolds-averaged Navier–Stokes problems in conjunction with a standard finite-volume discretisation on median-dual grids and a strong implementation of noslip walls, found in many unstructured general-purpose flow solvers. Going out from a base-line discretisation of force objective functionals which is independent of the boundary treatment in the flow solver, two improved flux-consistent schemes are presented; based on either body wall-defined or farfield-defined control-volumes they resolve the dual inconsistency. The behaviour of the schemes is investigated on a sequence of grids in 2D and 3D.« less

Effects of thermal vapor diffusion on seasonal dynamics of water in the unsaturated zone

USGS Publications Warehouse

Milly, Paul C.D.

1996-01-01

The response of water in the unsaturated zone to seasonal changes of temperature (T) is determined analytically using the theory of nonisothermal water transport in porous media, and the solutions are tested against field observations of moisture potential and bomb fallout isotopic (36Cl and 3H) concentrations. Seasonally varying land surface temperatures and the resulting subsurface temperature gradients induce thermal vapor diffusion. The annual mean vertical temperature gradient is close to zero; however, the annual mean thermal vapor flux is downward, because the temperature‐dependent vapor diffusion coefficient is larger, on average, during downward diffusion (occurring at high T) than during upward diffusion (low T). The annual mean thermal vapor flux is shown to decay exponentially with depth; the depth (about 1 m) at which it decays to e−1of its surface value is one half of the corresponding decay depth for the amplitude of seasonal temperature changes. This depth‐dependent annual mean flux is effectively a source of water, which must be balanced by a flux divergence associated with other transport processes. In a relatively humid environment the liquid fluxes greatly exceed the thermal vapor fluxes, so such a balance is readily achieved without measurable effect on the dynamics of water in the unsaturated zone. However, if the mean vertical water flux through the unsaturated zone is very small (<1 mm y−1), as it may be at many locations in a desert landscape, the thermal vapor flux must be balanced mostly by a matric‐potential‐induced upward flux of water. This return flux may include both vapor and liquid components. Below any near‐surface zone of weather‐related fluctuations of matric potential, maintenance of this upward flux requires an increase with depth in the annual mean matric potential; this theoretical prediction is supported by long‐term field measurements in the Chihuahuan Desert. The analysis also makes predictions, confirmed by the field observations, regarding the seasonal variations of matric potential at a given depth. The conceptual model of unsaturated zone water transport developed here implies the possibility of near‐surface trapping of any aqueous constituent introduced at the surface.

Ionomers for Ion-Conducting Energy Materials

NASA Astrophysics Data System (ADS)

Colby, Ralph

For ionic actuators and battery separators, it is vital to utilize single-ion conducting ionomers that avoid the detrimental polarization of other ions. Single-ion conducting ionomers are synthesized based on DFT calculations, with low glass transition temperatures (facile dynamics) to prepare ion-conducting membranes for battery separators that conduct Li+ or Na+. Characterization by X-ray scattering, dielectric spectroscopy, FTIR, NMR and linear viscoelasticity collectively develop a coherent picture of ionic aggregation and both counterion and polymer dynamics. 7Li NMR diffusion measurements find that diffusion is faster than expected by conductivity using the Nernst-Einstein equation, which means that the majority of Li diffusion occurs by ion pairs moving with the polymer segmental motion. Segmental motion only contributes to ionic conduction in the rare event that one of these ion pairs has an extra Li (a positive triple ion). This leads us to a new metric for ion-conducting soft materials, the product of the cation number density p0 and their diffusion coefficient D; p0D is the diffusive flux of lithium ions. This new metric has a maximum at intermediate ion content that corresponds to the overlap of ion pair polarizability volumes. At higher ion contents, the ion pairs interact strongly and form larger aggregation states that retard segmental motion of both mobile ion pairs and triple ions.

Tritium Fluxes through the Shallow Unsaturated Zone adjacent to a Radioactive Waste Disposal Facility in an Arid Environment

NASA Astrophysics Data System (ADS)

Maples, S.; Andraski, B. J.; Stonestrom, D. A.; Cooper, C. A.; Pohll, G.

2011-12-01

Studies at the U.S. Geological Survey's Amargosa Desert Research Site (ADRS) in southern Nevada have documented long-distance (>400-m) tritium (3H) transport adjacent to a commercial, low-level radioactive waste disposal facility. Transport at this scale is orders of magnitude greater than anticipated; however, lateral 3H fluxes through the shallow unsaturated zone (UZ) have not been investigated in detail. The objective of this study is to estimate and compare lateral and vertical tritiated water-vapor (3HHOg) fluxes in the shallow UZ and their relation to the observed plume migration. Previous studies have recognized two distinct plumes of 3H emanating from the facility. Shallow (0.5 and 1.5-m depth) soil-water vapor samples were collected yearly along 400-m long transects through both plumes from 2003-09. Within the south plume, 3H concentrations at 1.5-m depth have decreased by 44 ± 0.3% during this period, and plume advancement there has effectively ceased (i.e., rate of advance equals rate of decay). During the same period, the west plume showed a net decrease in concentration of 34 ± 0.9% within 100-m of the facility; however, plume advancement is observed at the leading edge of the plume, and concentrations 200-300-m from the facility show an increase in 3H concentration of 64 ± 28.4%. Lateral and vertical diffusive fluxes within both plumes were calculated using 3HHOg concentrations from 2006. Lateral 3HHOg diffusive fluxes within both plumes have been estimated 25-300-m from the facility at 1.5-m depth. Mean lateral 3HHOg diffusive fluxes are 10-14 g m-2 yr-1 within the south plume, and 10-13 g m-2 yr-1 within the west plume. Mean lateral fluxes in the south plume are an order of magnitude lower than in the west plume. This behavior corresponds with the observed relative immobility of the south plume, while the elevated west plume fluxes agree with the plume advancement seen there. Shallow, upward directed, mean vertical 3HHOg fluxes 25-300-m from the facility are estimated to be 10-12 g m-2 yr-1 in the south plume and 10-11 g m-2 yr-1 in the west plume. Within both plumes, mean vertical diffusive fluxes are two orders of magnitude greater than mean lateral diffusive fluxes. Lateral diffusive 3HHOg fluxes have been calculated similarly using 2001 south plume data and were compared to 2001 south plume vertical diffusive 3HHOg fluxes published by Andraski et al. (2005). Here, too, mean vertical fluxes dwarf mean lateral fluxes (10-11 g m-2 yr-1 vs. 10-14 g m-2 yr-1). This behavior highlights the importance of upward movement and release of 3H to the atmosphere. The potential role of advective lateral transport and its contribution to observed plume migration is also under investigation.

Technical Note: A simple calculation algorithm to separate high-resolution CH4 flux measurements into ebullition and diffusion-derived components

NASA Astrophysics Data System (ADS)

Hoffmann, M.; Schulz-Hanke, M.; Garcia Alba, J.; Jurisch, N.; Hagemann, U.; Sachs, T.; Sommer, M.; Augustin, J.

2015-08-01

Processes driving the production, transformation and transport of methane (CH4) in wetland ecosystems are highly complex. Thus, serious challenges are constitutes in terms of the mechanistic process understanding, the identification of potential environmental drivers and the calculation of reliable CH4 emission estimates. We present a simple calculation algorithm to separate open-water CH4 fluxes measured with automatic chambers into diffusion- and ebullition-derived components, which helps facilitating the identification of underlying dynamics and potential environmental drivers. Flux separation is based on ebullition related sudden concentration changes during single measurements. A variable ebullition filter is applied, using the lower and upper quartile and the interquartile range (IQR). Automation of data processing is achieved by using an established R-script, adjusted for the purpose of CH4 flux calculation. The algorithm was tested using flux measurement data (July to September 2013) from a former fen grassland site, converted into a shallow lake as a result of rewetting ebullition and diffusion contributed 46 and 55 %, respectively, to total CH4 emissions, which is comparable to those previously reported by literature. Moreover, the separation algorithm revealed a concealed shift in the diurnal trend of diffusive fluxes throughout the measurement period.

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.

Hemodiafiltration history, technology, and clinical results.

PubMed

Ronco, Claudio; Cruz, Dinna

2007-07-01

Hemodiafiltration (HDF) is an extracorporeal renal-replacement technique using a highly permeable membrane, in which diffusion and convection are conveniently combined to enhance solute removal in a wide spectrum of molecular weights. In this modality, ultrafiltration exceeds the desired fluid loss in the patient, and replacement fluid must be administered to achieve the target fluid balance. Over the years, various HDF variants have emerged, including acetate-free biofiltration, high-volume HDF, internal HDF, paired-filtration dialysis, middilution HDF, double high-flux HDF, push-pull HDF, and online HDF. Recent technology has allowed online production of large volumes of microbiologically ultrapure fluid for reinfusion, greatly simplifying the practice of HDF. Several advantages of HDF over purely diffusive hemodialysis techniques have been described in the literature, including a greater clearance of urea, phosphate, beta(2)-microglobulin and other larger solutes, reduction in dialysis hypotension, and improved anemia management. Although randomized controlled trials have failed to show a survival benefit of HDF, recent data from large observational studies suggest a positive effect of HDF on survival. This article provides a brief review of the history of HDF, the various HDF techniques, and summary of their clinical effects.

Understanding why the volume of suboxic waters does not increase over centuries of global warming in an Earth System Model

NASA Astrophysics Data System (ADS)

Gnanadesikan, A.; Dunne, J. P.; John, J.

2012-03-01

Global warming is expected to reduce oxygen solubility and vertical exchange in the ocean, changes which would be expected to result in an increase in the volume of hypoxic waters. A simulation made with a full Earth System model with dynamical atmosphere, ocean, sea ice and biogeochemical cycling (the Geophysical Fluid Dynamics Laboratory's Earth System Model 2.1) shows that this holds true if the condition for hypoxia is set relatively high. However, the volume of the most hypoxic (i.e., suboxic) waters does not increase under global warming, as these waters actually become more oxygenated. We show that the rise in dissolved oxygen in the tropical Pacific is associated with a drop in ventilation time. A term-by-term analysis within the least oxygenated waters shows an increased supply of dissolved oxygen due to lateral diffusion compensating an increase in remineralization within these highly hypoxic waters. This lateral diffusive flux is the result of an increase of ventilation along the Chilean coast, as a drying of the region under global warming opens up a region of wintertime convection in our model. The results highlight the potential sensitivity of suboxic waters to changes in subtropical ventilation as well as the importance of constraining lateral eddy transport of dissolved oxygen in such waters.

Comparison of the Radiative Two-Flux and Diffusion Approximations

NASA Technical Reports Server (NTRS)

Spuckler, Charles M.

2006-01-01

Approximate solutions are sometimes used to determine the heat transfer and temperatures in a semitransparent material in which conduction and thermal radiation are acting. A comparison of the Milne-Eddington two-flux approximation and the diffusion approximation for combined conduction and radiation heat transfer in a ceramic material was preformed to determine the accuracy of the diffusion solution. A plane gray semitransparent layer without a substrate and a non-gray semitransparent plane layer on an opaque substrate were considered. For the plane gray layer the material is semitransparent for all wavelengths and the scattering and absorption coefficients do not vary with wavelength. For the non-gray plane layer the material is semitransparent with constant absorption and scattering coefficients up to a specified wavelength. At higher wavelengths the non-gray plane layer is assumed to be opaque. The layers are heated on one side and cooled on the other by diffuse radiation and convection. The scattering and absorption coefficients were varied. The error in the diffusion approximation compared to the Milne-Eddington two flux approximation was obtained as a function of scattering coefficient and absorption coefficient. The percent difference in interface temperatures and heat flux through the layer obtained using the Milne-Eddington two-flux and diffusion approximations are presented as a function of scattering coefficient and absorption coefficient. The largest errors occur for high scattering and low absorption except for the back surface temperature of the plane gray layer where the error is also larger at low scattering and low absorption. It is shown that the accuracy of the diffusion approximation can be improved for some scattering and absorption conditions if a reflectance obtained from a Kubelka-Munk type two flux theory is used instead of a reflection obtained from the Fresnel equation. The Kubelka-Munk reflectance accounts for surface reflection and radiation scattered back by internal scattering sites while the Fresnel reflection only accounts for surface reflections.

Diffusion, Fluxes, Friction Forces, and Joule Heating in Two-Temperature Multicomponent Magnetohydrodynamics

NASA Technical Reports Server (NTRS)

Chang, C. H.

1999-01-01

The relationship between Joule heating, diffusion fluxes, and friction forces has been studied for both total and electron thermal energy equations, using general expressions for multicomponent diffusion in two-temperature plasmas with the velocity dependent Lorentz force acting on charged species in a magnetic field. It is shown that the derivation of Joule heating terms requires both diffusion fluxes and friction between species which represents the resistance experienced by the species moving at different relative velocities. It is also shown that the familiar Joule heating term in the electron thermal energy equation includes artificial effects produced by switching the convective velocity from the species velocity to the mass-weighted velocity, and thus should not be ignored even when there is no net energy dissipation.

Examination of factors dominating the sediment-water diffusion flux of DDT-related compounds measured by passive sampling in an urbanized estuarine bay.

PubMed

Feng, Yan; Wu, Chen-Chou; Bao, Lian-Jun; Shi, Lei; Song, Lin; Zeng, Eddy Y

2016-12-01

The fate of hydrophobic organic compounds in aquatic environment are largely determined by their exchange at sediment-water interface, which is highly dynamic and subject to rapidly evolving environmental conditions. In turn, environmental conditions may be governed by both physicochemical parameters and anthropogenic events. To examine the importance of various impact factors, passive sampling devices were deployed at the seafloor of Hailing Bay, an urbanized estuarine bay in Guangdong Province of South China to measure the sediment-water diffusion fluxes of several metabolites of dichlorodiphenyltrichloroethane (DDT), p,p'-DDE, p,p'-DDD and o,p'-DDD. The physicochemical properties of water (temperature, pH, salinity and dissolved oxygen) and surface sediment (sediment organic matter, physical composition, pH, water content, colony forming unit and catalase activity) were also measured. The results showed that the diffusion fluxes of o,p'-DDD, p,p'-DDD and p,p'-DDE at sites A1 and A2 near a fishing boat maintenance facility ranged from 0.42 to 4.73 ng m -2 d -1 (from sediment to overlying water), whereas those at offshore sites varied between -0.03 and -3.02 ng m -2 d -1 (from overlying water to sediment), implicating A1 and A2 as the sources of the target compounds. The distribution patterns of the diffusion fluxes of the target compounds were different from those of water and sediment parameters (water temperature, salinity, sediment texture, pH, colony forming unit and catalase activity) at six sampling sites. This finding suggested that none of these parameters were critical in dictating the sediment-water diffusion fluxes. Besides, decreases in the contents of kerogen and black carbon by 6.7% and 11% would enhance the diffusion fluxes of the target compounds by 11-14% and 12-23%, respectively, at site A1, indicating that kerogen and black carbon were the key factors in mediating the sediment-water diffusion fluxes of DDT-related compounds in field environments. Copyright © 2016 Elsevier Ltd. All rights reserved.

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

Internal loading of phosphorus in a sedimentation pond of a treatment wetland: effect of a phytoplankton crash.

PubMed

Palmer-Felgate, Elizabeth J; Mortimer, Robert J G; Krom, Michael D; Jarvie, Helen P; Williams, Richard J; Spraggs, Rachael E; Stratford, Charlie J

2011-05-01

Sedimentation ponds are widely believed to act as a primary removal process for phosphorus (P) in nutrient treatment wetlands. High frequency in-situ P, ammonium (NH(4)(+)) and dissolved oxygen measurements, alongside occasional water quality measurements, assessed changes in nutrient concentrations and productivity in the sedimentation pond of a treatment wetland between March and June. Diffusive equilibrium in thin films (DET) probes were used to measure in-situ nutrient and chemistry pore-water profiles. Diffusive fluxes across the sediment-water interface were calculated from the pore-water profiles, and dissolved oxygen was used to calculate rates of primary productivity and respiration. The sedimentation pond was a net sink for total P (TP), soluble reactive P (SRP) and NH(4)(+) in March, but became subject to a net internal loading of TP, SRP and NH(4)(+) in May, with SRP concentrations increasing by up to 41μM (1300μl(-1)). Reductions in chlorophyll a and dissolved oxygen concentrations also occurred at this time. The sediment changed from a small net sink of SRP in March (average diffusive flux: -8.2μmolm(-2)day(-1)) to a net source of SRP in June (average diffusive flux: +1324μmolm(-2)day(-1)). A diurnal pattern in water column P concentrations, with maxima in the early hours of the morning, and minima in the afternoon, occurred during May. The diurnal pattern and release of SRP from the sediment were attributed to microbial degradation of diatom biomass, causing reduction of the dissolved oxygen concentration and leading to redox-dependent release of P from the sediment. In June, 2.7mol-Pday(-1) were removed by photosynthesis and 23mol-Pday(-1) were supplied by respiration in the lake volume. SRP was also released through microbial respiration within the water column, including the decomposition of algal matter. It is imperative that consideration to internal recycling is given when maintaining sedimentation ponds, and before the installation of new ponds designed to treat nutrient waste. Copyright © 2011 Elsevier B.V. All rights reserved.

Computations of the three-dimensional flow and heat transfer within a coolant passage of a radial turbine blade

NASA Technical Reports Server (NTRS)

Shih, T. I.-P.; Roelke, R. J.; Steinthorsson, E.

1991-01-01

A numerical code is developed for computing three-dimensional, turbulent, compressible flow within coolant passages of turbine blades. The code is based on a formulation of the compressible Navier-Stokes equations in a rotating frame of reference in which the velocity dependent variable is specified with respect to the rotating frame instead of the inertial frame. The algorithm employed to obtain solutions to the governing equation is a finite-volume LU algorithm that allows convection, source, as well as diffusion terms to be treated implicitly. In this study, all convection terms are upwind differenced by using flux-vector splitting, and all diffusion terms are centrally differenced. This paper describes the formulation and algorithm employed in the code. Some computed solutions for the flow within a coolant passage of a radial turbine are also presented.

Gas transport in unsaturated porous media: the adequacy of Fick's law

USGS Publications Warehouse

Thorstenson, D.C.; Pollock, D.W.

1989-01-01

The increasing use of natural unsaturated zones as repositories for landfills and disposal sites for hazardous wastes (chemical and radioactive) requires a greater understanding of transport processes in the unsaturated zone. For volatile constituents an important potential transport mechanism is gaseous diffusion. Diffusion, however, cannot be treated as an independent isolated transport mechanism. A complete understanding of multicomponent gas transport in porous media (unsaturated zones) requires a knowledge of Knudsen transport, the molecular and nonequimolar components of diffusive flux, and viscous (pressure driven) flux. This review presents a brief discussion of the underlying principles and interrelationships among each of the above flux mechanisms. -from Authors

A microscale turbine driven by diffusive mass flux.

PubMed

Yang, Mingcheng; Liu, Rui; Ripoll, Marisol; Chen, Ke

2015-10-07

An external diffusive mass flux is shown to be able to generate a mechanical torque on a microscale object based on anisotropic diffusiophoresis. In light of this finding, we propose a theoretical prototype micro-turbine driven purely by diffusive mass flux, which is in strong contrast to conventional turbines driven by convective mass flows. The rotational velocity of the proposed turbine is determined by the external concentration gradient, the geometry and the diffusiophoretic properties of the turbine. This scenario is validated by performing computer simulations. Our finding thus provides a new type of chemo-mechanical response which could be used to exploit existing chemical energies at small scales.

Evaluation of the entropy consistent euler flux on 1D and 2D test problems

NASA Astrophysics Data System (ADS)

Roslan, Nur Khairunnisa Hanisah; Ismail, Farzad

2012-06-01

Perhaps most CFD simulations may yield good predictions of pressure and velocity when compared to experimental data. Unfortunately, these results will most likely not adhere to the second law of thermodynamics hence comprising the authenticity of predicted data. Currently, the test of a good CFD code is to check how much entropy is generated in a smooth flow and hope that the numerical entropy produced is of the correct sign when a shock is encountered. Herein, a shock capturing code written in C++ based on a recent entropy consistent Euler flux is developed to simulate 1D and 2D flows. Unlike other finite volume schemes in commercial CFD code, this entropy consistent flux (EC) function precisely satisfies the discrete second law of thermodynamics. This EC flux has an entropy-conserved part, preserving entropy for smooth flows and a numerical diffusion part that will accurately produce the proper amount of entropy, consistent with the second law. Several numerical simulations of the entropy consistent flux have been tested on two dimensional test cases. The first case is a Mach 3 flow over a forward facing step. The second case is a flow over a NACA 0012 airfoil while the third case is a hypersonic flow passing over a 2D cylinder. Local flow quantities such as velocity and pressure are analyzed and then compared with mainly the Roe flux. The results herein show that the EC flux does not capture the unphysical rarefaction shock unlike the Roe-flux and does not easily succumb to the carbuncle phenomenon. In addition, the EC flux maintains good performance in cases where the Roe flux is known to be superior.

Quantifying benthic nitrogen fluxes in Puget Sound, Washington: a review of available data

USGS Publications Warehouse

Sheibley, Richard W.; Paulson, Anthony J.

2014-01-01

Understanding benthic fluxes is important for understanding the fate of materials that settle to the Puget Sound, Washington, seafloor, as well as the impact these fluxes have on the chemical composition and biogeochemical cycles of marine waters. Existing approaches used to measure benthic nitrogen flux in Puget Sound and elsewhere were reviewed and summarized, and factors for considering each approach were evaluated. Factors for selecting an appropriate approach for gathering information about benthic flux include: availability of resources, objectives of projects, and determination of which processes each approach measures. An extensive search of literature was undertaken to summarize known benthic nitrogen fluxes in Puget Sound. A total of 138 individual flux chamber measurements and 38 sets of diffusive fluxes were compiled for this study. Of the diffusive fluxes, 35 new datasets were located, and new flux calculations are presented in this report. About 65 new diffusive flux calculations are provided across all nitrogen species (nitrate, NO3-; nitrite, NO2-; ammonium, NH4+). Data analysis of this newly compiled benthic flux dataset showed that fluxes beneath deep (greater than 50 meters) water tended to be lower than those beneath shallow (less than 50 meters) water. Additionally, variability in flux at the shallow depths was greater, possibly indicating a more dynamic interaction between the benthic and pelagic environments. The overall range of bottom temperatures from studies in the Puget Sound area were small (5–16 degrees Celsius), and only NH4+ flux showed any pattern with temperature. For NH4+, flux values and variability increased at greater than about 12 degrees Celsius. Collection of additional study site metadata about environmental factors (bottom temperature, depth, sediment porosity, sediment type, and sediment organic matter) will help with development of a broader regional understanding benthic nitrogen flux in the Puget Sound.

A new solar cycle model including meridional circulation

NASA Technical Reports Server (NTRS)

Wang, Y.-M.; Sheeley, N. R., Jr.; Nash, A. G.

1991-01-01

A kinematic model is presented for the solar cycle which includes not only the transport of magnetic flux by supergranular diffusion and a poleward bulk flow at the sun's surface, but also the effects of turbulent diffusion and an equatorward 'return flow' beneath the surface. As in the earlier models of Babcock and Leighton, the rotational shearing of a subsurface poloidal field generates toroidal flux that erupts at the surface in the form of bipolar magnetic regions. However, such eruptions do not result in any net loss of toroidal flux from the sun (as assumed by Babcock and Leighton); instead, the large-scale toroidal field is destroyed both by 'unwinding' as the local poloidal field reverses its polarity, and by diffusion as the toroidal flux is transported equatorward by the subsurface flow and merged with its opposite hemisphere counterpart. The inclusion of meridional circulation allows stable oscillations of the magnetic field, accompanied by the equatorward progression of flux eruptions, to be achieved even in the absence of a radial gradient in the angular velocity. An illustrative case in which a subsurface flow speed of order 1 m/s and subsurface diffusion rate of order 10 sq km/s yield 22-yr oscillations in qualitative agreement with observations.

Magnetic flux and heat losses by diffusive, advective, and Nernst effects in MagLIF-like plasma

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

Velikovich, A. L., E-mail: sasha.velikovich@nrl.navy.mil; Giuliani, J. L., E-mail: sasha.velikovich@nrl.navy.mil; Zalesak, S. T.

2014-12-15

The MagLIF approach to inertial confinement fusion involves subsonic/isobaric compression and heating of a DT plasma with frozen-in magnetic flux by a heavy cylindrical liner. The losses of heat and magnetic flux from the plasma to the liner are thereby determined by plasma advection and gradient-driven transport processes, such as thermal conductivity, magnetic field diffusion and thermomagnetic effects. Theoretical analysis based on obtaining exact self-similar solutions of the classical collisional Braginskii's plasma transport equations in one dimension demonstrates that the heat loss from the hot plasma to the cold liner is dominated by the transverse heat conduction and advection, andmore » the corresponding loss of magnetic flux is dominated by advection and the Nernst effect. For a large electron Hall parameter ω{sub e}τ{sub e} effective diffusion coefficients determining the losses of heat and magnetic flux are both shown to decrease with ω{sub e}τ{sub e} as does the Bohm diffusion coefficient, which is commonly associated with low collisionality and two-dimensional transport. This family of exact solutions can be used for verification of codes that model the MagLIF plasma dynamics.« less

An Analytical Diffusion–Expansion Model for Forbush Decreases Caused by Flux Ropes

NASA Astrophysics Data System (ADS)

Dumbović, Mateja; Heber, Bernd; Vršnak, Bojan; Temmer, Manuela; Kirin, Anamarija

2018-06-01

We present an analytical diffusion–expansion Forbush decrease (FD) model ForbMod, which is based on the widely used approach of an initially empty, closed magnetic structure (i.e., flux rope) that fills up slowly with particles by perpendicular diffusion. The model is restricted to explaining only the depression caused by the magnetic structure of the interplanetary coronal mass ejection (ICME). We use remote CME observations and a 3D reconstruction method (the graduated cylindrical shell method) to constrain initial boundary conditions of the FD model and take into account CME evolutionary properties by incorporating flux rope expansion. Several flux rope expansion modes are considered, which can lead to different FD characteristics. In general, the model is qualitatively in agreement with observations, whereas quantitative agreement depends on the diffusion coefficient and the expansion properties (interplay of the diffusion and expansion). A case study was performed to explain the FD observed on 2014 May 30. The observed FD was fitted quite well by ForbMod for all expansion modes using only the diffusion coefficient as a free parameter, where the diffusion parameter was found to correspond to an expected range of values. Our study shows that, in general, the model is able to explain the global properties of an FD caused by a flux rope and can thus be used to help understand the underlying physics in case studies.

Connection between encounter volume and diffusivity in geophysical flows

NASA Astrophysics Data System (ADS)

Rypina, Irina I.; Smith, Stefan G. Llewellyn; Pratt, Larry J.

2018-04-01

Trajectory encounter volume - the volume of fluid that passes close to a reference fluid parcel over some time interval - has been recently introduced as a measure of mixing potential of a flow. Diffusivity is the most commonly used characteristic of turbulent diffusion. We derive the analytical relationship between the encounter volume and diffusivity under the assumption of an isotropic random walk, i.e., diffusive motion, in one and two dimensions. We apply the derived formulas to produce maps of encounter volume and the corresponding diffusivity in the Gulf Stream region of the North Atlantic based on satellite altimetry, and discuss the mixing properties of Gulf Stream rings. Advantages offered by the derived formula for estimating diffusivity from oceanographic data are discussed, as well as applications to other disciplines.

NUMERICAL ANALYSES FOR TREATING DIFFUSION IN SINGLE-, TWO-, AND THREE-PHASE BINARY ALLOY SYSTEMS

NASA Technical Reports Server (NTRS)

Tenney, D. R.

1994-01-01

This package consists of a series of three computer programs for treating one-dimensional transient diffusion problems in single and multiple phase binary alloy systems. An accurate understanding of the diffusion process is important in the development and production of binary alloys. Previous solutions of the diffusion equations were highly restricted in their scope and application. The finite-difference solutions developed for this package are applicable for planar, cylindrical, and spherical geometries with any diffusion-zone size and any continuous variation of the diffusion coefficient with concentration. Special techniques were included to account for differences in modal volumes, initiation and growth of an intermediate phase, disappearance of a phase, and the presence of an initial composition profile in the specimen. In each analysis, an effort was made to achieve good accuracy while minimizing computation time. The solutions to the diffusion equations for single-, two-, and threephase binary alloy systems are numerically calculated by the three programs NAD1, NAD2, and NAD3. NAD1 treats the diffusion between pure metals which belong to a single-phase system. Diffusion in this system is described by a one-dimensional Fick's second law and will result in a continuous composition variation. For computational purposes, Fick's second law is expressed as an explicit second-order finite difference equation. Finite difference calculations are made by choosing the grid spacing small enough to give convergent solutions of acceptable accuracy. NAD2 treats diffusion between pure metals which form a two-phase system. Diffusion in the twophase system is described by two partial differential equations (a Fick's second law for each phase) and an interface-flux-balance equation which describes the location of the interface. Actual interface motion is obtained by a mass conservation procedure. To account for changes in the thicknesses of the two phases as diffusion progresses, a variable grid technique developed by Murray and Landis is employed. These equations are expressed in finite difference form and solved numerically. Program NAD3 treats diffusion between pure metals which form a two-phase system with an intermediate third phase. Diffusion in the three-phase system is described by three partial differential expressions of Fick's second law and two interface-flux-balance equations. As with the two-phase case, a variable grid finite difference is used to numerically solve the diffusion equations. Computation time is minimized without sacrificing solution accuracy by treating the three-phase problem as a two-phase problem when the thickness of the intermediate phase is less than a preset value. Comparisons between these programs and other solutions have shown excellent agreement. The programs are written in FORTRAN IV for batch execution on the CDC 6600 with a central memory requirement of approximately 51K (octal) 60 bit words.

Spatial variability of the Arctic Ocean's double-diffusive staircase

NASA Astrophysics Data System (ADS)

Shibley, N. C.; Timmermans, M.-L.; Carpenter, J. R.; Toole, J. M.

2017-02-01

The Arctic Ocean thermohaline stratification frequently exhibits a staircase structure overlying the Atlantic Water Layer that can be attributed to the diffusive form of double-diffusive convection. The staircase consists of multiple layers of O(1) m in thickness separated by sharp interfaces, across which temperature and salinity change abruptly. Through a detailed analysis of Ice-Tethered Profiler measurements from 2004 to 2013, the double-diffusive staircase structure is characterized across the entire Arctic Ocean. We demonstrate how the large-scale Arctic Ocean circulation influences the small-scale staircase properties. These staircase properties (layer thicknesses and temperature and salinity jumps across interfaces) are examined in relation to a bulk vertical density ratio spanning the staircase stratification. We show that the Lomonosov Ridge serves as an approximate boundary between regions of low density ratio (approximately 3-4) on the Eurasian side and higher density ratio (approximately 6-7) on the Canadian side. We find that the Eurasian Basin staircase is characterized by fewer, thinner layers than that in the Canadian Basin, although the margins of all basins are characterized by relatively thin layers and the absence of a well-defined staircase. A double-diffusive 4/3 flux law parametrization is used to estimate vertical heat fluxes in the Canadian Basin to be O(0.1) W m-2. It is shown that the 4/3 flux law may not be an appropriate representation of heat fluxes through the Eurasian Basin staircase. Here molecular heat fluxes are estimated to be between O(0.01) and O(0.1) W m-2. However, many uncertainties remain about the exact nature of these fluxes.

BOXER: Fine-flux Cross Section Condensation, 2D Few Group Diffusion and Transport Burnup Calculations

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

2010-02-01

Neutron transport, calculation of multiplication factor and neutron fluxes in 2-D configurations: cell calculations, 2-D diffusion and transport, and burnup. Preparation of a cross section library for the code BOXER from a basic library in ENDF/B format (ETOBOX).

Diffuse versus discrete venting at the Tour Eiffel vent site, Lucky Strike hydrothermal field

NASA Astrophysics Data System (ADS)

Mittelstaedt, E. L.; Escartin, J.; Gracias, N.; Olive, J. L.; Barreyre, T.; Davaille, A. B.; Cannat, M.

2010-12-01

Two styles of fluid flow at the seafloor are widely recognized: (1) localized outflows of high temperature (>300°C) fluids, often black or grey color in color (“black smokers”) and (2) diffuse, lower temperature (<100°C), fluids typically transparent and which escape through fractures, porous rock, and sediment. The partitioning of heat flux between these two types of hydrothermal venting is debated and estimates of the proportion of heat carried by diffuse flow at ridge axes range from 20% to 90% of the total axial heat flux. Here, we attempt to improve estimates of this partitioning by carefully characterizing the heat fluxes carried by diffuse and discrete flows at a single vent site, Tour Eiffel in the Lucky Strike hydrothermal field along the Mid-Atlantic Ridge. Fluid temperature and video data were acquired during the recent Bathyluck’09 cruise to the Lucky Strike hydrothermal field (September, 2009) by Victor aboard “Pourquoi Pas?” (IFREMER, France). Temperature measurements were made of fluid exiting discrete vents, of diffuse effluents immediately above the seafloor, and of vertical temperature gradients within discrete hydrothermal plumes. Video data allow us to calculate the fluid velocity field associated with these outflows: for diffuse fluids, Diffuse Flow Velocimetry tracks the displacement of refractive index anomalies through time; for individual hydrothermal plumes, Particle Image Velocimetry tracks eddies by cross-correlation of pixels intensities between subsequent images. Diffuse fluids exhibit temperatures of 8-60°C and fluid velocities of ~1-10 cm s-1. Discrete outflows at 204-300°C have velocities of ~1-2 m s-1. Combined fluid flow velocities, temperature measurements, and full image mosaics of the actively venting areas are used to estimate heat flux of both individual discrete vents and diffuse outflow. The total integrated heat flux and the partitioning between diffuse and discrete venting at Tour Eiffel, and its implications for the nature of hydrothermal activity across the Lucky Strike site are discussed along with the implications for crustal permeability, associated ecosystems, and mid-ocean ridge processes.

Infrasound Waveform Inversion and Mass Flux Validation from Sakurajima Volcano, Japan

NASA Astrophysics Data System (ADS)

Fee, D.; Kim, K.; Yokoo, A.; Izbekov, P. E.; Lopez, T. M.; Prata, F.; Ahonen, P.; Kazahaya, R.; Nakamichi, H.; Iguchi, M.

2015-12-01

Recent advances in numerical wave propagation modeling and station coverage have permitted robust inversion of infrasound data from volcanic explosions. Complex topography and crater morphology have been shown to substantially affect the infrasound waveform, suggesting that homogeneous acoustic propagation assumptions are invalid. Infrasound waveform inversion provides an exciting tool to accurately characterize emission volume and mass flux from both volcanic and non-volcanic explosions. Mass flux, arguably the most sought-after parameter from a volcanic eruption, can be determined from the volume flux using infrasound waveform inversion if the volcanic flow is well-characterized. Thus far, infrasound-based volume and mass flux estimates have yet to be validated. In February 2015 we deployed six infrasound stations around the explosive Sakurajima Volcano, Japan for 8 days. Here we present our full waveform inversion method and volume and mass flux estimates of numerous high amplitude explosions using a high resolution DEM and 3-D Finite Difference Time Domain modeling. Application of this technique to volcanic eruptions may produce realistic estimates of mass flux and plume height necessary for volcanic hazard mitigation. Several ground-based instruments and methods are used to independently determine the volume, composition, and mass flux of individual volcanic explosions. Specifically, we use ground-based ash sampling, multispectral infrared imagery, UV spectrometry, and multigas data to estimate the plume composition and flux. Unique tiltmeter data from underground tunnels at Sakurajima also provides a way to estimate the volume and mass of each explosion. In this presentation we compare the volume and mass flux estimates derived from the different methods and discuss sources of error and future improvements.

Large-scale patterns in summer diffusive CH4 fluxes across boreal lakes, and contribution to diffusive C emissions.

PubMed

Rasilo, Terhi; Prairie, Yves T; Del Giorgio, Paul A

2015-03-01

Lakes are a major component of boreal landscapes, and whereas lake CO2 emissions are recognized as a major component of regional C budgets, there is still much uncertainty associated to lake CH4 fluxes. Here, we present a large-scale study of the magnitude and regulation of boreal lake summer diffusive CH4 fluxes, and their contribution to total lake carbon (C) emissions, based on in situ measurements of concentration and fluxes of CH4 and CO2 in 224 lakes across a wide range of lake type and environmental gradients in Québec. The diffusive CH4 flux was highly variable (mean 11.6 ± 26.4 SD mg m(-2)  d(-1) ), and it was positively correlated with temperature and lake nutrient status, and negatively correlated with lake area and colored dissolved organic matter (CDOM). The relationship between CH4 and CO2 concentrations fluxes was weak, suggesting major differences in their respective sources and/or regulation. For example, increasing water temperature leads to higher CH4 flux but does not significantly affect CO2 flux, whereas increasing CDOM concentration leads to higher CO2 flux but lower CH4 flux. CH4 contributed to 8 ± 23% to the total lake C emissions (CH4  + CO2 ), but 18 ± 25% to the total flux in terms of atmospheric warming potential, expressed as CO2 -equivalents. The incorporation of ebullition and plant-mediated CH4 fluxes would further increase the importance of lake CH4 . The average Q10 of CH4 flux was 3.7, once other covarying factors were accounted for, but this apparent Q10 varied with lake morphometry and was higher for shallow lakes. We conclude that global climate change and the resulting shifts in temperature will strongly influence lake CH4 fluxes across the boreal biome, but these climate effects may be altered by regional patterns in lake morphometry, nutrient status, and browning. © 2014 John Wiley & Sons Ltd.

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

PubMed Central

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

2011-01-01

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. PMID:21177946

Double-diffusive instabilities in ancient seawater

NASA Astrophysics Data System (ADS)

Pawlowicz, Rich; Scheifele, Ben; Zaloga, Artem; Wuest, Alfred; Sommer, Tobias

2015-04-01

Powell Lake, British Columbia, Canada is a geothermally heated lake about 350m deep with a saline lower layer that was isolated from the ocean by coastal uplift about 11000 years ago, after the last ice age. Careful temperature and conductivity profiling measurements show consistent, stable, and spatially/temporally coherent steps resulting from double-diffusive processes in certain ranges of depth, vertically interspersed with other depth ranges where these signatures are not present. These features are quasi-stable for at least several years. Although molecular diffusion has removed about half the salt from the deepest waters and biogeochemical processes have slightly modified the water composition, the lack of tidal processes and shear-driven mixing, as well as an accurate estimate of heat flux from both sediment heat flux measurements and gradient measurements in a region not susceptible to diffusive instabilities, makes this a unique geophysical laboratory to study double diffusion. Here we present a detailed picture of the structure of Powell Lake and its double-diffusive stair cases, and suggest shortcomings with existing parameterizations for fluxes through such staircases.

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.