Generalized diffusion approximation for highly absorbing media and small source-detector separations
NASA Astrophysics Data System (ADS)
Zhu, Dan; Wu, Guiling; Luo, Qingming; Gong, Hui
2006-01-01
The diffusion approximation to the transport equation is commonly used in biomedical optical diagnostic techniques, but constrains its applicability to highly scattering system. The generalized diffusion approximation was developed by Venugopalan can be used to quantify optical properties of turbid media using small source detector separations and allow the measurement of media with highly absorption. Unfortunately, the simulated result from this theory was larger than the real value because δ-Eddington phase function contained too much forward scattering. Here a new independent control parameter is introduced to δ-Eddington phase function so as to modify the generalized diffusion approximation presented. The solution is presented in the stationary case for infinite media with a collimated source of finite size exhibiting spherical symmetry. The solution is compared to results given by the conventional diffusion theory, the generalized diffusion approximation as well as to the Monte-Carlo simulation in steady state diffusion equation for slab boundary condition. The simulation results show that the modified generalized diffusion approximation with an appropriate control parameter is more closed to Monte-Carlo simulation. The modified generalized formulation of diffusion theory presented here may enable the quantitative application of present optical diagnostic techniques to turbid systems which are more highly absorbing and allow these systems to be probed using smaller source-detector separations.
Multiple-source optical diffusion approximation for a multilayer scattering medium.
Hollmann, Joseph L; Wang, Lihong V
2007-08-10
A method for improving the accuracy of the optical diffusion theory for a multilayer scattering medium is presented. An infinitesimally narrow incident light beam is replaced by multiple isotropic point sources of different strengths that are placed in the scattering medium along the incident beam. The multiple sources are then used to develop a multilayer diffusion theory. Diffuse reflectance is then computed using the multilayer diffusion theory and compared with accurate data computed by the Monte Carlo method. This multisource method is found to be significantly more accurate than the previous single-source method. PMID:17694156
NASA Astrophysics Data System (ADS)
Pan, Kuo-Chuan; Liebendörfer, Matthias; Hempel, Matthias; Thielemann, Friedrich-Karl
2016-01-01
The neutrino mechanism of core-collapse supernova is investigated via non-relativistic, two-dimensional (2D), neutrino radiation-hydrodynamic simulations. For the transport of electron flavor neutrinos, we use the interaction rates defined by Bruenn and the isotropic diffusion source approximation (IDSA) scheme, which decomposes the transported particles into trapped-particle and streaming-particle components. Heavy neutrinos are described by a leakage scheme. Unlike the “ray-by-ray” approach in some other multidimensional supernova models, we use cylindrical coordinates and solve the trapped-particle component in multiple dimensions, improving the proto-neutron star resolution and the neutrino transport in angular and temporal directions. We provide an IDSA verification by performing one-dimensional (1D) and 2D simulations with 15 and 20 M⊙ progenitors from Woosley et al. and discuss the difference between our IDSA results and those existing in the literature. Additionally, we perform Newtonian 1D and 2D simulations from prebounce core collapse to several hundred milliseconds postbounce with 11, 15, 21, and 27 M⊙ progenitors from Woosley et al. with the HS(DD2) equation of state. General-relativistic effects are neglected. We obtain robust explosions with diagnostic energies Edia ≳ 0.1-0.5 B (1 B ≡ 1051 erg) for all considered 2D models within approximately 100-300 ms after bounce and find that explosions are mostly dominated by the neutrino-driven convection, although standing accretion shock instabilities are observed as well. We also find that the level of electron deleptonization during collapse dramatically affects the postbounce evolution, e.g., the neglect of neutrino-electron scattering during collapse will lead to a stronger explosion.
Approximate Solutions Of Equations Of Steady Diffusion
NASA Technical Reports Server (NTRS)
Edmonds, Larry D.
1992-01-01
Rigorous analysis yields reliable criteria for "best-fit" functions. Improved "curve-fitting" method yields approximate solutions to differential equations of steady-state diffusion. Method applies to problems in which rates of diffusion depend linearly or nonlinearly on concentrations of diffusants, approximate solutions analytic or numerical, and boundary conditions of Dirichlet type, of Neumann type, or mixture of both types. Applied to equations for diffusion of charge carriers in semiconductors in which mobilities and lifetimes of charge carriers depend on concentrations.
Approximate factorization with source terms
NASA Technical Reports Server (NTRS)
Shih, T. I.-P.; Chyu, W. J.
1991-01-01
A comparative evaluation is made of three methodologies with a view to that which offers the best approximate factorization error. While two of these methods are found to lead to more efficient algorithms in cases where factors which do not contain source terms can be diagonalized, the third method used generates the lowest approximate factorization error. This method may be preferred when the norms of source terms are large, and transient solutions are of interest.
Microfabricated diffusion source
Oborny, Michael C.; Frye-Mason, Gregory C.; Manginell, Ronald P.
2008-07-15
A microfabricated diffusion source to provide for a controlled diffusion rate of a vapor comprises a porous reservoir formed in a substrate that can be filled with a liquid, a headspace cavity for evaporation of the vapor therein, a diffusion channel to provide a controlled diffusion of the vapor, and an outlet to release the vapor into a gas stream. The microfabricated diffusion source can provide a calibration standard for a microanalytical system. The microanalytical system with an integral diffusion source can be fabricated with microelectromechanical systems technologies.
'LTE-diffusion approximation' for arc calculations
NASA Astrophysics Data System (ADS)
Lowke, J. J.; Tanaka, M.
2006-08-01
This paper proposes the use of the 'LTE-diffusion approximation' for predicting the properties of electric arcs. Under this approximation, local thermodynamic equilibrium (LTE) is assumed, with a particular mesh size near the electrodes chosen to be equal to the 'diffusion length', based on De/W, where De is the electron diffusion coefficient and W is the electron drift velocity. This approximation overcomes the problem that the equilibrium electrical conductivity in the arc near the electrodes is almost zero, which makes accurate calculations using LTE impossible in the limit of small mesh size, as then voltages would tend towards infinity. Use of the LTE-diffusion approximation for a 200 A arc with a thermionic cathode gives predictions of total arc voltage, electrode temperatures, arc temperatures and radial profiles of heat flux density and current density at the anode that are in approximate agreement with more accurate calculations which include an account of the diffusion of electric charges to the electrodes, and also with experimental results. Calculations, which include diffusion of charges, agree with experimental results of current and heat flux density as a function of radius if the Milne boundary condition is used at the anode surface rather than imposing zero charge density at the anode.
Diffusion approximation for modeling of 3-D radiation distributions
Zardecki, A.; Gerstl, S.A.W.; De Kinder, R.E. Jr.
1985-01-01
A three-dimensional transport code DIF3D, based on the diffusion approximation, is used to model the spatial distribution of radiation energy arising from volumetric isotropic sources. Future work will be concerned with the determination of irradiances and modeling of realistic scenarios, relevant to the battlefield conditions. 8 refs., 4 figs.
Extending the diffusion approximation to the boundary using an integrated diffusion model
NASA Astrophysics Data System (ADS)
Chen, Chen; Du, Zhidong; Pan, Liang
2015-06-01
The widely used diffusion approximation is inaccurate to describe the transport behaviors near surfaces and interfaces. To solve such stochastic processes, an integro-differential equation, such as the Boltzmann transport equation (BTE), is typically required. In this work, we show that it is possible to keep the simplicity of the diffusion approximation by introducing a nonlocal source term and a spatially varying diffusion coefficient. We apply the proposed integrated diffusion model (IDM) to a benchmark problem of heat conduction across a thin film to demonstrate its feasibility. We also validate the model when boundary reflections and uniform internal heat generation are present.
A diffuse-interface approximation for surface diffusion including adatoms
NASA Astrophysics Data System (ADS)
Rätz, Andreas; Voigt, A.
2007-01-01
We introduce a diffuse-interface approximation for solving partial differential equations on evolving surfaces. The model of interest is a fourth-order geometric evolution equation for a growing surface with an additional diffusive adatom density on the surface. Such models arise in the description of epitaxial growth, where the surface of interest is the solid-vapour interface. The model allows us to handle complex geometries in an implicit manner, by considering an evolution equation for a phase-field variable describing the surface and an evolution equation for an extended adatom concentration on a time-independent domain. Matched asymptotic analysis shows the formal convergence towards the sharp interface model and numerical results based on adaptive finite elements demonstrate the applicability of the approach.
An asymptotic homogenized neutron diffusion approximation. II. Numerical comparisons
Trahan, T. J.; Larsen, E. W.
2012-07-01
In a companion paper, a monoenergetic, homogenized, anisotropic diffusion equation is derived asymptotically for large, 3-D, multiplying systems with a periodic lattice structure [1]. In the present paper, this approximation is briefly compared to several other well known diffusion approximations. Although the derivation is different, the asymptotic diffusion approximation matches that proposed by Deniz and Gelbard, and is closely related to those proposed by Benoist. The focus of this paper, however, is a numerical comparison of the various methods for simple reactor analysis problems in 1-D. The comparisons show that the asymptotic diffusion approximation provides a more accurate estimate of the eigenvalue than the Benoist diffusion approximations. However, the Benoist diffusion approximations and the asymptotic diffusion approximation provide very similar estimates of the neutron flux. The asymptotic method and the Benoist methods both outperform the standard homogenized diffusion approximation, with flux weighted cross sections. (authors)
Multivariate Padé Approximations For Solving Nonlinear Diffusion Equations
NASA Astrophysics Data System (ADS)
Turut, V.
2015-11-01
In this paper, multivariate Padé approximation is applied to power series solutions of nonlinear diffusion equations. As it is seen from tables, multivariate Padé approximation (MPA) gives reliable solutions and numerical results.
Capturing correlations in chaotic diffusion by approximation methods.
Knight, Georgie; Klages, Rainer
2011-10-01
We investigate three different methods for systematically approximating the diffusion coefficient of a deterministic random walk on the line that contains dynamical correlations that change irregularly under parameter variation. Capturing these correlations by incorporating higher-order terms, all schemes converge to the analytically exact result. Two of these methods are based on expanding the Taylor-Green-Kubo formula for diffusion, while the third method approximates Markov partitions and transition matrices by using a slight variation of the escape rate theory of chaotic diffusion. We check the practicability of the different methods by working them out analytically and numerically for a simple one-dimensional map, study their convergence, and critically discuss their usefulness in identifying a possible fractal instability of parameter-dependent diffusion, in the case of dynamics where exact results for the diffusion coefficient are not available. PMID:22181115
Integral approximations to classical diffusion and smoothed particle hydrodynamics
Du, Qiang; Lehoucq, R. B.; Tartakovsky, A. M.
2014-12-31
The contribution of the paper is the approximation of a classical diffusion operator by an integral equation with a volume constraint. A particular focus is on classical diffusion problems associated with Neumann boundary conditions. By exploiting this approximation, we can also approximate other quantities such as the flux out of a domain. Our analysis of the model equation on the continuum level is closely related to the recent work on nonlocal diffusion and peridynamic mechanics. In particular, we elucidate the role of a volumetric constraint as an approximation to a classical Neumann boundary condition in the presence of physical boundary. The volume-constrained integral equation then provides the basis for accurate and robust discretization methods. As a result, an immediate application is to the understanding and improvement of the Smoothed Particle Hydrodynamics (SPH) method.
Integral approximations to classical diffusion and smoothed particle hydrodynamics
Du, Qiang; Lehoucq, R. B.; Tartakovsky, A. M.
2014-12-31
The contribution of the paper is the approximation of a classical diffusion operator by an integral equation with a volume constraint. A particular focus is on classical diffusion problems associated with Neumann boundary conditions. By exploiting this approximation, we can also approximate other quantities such as the flux out of a domain. Our analysis of the model equation on the continuum level is closely related to the recent work on nonlocal diffusion and peridynamic mechanics. In particular, we elucidate the role of a volumetric constraint as an approximation to a classical Neumann boundary condition in the presence of physical boundary.more » The volume-constrained integral equation then provides the basis for accurate and robust discretization methods. As a result, an immediate application is to the understanding and improvement of the Smoothed Particle Hydrodynamics (SPH) method.« less
Integral approximations to classical diffusion and smoothed particle hydrodynamics
Du, Q.; Lehoucq, Richard B.; Tartakovsky, Alexandre M.
2015-04-01
The contribution of the paper is the approximation of a classical diffusion operator by an integral equation with a volume constraint. A particular focus is on classical diffusion problems associated with Neumann boundary conditions. By exploiting this approximation, we can also approximate other quantities such as the flux out of a domain. Our analysis of the model equation on the continuum level is closely related to the recent work on nonlocal diffusion and peridynamic mechanics. In particular, we elucidate the role of a volumetric constraint as an approximation to a classical Neumann boundary condition in the presence of physical boundary. The volume-constrained integral equation then provides the basis for accurate and robust discretization methods. An immediate application is to the understanding and improvement of the Smoothed Particle Hydrodynamics (SPH) method.
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
Accuracy of the non-relativistic approximation for momentum diffusion
NASA Astrophysics Data System (ADS)
Liang, Shiuan-Ni; Lan, Boon Leong
2016-06-01
The accuracy of the non-relativistic approximation, which is calculated using the same parameter and the same initial ensemble of trajectories, to relativistic momentum diffusion at low speed is studied numerically for a prototypical nonlinear Hamiltonian system -the periodically delta-kicked particle. We find that if the initial ensemble is a non-localized semi-uniform ensemble, the non-relativistic approximation to the relativistic mean square momentum displacement is always accurate. However, if the initial ensemble is a localized Gaussian, the non-relativistic approximation may not always be accurate and the approximation can break down rapidly.
NASA Technical Reports Server (NTRS)
Burns, R. E.
1973-01-01
The problem with predicting pollutant diffusion from a line source of arbitrary geometry is treated. The concentration at the line source may be arbitrarily varied with time. Special attention is given to the meteorological inputs which act as boundary conditions for the problem, and a mixing layer of arbitrary depth is assumed. Numerical application of the derived theory indicates the combinations of meteorological parameters that may be expected to result in high pollution concentrations.
NASA Astrophysics Data System (ADS)
Yuan, Zhen; Zhang, Qizhi; Sobel, Eric; Jiang, Huabei
2009-09-01
In this study, a simplified spherical harmonics approximated higher order diffusion model is employed for 3-D diffuse optical tomography of osteoarthritis in the finger joints. We find that the use of a higher-order diffusion model in a stand-alone framework provides significant improvement in reconstruction accuracy over the diffusion approximation model. However, we also find that this is not the case in the image-guided setting when spatial prior knowledge from x-rays is incorporated. The results show that the reconstruction error between these two models is about 15 and 4%, respectively, for stand-alone and image-guided frameworks.
The diffusion approximation. An application to radiative transfer in clouds
NASA Technical Reports Server (NTRS)
Arduini, R. F.; Barkstrom, B. R.
1976-01-01
It is shown how the radiative transfer equation reduces to the diffusion equation. To keep the mathematics as simple as possible, the approximation is applied to a cylindrical cloud of radius R and height h. The diffusion equation separates in cylindrical coordinates and, in a sample calculation, the solution is evaluated for a range of cloud radii with cloud heights of 0.5 km and 1.0 km. The simplicity of the method and the speed with which solutions are obtained give it potential as a tool with which to study the effects of finite-sized clouds on the albedo of the earth-atmosphere system.
Lazarov, R D; Vassilevski, P S
1999-05-06
In this paper we introduce and study a least-squares finite element approximation for singularly perturbed convection-diffusion equations of second order. By introducing the flux (diffusive plus convective) as a new unknown, the problem is written in a mixed form as a first order system. Further, the flux is augmented by adding the lower order terms with a small parameter. The new first order system is approximated by the least-squares finite element method using the minus one norm approach of Bramble, Lazarov, and Pasciak [2]. Further, we estimate the error of the method and discuss its implementation and the numerical solution of some test problems.
The Approximate Number System Acuity Redefined: A Diffusion Model Approach
Park, Joonkoo; Starns, Jeffrey J.
2015-01-01
While all humans are capable of non-verbally representing numerical quantity using so-called the approximate number system (ANS), there exist considerable individual differences in its acuity. For example, in a non-symbolic number comparison task, some people find it easy to discriminate brief presentations of 14 dots from 16 dots while others do not. Quantifying individual ANS acuity from such a task has become an essential practice in the field, as individual differences in such a primitive number sense is thought to provide insights into individual differences in learned symbolic math abilities. However, the dominant method of characterizing ANS acuity—computing the Weber fraction (w)—only utilizes the accuracy data while ignoring response times (RT). Here, we offer a novel approach of quantifying ANS acuity by using the diffusion model, which accounts both accuracy and RT distributions. Specifically, the drift rate in the diffusion model, which indexes the quality of the stimulus information, is used to capture the precision of the internal quantity representation. Analysis of behavioral data shows that w is contaminated by speed-accuracy tradeoff, making it problematic as a measure of ANS acuity, while drift rate provides a measure more independent from speed-accuracy criterion settings. Furthermore, drift rate is a better predictor of symbolic math ability than w, suggesting a practical utility of the measure. These findings demonstrate critical limitations of the use of w and suggest clear advantages of using drift rate as a measure of primitive numerical competence. PMID:26733929
Extended source model for diffusive coupling.
González-Ochoa, Héctor O; Flores-Moreno, Roberto; Reyes, Luz M; Femat, Ricardo
2016-01-01
Motivated by the prevailing approach to diffusion coupling phenomena which considers point-like diffusing sources, we derived an analogous expression for the concentration rate of change of diffusively coupled extended containers. The proposed equation, together with expressions based on solutions to the diffusion equation, is intended to be applied to the numerical solution of systems exclusively composed of ordinary differential equations, however is able to account for effects due the finite size of the coupled sources. PMID:26802012
NASA Astrophysics Data System (ADS)
Yi, Xi; Wang, Bingyuan; Wan, Wenbo; Wang, Yihan; Zhang, Yanqi; Zhao, Huijuan; Gao, Feng
2015-05-01
Of the three measurement schemes established for diffuse fluorescence tomography (DFT), the time-domain scheme is well known to provide the richest information about the distribution of the targeting fluorophore in living tissues. However, the explicit use of the full time-resolved data usually leads to a considerably lengthy time for image reconstruction, limiting its applications to three-dimensional or small-volume imaging. To cope with the adversity, we propose herein a computationally efficient scheme for DFT image reconstruction where the time-dependent photon density is expanded to a Fourier-series and calculated by solving the independent frequency-domain diffusion equations at multiple sampling frequencies with the support of a combined multicore CPU-based coarse-grain and multithread GPU-based fine-grain parallelization strategy. With such a parallelized Fourier-series truncated diffusion approximation, both the time- and frequency-domain inversion procedures are developed and validated for their effectiveness and accuracy using simulative and phantom experiments. The results show that the proposed method can generate reconstructions comparable to the explicit time-domain scheme, with significantly reduced computational time.
Berkel, M. van; Tamura, N.; Ida, K.; Hogeweij, G. M. D.; Zwart, H. J.; Inagaki, S.; Baar, M. R. de
2014-11-15
In this paper, a number of new explicit approximations are introduced to estimate the perturbative diffusivity (χ), convectivity (V), and damping (τ) in cylindrical geometry. For this purpose, the harmonic components of heat waves induced by localized deposition of modulated power are used. The approximations are based on the heat equation in cylindrical geometry using the symmetry (Neumann) boundary condition at the plasma center. This means that the approximations derived here should be used only to estimate transport coefficients between the plasma center and the off-axis perturbative source. If the effect of cylindrical geometry is small, it is also possible to use semi-infinite domain approximations presented in Part I and Part II of this series. A number of new approximations are derived in this part, Part III, based upon continued fractions of the modified Bessel function of the first kind and the confluent hypergeometric function of the first kind. These approximations together with the approximations based on semi-infinite domains are compared for heat waves traveling towards the center. The relative error for the different derived approximations is presented for different values of the frequency, transport coefficients, and dimensionless radius. Moreover, it is shown how combinations of different explicit formulas can be used to estimate the transport coefficients over a large parameter range for cases without convection and damping, cases with damping only, and cases with convection and damping. The relative error between the approximation and its underlying model is below 2% for the case, where only diffusivity and damping are considered. If also convectivity is considered, the diffusivity can be estimated well in a large region, but there is also a large region in which no suitable approximation is found. This paper is the third part (Part III) of a series of three papers. In Part I, the semi-infinite slab approximations have been treated. In Part II
Diffusion approximation in turbulent two-particle dispersion
NASA Astrophysics Data System (ADS)
Eyink, Gregory L.; Benveniste, Damien
2013-10-01
We solve an inverse problem for fluid particle pair statistics: we show that a time sequence of probability density functions (PDFs) of separations can be exactly reproduced by solving the diffusion equation with a suitable time-dependent diffusivity. The diffusivity tensor is given by a time integral of a conditional Lagrangian velocity structure function, weighted by a ratio of PDFs. Physical hypotheses for hydrodynamic turbulence (sweeping, short memory, mean-field) yield simpler integral formulas, including one of Kraichnan and Lundgren (K-L). We evaluate the latter using a space-time database from a numerical Navier-Stokes solution for driven turbulence. The K-L formula reproduces PDFs well at root-mean-square separations, but growth rate of mean-square dispersion is overpredicted due to neglect of memory effects. More general applications of our approach are sketched.
An asymptotic homogenized neutron diffusion approximation. I. Theory
Trahan, T. J.; Larsen, E. W.
2012-07-01
A monoenergetic, homogenized, anisotropic diffusion equation is derived asymptotically for large, 3-D, multiplying systems with a periodic lattice structure. The primary assumption is that the system is slightly perturbed from an infinite, periodic lattice, and that the length scale of a lattice element is small relative to the total system size. The perturbed flux is slightly buckled, and the leading order term is the product of a slowly varying amplitude component, and a rapidly varying periodic component. The amplitude function is the solution to the homogenized diffusion equation, while the periodic component is the solution to the unperturbed, infinite system, and can be found using any high-order transport method. The first order term acts as a correction term, and makes it possible to obtain a zero flux extrapolation distance for the diffusion equation by applying the Marshak boundary condition. (authors)
Multiple Scattering in Random Mechanical Systems and Diffusion Approximation
NASA Astrophysics Data System (ADS)
Feres, Renato; Ng, Jasmine; Zhang, Hong-Kun
2013-10-01
This paper is concerned with stochastic processes that model multiple (or iterated) scattering in classical mechanical systems of billiard type, defined below. From a given (deterministic) system of billiard type, a random process with transition probabilities operator P is introduced by assuming that some of the dynamical variables are random with prescribed probability distributions. Of particular interest are systems with weak scattering, which are associated to parametric families of operators P h , depending on a geometric or mechanical parameter h, that approaches the identity as h goes to 0. It is shown that ( P h - I)/ h converges for small h to a second order elliptic differential operator on compactly supported functions and that the Markov chain process associated to P h converges to a diffusion with infinitesimal generator . Both P h and are self-adjoint (densely) defined on the space of square-integrable functions over the (lower) half-space in , where η is a stationary measure. This measure's density is either (post-collision) Maxwell-Boltzmann distribution or Knudsen cosine law, and the random processes with infinitesimal generator respectively correspond to what we call MB diffusion and (generalized) Legendre diffusion. Concrete examples of simple mechanical systems are given and illustrated by numerically simulating the random processes.
Linear method of fluorescent source reconstruction in a diffusion medium.
Janunts, Edgar; Pöschinger, Thomas; Brünner, Holger; Langenbucher, Achim
2008-01-01
A new method is described for obtaining a 2D reconstruction of a fluorescent source distribution inside a diffusion medium from planar measurements of the emission light at the surface after excitation by a plane wave. Point sources are implanted at known locations of a rectangular phantom. The forward model of the photon transport is based on the diffusion approximation of the radiative transport equation (RTE) for homogeneous media. This can be described by a hierarchical system of two time-independent RTE's, one for the excitation plane wave originating from the external light source to the medium and another one for the fluorescence emission originating from the fluorophore marker to the detector. A linear inverse source problem was solved for image reconstruction. The applicability of the theoretical method is demonstrated in some representative working examples. For an optimization of the problem we used least squares minimization technique. PMID:18826162
Berkel, M. van; Hogeweij, G. M. D.; Tamura, N.; Ida, K.; Zwart, H. J.; Inagaki, S.; Baar, M. R. de
2014-11-15
In this paper, a number of new explicit approximations are introduced to estimate the perturbative diffusivity (χ), convectivity (V), and damping (τ) in a cylindrical geometry. For this purpose, the harmonic components of heat waves induced by localized deposition of modulated power are used. The approximations are based upon the heat equation in a semi-infinite cylindrical domain. The approximations are based upon continued fractions, asymptotic expansions, and multiple harmonics. The relative error for the different derived approximations is presented for different values of frequency, transport coefficients, and dimensionless radius. Moreover, it is shown how combinations of different explicit formulas can yield good approximations over a wide parameter space for different cases, such as no convection and damping, only damping, and both convection and damping. This paper is the second part (Part II) of a series of three papers. In Part I, the semi-infinite slab approximations have been treated. In Part III, cylindrical approximations are treated for heat waves traveling towards the center of the plasma.
Berkel, M. van; Zwart, H. J.; Tamura, N.; Ida, K.; Hogeweij, G. M. D.; Inagaki, S.; Baar, M. R. de
2014-11-15
In this paper, a number of new approximations are introduced to estimate the perturbative diffusivity (χ), convectivity (V), and damping (τ) in cylindrical geometry. For this purpose, the harmonic components of heat waves induced by localized deposition of modulated power are used. The approximations are based on semi-infinite slab approximations of the heat equation. The main result is the approximation of χ under the influence of V and τ based on the phase of two harmonics making the estimate less sensitive to calibration errors. To understand why the slab approximations can estimate χ well in cylindrical geometry, the relationships between heat transport models in slab and cylindrical geometry are studied. In addition, the relationship between amplitude and phase with respect to their derivatives, used to estimate χ, is discussed. The results are presented in terms of the relative error for the different derived approximations for different values of frequency, transport coefficients, and dimensionless radius. The approximations show a significant region in which χ, V, and τ can be estimated well, but also regions in which the error is large. Also, it is shown that some compensation is necessary to estimate V and τ in a cylindrical geometry. On the other hand, errors resulting from the simplified assumptions are also discussed showing that estimating realistic values for V and τ based on infinite domains will be difficult in practice. This paper is the first part (Part I) of a series of three papers. In Part II and Part III, cylindrical approximations based directly on semi-infinite cylindrical domain (outward propagating heat pulses) and inward propagating heat pulses in a cylindrical domain, respectively, will be treated.
Gill, Michelle L; Palmer, Arthur G
2014-09-25
The present work demonstrates that NMR spin relaxation rate constants for molecules interconverting between states with different diffusion tensors can be modeled theoretically by combining orientational correlation functions for exchanging spherical molecules with locally isotropic approximations for the diffusion anisotropic tensors. The resulting expressions are validated by comparison with correlation functions obtained by Monte Carlo simulations and are accurate for moderate degrees of diffusion anisotropy typically encountered in investigations of globular proteins. The results are complementary to an elegant, but more complex, formalism that is accurate for all degrees of diffusion anisotropy [Ryabov, Y.; Clore, G. M.; Schwieters, C. D. J. Chem. Phys. 2012, 136, 034108]. PMID:25167331
NASA Astrophysics Data System (ADS)
Horowitz, Jordan M.
2015-07-01
The stochastic thermodynamics of a dilute, well-stirred mixture of chemically reacting species is built on the stochastic trajectories of reaction events obtained from the chemical master equation. However, when the molecular populations are large, the discrete chemical master equation can be approximated with a continuous diffusion process, like the chemical Langevin equation or low noise approximation. In this paper, we investigate to what extent these diffusion approximations inherit the stochastic thermodynamics of the chemical master equation. We find that a stochastic-thermodynamic description is only valid at a detailed-balanced, equilibrium steady state. Away from equilibrium, where there is no consistent stochastic thermodynamics, we show that one can still use the diffusive solutions to approximate the underlying thermodynamics of the chemical master equation.
Horowitz, Jordan M.
2015-07-28
The stochastic thermodynamics of a dilute, well-stirred mixture of chemically reacting species is built on the stochastic trajectories of reaction events obtained from the chemical master equation. However, when the molecular populations are large, the discrete chemical master equation can be approximated with a continuous diffusion process, like the chemical Langevin equation or low noise approximation. In this paper, we investigate to what extent these diffusion approximations inherit the stochastic thermodynamics of the chemical master equation. We find that a stochastic-thermodynamic description is only valid at a detailed-balanced, equilibrium steady state. Away from equilibrium, where there is no consistent stochastic thermodynamics, we show that one can still use the diffusive solutions to approximate the underlying thermodynamics of the chemical master equation.
S{sub N} Schemes, Linear Infinite-Medium Solutions, and the Diffusion Approximation
Larsen, E.W.
2001-06-17
It is standard practice to require an S{sub N} spatial discretization scheme to preserve the ''flat infinite-medium'' solution of the transport equation. This solution consists of a spatially independent source that gives rise to a spatially independent flux. However, there exist many other exact solutions of the transport equation that are typically not preserved by approximation schemes. Here, we discuss one of these: a source that is linear in space giving rise to an angular flux that is linear in space and angle. For one-group, planar-geometry S{sub N} problems, we show that (a) among the class of weighted-diamond schemes, only one - the diamond-difference scheme - preserves this exact ''linear'' solution; (b) consequently, only the diamond scheme preserves the correct Fick's Law; and (c) as a further consequence, nondiamond schemes can produce significant errors (not observed in the diamond solution) for diffusive problems with spatial cells that are not optically thin. These results demonstrate that it is advantageous for S{sub N} discretization schemes to preserve the ''flat'' and ''linear'' infinite-medium solutions.
The Heavy-Ion Approximation for Ambipolar Diffusion Calcuations for Weakly Ionized Plasmas
Li, P; McKee, C; Klein, R
2006-07-27
Ambipolar diffusion redistributes magnetic flux in weakly ionized plasmas and plays a critical role in star formation. Simulations of ambipolar diffusion using explicit MHD codes are prohibitively expensive for the level of ionization observed in molecular clouds ({approx}< 10{sup -6}) since an enormous number of time steps is required to represent the dynamics of the dominant neutral component with a time step determined by the trace ion component. Here we show that ambipolar diffusion calculations can be significantly accelerated by the 'heavy-ion approximation', in which the mass density of the ions is increased and the collisional coupling constant with the neutrals decreased such that the product remains constant. In this approximation, the ambipolar diffusion time and the ambipolar magnetic Reynolds number remain unchanged. We present three tests of the heavy-ion approximation: C-type shocks, the Wardle instability, and the 1D collapse of a magnetized slab. We show that this approximation is quite accurate provided that (1) the square of the Alfven Mach number is small compared to the ambipolar diffusion Reynolds number for dynamical problems, and that (2) the ion mass density is negligible for quasi-static problems; a specific criterion is given for the magnetized slab problem. The first condition can be very stringent for turbulent flows with large density fluctuations.
Gorpas, Dimitris; Andersson-Engels, Stefan
2012-12-01
The solution of the forward problem in fluorescence molecular imaging strongly influences the successful convergence of the fluorophore reconstruction. The most common approach to meeting this problem has been to apply the diffusion approximation. However, this model is a first-order angular approximation of the radiative transfer equation, and thus is subject to some well-known limitations. This manuscript proposes a methodology that confronts these limitations by applying the radiative transfer equation in spatial regions in which the diffusion approximation gives decreased accuracy. The explicit integro differential equations that formulate this model were solved by applying the Galerkin finite element approximation. The required spatial discretization of the investigated domain was implemented through the Delaunay triangulation, while the azimuthal discretization scheme was used for the angular space. This model has been evaluated on two simulation geometries and the results were compared with results from an independent Monte Carlo method and the radiative transfer equation by calculating the absolute values of the relative errors between these models. The results show that the proposed forward solver can approximate the radiative transfer equation and the Monte Carlo method with better than 95% accuracy, while the accuracy of the diffusion approximation is approximately 10% lower. PMID:23208221
Test particle propagation in magnetostatic turbulence. 1. Failure of the diffusion approximation
NASA Technical Reports Server (NTRS)
Klimas, A. J.; Sandri, G.; Scudder, J. D.; Howell, D. R.
1976-01-01
The equation which governs the quasi-linear approximation to the ensemble and gyro-phase averaged one-body probability distribution function is constructed from first principles. This derived equation is subjected to a thorough investigation in order to calculate the possible limitations of the quasi-linear approximation. It is shown that the reduction of this equation to a standard diffusion equation in the Markovian limit can be accomplished through the application of the adiabatic approximation. A numerical solution of the standard diffusion equation in the Markovian limit is obtained for the narrow parallel beam injection. Comparison of the diabatic and adiabatic results explicitly demonstrates the failure of the Markovian description of the probability distribution function. Through the use of a linear time-scale extension the failure of the adiabatic approximation, which leads to the Markovian limit, is shown to be due to mixing of the relaxation and interaction time scales in the presence of the strong mean field.
NASA Astrophysics Data System (ADS)
Sweilam, N. H.; Abou Hasan, M. M.
2016-08-01
This paper reports a new spectral algorithm for obtaining an approximate solution for the Lévy-Feller diffusion equation depending on Legendre polynomials and Chebyshev collocation points. The Lévy-Feller diffusion equation is obtained from the standard diffusion equation by replacing the second-order space derivative with a Riesz-Feller derivative. A new formula expressing explicitly any fractional-order derivatives, in the sense of Riesz-Feller operator, of Legendre polynomials of any degree in terms of Jacobi polynomials is proved. Moreover, the Chebyshev-Legendre collocation method together with the implicit Euler method are used to reduce these types of differential equations to a system of algebraic equations which can be solved numerically. Numerical results with comparisons are given to confirm the reliability of the proposed method for the Lévy-Feller diffusion equation.
NASA Astrophysics Data System (ADS)
Fontaine, G.; Brassard, P.; Dufour, P.; Tremblay, P.-E.
2015-06-01
The accretion-diffusion picture is the model par excellence for describing the presence of planetary debris polluting the atmospheres of relatively cool white dwarfs. Some important insights into the process may be derived using an approximate approach which combines static stellar models with estimates of diffusion timescales at the base of the outer convection zone or, in its absence, at the photosphere. Until recently, and to our knowledge, values of diffusion timescales in white dwarfs have all been obtained on the basis of the same physics as that developed initially by Paquette et al., including their diffusion coefficients and thermal diffusion coefficients. In view of the recent exciting discoveries of a plethora of metals (including some never seen before) polluting the atmospheres of an increasing number of cool white dwarfs, we felt that a new look at the estimates of settling timescales would be worthwhile. We thus provide improved estimates of diffusion timescales for all 27 elements from Li to Cu in the periodic table in a wide range of the surface gravity-effective temperature domain and for both DA and non-DA stars.
Barth, Andrea Lang, Annika
2012-12-15
In this paper, the strong approximation of a stochastic partial differential equation, whose differential operator is of advection-diffusion type and which is driven by a multiplicative, infinite dimensional, cadlag, square integrable martingale, is presented. A finite dimensional projection of the infinite dimensional equation, for example a Galerkin projection, with nonequidistant time stepping is used. Error estimates for the discretized equation are derived in L{sup 2} and almost sure senses. Besides space and time discretizations, noise approximations are also provided, where the Milstein double stochastic integral is approximated in such a way that the overall complexity is not increased compared to an Euler-Maruyama approximation. Finally, simulations complete the paper.
Shirdel-Havar, A. H. Masoudian Saadabad, R.
2015-03-21
Based on ballistic-diffusive approximation, a method is presented to model heat transfer in nanocomposites containing metal nanoparticles. This method provides analytical expression for the temperature dynamics of metallic nanoparticles embedded in a dielectric medium. In this study, nanoparticles are considered as spherical shells, so that Boltzmann equation is solved using ballistic-diffusive approximation to calculate the electron and lattice thermal dynamics in gold nanoparticles, while thermal exchange between the particles is taken into account. The model was used to investigate the influence of particle size and metal concentration of the medium on the electron and lattice thermal dynamics. It is shown that these two parameters are crucial in determining the nanocomposite thermal behavior. Our results showed that the heat transfer rate from nanoparticles to the matrix decreases as the nanoparticle size increases. On the other hand, increasing the metal concentration of the medium can also decrease the heat transfer rate.
The effects of hyper-spherical approximation of Yukawa potential to diffusion properties
NASA Astrophysics Data System (ADS)
Kim, In Gee; Murillo, Michael S.
2015-11-01
The effects of Yukawa potential to the diffusion properties of binary ionic mixtures are investigated in terms of both the classical molecular dynamics and the kinetic theory. The Yukawa interatomic potential is treated by means of the hyper-spherical approximation, which replaces the Ewald summation by a multiple of the hyperbolic trigonometric functions and the lattice summation of screening. The influence of the hyper-spherical approximation of the Yukawa potential is able to be understood through the calculations of transport coefficients with the relationship to Coulomb logarithm. Numerical studies over a various range of the Debye-Hückel screening parameter and of the plasma coupling parameter to binary ionic mixtures will be provided. We consider primarily the interdiffusion coefficients and then discuss about the mixing properties of self-diffusion coefficients.
Spectral approximation to advection-diffusion problems by the fictitious interface method
NASA Astrophysics Data System (ADS)
Frati, A.; Pasquarelli, F.; Quarteroni, A.
1993-08-01
The algorithmic aspects of the 'fictitious interface' method used in numerical approximations of convection-dominated flows are discussed. The solution algorithm presented alternates the advection-equation solution with that of the advection-diffusion equation within complementary subdomains. For the problems presently considered, spatial discretization is obtained by the spectral collocation method via Legendre-Gaussian modes. Attention is given to the the fictitious interface method's application to the Burgers equation.
Approximations of the operator exponential in a periodic diffusion problem with drift
Pastukhova, Svetlana E
2013-02-28
A Cauchy problem for a parabolic diffusion equation with 1-periodic coefficients containing first order terms is studied. For the corresponding semigroup we construct approximations in the L{sup 2}-operator norm on sections t=const of order O(t{sup -m/2}) as t{yields}{infinity} for m=1 or m=2. The spectral method based on the Bloch representation of an operator with periodic coefficients is used. Bibliography: 25 titles.
Thygesen, Uffe Høgsbro
2016-03-01
We consider organisms which use a renewal strategy such as run-tumble when moving in space, for example to perform chemotaxis in chemical gradients. We derive a diffusion approximation for the motion, applying a central limit theorem due to Anscombe for renewal-reward processes; this theorem has not previously been applied in this context. Our results extend previous work, which has established the mean drift but not the diffusivity. For a classical model of tumble rates applied to chemotaxis, we find that the resulting chemotactic drift saturates to the swimming velocity of the organism when the chemical gradients grow increasingly steep. The dispersal becomes anisotropic in steep gradients, with larger dispersal across the gradient than along the gradient. In contrast to one-dimensional settings, strong bias increases dispersal. We next include Brownian rotation in the model and find that, in limit of high chemotactic sensitivity, the chemotactic drift is 64% of the swimming velocity, independent of the magnitude of the Brownian rotation. We finally derive characteristic timescales of the motion that can be used to assess whether the diffusion limit is justified in a given situation. The proposed technique for obtaining diffusion approximations is conceptually and computationally simple, and applicable also when statistics of the motion is obtained empirically or through Monte Carlo simulation of the motion. PMID:27012850
Neutron diffusion in a randomly inhomogeneous multiplying medium with random phase approximation
Imre, Kaya; Akcasu, A. Ziya
2012-06-15
Neutron diffusion in a randomly inhomogeneous multiplying medium is studied. By making use of a random phase assumption we show that the average neutron density approximately satisfies an integral equation in Fourier space, which is solved using Kummer functions. We used multi-dimensional formulation. In the case of one dimension, we obtain the result of Rosenbluth and Tao for the mean total density for large t. In the three-dimensional case, a closed form of solution is derived for the mean total neutron density. Its asymptotic behavior is also investigated for large t.
Neutron diffusion in a randomly inhomogeneous multiplying medium with random phase approximation
NASA Astrophysics Data System (ADS)
Imre, Kaya; Akcasu, A. Ziya
2012-06-01
Neutron diffusion in a randomly inhomogeneous multiplying medium is studied. By making use of a random phase assumption we show that the average neutron density approximately satisfies an integral equation in Fourier space, which is solved using Kummer functions. We used multi-dimensional formulation. In the case of one dimension, we obtain the result of Rosenbluth and Tao for the mean total density for large t. In the three-dimensional case, a closed form of solution is derived for the mean total neutron density. Its asymptotic behavior is also investigated for large t.
A Diffusion Approximation and Numerical Methods for Adaptive Neuron Models with Stochastic Inputs.
Rosenbaum, Robert
2016-01-01
Characterizing the spiking statistics of neurons receiving noisy synaptic input is a central problem in computational neuroscience. Monte Carlo approaches to this problem are computationally expensive and often fail to provide mechanistic insight. Thus, the field has seen the development of mathematical and numerical approaches, often relying on a Fokker-Planck formalism. These approaches force a compromise between biological realism, accuracy and computational efficiency. In this article we develop an extension of existing diffusion approximations to more accurately approximate the response of neurons with adaptation currents and noisy synaptic currents. The implementation refines existing numerical schemes for solving the associated Fokker-Planck equations to improve computationally efficiency and accuracy. Computer code implementing the developed algorithms is made available to the public. PMID:27148036
A Diffusion Approximation and Numerical Methods for Adaptive Neuron Models with Stochastic Inputs
Rosenbaum, Robert
2016-01-01
Characterizing the spiking statistics of neurons receiving noisy synaptic input is a central problem in computational neuroscience. Monte Carlo approaches to this problem are computationally expensive and often fail to provide mechanistic insight. Thus, the field has seen the development of mathematical and numerical approaches, often relying on a Fokker-Planck formalism. These approaches force a compromise between biological realism, accuracy and computational efficiency. In this article we develop an extension of existing diffusion approximations to more accurately approximate the response of neurons with adaptation currents and noisy synaptic currents. The implementation refines existing numerical schemes for solving the associated Fokker-Planck equations to improve computationally efficiency and accuracy. Computer code implementing the developed algorithms is made available to the public. PMID:27148036
Dworak, D; Loskiewicz, J; Janik, M
2001-05-01
The diffusion approximation solution for neutron transport has been used in well-logging geophysics for calculating tool responses in boreholes, sometimes with success. The problem of the dimension of different materials to which it can be applied with success is important for the borehole environment. The results obtained show that the diffusion approximation can be used for distances greater than a few millimetre in some rock types. For iron, barium, and other highly absorbing media the use of the diffusion approximation is inappropriate even for large distances. PMID:11258535
Approximating edges of source bodies from magnetic or gravity anomalies.
Blakely, R.J.; Simpson, R.W.
1986-01-01
Cordell and Grauch (1982, 1985) discussed a technique to estimate the location of abrupt lateral changes in magnetization or mass density of upper crustal rocks. The final step of their procedure is to identify maxima on a contoured map of horizontal gradient magnitudes. Attempts to automate their final step. The method begins with gridded magnetic or gravity anomaly data and produces a plan view of inferred boundaries of magnetic or gravity sources. The method applies to both local surveys and to continent-wide compilations of magnetic and gravity data.-from Authors
NASA Astrophysics Data System (ADS)
Sarwar, S.; Rashidi, M. M.
2016-07-01
This paper deals with the investigation of the analytical approximate solutions for two-term fractional-order diffusion, wave-diffusion, and telegraph equations. The fractional derivatives are defined in the Caputo sense, whose orders belong to the intervals [0,1], (1,2), and [1,2], respectively. In this paper, we extended optimal homotopy asymptotic method (OHAM) for two-term fractional-order wave-diffusion equations. Highly approximate solution is obtained in series form using this extended method. Approximate solution obtained by OHAM is compared with the exact solution. It is observed that OHAM is a prevailing and convergent method for the solutions of nonlinear-fractional-order time-dependent partial differential problems. The numerical results rendering that the applied method is explicit, effective, and easy to use, for handling more general fractional-order wave diffusion, diffusion, and telegraph problems.
Diffusion filter eliminates fringe effects of coherent laser light source
NASA Technical Reports Server (NTRS)
Olsasky, M. J.
1970-01-01
Diffusion filter comprised of small particles in colloidal suspension reduces the coherence of a laser beam used as a photographic light source. Interference patterns which obscure details in photographic film are eliminated, the intensity and collimation are moderately affected.
Nonlinear Solver Approaches for the Diffusive Wave Approximation to the Shallow Water Equations
NASA Astrophysics Data System (ADS)
Collier, N.; Knepley, M.
2015-12-01
The diffusive wave approximation to the shallow water equations (DSW) is a doubly-degenerate, nonlinear, parabolic partial differential equation used to model overland flows. Despite its challenges, the DSW equation has been extensively used to model the overland flow component of various integrated surface/subsurface models. The equation's complications become increasingly problematic when ponding occurs, a feature which becomes pervasive when solving on large domains with realistic terrain. In this talk I discuss the various forms and regularizations of the DSW equation and highlight their effect on the solvability of the nonlinear system. In addition to this analysis, I present results of a numerical study which tests the applicability of a class of composable nonlinear algebraic solvers recently added to the Portable, Extensible, Toolkit for Scientific Computation (PETSc).
NASA Astrophysics Data System (ADS)
Kalenkov, Mikhail S.; Zaikin, Andrei D.
2015-07-01
We extend the quasiclassical formalism for diffusive superconductors by deriving anisotropic (gradient) corrections to the Usadel equation. We demonstrate that in a number of physical situations such corrections may play a crucial role, being responsible for the effects which cannot be recovered within the standard Usadel approximation. One of them is the so-called photoelectric effect in superconductors and superconducting-normal (SN) hybrid structures. Provided a superconducting part of the system is irradiated by an external ac electromagnetic field, the charge imbalance develops and a nonvanishing dc voltage is induced across the SN interface. In the presence of magnetic impurities in a superconductor the magnitude of this effect becomes large and can easily be detected in modern experiments.
Svyatskiy, Daniil; Shashkov, Mikhail; Kuzmin, D
2008-01-01
A new approach to the design of constrained finite element approximations to second-order elliptic problems is introduced. This approach guarantees that the finite element solution satisfies the discrete maximum principle (DMP). To enforce these monotonicity constrains the sufficient conditions for elements of the stiffness matrix are formulated. An algebraic splitting of the stiffness matrix is employed to separate the contributions of diffusive and antidiffusive numerical fluxes, respectively. In order to prevent the formation of spurious undershoots and overshoots, a symmetric slope limiter is designed for the antidiffusive part. The corresponding upper and lower bounds are defined using an estimate of the steepest gradient in terms of the maximum and minimum solution values at surrounding nodes. The recovery of nodal gradients is performed by means of a lumped-mass L{sub 2} projection. The proposed slope limiting strategy preserves the consistency of the underlying discrete problem and the structure of the stiffness matrix (symmetry, zero row and column sums). A positivity-preserving defect correction scheme is devised for the nonlinear algebraic system to be solved. Numerical results and a grid convergence study are presented for a number of anisotropic diffusion problems in two space dimensions.
A diffusion approximation model of light transport in multilayered skin tissue
NASA Astrophysics Data System (ADS)
Makropoulou, M.; Kaselouris, E.; Drakaki, E.; Serafetinides, A. A.; Sianoudis, J. A.
2007-07-01
In dermatology, biophotonic methods offer high sensitivity and non-invasive measurements of skin tissue optical properties, in various physiological and pathological conditions. There are numerous skin processes, which can be examined and characterized using diagnostic optical spectroscopy, as the monitoring of skin aging, diagnosis of benign and malignant cutaneous lesions, dosimetry in photodynamic therapy (PDT), etc. Several mathematical models have been used to calculate the tissue optical properties from experimental measurements and to predict the light propagation in soft tissues, like skin, based on transport theory or Monte Carlo modeling. This work analyses the phenomena which are observed experimentally during the irradiation of skin, such as the absorption, reflectance, scattering, fluorescence and transmission of laser light. The study was carried out on animal skin samples, extracted post-mortem. In this work we also tried to evaluate the utility of diffusion approximation modeling for measuring the light intensity distribution in the skin samples with cw visible laser beam (λ=632.8 nm). The diffusion theory model was tested for the simulation results of the spatial light distribution within a five-layer model of animal skin tissue. We have studied the dependence towards the depth and the radial distance of the photon density of the incident radiation.
Linaro, Daniele; Storace, Marco; Giugliano, Michele
2011-01-01
Stochastic channel gating is the major source of intrinsic neuronal noise whose functional consequences at the microcircuit- and network-levels have been only partly explored. A systematic study of this channel noise in large ensembles of biophysically detailed model neurons calls for the availability of fast numerical methods. In fact, exact techniques employ the microscopic simulation of the random opening and closing of individual ion channels, usually based on Markov models, whose computational loads are prohibitive for next generation massive computer models of the brain. In this work, we operatively define a procedure for translating any Markov model describing voltage- or ligand-gated membrane ion-conductances into an effective stochastic version, whose computer simulation is efficient, without compromising accuracy. Our approximation is based on an improved Langevin-like approach, which employs stochastic differential equations and no Montecarlo methods. As opposed to an earlier proposal recently debated in the literature, our approximation reproduces accurately the statistical properties of the exact microscopic simulations, under a variety of conditions, from spontaneous to evoked response features. In addition, our method is not restricted to the Hodgkin-Huxley sodium and potassium currents and is general for a variety of voltage- and ligand-gated ion currents. As a by-product, the analysis of the properties emerging in exact Markov schemes by standard probability calculus enables us for the first time to analytically identify the sources of inaccuracy of the previous proposal, while providing solid ground for its modification and improvement we present here. PMID:21423712
Kotlarchyk, M; Chen, S H; Asano, S
1979-07-15
The quasi-elastic light scattering has become an established technique for a rapid and quantitative characterization of an average motility pattern of motile bacteria in suspensions. Essentially all interpretations of the measured light scattering intensities and spectra so far are based on the Rayleigh-Gans-Debye (RGD) approximation. Since the range of sizes of bacteria of interest is generally larger than the wavelength of light used in the measurement, one is not certain of the justification for the use of the RGD approximation. In this paper we formulate a method by which both the scattering intensity and the quasi-elastic light scattering spectra can be calculated from a rigorous scattering theory. For a specific application we study the case of bacteria Escherichia coli (about 1 microm in size) by using numerical solutions of the scattering field amplitudes from a prolate spheroid, which is known to simulate optical properties of the bacteria well. We have computed (1) polarized scattered light intensity vs scattering angle for a randomly oriented bacteria population; (2) polarized scattered field correlation functions for both a freely diffusing bacterium and for a bacterium undergoing a straight line motion in random directions and with a Maxwellian speed distribution; and (3) the corresponding depolarized scattered intensity and field correlation functions. In each case sensitivity of the result to variations of the index of refraction and size of the bacterium is investigated. The conclusion is that within a reasonable range of parameters applicable to E. coli, the accuracy of the RGD is good to within 10% at all angles for the properties (1) and (2), and the depolarized contributions in (3) are generally very small. PMID:20212685
Gjini, Erida; Haydon, Daniel T; David Barry, J; Cobbold, Christina A
2014-01-21
Genetic diversity in multigene families is shaped by multiple processes, including gene conversion and point mutation. Because multi-gene families are involved in crucial traits of organisms, quantifying the rates of their genetic diversification is important. With increasing availability of genomic data, there is a growing need for quantitative approaches that integrate the molecular evolution of gene families with their higher-scale function. In this study, we integrate a stochastic simulation framework with population genetics theory, namely the diffusion approximation, to investigate the dynamics of genetic diversification in a gene family. Duplicated genes can diverge and encode new functions as a result of point mutation, and become more similar through gene conversion. To model the evolution of pairwise identity in a multigene family, we first consider all conversion and mutation events in a discrete manner, keeping track of their details and times of occurrence; second we consider only the infinitesimal effect of these processes on pairwise identity accounting for random sampling of genes and positions. The purely stochastic approach is closer to biological reality and is based on many explicit parameters, such as conversion tract length and family size, but is more challenging analytically. The population genetics approach is an approximation accounting implicitly for point mutation and gene conversion, only in terms of per-site average probabilities. Comparison of these two approaches across a range of parameter combinations reveals that they are not entirely equivalent, but that for certain relevant regimes they do match. As an application of this modelling framework, we consider the distribution of nucleotide identity among VSG genes of African trypanosomes, representing the most prominent example of a multi-gene family mediating parasite antigenic variation and within-host immune evasion. PMID:24120993
Mean field approximation for biased diffusion on Japanese inter-firm trading network.
Watanabe, Hayafumi
2014-01-01
By analysing the financial data of firms across Japan, a nonlinear power law with an exponent of 1.3 was observed between the number of business partners (i.e. the degree of the inter-firm trading network) and sales. In a previous study using numerical simulations, we found that this scaling can be explained by both the money-transport model, where a firm (i.e. customer) distributes money to its out-edges (suppliers) in proportion to the in-degree of destinations, and by the correlations among the Japanese inter-firm trading network. However, in this previous study, we could not specifically identify what types of structure properties (or correlations) of the network determine the 1.3 exponent. In the present study, we more clearly elucidate the relationship between this nonlinear scaling and the network structure by applying mean-field approximation of the diffusion in a complex network to this money-transport model. Using theoretical analysis, we obtained the mean-field solution of the model and found that, in the case of the Japanese firms, the scaling exponent of 1.3 can be determined from the power law of the average degree of the nearest neighbours of the network with an exponent of -0.7. PMID:24626149
Mean Field Approximation for Biased Diffusion on Japanese Inter-Firm Trading Network
Watanabe, Hayafumi
2014-01-01
By analysing the financial data of firms across Japan, a nonlinear power law with an exponent of 1.3 was observed between the number of business partners (i.e. the degree of the inter-firm trading network) and sales. In a previous study using numerical simulations, we found that this scaling can be explained by both the money-transport model, where a firm (i.e. customer) distributes money to its out-edges (suppliers) in proportion to the in-degree of destinations, and by the correlations among the Japanese inter-firm trading network. However, in this previous study, we could not specifically identify what types of structure properties (or correlations) of the network determine the 1.3 exponent. In the present study, we more clearly elucidate the relationship between this nonlinear scaling and the network structure by applying mean-field approximation of the diffusion in a complex network to this money-transport model. Using theoretical analysis, we obtained the mean-field solution of the model and found that, in the case of the Japanese firms, the scaling exponent of 1.3 can be determined from the power law of the average degree of the nearest neighbours of the network with an exponent of −0.7. PMID:24626149
de Boissieu, M; Francoual, S; Kaneko, Y; Ishimasa, T
2005-09-01
We report on the absolute scale measurement of the x-ray diffuse scattering in the ZnMgSc icosahedral quasicrystal and its periodic approximant. Whereas the diffuse scattering in the approximant is purely accounted for by thermal diffuse scattering, an additional signal is observed in the quasicrystal. It is related to phason fluctuations as indicated by its Q(2)(per) dependence. Moreover, when compared to previous measurements carried out on the i-AlPdMn phase, we find that the amount of diffuse scattering is smaller in the i-ZnMgSc phase, in agreement with larger phason elastic constants in this phase. This is confirmed by the observation of a large number of weak Bragg peaks having a high Q(per) reciprocal space component. PMID:16196940
Extracting spatial information from large aperture exposures of diffuse sources
NASA Technical Reports Server (NTRS)
Clarke, J. T.; Moos, H. W.
1981-01-01
The spatial properties of large aperture exposures of diffuse emission can be used both to investigate spatial variations in the emission and to filter out camera noise in exposures of weak emission sources. Spatial imaging can be accomplished both parallel and perpendicular to dispersion with a resolution of 5-6 arc sec, and a narrow median filter running perpendicular to dispersion across a diffuse image selectively filters out point source features, such as reseaux marks and fast particle hits. Spatial information derived from observations of solar system objects is presented.
An approximate algorithm for the flux from a rectangular volume source
Wallace, O.J.
1994-11-09
An exact semi-analytic formula for the flux from a rectangular surface source with a slab shield has been derived and the required function table has been calculated. This formula is the basis for an algorithm which gives a good approximation for the flux from a rectangular volume source. No other hand calculation method for this source geometry is available in the literature.
NASA Technical Reports Server (NTRS)
Mair, R. W.; Sen, P. N.; Hurlimann, M. D.; Patz, S.; Cory, D. G.; Walsworth, R. L.
2002-01-01
We report a systematic study of xenon gas diffusion NMR in simple model porous media, random packs of mono-sized glass beads, and focus on three specific areas peculiar to gas-phase diffusion. These topics are: (i) diffusion of spins on the order of the pore dimensions during the application of the diffusion encoding gradient pulses in a PGSE experiment (breakdown of the narrow pulse approximation and imperfect background gradient cancellation), (ii) the ability to derive long length scale structural information, and (iii) effects of finite sample size. We find that the time-dependent diffusion coefficient, D(t), of the imbibed xenon gas at short diffusion times in small beads is significantly affected by the gas pressure. In particular, as expected, we find smaller deviations between measured D(t) and theoretical predictions as the gas pressure is increased, resulting from reduced diffusion during the application of the gradient pulse. The deviations are then completely removed when water D(t) is observed in the same samples. The use of gas also allows us to probe D(t) over a wide range of length scales and observe the long time asymptotic limit which is proportional to the inverse tortuosity of the sample, as well as the diffusion distance where this limit takes effect (approximately 1-1.5 bead diameters). The Pade approximation can be used as a reference for expected xenon D(t) data between the short and the long time limits, allowing us to explore deviations from the expected behavior at intermediate times as a result of finite sample size effects. Finally, the application of the Pade interpolation between the long and the short time asymptotic limits yields a fitted length scale (the Pade length), which is found to be approximately 0.13b for all bead packs, where b is the bead diameter. c. 2002 Elsevier Sciences (USA).
Mair, R W; Sen, P N; Hürlimann, M D; Patz, S; Cory, D G; Walsworth, R L
2002-06-01
We report a systematic study of xenon gas diffusion NMR in simple model porous media, random packs of mono-sized glass beads, and focus on three specific areas peculiar to gas-phase diffusion. These topics are: (i) diffusion of spins on the order of the pore dimensions during the application of the diffusion encoding gradient pulses in a PGSE experiment (breakdown of the narrow pulse approximation and imperfect background gradient cancellation), (ii) the ability to derive long length scale structural information, and (iii) effects of finite sample size. We find that the time-dependent diffusion coefficient, D(t), of the imbibed xenon gas at short diffusion times in small beads is significantly affected by the gas pressure. In particular, as expected, we find smaller deviations between measured D(t) and theoretical predictions as the gas pressure is increased, resulting from reduced diffusion during the application of the gradient pulse. The deviations are then completely removed when water D(t) is observed in the same samples. The use of gas also allows us to probe D(t) over a wide range of length scales and observe the long time asymptotic limit which is proportional to the inverse tortuosity of the sample, as well as the diffusion distance where this limit takes effect (approximately 1-1.5 bead diameters). The Padé approximation can be used as a reference for expected xenon D(t) data between the short and the long time limits, allowing us to explore deviations from the expected behavior at intermediate times as a result of finite sample size effects. Finally, the application of the Padé interpolation between the long and the short time asymptotic limits yields a fitted length scale (the Padé length), which is found to be approximately 0.13b for all bead packs, where b is the bead diameter. PMID:12165255
Singh, Brajesh K.; Srivastava, Vineet K.
2015-01-01
The main goal of this paper is to present a new approximate series solution of the multi-dimensional (heat-like) diffusion equation with time-fractional derivative in Caputo form using a semi-analytical approach: fractional-order reduced differential transform method (FRDTM). The efficiency of FRDTM is confirmed by considering four test problems of the multi-dimensional time fractional-order diffusion equation. FRDTM is a very efficient, effective and powerful mathematical tool which provides exact or very close approximate solutions for a wide range of real-world problems arising in engineering and natural sciences, modelled in terms of differential equations. PMID:26064639
NASA Astrophysics Data System (ADS)
Ghanbarian, Behzad; Daigle, Hugh; Hunt, Allen G.; Ewing, Robert P.; Sahimi, Muhammad
2015-01-01
Understanding and accurate prediction of gas or liquid phase (solute) diffusion are essential to accurate prediction of contaminant transport in partially saturated porous media. In this study, we propose analytical equations, using concepts from percolation theory and the Effective Medium Approximation (EMA) to model the saturation dependence of both gas and solute diffusion in porous media. The predictions of our theoretical approach agree well with the results of nine lattice Boltzmann simulations. We find that the universal quadratic scaling predicted by percolation theory, combined with the universal linear scaling predicted by the EMA, describes diffusion in porous media with both relatively broad and extremely narrow pore size distributions.
Multi-source information diffusion in online social networks
NASA Astrophysics Data System (ADS)
Xiong, Fei; Liu, Yun; Zhang, Hai-Feng
2015-07-01
Individual spreading behavior in online social networks is closely related to user activity, tie strength, and other user and network features. The results concentrate on personal spreading decisions; however, whether these features promote the global information diffusion and increase the macroscopic density of infected agents, remains unclear to us. In this paper, we propose a multi-source diffusion model in which agents may create new messages and spread other agents’ messages. Agents receive many messages, and each time they select a certain message preferentially to spread in consideration of different features. Simulation results show the density of infected agents for different messages follows a power-law distribution in both scale-free and small-world networks. Selecting the largest author degree, author activity and tie strength preferentially can advance the overall diffusion process. Weak tie bias is the least effective feature for multiple information diffusion, but it helps to diffuse a single message. Unexpectedly, the bias of interest similarity does not have an apparent effect. Integrated with the influence on individual diffusion behavior, strong tie bias is a significant feature both for local and global diffusion.
Application of diffusion approximation for risk assessments of sea turtle populations.
Snover, Melissa L; Heppell, Selina S
2009-04-01
Population viability analysis (PVA) to forecast extinction risk is a commonly used tool in decision- and policy-making processes of governments and conservation organizations. A drawback to PVA is the high degree of uncertainty in these forecasts due to both population stochasticity and parameter estimation uncertainty. With sparse or noisy data, extinction probabilities frequently have 95% confidence intervals ranging from 0 to 1. To make stochastic simulation results more interpretable, we present a new metric, susceptibility to quasi-extinction (SQE), to assess whether or not a population is at risk of declining to a prespecified level (quasi-extinction). Following standard methods for diffusion approximation of extinction risk, we use a parametric bootstrap to determine the 95% CI for the probability of quasi-extinction. SQE is the proportion of this parametric bootstrap that indicates a high (defined as > or = 0.90) probability of quasi-extinction, resulting in a point estimate that integrates both parameter uncertainty and stochasticity in extinction forecasting. We demonstrate the application of the metric with sea turtle nest census data, which have a high degree of year-to-year variance and represent only a small fraction of the total population. Using population simulations, we found that for these types of data a critical SQE value of 0.40 corresponds to populations that have a true risk of quasi-extinction. The metric has an accuracy of > 80%, which can be increased further by lowering the 0.40 threshold and trading off Type I error (considering a population to be "not at risk" when it actually is) and Type II error (considering a population to be "at risk" when it actually is not), giving managers a flexible and quantitative tool for assessments of population status. PMID:19425438
Diffusion approximation-based simulation of stochastic ion channels: which method to use?
Pezo, Danilo; Soudry, Daniel; Orio, Patricio
2014-01-01
To study the effects of stochastic ion channel fluctuations on neural dynamics, several numerical implementation methods have been proposed. Gillespie's method for Markov Chains (MC) simulation is highly accurate, yet it becomes computationally intensive in the regime of a high number of channels. Many recent works aim to speed simulation time using the Langevin-based Diffusion Approximation (DA). Under this common theoretical approach, each implementation differs in how it handles various numerical difficulties—such as bounding of state variables to [0,1]. Here we review and test a set of the most recently published DA implementations (Goldwyn et al., 2011; Linaro et al., 2011; Dangerfield et al., 2012; Orio and Soudry, 2012; Schmandt and Galán, 2012; Güler, 2013; Huang et al., 2013a), comparing all of them in a set of numerical simulations that assess numerical accuracy and computational efficiency on three different models: (1) the original Hodgkin and Huxley model, (2) a model with faster sodium channels, and (3) a multi-compartmental model inspired in granular cells. We conclude that for a low number of channels (usually below 1000 per simulated compartment) one should use MC—which is the fastest and most accurate method. For a high number of channels, we recommend using the method by Orio and Soudry (2012), possibly combined with the method by Schmandt and Galán (2012) for increased speed and slightly reduced accuracy. Consequently, MC modeling may be the best method for detailed multicompartment neuron models—in which a model neuron with many thousands of channels is segmented into many compartments with a few hundred channels. PMID:25404914
The discovery of diffuse steep spectrum sources in Abell 2256
NASA Astrophysics Data System (ADS)
van Weeren, R. J.; Intema, H. T.; Oonk, J. B. R.; Röttgering, H. J. A.; Clarke, T. E.
2009-12-01
Context: Hierarchical galaxy formation models indicate that during their lifetime galaxy clusters undergo several mergers. An example of such a merging cluster is Abell 2256. Here we report on the discovery of three diffuse radio sources in the periphery of Abell 2256, using the Giant Metrewave Radio Telescope (GMRT). Aims: The aim of the observations was to search for diffuse ultra-steep spectrum radio sources within the galaxy cluster Abell 2256. Methods: We have carried out GMRT 325 MHz radio continuum observations of Abell 2256. V, R and I band images of the cluster were taken with the 4.2 m William Herschel Telescope (WHT). Results: We have discovered three diffuse elongated radio sources located about 1 Mpc from the cluster center. Two are located to the west of the cluster center, and one to the southeast. The sources have a measured physical extent of 170, 140 and 240 kpc, respectively. The two western sources are also visible in deep low-resolution 115-165 MHz Westerbork Synthesis Radio Telescope (WSRT) images, although they are blended into a single source. For the combined emission of the blended source we find an extreme spectral index (α) of -2.05 ± 0.14 between 140 and 351 MHz. The extremely steep spectral index suggests these two sources are most likely the result of adiabatic compression of fossil radio plasma due to merger shocks. For the source to the southeast, we find that {α < -1.45} between 1369 and 325 MHz. We did not find any clear optical counterparts to the radio sources in the WHT images. Conclusions: The discovery of the steep spectrum sources implies the existence of a population of faint diffuse radio sources in (merging) clusters with such steep spectra that they have gone unnoticed in higher frequency (⪆1 GHz) observations. Simply considering the timescales related to the AGN activity, synchrotron losses, and the presence of shocks, we find that most massive clusters should possess similar sources. An exciting possibility
Cong, Wenxiang; Wang, Lihong V; Wang, Ge
2004-01-01
Background The bioluminescent enzyme firefly luciferase (Luc) or variants of green fluorescent protein (GFP) in transformed cells can be effectively used to reveal molecular and cellular features of neoplasia in vivo. Tumor cell growth and regression in response to various therapies can be evaluated by using bioluminescent imaging. In bioluminescent imaging, light propagates in highly scattering tissue, and the diffusion approximation is sufficiently accurate to predict the imaging signal around the biological tissue. The numerical solutions to the diffusion equation take large amounts of computational time, and the studies for its analytic solutions have attracted more attention in biomedical engineering applications. Methods Biological tissue is a turbid medium that both scatters and absorbs photons. An accurate model for the propagation of photons through tissue can be adopted from transport theory, and its diffusion approximation is applied to predict the imaging signal around the biological tissue. The solution to the diffusion equation is formulated by the convolution between its Green's function and source term. The formulation of photon diffusion from spherical bioluminescent sources in an infinite homogeneous medium can be obtained to accelerate the forward simulation of bioluminescent phenomena. Results The closed form solutions have been derived for the time-dependent diffusion equation and the steady-state diffusion equation with solid and hollow spherical sources in a homogeneous medium, respectively. Meanwhile, the relationship between solutions with a solid sphere source and ones with a surface sphere source is obtained. Conclusion We have formulated solutions for the diffusion equation with solid and hollow spherical sources in an infinite homogeneous medium. These solutions have been verified by Monte Carlo simulation for use in biomedical optical imaging studies. The closed form solution is highly accurate and more computationally efficient in
NASA Astrophysics Data System (ADS)
Martín, Jonatan; Torres, Diego F.; Rea, Nanda
2012-11-01
In this work, we present a leptonic, time-dependent model of pulsar wind nebulae (PWNe). The model seeks a solution for the lepton distribution function considering the full time-energy-dependent diffusion-loss equation. The time-dependent lepton population is balanced by injection, energy losses and escape. We include synchrotron, inverse-Compton (IC; with the cosmic-microwave background as well as with IR/optical photon fields), self-synchrotron Compton, and bremsstrahlung processes, all devoid of any radiative approximations. With this model in place we focus on the Crab nebula as an example and present its time-dependent evolution. Afterwards, we analyse the impact of different approximations made at the level of the diffusion-loss equation, as can be found in the literature. Whereas previous models ignored the escape term, e.g. with the diffusion-loss equation becoming advective, others approximated the losses as catastrophic, so that the equation has only time derivatives. Additional approximations are also described and computed. We study what the impact of these approaches is on the determination of the PWN evolution. In particular, we find the time-dependent deviation of the multi-wavelength spectrum and the best-fitting parameters obtained with the complete and the approximate models.
Traytak, Sergey D.
2014-06-14
The anisotropic 3D equation describing the pointlike particles diffusion in slender impermeable tubes of revolution with cross section smoothly depending on the longitudinal coordinate is the object of our study. We use singular perturbations approach to find the rigorous asymptotic expression for the local particles concentration as an expansion in the ratio of the characteristic transversal and longitudinal diffusion relaxation times. The corresponding leading-term approximation is a generalization of well-known Fick-Jacobs approximation. This result allowed us to delineate the conditions on temporal and spatial scales under which the Fick-Jacobs approximation is valid. A striking analogy between solution of our problem and the method of inner-outer expansions for low Knudsen numbers gas kinetic theory is established. With the aid of this analogy we clarify the physical and mathematical meaning of the obtained results.
NASA Astrophysics Data System (ADS)
Shendeleva, Margarita L.
2004-03-01
Diffuse photon-density waves generated by an instantaneous line source that is parallel to the interface between two semi-infinite turbid media are studied by use of the diffusion approximation. For two nonabsorbing media the Green functions for diffuse light are obtained based on the Green functions for temperature fields that were derived with the Cagniard-de Hoop method. The boundary conditions for diffuse light take into account the discontinuity in the specific intensity at the interface between two media with different refractive indices. The results of the calculations of the specific intensities and the gradient lines for different sets of parameters are presented.
NASA Astrophysics Data System (ADS)
Armstrong, N. M. R.; Mortimer, K. D.; Kong, T.; Bud'ko, S. L.; Canfield, P. C.; Basov, D. N.; Timusk, T.
2016-04-01
Icosahedral quasicrystals are characterised by the absence of a distinct Drude peak in their low-frequency optical conductivity and the same is true of their crystalline approximants. We have measured the optical conductivity of i-GdCd?, an icosahedral quasicrystal, and two approximants, GdCd? and YCd?. We find that there is a significant difference in the optical properties of these compounds. The approximants have a zero frequency peak, characteristic of a metal, whereas the quasicrystal has a striking minimum. This is the first example where the transport properties of a quasicrystal and its approximant differ in such a fundamental way. Using a generalised Drude model introduced by Mayou, we find that our data are well described by this model. It implies that the quantum diffusion of electron wave packets through the periodic and quasiperiodic lattices is responsible for these dramatic differences: in the approximants, the transport is superdiffusive, whereas the quasicrystals show subdiffusive motion of the electrons.
Diffusion of Heat from a Line Source in Isotropic Turbulence
NASA Technical Reports Server (NTRS)
Uberoi, Mahinder S; Corrsin, Stanley
1953-01-01
An experimental and analytical study has been made of some features of the turbulent heat diffusion behind a line heated wire stretched perpendicular to a flowing isotropic turbulence. The mean temperature distributions have been measured with systematic variations in wind speed, size of turbulence-producing grid, and downstream location of heat source. The nature of the temperature fluctuation field has been studied. A comparison of Lagrangian and Eulerian analyses for diffusion in a nondecaying turbulence yields an expression for turbulent-heat-transfer coefficient in terms of turbulence velocity and a Lagrangian "scale." the ratio of Eulerian to Lagrangian microscale has been determined theoretically by generalization of a result of Heisenberg and with arbitrary constants taken from independent sources, shows rough agreement with experimental results. A convenient form has been deduced for the criterion of interchangeability of instantaneous space and time derivatives in a flowing turbulence.
NASA Astrophysics Data System (ADS)
Gorpas, Dimitris; Andersson-Engels, Stefan
2012-03-01
The solution of the forward problem in fluorescence molecular imaging is among the most important premises for the successful confrontation of the inverse reconstruction problem. To date, the most typical approach has been the application of the diffusion approximation as the forward model. This model is basically a first order angular approximation for the radiative transfer equation, and thus it presents certain limitations. The scope of this manuscript is to present the dual coupled radiative transfer equation and diffusion approximation model for the solution of the forward problem in fluorescence molecular imaging. The integro-differential equations of its weak formalism were solved via the finite elements method. Algorithmic blocks with cubature rules and analytical solutions of the multiple integrals have been constructed for the solution. Furthermore, specialized mapping matrices have been developed to assembly the finite elements matrix. As a radiative transfer equation based model, the integration over the angular discretization was implemented analytically, while quadrature rules were applied whenever required. Finally, this model was evaluated on numerous virtual phantoms and its relative accuracy, with respect to the radiative transfer equation, was over 95%, when the widely applied diffusion approximation presented almost 85% corresponding relative accuracy for the fluorescence emission.
Huang, Jianping; Wang, Lihui; Chu, Chunyu; Zhang, Yanli; Liu, Wanyu; Zhu, Yuemin
2016-04-29
Diffusion tensor magnetic resonance (DTMR) imaging and diffusion tensor imaging (DTI) have been widely used to probe noninvasively biological tissue structures. However, DTI suffers from long acquisition times, which limit its practical and clinical applications. This paper proposes a new Compressed Sensing (CS) reconstruction method that employs joint sparsity and rank deficiency to reconstruct cardiac DTMR images from undersampled k-space data. Diffusion-weighted images acquired in different diffusion directions were firstly stacked as columns to form the matrix. The matrix was row sparse in the transform domain and had a low rank. These two properties were then incorporated into the CS reconstruction framework. The underlying constrained optimization problem was finally solved by the first-order fast method. Experiments were carried out on both simulation and real human cardiac DTMR images. The results demonstrated that the proposed approach had lower reconstruction errors for DTI indices, including fractional anisotropy (FA) and mean diffusivities (MD), compared to the existing CS-DTMR image reconstruction techniques. PMID:27163322
Diffusion and reaction for a spherical source and sink
NASA Astrophysics Data System (ADS)
McDonald, Nyrée; Strieder, William
2003-03-01
Two chemically active spheres in an infinite medium, one a zeroth-order reactant source and the other a first-order sink, are considered for various sphere size ratios, center-to-center distances, and sink strengths from chemical to diffusion controlled conditions. This source-sink model simulates some aspects of biological mutualism interactions between different cells. Infinite series expansions in a single index n are obtained for the sink reaction rate and reactant concentration profiles using the bispherical expansion. Each of the coefficients, generated exactly by a matrix elimination method, is expressed in terms of nested, continued fractions easily evaluated for the given n. At intermediate and larger sink-source separation distances the sink reaction rate decays harmonically. For smaller sink-source separations with a highly reactive small sink, a local maximum in the sink reaction rate is found.
NASA Astrophysics Data System (ADS)
Bologna, Mauro; Svenkeson, Adam; West, Bruce J.; Grigolini, Paolo
2015-07-01
Diffusion processes in heterogeneous media, and biological systems in particular, are riddled with the difficult theoretical issue of whether the true origin of anomalous behavior is renewal or memory, or a special combination of the two. Accounting for the possible mixture of renewal and memory sources of subdiffusion is challenging from a computational point of view as well. This problem is exacerbated by the limited number of techniques available for solving fractional diffusion equations with time-dependent coefficients. We propose an iterative scheme for solving fractional differential equations with time-dependent coefficients that is based on a parametric expansion in the fractional index. We demonstrate how this method can be used to predict the long-time behavior of nonautonomous fractional differential equations by studying the anomalous diffusion process arising from a mixture of renewal and memory sources.
Super-resolution image reconstruction using diffuse source models.
Ellis, Michael A; Viola, Francesco; Walker, William F
2010-06-01
Image reconstruction is central to many scientific fields, from medical ultrasound and sonar to computed tomography and computer vision. Although lenses play a critical reconstruction role in these fields, digital sensors enable more sophisticated computational approaches. A variety of computational methods have thus been developed, with the common goal of increasing contrast and resolution to extract the greatest possible information from raw data. This paper describes a new image reconstruction method named the Diffuse Time-domain Optimized Near-field Estimator (dTONE). dTONE represents each hypothetical target in the system model as a diffuse region of targets rather than a single discrete target, which more accurately represents the experimental data that arise from signal sources in continuous space, with no additional computational requirements at the time of image reconstruction. Simulation and experimental ultrasound images of animal tissues show that dTONE achieves image resolution and contrast far superior to those of conventional image reconstruction methods. We also demonstrate the increased robustness of the diffuse target model to major sources of image degradation through the addition of electronic noise, phase aberration and magnitude aberration to ultrasound simulations. Using experimental ultrasound data from a tissue-mimicking phantom containing a 3-mm-diameter anechoic cyst, the conventionally reconstructed image has a cystic contrast of -6.3 dB, whereas the dTONE image has a cystic contrast of -14.4 dB. PMID:20447760
Super-Resolution Image Reconstruction Using Diffuse Source Models
Ellis, Michael A.; Viola, Francesco; Walker, William F.
2010-01-01
Image reconstruction is central to many scientific fields, from medical ultrasound and sonar to computed tomography and computer vision. While lenses play a critical reconstruction role in these fields, digital sensors enable more sophisticated computational approaches. A variety of computational methods have thus been developed, with the common goal of increasing contrast and resolution to extract the greatest possible information from raw data. This paper describes a new image reconstruction method named the Diffuse Time-domain Optimized Near-field Estimator (dTONE). dTONE represents each hypothetical target in the system model as a diffuse region of targets rather than a single discrete target, which more accurately represents the experimental data that arise from signal sources in continuous space, with no additional computational requirements at the time of image reconstruction. Simulation and experimental ultrasound images of animal tissues show that dTONE achieves image resolution and contrast far superior to those of conventional image reconstruction methods. We also demonstrate the increased robustness of the diffuse target model to major sources of image degradation, through the addition of electronic noise, phase aberration, and magnitude aberration to ultrasound simulations. Using experimental ultrasound data from a tissue-mimicking phantom containing a 3 mm diameter anechoic cyst, the conventionally reconstructed image has a cystic contrast of −6.3 dB whereas the dTONE image has a cystic contrast of −14.4 dB. PMID:20447760
An Approximate Markov Model for the Wright-Fisher Diffusion and Its Application to Time Series Data.
Ferrer-Admetlla, Anna; Leuenberger, Christoph; Jensen, Jeffrey D; Wegmann, Daniel
2016-06-01
The joint and accurate inference of selection and demography from genetic data is considered a particularly challenging question in population genetics, since both process may lead to very similar patterns of genetic diversity. However, additional information for disentangling these effects may be obtained by observing changes in allele frequencies over multiple time points. Such data are common in experimental evolution studies, as well as in the comparison of ancient and contemporary samples. Leveraging this information, however, has been computationally challenging, particularly when considering multilocus data sets. To overcome these issues, we introduce a novel, discrete approximation for diffusion processes, termed mean transition time approximation, which preserves the long-term behavior of the underlying continuous diffusion process. We then derive this approximation for the particular case of inferring selection and demography from time series data under the classic Wright-Fisher model and demonstrate that our approximation is well suited to describe allele trajectories through time, even when only a few states are used. We then develop a Bayesian inference approach to jointly infer the population size and locus-specific selection coefficients with high accuracy and further extend this model to also infer the rates of sequencing errors and mutations. We finally apply our approach to recent experimental data on the evolution of drug resistance in influenza virus, identifying likely targets of selection and finding evidence for much larger viral population sizes than previously reported. PMID:27038112
NASA Astrophysics Data System (ADS)
Seyedabbasi, Mir Ahmad; Newell, Charles J.; Adamson, David T.; Sale, Thomas C.
2012-06-01
The relative contribution of dense non-aqueous phase liquid (DNAPL) dissolution versus matrix diffusion processes to the longevity of chlorinated source zones was investigated. Matrix diffusion is being increasingly recognized as an important non-DNAPL component of source behavior over time, and understanding the persistence of contaminants that have diffused into lower permeability units can impact remedial decision-making. In this study, a hypothetical DNAPL source zone architecture consisting of several different sized pools and fingers originally developed by Anderson et al. (1992) was adapted to include defined low permeability layers. A coupled dissolution-diffusion model was developed to allow diffusion into these layers while in contact with DNAPL, followed by diffusion out of these same layers after complete DNAPL dissolution. This exercise was performed for releases of equivalent masses (675 kg) of three different compounds, including chlorinated solvents with solubilities ranging from low (tetrachloroethene (PCE)), moderate (trichloroethene (TCE)) to high (dichloromethane (DCM)). The results of this simple modeling exercise demonstrate that matrix diffusion can be a critical component of source zone longevity and may represent a longer-term contributor to source longevity (i.e., longer time maintaining concentrations above MCLs) than DNAPL dissolution alone at many sites. For the hypothetical TCE release, the simulation indicated that dissolution of DNAPL would take approximately 38 years, while the back diffusion from low permeability zones could maintain the source for an additional 83 years. This effect was even more dramatic for the higher solubility DCM (97% of longevity due to matrix diffusion), while the lower solubility PCE showed a more equal contribution from DNAPL dissolution vs. matrix diffusion. Several methods were used to describe the resulting source attenuation curves, including a first-order decay model which showed that half-life of
A New 2D-Transport, 1D-Diffusion Approximation of the Boltzmann Transport equation
Larsen, Edward
2013-06-17
The work performed in this project consisted of the derivation, implementation, and testing of a new, computationally advantageous approximation to the 3D Boltz- mann transport equation. The solution of the Boltzmann equation is the neutron flux in nuclear reactor cores and shields, but solving this equation is difficult and costly. The new “2D/1D” approximation takes advantage of a special geometric feature of typical 3D reactors to approximate the neutron transport physics in a specific (ax- ial) direction, but not in the other two (radial) directions. The resulting equation is much less expensive to solve computationally, and its solutions are expected to be sufficiently accurate for many practical problems. In this project we formulated the new equation, discretized it using standard methods, developed a stable itera- tion scheme for solving the equation, implemented the new numerical scheme in the MPACT code, and tested the method on several realistic problems. All the hoped- for features of this new approximation were seen. For large, difficult problems, the resulting 2D/1D solution is highly accurate, and is calculated about 100 times faster than a 3D discrete ordinates simulation.
Far-Infrared sources and diffuse emission in M31
NASA Technical Reports Server (NTRS)
Xu, Cong; Helou, George
1994-01-01
A study on the far-infrared (FIR) emission of M31 has been carried out with the High Resolution (HiRes) maps (approx. 1 min) derived from IRAS data. Sixty-eight FIR sources are detected in M31, which in general coincide with optical HII regions, and contribute 15, 23, 29, and 23 percent to the fluxes in 12, 25, 60, and 100 micron bands, respectively. The remaining diffuse emission, which dominates the FIR emission of M31, is studied using a dust heating model which utilizes the UV and optical photometry maps and the HI maps available in the literature. It is found that the global dust-to-gas ratio in M31 disk is 6.5 10(exp -3), very close to the dust-to-gas ratio in the solar neighborhood. There is a significant galactocentric gradient of the dust-to-HI-gas ratio, with an e-folding scale length of 9 kpc. The diffuse dust correlates tightly with the HI gas. The model indicates that the non-ionizing UV (913-4000A) radiation from massive and intermediate massive stars contributes only about 30 percent of the heating of the diffuse dust, while the optical-NIR (4000-9000A) radiation from the old stellar population is responsible for the most of the heating.
NASA Astrophysics Data System (ADS)
Zhang, Tianhe C.; Grill, Warren M.
2010-12-01
Deep brain stimulation (DBS) has emerged as an effective treatment for movement disorders; however, the fundamental mechanisms by which DBS works are not well understood. Computational models of DBS can provide insights into these fundamental mechanisms and typically require two steps: calculation of the electrical potentials generated by DBS and, subsequently, determination of the effects of the extracellular potentials on neurons. The objective of this study was to assess the validity of using a point source electrode to approximate the DBS electrode when calculating the thresholds and spatial distribution of activation of a surrounding population of model neurons in response to monopolar DBS. Extracellular potentials in a homogenous isotropic volume conductor were calculated using either a point current source or a geometrically accurate finite element model of the Medtronic DBS 3389 lead. These extracellular potentials were coupled to populations of model axons, and thresholds and spatial distributions were determined for different electrode geometries and axon orientations. Median threshold differences between DBS and point source electrodes for individual axons varied between -20.5% and 9.5% across all orientations, monopolar polarities and electrode geometries utilizing the DBS 3389 electrode. Differences in the percentage of axons activated at a given amplitude by the point source electrode and the DBS electrode were between -9.0% and 12.6% across all monopolar configurations tested. The differences in activation between the DBS and point source electrodes occurred primarily in regions close to conductor-insulator interfaces and around the insulating tip of the DBS electrode. The robustness of the point source approximation in modeling several special cases—tissue anisotropy, a long active electrode and bipolar stimulation—was also examined. Under the conditions considered, the point source was shown to be a valid approximation for predicting excitation
NASA Astrophysics Data System (ADS)
Carrasco-Hernandez, Roberto; Smedley, Andrew R. D.; Webb, Ann R.
2016-05-01
Two radiative transfer models are presented that simplify calculations of street canyon spectral irradiances with minimum data input requirements, allowing better assessment of urban exposures than can be provided by standard unobstructed radiation measurements alone. Fast calculations improve the computational performance of radiation models, when numerous repetitions are required in time and location. The core of the models is the calculation of the spectral diffuse-to-global ratios (DGR) from an unobstructed global spectral measurement. The models are based on, and have been tested against, outcomes of the SMARTS2 algorithm (i.e. Simple Model of the Atmospheric Radiative Transfer of Sunshine). The modelled DGRs can then be used to partition global spectral irradiance values into their direct and diffuse components for different solar zenith angles. Finally, the effects of canyon obstructions can be evaluated independently on the direct and diffuse components, which are then recombined to give the total canyon irradiance. The first model allows ozone and aerosol inputs, while the second provides a further simplification, restricted to average ozone and aerosol contents but specifically designed for faster calculations. To assess the effect of obstructions and validate the calculations, a set of experiments with simulated obstructions (simulated canyons) were performed. The greatest source of uncertainty in the simplified calculations is in the treatment of diffuse radiation. The measurement-model agreement is therefore dependent on the region of the sky obscured and ranges from <5 % at all wavelengths to 20-40 % (wavelength dependent) when diffuse sky only is visible from the canyon.
Nonlocal diffusion problems that approximate a parabolic equation with spatial dependence
NASA Astrophysics Data System (ADS)
Molino, Alexis; Rossi, Julio D.
2016-06-01
In this paper, we show that smooth solutions to the Dirichlet problem for the parabolic equation v_t(x,t)=sum_{i,j=1}N a_{ij}(x)partial2v(x,t)/partial{xipartial{x}j} + sum_{i =1}N bi(x)partial{v}(x,t)/partial{x_i} qquad x in Ω, with v( x, t) = g( x, t), {x in partial Ω,} can be approximated uniformly by solutions of nonlocal problems of the form ut^{\\varepsilon}(x,t)=int_{mathbb{R}n} K_{\\varepsilon}(x,y)(u^{\\varepsilon}(y,t)-u^{\\varepsilon}(x,t))dy, quad x in Ω, with {u^{\\varepsilon}(x,t)=g(x,t)}, {x notin Ω}, as {\\varepsilon to 0}, for an appropriate rescaled kernel {K_{\\varepsilon}}. In this way, we show that the usual local evolution problems with spatial dependence can be approximated by nonlocal ones. In the case of an equation in divergence form, we can obtain an approximation with symmetric kernels, that is, {K_{\\varepsilon}(x,y) = K_{\\varepsilon}(y,x)}.
Locating the Source of Diffusion in Large-Scale Networks
NASA Astrophysics Data System (ADS)
Pinto, Pedro C.; Thiran, Patrick; Vetterli, Martin
2012-08-01
How can we localize the source of diffusion in a complex network? Because of the tremendous size of many real networks—such as the internet or the human social graph—it is usually unfeasible to observe the state of all nodes in a network. We show that it is fundamentally possible to estimate the location of the source from measurements collected by sparsely placed observers. We present a strategy that is optimal for arbitrary trees, achieving maximum probability of correct localization. We describe efficient implementations with complexity O(Nα), where α=1 for arbitrary trees and α=3 for arbitrary graphs. In the context of several case studies, we determine how localization accuracy is affected by various system parameters, including the structure of the network, the density of observers, and the number of observed cascades.
Approximately a Thousand Ultra-diffuse Galaxies in the Coma Cluster
NASA Astrophysics Data System (ADS)
Koda, Jin; Yagi, Masafumi; Yamanoi, Hitomi; Komiyama, Yutaka
2015-07-01
We report the discovery of 854 ultra-diffuse galaxies (UDGs) in the Coma cluster using deep R band images, with partial B, i, and Hα band coverage, obtained with the Subaru telescope. Many of them (332) are Milky Way (MW) sized with very large effective radii of {r}{e}\\gt 1.5 {kpc}. This study was motivated by the recent discovery of 47 UDGs by Dokkum et al.; our discovery suggests \\gt 1000 UDGs after accounting for the smaller Subaru field (4.1 {{degree}}2; about one-half of Dragonfly). The new Subaru UDGs show a distribution concentrated around the cluster center, strongly suggesting that the great majority are (likely longtime) cluster members. They are a passively evolving population, lying along the red sequence in the color-magnitude diagram with no signature of Hα emission. Star formation was, therefore, quenched in the past. They have exponential light profiles, effective radii {r}{e}˜ 800\\{pc}-5 {kpc}, effective surface brightnesses {μ }{e}({\\text{}}R) = 25-28 mag arcsec-2, and stellar masses ˜ 1× {10}7{\\text{}}{M}⊙ -5× {10}8{\\text{}}{M}⊙ . There is also a population of nucleated UDGs. Some MW-sized UDGs appear closer to the cluster center than previously reported; their survival in the strong tidal field, despite their large sizes, possibly indicates a large dark matter fraction protecting the diffuse stellar component. The indicated baryon fraction ≲ 1% is less than the cosmic average, and thus the gas must have been removed (from the possibly massive dark halo). The UDG population is elevated in the Coma cluster compared to the field, indicating that the gas removal mechanism is related primarily to the cluster environment.
NASA Technical Reports Server (NTRS)
Kottarchyk, M.; Chen, S.-H.; Asano, S.
1979-01-01
The study tests the accuracy of the Rayleigh-Gans-Debye (RGD) approximation against a rigorous scattering theory calculation for a simplified model of E. coli (about 1 micron in size) - a solid spheroid. A general procedure is formulated whereby the scattered field amplitude correlation function, for both polarized and depolarized contributions, can be computed for a collection of particles. An explicit formula is presented for the scattered intensity, both polarized and depolarized, for a collection of randomly diffusing or moving particles. Two specific cases for the intermediate scattering functions are considered: diffusing particles and freely moving particles with a Maxwellian speed distribution. The formalism is applied to microorganisms suspended in a liquid medium. Sensitivity studies revealed that for values of the relative index of refraction greater than 1.03, RGD could be in serious error in computing the intensity as well as correlation functions.
Continued Study on Hohlraum Radiation Source with Approximately Constant Radiation Temperature
NASA Astrophysics Data System (ADS)
Song, Tianming; Yang, Jiamin; Zhu, Tuo; Li, Zhichao; Huang, Chengwu
2016-04-01
An experiment was performed on the Shenguang III prototype laser facility to continue the study on hohlraum radiation source with approximately constant radiation temperature using a continuously shaped laser pulse. A radiation source with a flattop temperature of about 130 eV that lasted about 5 ns was obtained. The previous analytical iteration method based on power balance and self-similar solution of ablation was modified taking into account the plasma movements and it was used to design the laser pulse shape for experiment. A comparison between experimental results and simulation is presented and better agreement was achieved using the modified method. Further improvements are discussed.
Contribution to the diffuse radio background from extragalactic radio sources
NASA Astrophysics Data System (ADS)
Vernstrom, T.; Scott, Douglas; Wall, J. V.
2011-08-01
We examine the brightness of the cosmic radio background (CRB) by comparing the contribution from individual source counts to absolute measurements. We use a compilation of radio counts to estimate the contribution of detected sources to the CRB in several different frequency bands. Using a Monte Carlo Markov chain technique, we estimate the brightness values and uncertainties, paying attention to various sources of systematic error. At ν= 150, 325, 408, 610, 1.4, 4.8 and 8.4 GHz, our calculated contributions to the background sky temperature are 18, 2.8, 1.6, 0.71, 0.11, 0.0032 and 0.0059 K, respectively. We then compare our results to absolute measurements from the Absolute Radiometer for Cosmology, Astrophysics and Diffuse Emission (ARCADE 2) experiment. If the ARCADE 2 measurements are correct and come from sources, then there must be an additional population of radio galaxies, fainter than where current data are probing. More specifically, the Euclidean-normalized counts at 1.4 GHz have to have an additional bump below about 10 μJy.
Verdu, G.; Capilla, M.; Talavera, C. F.; Ginestar, D.
2012-07-01
PL equations are classical high order approximations to the transport equations which are based on the expansion of the angular dependence of the angular neutron flux and the nuclear cross sections in terms of spherical harmonics. A nodal collocation method is used to discretize the PL equations associated with a neutron source transport problem. The performance of the method is tested solving two 1D problems with analytical solution for the transport equation and a classical 2D problem. (authors)
NASA Astrophysics Data System (ADS)
Nazarov, R.; Shulenburger, L.; Morales, M.; Hood, Randolph Q.
2016-03-01
We performed diffusion Monte Carlo (DMC) calculations of the spectroscopic properties of a large set of molecules, assessing the effect of different approximations. In systems containing elements with large atomic numbers, we show that the errors associated with the use of nonlocal mean-field-based pseudopotentials in DMC calculations can be significant and may surpass the fixed-node error. We suggest practical guidelines for reducing these pseudopotential errors, which allow us to obtain DMC-computed spectroscopic parameters of molecules and equation of state properties of solids in excellent agreement with experiment.
Nazarov, Roman; Shulenburger, Luke; Morales, Miguel A.; Hood, Randolph Q.
2016-03-28
Diffusion Monte Carlo (DMC) calculations of the spectroscopic properties of a large set of molecules were performed, assessing the effect of different approximations. In systems containing elements with large atomic numbers, we show that the errors associated with the use of nonlocal mean-field-based pseudopotentials in DMC calculations can be significant and may surpass the fixed-node error. We suggest practical guidelines for reducing these pseudopotential errors, which allow us to obtain DMC-computed spectroscopic parameters of molecules and equation of state properties of solids in excellent agreement with experiment.
NASA Astrophysics Data System (ADS)
Mourenas, D.; Artemyev, A. V.; Agapitov, O. V.
2015-09-01
Modeling the spatiotemporal evolution of relativistic electron fluxes trapped in the Earth's radiation belts in the presence of radial diffusion coupled with wave-induced losses should address one important question: how deep can relativistic electrons penetrate into the inner magnetosphere? However, a full modeling requires extensive numerical simulations solving the comprehensive quasi-linear equations describing pitch angle and radial diffusion of the electron distribution, making it rather difficult to perform parametric studies of the flux behavior. Here we consider the particular situation where a localized flux peak (or storage ring) has been produced at low L < 4 during a period of strong disturbances, through a combination of chorus-induced energy diffusion (or direct injection) at low L together with enhanced wave-induced losses and outward radial transport at higher L. Assuming that radial diffusion can be further described as the spatial broadening within the plasmasphere of this preexisting flux peak, simple approximate analytical solutions for the distribution of trapped relativistic electrons are derived. Such a simplified formalism provides a convenient means for easily determining whether radial diffusion actually prevails over atmospheric losses at any particular time for given electron energy E and location L. It is further used to infer favorable conditions for relativistic electron access to the inner belt, providing an explanation for the relative scarcity of such a feat under most circumstances. Comparisons with electron flux measurements on board the Van Allen Probes show a reasonable agreement between a few weeks and 4 months after the formation of a flux peak.
An inverse time-dependent source problem for a time-fractional diffusion equation
NASA Astrophysics Data System (ADS)
Wei, T.; Li, X. L.; Li, Y. S.
2016-08-01
This paper is devoted to identifying a time-dependent source term in a multi-dimensional time-fractional diffusion equation from boundary Cauchy data. The existence and uniqueness of a strong solution for the corresponding direct problem with homogeneous Neumann boundary condition are firstly proved. We provide the uniqueness and a stability estimate for the inverse time-dependent source problem. Then we use the Tikhonov regularization method to solve the inverse source problem and propose a conjugate gradient algorithm to find a good approximation to the minimizer of the Tikhonov regularization functional. Numerical examples in one-dimensional and two-dimensional cases are provided to show the effectiveness of the proposed method. This paper was supported by the NSF of China (11371181) and the Fundamental Research Funds for the Central Universities (lzujbky-2013-k02).
NASA Astrophysics Data System (ADS)
Chen, Hu; Lü, Shujuan; Chen, Wenping
2016-06-01
The numerical approximation of the distributed order time fractional reaction-diffusion equation on a semi-infinite spatial domain is discussed in this paper. A fully discrete scheme based on finite difference method in time and spectral approximation using Laguerre functions in space is proposed. The scheme is unconditionally stable and convergent with order O (τ2 + Δα2 +N (1 - m) / 2), where τ, Δα, N, and m are the time-step size, step size in distributed-order variable, polynomial degree, and regularity in the space variable of the exact solution, respectively. A pseudospectral scheme is also proposed and analyzed. Some numerical examples are presented to demonstrate the efficiency of the proposed scheme.
Tramline wheelings as diffuse pollution sources: integrating across spatial scales
NASA Astrophysics Data System (ADS)
Silgram, Martyn; Collins, Adrian
2010-05-01
Identifying the sources of diffuse pollution, characterising transport pathways, and understanding the factors controlling delivery to water bodies, are generic challenges faced by scientists and policymakers supporting the implementation of the EU Water Framework Directive and other water quality legislation. Land management is intrinsically field scale, which is the scale at which much fieldwork is focused. However, processes identified at plot or field scale may not be significant in terms of water quality at catchment scale, which poses the problem of how best to upscale relevant field-scale characterisations. A consistent framework is urgently needed which is capable of bringing together experimental techniques at these disparate spatial scales (plot, field and river catchment scale), to provide high resolution evidence to support catchment characterisations (which identify "hot spots") and the spatial targeting of mitigation measures aimed at diffuse pollution control. In the absence of such frameworks, comparing existing data from disparate projects represents a starting point for assessing scaling issues. The emerging evidence base is increasingly underscoring the importance of tramline wheelings as sources and pathways for diffuse pollution loss to the edge of cereal fields. Recent UK research has demonstrated that these compacted, unvegetated areas (which are used as markers in spraying operations) can account for over 80% of surface losses of runoff, sediment and phosphorus to edge-of-field. Current government and industry-funded projects are now exploring a range of novel practical land management techniques to minimise these field-scale losses. However, given the mitigating role of field headlands and boundaries, and riparian buffer strips in limiting connectivity and delivery to water bodies, it has not been clear whether such losses represent an important element of the net pollutant load reaching catchment outlets. Such scaling issues represent an
NASA Astrophysics Data System (ADS)
Barrett, Steven R. H.; Britter, Rex E.
Predicting long-term mean pollutant concentrations in the vicinity of airports, roads and other industrial sources are frequently of concern in regulatory and public health contexts. Many emissions are represented geometrically as ground-level line or area sources. Well developed modelling tools such as AERMOD and ADMS are able to model dispersion from finite (i.e. non-point) sources with considerable accuracy, drawing upon an up-to-date understanding of boundary layer behaviour. Due to mathematical difficulties associated with line and area sources, computationally expensive numerical integration schemes have been developed. For example, some models decompose area sources into a large number of line sources orthogonal to the mean wind direction, for which an analytical (Gaussian) solution exists. Models also employ a time-series approach, which involves computing mean pollutant concentrations for every hour over one or more years of meteorological data. This can give rise to computer runtimes of several days for assessment of a site. While this may be acceptable for assessment of a single industrial complex, airport, etc., this level of computational cost precludes national or international policy assessments at the level of detail available with dispersion modelling. In this paper, we extend previous work [S.R.H. Barrett, R.E. Britter, 2008. Development of algorithms and approximations for rapid operational air quality modelling. Atmospheric Environment 42 (2008) 8105-8111] to line and area sources. We introduce approximations which allow for the development of new analytical solutions for long-term mean dispersion from line and area sources, based on hypergeometric functions. We describe how these solutions can be parameterized from a single point source run from an existing advanced dispersion model, thereby accounting for all processes modelled in the more costly algorithms. The parameterization method combined with the analytical solutions for long-term mean
Diffuse optical tomography using multi-directional sources and detectors
Shimokawa, Takeaki; Ishii, Toshihiro; Takahashi, Yoichiro; Sugawara, Satoru; Sato, Masa-aki; Yamashita, Okito
2016-01-01
Diffuse optical tomography (DOT) is an advanced imaging method used to visualize the internal state of biological tissues as 3D images. However, current continuous-wave DOT requires high-density probe arrays for measurement (less than 15-mm interval) to gather enough information for 3D image reconstruction, which makes the experiment time-consuming. In this paper, we propose a novel DOT measurement system using multi-directional light sources and multi-directional photodetectors instead of high-density probe arrays. We evaluated this system’s multi-directional DOT through computer simulation and a phantom experiment. From the results, we achieved DOT with less than 5-mm localization error up to a 15-mm depth with low-density probe arrays (30-mm interval), indicating that the multi-directional measurement approach allows DOT without requiring high-density measurement. PMID:27446694
Diffuse optical tomography using multi-directional sources and detectors.
Shimokawa, Takeaki; Ishii, Toshihiro; Takahashi, Yoichiro; Sugawara, Satoru; Sato, Masa-Aki; Yamashita, Okito
2016-07-01
Diffuse optical tomography (DOT) is an advanced imaging method used to visualize the internal state of biological tissues as 3D images. However, current continuous-wave DOT requires high-density probe arrays for measurement (less than 15-mm interval) to gather enough information for 3D image reconstruction, which makes the experiment time-consuming. In this paper, we propose a novel DOT measurement system using multi-directional light sources and multi-directional photodetectors instead of high-density probe arrays. We evaluated this system's multi-directional DOT through computer simulation and a phantom experiment. From the results, we achieved DOT with less than 5-mm localization error up to a 15-mm depth with low-density probe arrays (30-mm interval), indicating that the multi-directional measurement approach allows DOT without requiring high-density measurement. PMID:27446694
Scott, David J; Harding, Stephen E; Winzor, Donald J
2015-12-01
This investigation examined the feasibility of manipulating the rotor speed in sedimentation velocity experiments to spontaneously generate an approximate steady-state condition where the extent of diffusional spreading is matched exactly by the boundary sharpening arising from negative s-c dependence. Simulated sedimentation velocity distributions based on the sedimentation characteristics for a purified mucin preparation were used to illustrate a simple procedure for determining the diffusion coefficient from such steady-state distributions in situations where the concentration dependence of the sedimentation coefficient, s = s(0)/(1 + Kc), was quantified in terms of the limiting sedimentation coefficient as c → 0 (s(0)) and the concentration coefficient (K). Those simulations established that spontaneous generation of the approximate steady state could well be a feature of sedimentation velocity distributions for many unstructured polymer systems because the requirement that Kcoω(2)s(0)/D be between 46 and 183 cm(-2) is not unduly restrictive. Although spontaneous generation of the approximate steady state is also a theoretical prediction for structured macromolecular solutes exhibiting linear concentration dependence of the sedimentation coefficient, s = s(0)(1 - kc), the required value of k is far too large for any practical advantage to be taken of this approach with globular proteins. PMID:26321223
Orio, Patricio; Soudry, Daniel
2012-01-01
Background The phenomena that emerge from the interaction of the stochastic opening and closing of ion channels (channel noise) with the non-linear neural dynamics are essential to our understanding of the operation of the nervous system. The effects that channel noise can have on neural dynamics are generally studied using numerical simulations of stochastic models. Algorithms based on discrete Markov Chains (MC) seem to be the most reliable and trustworthy, but even optimized algorithms come with a non-negligible computational cost. Diffusion Approximation (DA) methods use Stochastic Differential Equations (SDE) to approximate the behavior of a number of MCs, considerably speeding up simulation times. However, model comparisons have suggested that DA methods did not lead to the same results as in MC modeling in terms of channel noise statistics and effects on excitability. Recently, it was shown that the difference arose because MCs were modeled with coupled gating particles, while the DA was modeled using uncoupled gating particles. Implementations of DA with coupled particles, in the context of a specific kinetic scheme, yielded similar results to MC. However, it remained unclear how to generalize these implementations to different kinetic schemes, or whether they were faster than MC algorithms. Additionally, a steady state approximation was used for the stochastic terms, which, as we show here, can introduce significant inaccuracies. Main Contributions We derived the SDE explicitly for any given ion channel kinetic scheme. The resulting generic equations were surprisingly simple and interpretable – allowing an easy, transparent and efficient DA implementation, avoiding unnecessary approximations. The algorithm was tested in a voltage clamp simulation and in two different current clamp simulations, yielding the same results as MC modeling. Also, the simulation efficiency of this DA method demonstrated considerable superiority over MC methods, except when
Lomov, I; Liu, B
2005-09-29
Multimaterial Eulerian and Arbitrary Lagragian-Eulerian (ALE) codes usually use volume fractions of materials to track individual components in mixed cells. Material advection usually is calculated either by interface capturing, where a high-order van Leer-like slope reconstruction technique is applied, or interface tracking, where a normal reconstruction technique is applied. The former approach is more appropriate for gas-like substances, and the latter is ideal for solids and liquids, since it does not smear out material interfaces. A wide range of problems involves both diffuse and sharp interfaces between substances and demands a combination of these techniques. It is possible to treat all substances that can diffuse into each other as a single material and only keep mass fractions of the individual components of the mixture. The material response can be determined based on the assumption of pressure and temperature equilibrium between components of the mixture. Unfortunately, it is extremely difficult to solve the corresponding system of equations. In order to avoid these problems one can introduce an effective gamma and employ the ideal gas approximation to calculate mixture response. This method provides reliable results, is able to compute strong shock waves, and deals with complex equations of state. Results from a number of simulations using this scheme are presented.
The effects of As overpressure and diffusion source on the diffusion of Mn in GaAs
NASA Astrophysics Data System (ADS)
Wu, C. H.; Hsieh, K. C.
1992-12-01
Data are presented to show the effect of As overpressure on the diffusion of Mn in GaAs using four different Mn sources. These sources include solid Mn thin film deposited directly on the GaAs substrate and Mn vapors from pure Mn, MnAs, and Mn3As solids. In the circumstance for which a solid Mn film is used as the diffusion source, a nonuniform doping distribution and poor surface morphology is obtained due to a reaction between the Mn film and the GaAs matrix. The degraded surface consists of a layer of polycrystalline cubic alloy having a lattice constant of nearly 8.4 Å and a composition close to MnGa2 with a small amount of As. Of the remaining diffusion sources (Mn, MnAs, and Mn3As), only MnAs consistently produces a uniform doping distribution and smooth surface morphology. For diffusions at 800 °C, a uniform surface hole carrier concentration as high as 1020/cm3 can be obtained using MnAs as the source. The As overpressure is found to drastically alter the Mn diffusion profile, and Mn, like Zn, may diffuse in GaAs interstitial-substitutionally. Vapor from both the Mn and Mn3As solids degrade the GaAs surface. Mn3As, however, uncharacteristically degrades the surface more rapidly although the details of such are not well understood. With the presence of a high As overpressure, however, both surfaces of the Mn and Mn3As sources are converted to (Mn,As) compounds, the compositions being close to MnAs. High enough As overpressures are shown to completely suppress the GaAs surface degradation which is evident when Mn3As alone is used as the diffusion source.
Yang, Defu; Chen, Xueli; Cao, Xu; Wang, Jing; Liang, Jimin; Tian, Jie
2015-09-01
The third-order simplified harmonic spherical approximation (SP3) and diffusion approximation (DA) equations have been widely used in the three-dimensional (3D) whole-body optical imaging of small animals. With different types of tissues, which were classified by the ratio of µ s'/µ ɑ, the two equations have their own application scopes. However, the classification criterion was blurring and unreasonable, and the scope has not been systematically investigated until now. In this study, a new criterion for classifying tissues was established based on the absolute value of absorption and reduced scattering coefficients. Using the newly defined classification criterion, the performance and applicability of the SP3 and DA equations were evaluated with a series of investigation experiments. Extensive investigation results showed that the SP3 equation exhibited a better performance and wider applicability than the DA one in most of the observed cases, especially in tissues of low-scattering-low-absorption and low-scattering-high-absorption range. For the case of tissues with the high-scattering-low-absorption properties, a similar performance was observed for both the SP3 and the DA equations, in which case the DA was the preferred option for 3D whole-body optical imaging. Results of this study would provide significant reference for the study of hybrid light transport models. PMID:25850985
NASA Astrophysics Data System (ADS)
Hofman, Radek; Seibert, Petra; Kovalets, Ivan; Andronopoulos, Spyros
2015-04-01
We are concerned with source term retrieval in the case of an accident in a nuclear power with off-site consequences. The goal is to optimize atmospheric dispersion model inputs using inverse modeling of gamma dose rate measurements (instantaneous or time-integrated). These are the most abundant type of measurements provided by various radiation monitoring networks across Europe and available continuously in near-real time. Usually, a source term of an accidental release comprises of a mixture of nuclides. Unfortunately, gamma dose rate measurements do not provide a direct information on the source term composition; however, physical properties of respective nuclides (deposition properties, decay half-life) can yield some insight. In the method presented, we assume that nuclide ratios are known at least approximately, e.g. from nuclide specific observations or reactor inventory and assumptions on the accident type. The source term can be in multiple phases, each being characterized by constant nuclide ratios. The method is an extension of a well-established source term inversion approach based on the optimization of an objective function (minimization of a cost function). This function has two quadratic terms: mismatch between model and measurements weighted by an observation error covariance matrix and the deviation of the solution from a first guess weighted by the first-guess error covariance matrix. For simplicity, both error covariance matrices are approximated as diagonal. Analytical minimization of the cost function leads to a liner system of equations. Possible negative parts of the solution are iteratively removed by the means of first guess error variance reduction. Nuclide ratios enter the problem in the form of additional linear equations, where the deviations from prescribed ratios are weighted by factors; the corresponding error variance allows us to control how strongly we want to impose the prescribed ratios. This introduces some freedom into the
Skyglow: a retrieval of the approximate radiant intensity function of ground-based light sources
NASA Astrophysics Data System (ADS)
Kocifaj, M.; Solano Lamphar, H. A.
2014-04-01
The angular distribution of the light emitted from a city is an important source of information about public lighting systems and it also plays a key role in modelling the skyglow. Usually, the upwardly directed radiation is characterized through a parametrized emission function - a semi-empirical approach as a reasonable approximation that allows for fast computations. However, theoretical or experimental retrievals of emission characteristics are extremely difficult to obtain because of both the complexity of radiative transfer methods and/or the lack of highly specialized measuring devices. Our research has been conducted with the specific objective to identify an efficient theoretical technique for retrieval of the emission pattern of ground-based light sources in order to determine the optimum values of the scaling parameters of the Garstang function. In particular, the input data involve the zenith luminance or radiance with horizontal illuminance or irradiance. Theoretical ratios of zenith luminance LV(0) to horizontal illuminance DV are calculated for a set of distances d that separate a hypothetical observer from the light source (a city or town). This approach is advantageous because inexpensive traditional equipment can be used to obtain the mean values of the Garstang parameters. Furthermore, it can also be applied to other parametrizable emission functions and to any measuring site, even one with a masked horizon.
Source-position transformation: an approximate invariance in strong gravitational lensing
NASA Astrophysics Data System (ADS)
Schneider, Peter; Sluse, Dominique
2014-04-01
The main obstacle that gravitational lensing has in determining accurate masses of deflectors, or in determining precise estimates for the Hubble constant, is the degeneracy of lensing observables with respect to the mass-sheet transformation (MST). The MST is a global modification of the mass distribution which leaves all image positions, shapes, and flux ratios invariant, but which changes the time delay. Here we show that another global transformation of lensing mass distributions exists which leaves image positions and flux ratios almost invariant, and of which the MST is a special case. As is the case for the MST, this new transformation only applies if one considers only those source components that are at the same distance from us. Whereas for axi-symmetric lenses this source position transformation exactly reproduces all strong lensing observables, it does so only approximately for more general lens situations. We provide crude estimates for the accuracy with which the transformed mass distribution can reproduce the same image positions as the original lens model, and present an illustrative example of its performance. This new invariance transformation is most likely the reason why the same strong lensing information can be accounted for with rather different mass models.
NASA Astrophysics Data System (ADS)
Kopka, P.; Wawrzynczak, A.; Borysiewicz, M.
2015-09-01
In many areas of application, a central problem is a solution to the inverse problem, especially estimation of the unknown model parameters to model the underlying dynamics of a physical system precisely. In this situation, the Bayesian inference is a powerful tool to combine observed data with prior knowledge to gain the probability distribution of searched parameters. We have applied the modern methodology named Sequential Approximate Bayesian Computation (S-ABC) to the problem of tracing the atmospheric contaminant source. The ABC is technique commonly used in the Bayesian analysis of complex models and dynamic system. Sequential methods can significantly increase the efficiency of the ABC. In the presented algorithm, the input data are the on-line arriving concentrations of released substance registered by distributed sensor network from OVER-LAND ATMOSPHERIC DISPERSION (OLAD) experiment. The algorithm output are the probability distributions of a contamination source parameters i.e. its particular location, release rate, speed and direction of the movement, start time and duration. The stochastic approach presented in this paper is completely general and can be used in other fields where the parameters of the model bet fitted to the observable data should be found.
Boitard, Simon; Loisel, Patrice
2007-05-01
The probability distribution of haplotype frequencies in a population, and the way it is influenced by genetical forces such as recombination, selection, random drift ...is a question of fundamental interest in population genetics. For large populations, the distribution of haplotype frequencies for two linked loci under the classical Wright-Fisher model is almost impossible to compute because of numerical reasons. However the Wright-Fisher process can in such cases be approximated by a diffusion process and the transition density can then be deduced from the Kolmogorov equations. As no exact solution has been found for these equations, we developed a numerical method based on finite differences to solve them. It applies to transient states and models including selection or mutations. We show by several tests that this method is accurate for computing the conditional joint density of haplotype frequencies given that no haplotype has been lost. We also prove that it is far less time consuming than other methods such as Monte Carlo simulations. PMID:17316725
NASA Astrophysics Data System (ADS)
Mourenas, D.; Artemyev, A. V.; Agapitov, O. V.; Krasnoselskikh, V.; Li, W.
2014-12-01
The distribution of trapped energetic electrons inside the Earth's radiation belts is the focus of intense studies aiming at better describing the evolution of the space environment in the presence of various disturbances induced by the solar wind or by an enhanced lightning activity. Such studies are usually performed by means of comparisons with full numerical simulations solving the Fokker-Planck quasi-linear diffusion equation for the particle distribution function. Here we present for the first time approximate but realistic analytical solutions for the electron distribution, which are shown to be in good agreement with exact numerical solutions in situations where resonant scattering of energetic electrons by whistler mode hiss, lightning-generated or chorus waves, is the dominant process. Quiet time distributions are well recovered, as well as the evolution of energized relativistic electron distributions during disturbed geomagnetic conditions. It is further shown that careful comparisons between the analytical solutions and measured distributions may allow to infer important bounce- and drift-averaged wave characteristics (such as wave amplitude). It could also help to improve the global understanding of underlying physical phenomena.
Location of the effective diffusing-photon source in a strongly scattering medium.
Kostko, A F; Pavlov, V A
1997-10-20
When a narrow laser beam illuminates a strongly scattering medium, the effective pointlike source of diffusing photons appears inside the medium. By the method worked out, which is based on measurements of the diffusive intensity of light emerging from a turbid spherical sample, the depth of this source site (the penetration depth) is determined relatively to the sample diameter, which is known accurately. By using this method of locating the effective source, we have discovered that its position inside the medium is unexpectedly deep. We obtained the penetration depth D(0) = 4.6 l* +/- 0.7 l* instead of one transport mean free path, where l* is the value of D(0) in the standard diffusion theory. Information about this source dipping is useful in diffusing-photon correlation spectroscopy because of its influence on the geometric factor calculated from the diffusion equation. PMID:18264271
Thermal diffusivity measurements of particulates using the differentiated line source
NASA Technical Reports Server (NTRS)
West, E. A.; Fountain, J. A.
1975-01-01
A method is described in which thermal diffusivity measurements can be made on particulate materials using the identical instrumentation as described in previous papers for measuring thermal conductivity. The measurements for the two properties can be made simultaneously, thus eliminating the changes in conditions when they are made separately. This system has particular application for studies of simulated lunar, planetary, or asteroidal surfaces in which laboratory measurements can be correlated with astronomical observations of thermal inertia. A representative set of data is shown which gives thermal diffusivity and thermal conductivity measurements on a particulate basalt in which the values for each temperature were taken simultaneously. A value for specific heat is calculated for measurements taken at each temperature.
Trivett, D.A.
1994-09-01
This paper focuses on modeling the fate of effluent from diffuse seafloor hydrothermal activity after it has been vented into the water column. The model was formulated using a number of simplifying assumptions which permit direct application of this model to field measurements. I have limited the configurations to those where the hydrothermal outflow velocities are smaller than horizontal current. I assume that the entrainment of ambient seawater into the plume is constant over the length of the plume. This permits formulation of a first-order relation for the rise height and dilution in a diffuse hydrothermal plume as a function of downstream distance. The analytic model is compared with a simple laboratory simulation of the hydrothermal flow. The results suggest that diffuse hydrothermal effluent will penetrate to a height in the water column that is proportional to the overall dimension of the diffuse vent patch, multiplied by a dimensionless plume intensity parameter. I also ahow relations for plume dilution which will be compared with field data in part 2 of this work.
NASA Astrophysics Data System (ADS)
Beatty, Thomas G.; Gaudi, B. Scott
2015-12-01
We investigate various astrophysical contributions to the statistical uncertainty of precision radial velocity measurements of stellar spectra. We first analytically determine the intrinsic uncertainty in centroiding isolated spectral lines broadened by Gaussian, Lorentzian, Voigt, and rotational profiles, finding that for all cases and assuming weak lines, the uncertainty in the line centroid is σV ≈ C\\Theta3/2/(WI1/20), where Θ is the full-width at half-maximum of the line, W is the equivalent width, and I0 is the continuum signal-to-noise ratio, with C a constant of order unity that depends on the specific line profile. We use this result to motivate approximate analytic expressions to the total radial velocity uncertainty for a stellar spectrum with a given photon noise, resolution, wavelength, effective temperature, surface gravity, metallicity, macroturbulence, and stellar rotation. We use these relations to determine the dominant contributions to the statistical uncertainties in precision radial velocity measurements as a function of effective temperature and mass for main-sequence stars. For stars more massive than ~1.1 Msolar we find that stellar rotation dominates the velocity uncertainties for moderate and high-resolution spectra (R gsim 30,000). For less-massive stars, a variety of sources contribute depending on the spectral resolution and wavelength, with photon noise due to decreasing bolometric luminosity generally becoming increasingly important for low-mass stars at fixed exposure time and distance. In most cases, resolutions greater than 60,000 provide little benefit in terms of statistical precision, although higher resolutions would likely allow for better control of systematic uncertainties. We find that the spectra of cooler stars and stars with higher metallicity are intrinsically richer in velocity information, as expected. We determine the optimal wavelength range for stars of various spectral types, finding that the optimal region
Modeling diffuse sources of surface water contamination with plant protection products
NASA Astrophysics Data System (ADS)
Wendland, Sandra; Bock, Michael; Böhner, Jürgen; Lembrich, David
2015-04-01
Entries of chemical pollutants in surface waters are a serious environmental problem. Among water pollutants plant protection products (ppp) from farming practice are of major concern not only for water suppliers and environmental agencies, but also for farmers and industrial manufacturers. Lost chemicals no longer fulfill their original purpose on the field, but lead to severe damage of the environment and surface waters. Besides point-source inputs of chemical pollutants, the diffuse-source inputs from agricultural procedures play an important and not yet sufficiently studied role concerning water quality. The two most important factors for diffuse inputs are erosion and runoff. The latter usually occurs before erosion begins, and is thus often not visible in hindsight. Only if it has come to erosion, it is obvious to expect runoff in foresight at this area, too. In addition to numerous erosion models, there are also few applications to model runoff processes available. However, these conventional models utilize approximations of catchment parameters based on long-term average values or theoretically calculated concentration peaks which can only provide indications to relative amounts. Our study aims to develop and validate a simplified spatially-explicit dynamic model with high spatiotemporal resolution that enables to measure current and forecast runoff potential not only at catchment scale but field-differentiated. This method allows very precise estimations of runoff risks and supports risk reduction measures to be targeted before fields are treated. By focusing on water pathways occurring on arable land, targeted risk reduction measures like buffer strips at certain points and adapted ppp use can be taken early and pollution of rivers and other surface waters through transported pesticides, fertilizers and their products could be nearly avoided or largely minimized. Using a SAGA-based physical-parametric modeling approach, major factors influencing runoff
Mimicking diffuse supernova antineutrinos with the sun as a source
Raffelt, G. G.; Rashba, T. I.
2010-04-15
Measuring the {nu}-bar{sub e} component of the cosmic diffuse supernova neutrino background (DSNB) is the next ambitious goal for low-energy neutrino astronomy. The largest flux is expected in the lowest accessible energy bin. However, for E {<=} 15 MeV a possible signal can be mimicked by a solar {nu}-bar{sub e} flux that originates from the usual {sup 8}B neutrinos by spin-flavor oscillations. We show that such an interpretation is possible within the allowed range of neutrino electromagnetic transition moments and solar turbulent field strengths and distributions. Therefore, an unambiguous detection of the DSNB requires a significant number of events at E {>=} 15 MeV.
NASA Astrophysics Data System (ADS)
Moura, R. C.; Sherwin, S. J.; Peiró, J.
2016-02-01
This study addresses linear dispersion-diffusion analysis for the spectral/hp continuous Galerkin (CG) formulation in one dimension. First, numerical dispersion and diffusion curves are obtained for the advection-diffusion problem and the role of multiple eigencurves peculiar to spectral/hp methods is discussed. From the eigencurves' behaviour, we observe that CG might feature potentially undesirable non-smooth dispersion/diffusion characteristics for under-resolved simulations of problems strongly dominated by either convection or diffusion. Subsequently, the linear advection equation augmented with spectral vanishing viscosity (SVV) is analysed. Dispersion and diffusion characteristics of CG with SVV-based stabilization are verified to display similar non-smooth features in flow regions where convection is much stronger than dissipation or vice-versa, owing to a dependency of the standard SVV operator on a local Péclet number. First a modification is proposed to the traditional SVV scaling that enforces a globally constant Péclet number so as to avoid the previous issues. In addition, a new SVV kernel function is suggested and shown to provide a more regular behaviour for the eigencurves along with a consistent increase in resolution power for higher-order discretizations, as measured by the extent of the wavenumber range where numerical errors are negligible. The dissipation characteristics of CG with the SVV modifications suggested are then verified to be broadly equivalent to those obtained through upwinding in the discontinuous Galerkin (DG) scheme. Nevertheless, for the kernel function proposed, the full upwind DG scheme is found to have a slightly higher resolution power for the same dissipation levels. These results show that improved CG-SVV characteristics can be pursued via different kernel functions with the aid of optimization algorithms.
Drabina, A; Dwora, D; Loskiewicz, J; Zorski, T
2001-05-01
Czubek has obtained a practically successful solution of a two-layer borehole geometry using the neutron diffusion approximation. Czubek and Woznicka solved the more realistic three-layer borehole model, which includes a middle layer consisting of mud-cake or steel tubing, in this approximation. The comparison with experimental data for 15-mm thick steel tubes shows some discrepancies. By using an exact Monte Carlo solution of the Boltzmann equation for calculating the migration length, measurements agree more closely with the "general calibration curve". This opens opportunities for cased borehole logging. PMID:11258536
COS-B gamma ray sources beyond the predicted diffuse emission
NASA Technical Reports Server (NTRS)
Mayer-Hasselwander, H. A.; Simpson, G.
1990-01-01
COS-B data were reanalyzed using for background subtraction the modeled galactic diffuse gamma-ray emission based on HI- and CO-line surveys and the gamma-ray data itself. A methodology was developed for this purpose with the following three features: automatic generation of source catalogs using correlation analysis, simulation of trials to derive significance thresholds for source detection, and bootstrap sampling to drive error boxes and confidence intervals for source parameters. The analysis shows that about half of the 2CG sources are explained by concentrations in the distribution of molecular hydrogen. Indication for a few weak new sources is also obtained.
Broadband signal generator for the approximation of a magnetotelluric source for indoor testing
NASA Astrophysics Data System (ADS)
Ge, Shuang-chao; Deng, Ming; Chen, Kai; Shi, Xin-yu
2016-08-01
To test the frequency response of a magnetotelluric (MT) receiver, a broadband source, especially white noise is more efficient and intuitive than single frequency signals. In view of the absence of an appropriate source generator for MT receiver indoor testing, we designed a broadband signal generator based on a pseudo-random binary sequence (PRBS). Firstly, we divided the whole MT band into two segments to avoid data redundancy and simplify calculation in data processing and designed a generator composed of several modules: a clock module, a PRBS logic module, and a voltage level conversion module. We conducted a detailed analysis of the optimal parameter selection methods for each module, and key parameters including clock frequency, order, the primitive polynomial and the original states of the linear registers were determined. The generator provides four-channel PRBS signals with two effective bandwidths of 5 × 10‑4–714 Hz and 0.1 Hz–14 kHz which are broad enough to cover the frequency range for different MT methods. These four-channel signals were used to simulate two modes of sources (xy and yx) with strong auto-correlation and weak cross-correlation. The power spectral density is quite stable in the whole passband. The new generator is characterized by broadband output in low-frequency bands, low power consumption, simple operation and reliable performance. Indoor and field tests indicated that the generator can provide an analog MT source and is a practical tool for MT receiver indoor testing.
Diffuse and point sources of silica in the Seine River watershed.
Sferratore, Agata; Garnier, Josette; Billen, Gilles; Conley, Daniel I; Pinault, Séverine
2006-11-01
Dissolved silica (DSi) is believed to enter aquatic ecosystems primarily through diffuse sources by weathering. Point sources have generally been considered negligible, although recent reports of DSi inputs from domestic and industrial sources suggest otherwise. In addition, particulate amorphous silica (ASi) inputs from terrestrial ecosystems during soil erosion and in vegetation can dissolve and also be a significant source of DSi. We quantify here both point and diffuse sources of DSi and particulate ASi to the Seine River watershed. The total per capita point source inputs of Si (DSi + ASi) were found to be 1.0 and 0.8 g Si inhabitant(-1) d(-1) in raw and treated waters of the Achères wastewater treatment plant, in agreement with calculations based on average food intake and silica-containing washing products consumption. A mass balance of Si inputs and outputs for the Seine drainage network was established for wet and dry hydrological conditions (2001 and 2003, respectively). Diffuse sources of Si are of 1775 kg Si km(-2) y(-1) in wet conditions and 762 kg Si km(-2) y(-1) in dry conditions, with the proportion of ASi around 6%. Point sources of Si from urban discharge can contribute to more than 8% of the total Si inputs at the basin scale in hydrologically dry years. An in-stream retention of 6% of total inputs in dry conditions and 12% in wet conditions is inferred from the budget. PMID:17144288
NASA Astrophysics Data System (ADS)
Baran, Timothy M.
2016-07-01
We demonstrate recovery of optical properties using arrays of interstitial cylindrical diffusing fibers as sources and detectors. A single 1-cm diffuser delivered laser illumination at 665 nm, while seven 1- and 2-cm diffusers at 1-cm grid spacing acted as detectors. Extraction of optical properties from these measurements was based upon a diffusion model of emission and detection distributions for these diffuser fibers, informed by previous measurements of heterogeneous axial detection. Verification of the technique was performed in 15 liquid tissue-simulating phantoms consisting of deionized water, India ink as absorber, and Intralipid 20% as scatterer. For the range of optical properties tested, mean errors were 4.4% for effective attenuation coefficient, 12.6% for absorption coefficient, and 7.6% for reduced scattering coefficient. Error in recovery tended to increase with decreasing transport albedo. For therapeutic techniques involving the delivery of light to locations deep within the body, such as interstitial photodynamic and photothermal therapies, the methods described here would allow the treatment diffuser fibers also to be used as sources and detectors for recovery of optical properties. This would eliminate the need for separately inserted fibers for spectroscopy, reducing clinical complexity and improving the accuracy of treatment planning.
Roebeling, P C; Cunha, M C; Arroja, L; van Grieken, M E
2015-01-01
Marine ecosystems are affected by water pollution originating from coastal catchments. The delivery of water pollutants can be reduced through water pollution abatement as well as water pollution treatment. Hence, sustainable economic development of coastal regions requires balancing of the marginal costs from water pollution abatement and/or treatment and the associated marginal benefits from marine resource appreciation. Water pollution delivery reduction costs are, however, not equal across abatement and treatment options. In this paper, an optimal control approach is developed and applied to explore welfare maximizing rates of water pollution abatement and/or treatment for efficient diffuse source water pollution management in terrestrial-marine systems. For the case of diffuse source dissolved inorganic nitrogen water pollution in the Tully-Murray region, Queensland, Australia, (agricultural) water pollution abatement cost, (wetland) water pollution treatment cost and marine benefit functions are determined to explore welfare maximizing rates of water pollution abatement and/or treatment. Considering partial (wetland) treatment costs and positive water quality improvement benefits, results show that welfare gains can be obtained, primarily, through diffuse source water pollution abatement (improved agricultural management practices) and, to a minor extent, through diffuse source water pollution treatment (wetland restoration). PMID:26287831
Required distribution of noise sources for Green's function recovery in diffusive fields
NASA Astrophysics Data System (ADS)
Shamsalsadati, S.; Weiss, C. J.
2011-12-01
In the most general sense, noise is the part of the signal of little or no interest, due to a multitude of reasons such as operator error, imperfect instrumentation, experiment design, or inescapable background interference. Considering the latter, it has been shown that Green's function can be extracted from cross-correlation of the ambient, diffusive wavefields arising from background random noise sources. Pore pressure and low-frequency electromagnetic induction are two such examples of diffusive fields. In theory, applying Green's function method in geophysical exploration requires infinity of volumetrically distributed sources; however, in the real world the number of noise sources in an area is limited, and furthermore, unevenly distributed in time, space and spectral content. Hence, quantification of the requisite noise sources that enable us to calculate Green's function acceptably well remains an open research question. The purpose of this study is to find the area of noise sources that contribute most to the Green's function estimation in diffusive systems. We call such a region the Volume of Relevance (VoR). Our analysis builds upon recent work in 1D homogeneous system where it was shown that sources located between two receivers positions are the most important ones for the purpose of Green's function recovery. Our results confirm the previous finding but we also examine the effect of heterogeneity, dimensionality and receiver location in both 1D and 2D at a fixed frequency. We demonstrate that for receivers located symmetrically across an interface between regions of contrasting diffusivity, the VoR rapidly shifts from one side of the interface to the other, and back again, as receiver separation increases. We also demonstrate that where the receiver pair is located on the interface itself, the shifting is less rapid, and for moderate to high diffusivity contrasts, the VoR remains entirely on the more diffusive side. In addition, because classical
Three-dimensional analytical solutions of the atmospheric diffusion equation with multiple sources and height-dependent wind speed and eddy diffusivities are derived in a systematic fashion. For homogeneous Neumann (total reflection), Dirichlet (total adsorpti...
Three-dimensional analytical solutions of the atmospheric diffusion equation with multiple sources and height-dependent wind speed and eddy diffusivities are derived in a systematic fashion. For homogeneous Neumann (total reflection), Dirichlet (total adsorpti...
Diffusion of alkali species in porous tungsten substrates used in contact ionization sources
Chacon-Golcher, E.; Kwan, J.; Morse, E. C.
2002-01-01
Contact ionization (doped) sources used in current Heavy Ion Fusion experiments consist of a porous tungsten substrate doped with an alkali carbonate. In the early stages of the heating cycle (T - 600 C), the carbonate breaks down and releases the alkali atoms that then diffuse through the substrate. At the emitter surface there is a balance between the fast desorption rate of the alkali atoms (mostly as neutrals) and the slower replenishment rate from the substrate by diffusion. Time-resolved measurements of neutral particle evaporation rates at the emitter surface have been used to estimate the effective diffusion coefficient (D) that characterizes the migration of alkali species in the substrate. These estimates are consistent with the observed source lifetimes (tens of hrs.) and establish the alkali migration in the bulk as a diffusion-limited process. The measurements suggest that the faster migration rates (D {approx} 10{sup -5} - 10{sup -6} cm2/s) occur early during the heating cycle when the dominant species are the neutral alkali atoms. At operating temperatures there is a slower migration rate (D x 10-7cm2/sd) ue to the dominance of ions, which difise by a slower surface diffusion process.
Locating the source of diffusion in complex networks by time-reversal backward spreading
NASA Astrophysics Data System (ADS)
Shen, Zhesi; Cao, Shinan; Wang, Wen-Xu; Di, Zengru; Stanley, H. Eugene
2016-03-01
Locating the source that triggers a dynamical process is a fundamental but challenging problem in complex networks, ranging from epidemic spreading in society and on the Internet to cancer metastasis in the human body. An accurate localization of the source is inherently limited by our ability to simultaneously access the information of all nodes in a large-scale complex network. This thus raises two critical questions: how do we locate the source from incomplete information and can we achieve full localization of sources at any possible location from a given set of observable nodes. Here we develop a time-reversal backward spreading algorithm to locate the source of a diffusion-like process efficiently and propose a general locatability condition. We test the algorithm by employing epidemic spreading and consensus dynamics as typical dynamical processes and apply it to the H1N1 pandemic in China. We find that the sources can be precisely located in arbitrary networks insofar as the locatability condition is assured. Our tools greatly improve our ability to locate the source of diffusion in complex networks based on limited accessibility of nodal information. Moreover, they have implications for controlling a variety of dynamical processes taking place on complex networks, such as inhibiting epidemics, slowing the spread of rumors, pollution control, and environmental protection.
Locating the source of diffusion in complex networks by time-reversal backward spreading.
Shen, Zhesi; Cao, Shinan; Wang, Wen-Xu; Di, Zengru; Stanley, H Eugene
2016-03-01
Locating the source that triggers a dynamical process is a fundamental but challenging problem in complex networks, ranging from epidemic spreading in society and on the Internet to cancer metastasis in the human body. An accurate localization of the source is inherently limited by our ability to simultaneously access the information of all nodes in a large-scale complex network. This thus raises two critical questions: how do we locate the source from incomplete information and can we achieve full localization of sources at any possible location from a given set of observable nodes. Here we develop a time-reversal backward spreading algorithm to locate the source of a diffusion-like process efficiently and propose a general locatability condition. We test the algorithm by employing epidemic spreading and consensus dynamics as typical dynamical processes and apply it to the H1N1 pandemic in China. We find that the sources can be precisely located in arbitrary networks insofar as the locatability condition is assured. Our tools greatly improve our ability to locate the source of diffusion in complex networks based on limited accessibility of nodal information. Moreover, they have implications for controlling a variety of dynamical processes taking place on complex networks, such as inhibiting epidemics, slowing the spread of rumors, pollution control, and environmental protection. PMID:27078360
Baran, Timothy M. Foster, Thomas H.
2014-02-15
Purpose: For interstitial photodynamic therapy (iPDT) of bulky tumors, careful treatment planning is required in order to ensure that a therapeutic dose is delivered to the tumor, while minimizing damage to surrounding normal tissue. In clinical contexts, iPDT has typically been performed with either flat cleaved or cylindrical diffusing optical fibers as light sources. Here, the authors directly compare these two source geometries in terms of the number of fibers and duration of treatment required to deliver a prescribed light dose to a tumor volume. Methods: Treatment planning software for iPDT was developed based on graphics processing unit enhanced Monte Carlo simulations. This software was used to optimize the number of fibers, total energy delivered by each fiber, and the position of individual fibers in order to deliver a target light dose (D{sub 90}) to 90% of the tumor volume. Treatment plans were developed using both flat cleaved and cylindrical diffusing fibers, based on tissue volumes derived from CT data from a head and neck cancer patient. Plans were created for four cases: fixed energy per fiber, fixed number of fibers, and in cases where both or neither of these factors were fixed. Results: When the number of source fibers was fixed at eight, treatment plans based on flat cleaved fibers required each to deliver 7180–8080 J in order to deposit 90 J/cm{sup 2} in 90% of the tumor volume. For diffusers, each fiber was required to deliver 2270–2350 J (333–1178 J/cm) in order to achieve this same result. For the case of fibers delivering a fixed 900 J, 13 diffusers or 19 flat cleaved fibers at a spacing of 1 cm were required to deliver the desired dose. With energy per fiber fixed at 2400 J and the number of fibers fixed at eight, diffuser fibers delivered the desired dose to 93% of the tumor volume, while flat cleaved fibers delivered this dose to 79%. With both energy and number of fibers allowed to vary, six diffusers delivering 3485–3600 J
Turbulent diffusion from a heated line source in non-equilibrium grid turbulence
NASA Astrophysics Data System (ADS)
Nedic, Jovan; Tavoularis, Stavros
2015-11-01
We have investigated turbulent diffusion of heat injected passively from a line source in equilibrium and non-equilibrium grid-generated turbulence, which are, respectively, flows in which the value of the non-dimensional rate of kinetic energy dissipation is constant or changes with streamwise distance from the grid. We used three grids with uniform square meshes and one fractal square grid (FSG), all of the same solidity, to generate non-equilibrium and equilibrium turbulence in a wind-tunnel. The regular grids have mesh sizes that are comparable to the first (RG160), second (RG80) and fourth (RG18) iterations of the fractal grid. The heated line source was inserted on the centre-plane of the grids at either of two downstream locations or an upstream one and it spanned the entire width of the wind-tunnel. We found that RG160 produced the greatest heat diffusion, followed by FSG, RG80 and RG18, in this order. The apparent turbulent diffusivity produced by the four grids also decreased in the same order. These findings conform with Taylor's theory of diffusion by continuous movements. Moreover, the present study demonstrates that the fractal space-scale unfolding (SSU) mechanism does not apply to grids with the same solidity but different effective mesh sizes. Supported by NSERC.
Aldaba, D; García-Gutiérrez, M; Rigol, A; Vidal, M
2010-11-01
A planar-source method, initially designed to obtain diffusion coefficients in compacted clay, is adapted here to determine the apparent diffusion coefficient (D(a)) of radiostrontium in soils representative of the Spanish territory. Experiments were carried out by varying the moisture content (F(moist)), and bulk dry density (ρ(bulk)) of the soil samples, in order to study the influence of these soil packing parameters on D(a) values. The moisture in the soil samples was established as the percentage of occupancy of each soil's field capacity (OFC). For a similar OFC, D(a) values in the examined soils ranged by approximately one order of magnitude (e.g. from 6.2 × 10(-)(11) to 6.5 × 10(-)(12)m(2)s(-)(1), at 100% of OFC; from 3.0 × 10(-)(11) to 3.8 × 10(-)(12)m(2)s(-)(1), at 60% of OFC). For a given soil, D(a) values increased when water content was increased. F(moist), and tortuosity (τ) explained D(a) variability, with R(2) values usually over 0.9. However, no good simple or multiple regressions between the soil packing parameters and D(a) were obtained with the whole dataset of all soils, which indicated that soil sorption capacity affects the diffusion of reactive radionuclides in soils. The inclusion of calculated K(d) values in the multiple regressions improved the correlations in all cases. Finally, D(a) values were compared with those obtained by the application of a half-cell method. The values of D(a) obtained by the planar-source methods were systematically lower than the half-cell ones, with a good correlation between the D(a) derived from both methods (R(2)=0.98). PMID:20850167
Anisotropies of ultrahigh energy cosmic ray nuclei diffusing from extragalactic sources
NASA Astrophysics Data System (ADS)
Harari, Diego; Mollerach, Silvia; Roulet, Esteban
2015-09-01
We obtain the dipolar anisotropies in the arrival directions of ultrahigh energy cosmic ray nuclei diffusing from nearby extragalactic sources. We consider mixed-composition scenarios in which different cosmic ray nuclei are accelerated up to the same maximum rigidity, so that E
On the source function of the soft X-ray diffuse background
NASA Technical Reports Server (NTRS)
Burrows, David N.; Kraft, Ralph P.
1993-01-01
Radiation transfer theory has been used recently to derive the source function of the soft X-ray diffuse background, resulting in the claim of evidence for 10 exp 6 K gas in the Galactic halo. We show that this analysis has several errors that invalidate its conclusions. We argue that the case for an extensive hot halo remains open, pending further work, but may be settled by the continuing series of Rosat observations of high-latitude soft X-ray shadows.
Nutrient pollution in central California rivers: a comparison of point and diffuse sources
NASA Astrophysics Data System (ADS)
Sobota, D. J.; Harrison, J. A.
2009-12-01
Eutrophication of aquatic ecosystems resulting from nutrient pollution is a persistent problem in regions where agriculture and urban development are significant. An important question for ecosystem management in these areas is whether nutrient pollution originates from point sources (e.g., sewage effluent) or diffuse sources (e.g., runoff of fertilizers from agricultural fields). In this poster, we present spatially explicit estimates for the input of two key constituents of nutrient pollution, nitrogen (N) and phosphorus (P), to 22 individual river basins in central California for the early 2000s using information on land cover, agricultural practices, and sewage treatment facilities. Net anthropogenic input (agricultural sources + urban sources - crop harvest) of N ranged from 194 to 11094 kg N km-2 yr-1 while net anthropogenic P inputs ranged from -22 to 346 kg P km-2 yr-1. We estimated that diffuse inputs of inorganic fertilizer and manure were the largest fraction of anthropogenic N and P inputs to these basins. Point source inputs from sewage effluent was nonexistent in 17 of the basins and contributed up to 2 and 15% of N and P inputs in the remaining 5 basins. Using measured data on riverine N and P export from these basins during the 2000 - 2003 water years, we estimate that these basins exported a median of 11 and 12% of their annual net anthropogenic N and P inputs. Among basins, export of net anthropogenic P inputs was on average 53% larger than export of net anthropogenic N inputs. This suggests either that N was retained more efficiently than P or it was lost from basins in gaseous form (e.g., N2 and N2O). Our results suggest that ecosystem management in these central California river basins should focus on diffuse agricultural sources of both N and P to reduce nutrient pollution of downstream aquatic ecosystems.
NASA Astrophysics Data System (ADS)
Strahm, Ivo; Munz, Nicole; Braun, Christian; Gälli, René; Leu, Christian; Stamm, Christian
2014-05-01
Water quality in the Swiss river network is affected by many micropollutants from a variety of diffuse sources. This study compares, for the first time, in a comprehensive manner the diffuse sources and the substance groups that contribute the most to water contamination in Swiss streams and highlights the major regions for water pollution. For this a simple but comprehensive model was developed to estimate emission from diffuse sources for the entire Swiss river network of 65 000 km. Based on emission factors the model calculates catchment specific losses to streams for more than 15 diffuse sources (such as crop lands, grassland, vineyards, fruit orchards, roads, railways, facades, roofs, green space in urban areas, landfills, etc.) and more than 130 different substances from 5 different substance groups (pesticides, biocides, heavy metals, human drugs, animal drugs). For more than 180 000 stream sections estimates of mean annual pollutant loads and mean annual concentration levels were modeled. This data was validated with a set of monitoring data and evaluated based on annual average environmental quality standards (AA-EQS). Model validation showed that the estimated mean annual concentration levels are within the range of measured data. Therefore simulations were considered as adequately robust for identifying the major sources of diffuse pollution. The analysis depicted that in Switzerland widespread pollution of streams can be expected. Along more than 18 000 km of the river network one or more simulated substances has a concentration exceeding the AA-EQS. In single stream sections it could be more than 50 different substances. Moreover, the simulations showed that in two-thirds of small streams (Strahler order 1 and 2) at least one AA-EQS is always exceeded. The highest number of substances exceeding the AA-EQS are in areas with large fractions of arable cropping, vineyards and fruit orchards. Urban areas are also of concern even without considering
Discrete-continuous reaction-diffusion model with mobile point-like sources and sinks.
Kondrat, Svyatoslav; Zimmermann, Olav; Wiechert, Wolfgang; von Lieres, Eric
2016-01-01
In many applications in soft and biological physics, there are multiple time and length scales involved but often with a distinct separation between them. For instance, in enzyme kinetics, enzymes are relatively large, move slowly and their copy numbers are typically small, while the metabolites (being transformed by these enzymes) are often present in abundance, are small in size and diffuse fast. It seems thus natural to apply different techniques to different time and length levels and couple them. Here we explore this possibility by constructing a stochastic-deterministic discrete-continuous reaction-diffusion model with mobile sources and sinks. Such an approach allows in particular to separate different sources of stochasticity. We demonstrate its application by modelling enzyme-catalysed reactions with freely diffusing enzymes and a heterogeneous source of metabolites. Our calculations suggest that using a higher amount of less active enzymes, as compared to fewer more active enzymes, reduces the metabolite pool size and correspondingly the lag time, giving rise to a faster response to external stimuli. The methodology presented can be extended to more complex systems and offers exciting possibilities for studying problems where spatial heterogeneities, stochasticity or discreteness play a role. PMID:26830760
Shumaker, D E; Woodward, C S
2005-05-03
In this paper, the authors investigate performance of a fully implicit formulation and solution method of a diffusion-reaction system modeling radiation diffusion with material energy transfer and a fusion fuel source. In certain parameter regimes this system can lead to a rapid conversion of potential energy into material energy. Accuracy in time integration is essential for a good solution since a major fraction of the fuel can be depleted in a very short time. Such systems arise in a number of application areas including evolution of a star and inertial confinement fusion. Previous work has addressed implicit solution of radiation diffusion problems. Recently Shadid and coauthors have looked at implicit and semi-implicit solution of reaction-diffusion systems. In general they have found that fully implicit is the most accurate method for difficult coupled nonlinear equations. In previous work, they have demonstrated that a method of lines approach coupled with a BDF time integrator and a Newton-Krylov nonlinear solver could efficiently and accurately solve a large-scale, implicit radiation diffusion problem. In this paper, they extend that work to include an additional heating term in the material energy equation and an equation to model the evolution of the reactive fuel density. This system now consists of three coupled equations for radiation energy, material energy, and fuel density. The radiation energy equation includes diffusion and energy exchange with material energy. The material energy equation includes reaction heating and exchange with radiation energy, and the fuel density equation includes its depletion due to the fuel consumption.
NASA Astrophysics Data System (ADS)
Kim, Eugene D.; Murgai, Vijay; Menzel, Reinhard W.
2012-09-01
The Visible Infrared Imaging Radiometer Suite (VIIRS) sensor on the Joint Polar-orbiting Satellite System (JPSS) mission has a solar diffuser as a reflective band calibrator. Due to UV solarization of the solar diffuser, the Solar Diffuser Stability Monitor (SDSM) is on-board to track the reflectance change of the solar diffuser in visible to near IR wavelengths. A 100 cm Sphere Integrating Source (SIS) has been used to configure and test the SDSM on the ground since MODerate resolution Imaging Spectroradiometer (MODIS) programs. Recent upgrades of the radiance transfer and BRDF measurement instruments in Raytheon have enabled more spectral data and faster measurement time with comparable uncertainty to the previous methods. The SIS has a Radiance Monitor, which has been mainly used as a SIS real-time health checker. It has been observed that the Radiance Monitor response is sufficiently linear and stable thus the Radiance Monitor can be used as a calibrator for ground tests. This paper describes the upgraded SIS calibration instruments, and the changes in the calibration philosophy of the SIS for the SDSM bands.
Gamma-ray observations of Ophiuchus with EGRET: The diffuse emission and point sources
NASA Technical Reports Server (NTRS)
Hunter, S. D.; Digel, S. W.; De Geus, E. J.; Kanbach, G.
1994-01-01
Observations of the Ophiuchus region made with the Energetic Gamma-Ray Experiment Telescope (EGRET) on the Compton Gamma Ray Observatory (CGRO) during the first 2 1/2 years of operation show the diffuse emission from the interstellar gas in Ophiuchus as well as variable emission from two point sources. The gamma-ray emission is modeled in terms of cosmic-ray interactions with atomic and molecular hydrogen in Ophiuchus and with low-energy photons along the line of sight. The model also includes the flux from the two point sources and an isotropic diffuse contribution. The cosmic-ray density is assumed to be uniform. The derived ratio of molecular hydrogen column density to integrated CO intensity is (1.1 +/- 0.2) x 10(exp 20) H-mols/sq cm (K km/s)(exp -1). At the sensitivity and resolution of the gamma-ray data, no variation of this ratio over the modeled region is discernible, nor are any regions of enhanced cosmic-ray density apparent. The model was fitted to seven narrow energy bands to obtain the energy depedence of the gamma-ray production function and the spectra of the point sources. The derived production function is in good agreement with theoretical calculations and the local cosmic-ray electron and proton spectra. The positions of the point sources were determined from maximum likelihood analysis of the gamma-ray emission observed in excess of the diffuse model. We identify one point source with the quasar PKS 1622-253, which has an average flux, E greater than 100 MeV, of (2.5 +/- 0.5) x 10(exp -7) photons/sq cm/s and photon spectral index -1.9 +/- 0.3. The other source, denoted GRO J1631-27, has not yet been identified at other wavelengths. Its average flux, E greater than 100 MeV, is (1.1 +/- 0.4) x 10(exp -7) photons/sq cm/s; however, its spectral index is poorly determined. The spectral index and intensity of the isotropic contribution to the model agree well with the extragalactic diffuse emission derived from the SAS 2 data.
Si diffusion and intermixing in AlGaAs/GaAs structures using buried impurity sources
Beernink, K.J.; Thornton, R.L.; Anderson, G.B.; Emanuel, M.A.
1995-05-08
Si diffusion and impurity-induced layer intermixing from a buried impurity source have been studied by transmission electron microscopy and secondary ion mass spectroscopy of isolated, Si-doped GaAs layers in an undoped Al{sub 0.4}Ga{sub 0.6}As/GaAs superlattice and by photoluminescence measurements on Si-doped GaAs quantum wells with undoped Al{sub 0.4}Ga{sub 0.6}As barriers. In annealed samples, the Si profile suggests a Si diffusion process involving multiply ionized column III vacancies. The width of the resulting Si profile and the spatial extent and completeness of intermixing strongly depend on the initial Si concentration in the doped layer.
Analytic solutions of the time-dependent quasilinear diffusion equation with source and loss terms
Hassan, M.H.A. ); Hamza, E.A. )
1993-08-01
A simplified one-dimensional quasilinear diffusion equation describing the time evolution of collisionless ions in the presence of ion-cyclotron-resonance heating, sources, and losses is solved analytically for all harmonics of the ion cyclotron frequency. Simple time-dependent distribution functions which are initially Maxwellian and vanish at high energies are obtained and calculated numerically for the first four harmonics of resonance heating. It is found that the strongest ion tail of the resulting anisotropic distribution function is driven by heating at the second harmonic followed by heating at the fundamental frequency.
Impacts by point and diffuse micropollutant sources on the stream water quality at catchment scale
NASA Astrophysics Data System (ADS)
Petersen, M. F.; Eriksson, E.; Binning, P. J.; Bjerg, P. L.
2012-04-01
The water quality of surface waters is threatened by multiple anthropogenic pollutants and the large variety of pollutants challenges the monitoring and assessment of the water quality. The aim of this study was to characterize and quantify both point and diffuse sources of micropollutants impacting the water quality of a stream at catchment scale. Grindsted stream in western Jutland, Denmark was used as a study site. The stream passes both urban and agricultural areas and is impacted by severe groundwater contamination in Grindsted city. Along a 12 km reach of Grindsted stream, the potential pollution sources were identified including a pharmaceutical factory site with a contaminated old drainage ditch, two waste deposits, a wastewater treatment plant, overflow structures, fish farms, industrial discharges and diffuse agricultural and urban sources. Six water samples were collected along the stream and analyzed for general water quality parameters, inorganic constituents, pesticides, sulfonamides, chlorinated solvents, BTEXs, and paracetamol and ibuprofen. The latter two groups were not detected. The general water quality showed typical conditions for a stream in western Jutland. Minor impacts by releases of organic matter and nutrients were found after the fish farms and the waste water treatment plant. Nickel was found at concentrations 5.8 - 8.8 μg/l. Nine pesticides and metabolites of both agricultural and urban use were detected along the stream; among these were the two most frequently detected and some rarely detected pesticides in Danish water courses. The concentrations were generally consistent with other findings in Danish streams and in the range 0.01 - 0.09 μg/l; except for metribuzin-diketo that showed high concentrations up to 0.74 μg/l. The groundwater contamination at the pharmaceutical factory site, the drainage ditch and the waste deposits is similar in composition containing among others sulfonamides and chlorinated solvents (including vinyl
Modeling the reversible, diffusive sink effect in response to transient contaminant sources.
Zhao, D; Little, J C; Hodgson, A T
2002-09-01
A physically based diffusion model is used to evaluate the sink effect of diffusion-controlled indoor materials and to predict the transient contaminant concentration in indoor air in response to several time-varying contaminant sources. For simplicity, it is assumed the predominant indoor material is a homogeneous slab, initially free of contaminant, and the air within the room is well mixed. The model enables transient volatile organic compound (VOC) concentrations to be predicted based on the material/air partition coefficient (K) and the material-phase diffusion coefficient (D) of the sink. Model predictions are made for three scenarios, each mimicking a realistic situation in a building. Styrene, phenol, and naphthalene are used as representative VOCs. A styrene butadiene rubber (SBR) backed carpet, vinyl flooring (VF), and a polyurethane foam (PUF) carpet cushion are considered as typical indoor sinks. In scenarios involving a sinusoidal VOC input and a double exponential decaying input, the model predicts the sink has a modest impact for SBR/styrene, but the effect increases for VF/phenol and PUF/naphthalene. In contrast, for an episodic chemical spill, SBR is predicted to reduce the peak styrene concentration considerably. A parametric study reveals for systems involving a large equilibrium constant (K), the kinetic constant (D) will govern the shape of the resulting gasphase concentration profile. On the other hand, for systems with a relaxed mass transfer resistance, K will dominate the profile. PMID:12244748
Identifiability for the pointwise source detection in Fisher’s reaction-diffusion equation
NASA Astrophysics Data System (ADS)
Ben Belgacem, Faker
2012-06-01
We are interested in the detection of a pointwise source in a class of semi-linear advection-diffusion-reaction equations of Fisher type. The source is determined by its location, which may be steady or unsteady, and its time-dependent intensity. Observations recorded at a couple of points are the available data. One observing station is located upstream of the source and the other downstream. This is a severely ill-posed nonlinear inverse problem. In this paper, we pursue an identifiability result. The process we follow has been developed earlier for the linear model and may be sharpened to operate for the semi-linear equation. It is based on the uniqueness for a parabolic (semi-linear) sideways problem, which is obtained by a suitable unique continuation theorem. We state a maximum principle that turns out to be necessary for our proof. The identifiability is finally obtained for a stationary or a moving source. Many applications may be found in biology, chemical physiology or environmental science. The problem we deal with is the detection of pointwise organic pollution sources in rivers and channels. The basic equation to consider is the one-dimensional biochemical oxygen demand equation, with a nonlinear power growth inhibitor and/or the Michaelis-Menten reaction coefficient.
NASA Astrophysics Data System (ADS)
Tetsu, Hiroyuki; Nakamoto, Taishi
2016-03-01
Radiation is an important process of energy transport, a force, and a basis for synthetic observations, so radiation hydrodynamics (RHD) calculations have occupied an important place in astrophysics. However, although the progress in computational technology is remarkable, their high numerical cost is still a persistent problem. In this work, we compare the following schemes used to solve the nonlinear simultaneous equations of an RHD algorithm with the flux-limited diffusion approximation: the Newton-Raphson (NR) method, operator splitting, and linearization (LIN), from the perspective of the computational cost involved. For operator splitting, in addition to the traditional simple operator splitting (SOS) scheme, we examined the scheme developed by Douglas & Rachford (DROS). We solve three test problems (the thermal relaxation mode, the relaxation and the propagation of linear waves, and radiating shock) using these schemes and then compare their dependence on the time step size. As a result, we find the conditions of the time step size necessary for adopting each scheme. The LIN scheme is superior to other schemes if the ratio of radiation pressure to gas pressure is sufficiently low. On the other hand, DROS can be the most efficient scheme if the ratio is high. Although the NR scheme can be adopted independently of the regime, especially in a problem that involves optically thin regions, the convergence tends to be worse. In all cases, SOS is not practical.
NASA Astrophysics Data System (ADS)
Martin-Samos, Layla; Bussi, Giovanni
2009-08-01
We present here SaX (Self-energies and eXcitations), a plane-waves package aimed at electronic-structure and optical-properties calculations in the GW framework, namely using the GW approximation for quasi-particle properties and the Bethe-Salpeter equation for the excitonic effects. The code is mostly written in FORTRAN90 in a modern style, with extensive use of data abstraction (i.e. objects). SaX employs state of the art techniques and can treat large systems. The package is released with an open source license and can be also download from http://www.sax-project.org/. Program summaryProgram title: SaX (Self-energies and eXcitations) Catalogue identifier: AEDF_v1_0 Program summary URL:http://cpc.cs.qub.ac.uk/summaries/AEDF_v1_0.html Program obtainable from: CPC Program Library, Queen's University, Belfast, N. Ireland Licensing provisions: GNU General Public License No. of lines in distributed program, including test data, etc.: 779 771 No. of bytes in distributed program, including test data, etc.: 4 894 755 Distribution format: tar.gz Programming language: FORTRAN, plus some C utilities Computer: Linux PC, Linux clusters, IBM-SP5 Operating system: Linux, Aix Has the code been vectorised or parallelized?: Yes RAM: depending on the system complexity Classification: 7.3 External routines: Message-Passing Interface (MPI) to perform parallel computations. ESPRESSO ( http://www.quantum-espresso.org) Nature of problem: SaX is designed to calculate the electronic band-structure of semiconductors, including quasi-particle effects and optical properties including excitonic effects. Solution method: The electronic band-structure is calculated using the GW approximation for the self-energy operator. The optical properties are calculated solving the Bethe-Salpeter equation in the GW approximation. The wavefunctions are expanded on a plane-waves basis set, using norm-conserving pseudopotentials. Restrictions: Many objects are non-local matrix represented in plane wave basis
The ionization sources of the diffuse ionized gas in nearby disk galaxies
NASA Astrophysics Data System (ADS)
Voges, Erica Susan
Diffuse ionized gas (DIG) has been shown to be an important component of the interstellar medium (ISM), with its large filling factor (>= 20%) and a mass that makes it the most massive component of the Galactic ionized ISM. Given that it has been found to be ubiquitous in both the Galaxy and external disk galaxies, the energy source to create and maintain the DIG must necessarily be large. Massive OB stars are the only known sources with enough energy to power the DIG, and DIG is also linked morphologically to OB stars as it is brightest near bright star forming regions. However, the details of the location and spectral types of the ionizing stars, as well as the relevance of other ionizing mechanisms, are still not clear. I present the results of three different studies aimed at exploring the ionization sources of the DIG. Optical spectroscopy of DIG in M33 and NGC 891 using the Gemini-North telescope has been obtained to compare diagnostic emission line ratios with photoionization models. The first detection of (O I] l6300 was made in the DIG of M33. In M33, models in which ionizing photons leaking from H II regions are responsible for the ionization of the DIG best fit our observed line ratios. In NGC 891, we found evidence that shock ionization may need to be included along with photoionization in order to explain our observed emission line ratios. The diffuse Ha fraction in eight nearby galaxies was studied as a function of radius and star formation rate per unit area. We found no correlation with radius, but we did find that regions with higher star formation rates have lower diffuse fractions. Neither of these results had any dependence on galaxy type. These results have implications regarding the circumstances under which H II regions may be leaking ionizing photons and thus ionizing DIG. We also compared observed and predicted ionizing photon emission rates for 39 H II regions in the Large Magellanic Cloud. Our results indicate that five of the H II
Ding, Jingtao; Xi, Beidou; Gao, Rutai; He, Liansheng; Liu, Hongliang; Dai, Xuanli; Yu, Yijun
2014-06-15
Nitrate (NO3(-)) pollution is a severe problem in aquatic systems in Taihu Lake Basin in China. A dual isotope approach (δ(15)NNO3(-) and δ(18)ONO3(-)) was applied to identify diffused NO3(-) inputs in a stream in an agricultural field at the basin in 2013. The site-specific isotopic characteristics of five NO3(-) sources (atmospheric deposition, AD; NO3(-) derived from soil organic matter nitrification, NS; NO3(-) derived from chemical fertilizer nitrification, NF; groundwater, GW; and manure and sewage, M&S) were identified. NO3(-) concentrations in the stream during the rainy season [mean±standard deviation (SD)=2.5±0.4mg/L] were lower than those during the dry season (mean±SD=4.0±0.5mg/L), whereas the δ(18)ONO3(-) values during the rainy season (mean±SD=+12.3±3.6‰) were higher than those during the dry season (mean±SD=+0.9±1.9‰). Both chemical and isotopic characteristics indicated that mixing with atmospheric NO3(-) resulted in the high δ(18)O values during the rainy season, whereas NS and M&S were the dominant NO3(-) sources during the dry season. A Bayesian model was used to determine the contribution of each NO3(-) source to total stream NO3(-). Results showed that reduced N nitrification in soil zones (including soil organic matter and fertilizer) was the main NO3(-) source throughout the year. M&S contributed more NO3(-) during the dry season (22.4%) than during the rainy season (17.8%). AD generated substantial amounts of NO3(-) in May (18.4%), June (29.8%), and July (24.5%). With the assessment of temporal variation of diffused NO3(-) sources in agricultural field, improved agricultural management practices can be implemented to protect the water resource and avoid further water quality deterioration in Taihu Lake Basin. PMID:24686140
NASA Astrophysics Data System (ADS)
Jaynes, A. N.; Lessard, M.; Rodriguez, J. V.; Rychert, K. M.; Donovan, E.; Michell, R. G.; Samara, M.
2012-12-01
Pulsating aurora (PA) is a common ionospheric phenomenon and as such offers a unique opportunity to study the source of the precipitating particle populations. While the suggestion that the source of these electrons originates near the equator was made decades ago, this conjecture has not been confirmed until now. The dominant source of loss-cone scattering for energetic equatorial electrons, which can then precipitate as PA, has been explored, but not yet clearly identified. In this study, we compare the frequencies of equatorial electron flux pulsations out in space and pulsating aurora luminosity fluctuations in the corresponding magnetic footprint within the ionosphere. We use simultaneous satellite- and ground-based data from GOES 13 and THEMIS instrumentation to show that there is a direct correlation between diffuse luminosity fluctuations near the ground and particle pulsations in equatorial space. To study the local morphology of such instances, observations of a dynamic pulsating aurora event were taken with a pair of colocated allsky imagers at Poker Flat, one filtered at 4278 (blue) and one at 5577 (green). Here we show preliminary results of differences in the energy channels and the structure that emerges and disappears as pulsating starts and stops.;
Huang, Chong; Lin, Yu; He, Lian; Irwin, Daniel; Szabunio, Margaret M.; Yu, Guoqiang
2016-01-01
Noncontact diffuse correlation tomography (ncDCT) is an emerging technology for 3D imaging of deep tissue blood flow distribution without distorting hemodynamic properties. To adapt the ncDCT for imaging in vivo breast tumors, we designed a motorized ncDCT probe to scan over the breast surface. A computer-aided design (CAD)-based approach was proposed to create solid volume mesh from arbitrary breast surface obtained by a commercial 3D camera. The sources and detectors of ncDCT were aligned on the breast surface through ray tracing to mimic the ncDCT scanning with CAD software. The generated breast volume mesh along with the boundary data of ncDCT at the aligned source and detector pairs were used for finite-element-method-based flow image reconstruction. We evaluated the accuracy of source alignments on mannequin and human breasts; largest alignment errors were less than 10% in both tangential and radial directions of scanning. The impact of alignment errors (assigned 10%) on the tumor reconstruction was estimated using computer simulations. The deviations of simulated tumor location and blood flow contrast resulted from the alignment errors were 0.77 mm (less than the node distance of 1 mm) and 1%, respectively, which result in minor impact on flow image reconstruction. Finally, a case study on a human breast tumor was conducted and a tumor-to-normal flow contrast was reconstructed, demonstrating the feasibility of ncDCT in clinical application. PMID:26479823
Qin, G.; He, H.-Q.; Zhang, M. E-mail: hqhe@spaceweather.ac.cn
2011-09-01
Recently, Tan and coworkers studied the 2001 September 24 solar energetic particle (SEP) event observed by the Wind spacecraft at 1 AU and found that there is a counter-streaming particle beam with a deep depression of flux at 90{sup 0} pitch angle during the beginning of the event. They suggested that it is a result of a reflecting boundary at some distance outside of 1 AU. While this scenario could be true under some specific configuration of an interplanetary magnetic field, in this paper we offer another possible explanation. We simulated the SEP event by solving the five-dimensional focused transport equation numerically for 40 keV electrons with perpendicular diffusion. We find that a counter-streaming particle beam with deep depression at 90{sup 0} pitch angle can form on Parker magnetic field lines that do not directly connect to the main particle source on the Sun in the beginning of an SEP event. It can happen when a significant number of observed particles come from adjacent field lines through parallel transport to large radial distance first, hopping across field lines through perpendicular diffusion, and then getting scattered back to 1 AU, where they combine with the particles directly coming from the Sun to form a counter-streaming beam.
Non-contact time-resolved diffuse reflectance imaging at null source-detector separation.
Mazurenka, M; Jelzow, A; Wabnitz, H; Contini, D; Spinelli, L; Pifferi, A; Cubeddu, R; Mora, A Dalla; Tosi, A; Zappa, F; Macdonald, R
2012-01-01
We report results of the proof-of-principle tests of a novel non-contact tissue imaging system. The system utilizes a quasi-null source-detector separation approach for time-domain near-infrared spectroscopy, taking advantage of an innovative state-of-the-art fast-gated single photon counting detector. Measurements on phantoms demonstrate the feasibility of the non-contact approach for the detection of optically absorbing perturbations buried up to a few centimeters beneath the surface of a tissue-like turbid medium. The measured depth sensitivity and spatial resolution of the new system are close to the values predicted by Monte Carlo simulations for the inhomogeneous medium and an ideal fast-gated detector, thus proving the feasibility of the non-contact approach for high density diffuse reflectance measurements on tissue. Potential applications of the system are also discussed. PMID:22274351
Electrosphere of macroscopic 'quark nuclei': A source for diffuse MeV emissions from dark matter
Forbes, Michael McNeil; Lawson, Kyle; Zhitnitsky, Ariel R.
2010-10-15
Using a Thomas-Fermi model, we calculate the structure of the electrosphere of the quark antimatter nuggets postulated to comprise much of the dark matter. This provides a single self-consistent density profile from ultrarelativistic densities to the nonrelativistic Boltzmann regime that we use to present microscopically justified calculations of several properties of the nuggets, including their net charge, and the ratio of MeV to 511 keV emissions from electron annihilation. We find that the calculated parameters agree with previous phenomenological estimates based on the observational supposition that the nuggets are a source of several unexplained diffuse emissions from the Galaxy. As no phenomenological parameters are required to describe these observations, the calculation provides another nontrivial verification of the dark-matter proposal. The structure of the electrosphere is quite general and will also be valid at the surface of strange-quark stars, should they exist.
Introducing Carbon Diffusion Barriers for Uniform, High-Quality Graphene Growth from Solid Sources
2013-01-01
Carbon diffusion barriers are introduced as a general and simple method to prevent premature carbon dissolution and thereby to significantly improve graphene formation from the catalytic transformation of solid carbon sources. A thin Al2O3 barrier inserted into an amorphous-C/Ni bilayer stack is demonstrated to enable growth of uniform monolayer graphene at 600 °C with domain sizes exceeding 50 μm, and an average Raman D/G ratio of <0.07. A detailed growth rationale is established via in situ measurements, relevant to solid-state growth of a wide range of layered materials, as well as layer-by-layer control in these systems. PMID:24024736
Dust and gas jets: Evidence for a diffuse source in Halley's coma
NASA Technical Reports Server (NTRS)
Clairemidi, Jacques; Rousselot, Philippe; Vernotte, F.; Moreels, Guy
1992-01-01
The distribution of dust-scattered intensity in Halley's inner coma is measured with the Vega three-channel spectrometer at three selected wavelengths: 377, 482, and 607 nm. The variation along a cometo-centric radius may be described by a p(sup -s) law where p is the distance between nucleus and optical axis and s is an exponent which is equal to 1 except in an intermediate 3000 less than p less than 7000 km region where s = 1.5. The shape of the radial distribution may be explained with a model including solar radiation pressure effect and quantum scattering efficiencies calculated from Mie theory. Monochromatic images inside an angular sector having its apex at the nucleus show evidence of two dust jets which extend to 40,000 Km. The pixel-to-pixel ratio of two images of dust intensity at 377 and 482 nm shows that the scattered intensity presents an excess of blue coloration in a zone located around the jets between 10,000 and 25,000 km. This coloration is interpreted as being due to a population of sub-micronic grains which result of the fragmentation of dust particles transported in the jets. It is suggested that the diffuse source where an additional quantity of CO was detected might be connected with the presence of a dust jet. In the present scheme, grain particles with a size of several micron or 10 micron would be transported inside a dust jet to distances of several 10,000 km where they would suffer fragmentation and produce sub-micronic particles and a release of gas which would be at the origin of the diffuse source.
Ceschin, Rafael; Panigrahy, Ashok; Gopalakrishnan, Vanathi
2015-01-01
A major challenge in the diagnosis and treatment of brain tumors is tissue heterogeneity leading to mixed treatment response. Additionally, they are often difficult or at very high risk for biopsy, further hindering the clinical management process. To overcome this, novel advanced imaging methods are increasingly being adapted clinically to identify useful noninvasive biomarkers capable of disease stage characterization and treatment response prediction. One promising technique is called functional diffusion mapping (fDM), which uses diffusion-weighted imaging (DWI) to generate parametric maps between two imaging time points in order to identify significant voxel-wise changes in water diffusion within the tumor tissue. Here we introduce serial functional diffusion mapping (sfDM), an extension of existing fDM methods, to analyze the entire tumor diffusion profile along the temporal course of the disease. sfDM provides the tools necessary to analyze a tumor data set in the context of spatiotemporal parametric mapping: the image registration pipeline, biomarker extraction, and visualization tools. We present the general workflow of the pipeline, along with a typical use case for the software. sfDM is written in Python and is freely available as an open-source package under the Berkley Software Distribution (BSD) license to promote transparency and reproducibility. PMID:25673970
Curotto, E; Mella, Massimo
2015-03-21
We test the second order Milstein method adapted to simulate diffusion in general compact Riemann manifolds on a number of systems characterized by nonconfining potential energy surfaces of increasing complexity. For the 2-sphere and more complex spaces derived from it, we compare the Milstein method with a number of other first and second order approaches. In each case tested, we find evidence that demonstrate the versatility and relative ease of implementation of the Milstein method derived in Part I. PMID:25796235
NASA Astrophysics Data System (ADS)
Legout, Cédric; Poulenard, Jérôme; Nemery, Julien; Navratil, Oldrich; Grangeon, Thomas; Evrard, Olivier; Esteves, Michel
2013-04-01
Increasing the understanding of the hydro-sedimentary dynamics at the catchment scale requires data on the origin of suspended sediments in addition to the classical measurements of suspended sediment concentrations and discharge. In mountainous environments the extremely high spatial and temporal variability of suspended sediment fluxes suggests that the proportions of suspended sediment sources transiting at outlets may also exhibit strong variations during flood events. However, conventional fingerprinting techniques based on geochemical and radionuclide measurements are time-consuming and rather expensive. They constitute a major limitation to conduct routine characterisation of the source of suspended sediment transiting at outlets that would require the analysis of a large number of samples. We investigated how visible or infrared diffuse reflectance spectroscopy coupled with partial least squares chemometrics techniques could be used to predict the proportion of source material in suspended sediment samples and how it could help understanding the hydro-sedimentary processes occurring during floods. The studied catchment is located in the southern French Alps, draining an area of 22-km². It is composed of black marls, limestones, molasses, undifferentiated deposits and gypsum. Forty-eight source material samples were collected in badlands areas and 328 suspended sediment samples were collected at the outlet during 23 flood events. Spectroscopic measurements were carried out on dried samples. Given that the erosion processes are particle size selective, five size fractions of source material were measured in order to assess the potential alteration of the signatures. As the biogeochemical processes occurring in the river such as iron oxidation could also affect the signatures, source materials that were immersed in the river for durations ranging from 1 day to 9 weeks were analysed. Finally, partial least-squares regression models were constructed on 81
Glass diffusion source for constraining BSF region of a solar cell
Lesk, I.A.; Pryor, R.A.; Coleman, M.G.
1982-08-27
The present invention is directed to a method of fabricating a solar cell comprising simultaneous diffusion of the p and n dopant materials into the solar cell substrate. The simultaneous diffusion process is preceded by deposition of a capping layer impervious to doping by thermal diffusion processes.
Pérez-Calatayud, J; Granero, D; Ballester, F; Puchades, V; Casal, E
2004-03-01
The purpose of this study is to calculate the dose rate distribution for the Amersham Cs-137 pellet source used in brachytherapy with the Selectron low-dose-rate remote afterloading system in gynaecological applications using the Monte Carlo code GEANT4. The absolute dose rate distribution for the pellet source was obtained and presented as a one-dimensional absolute dose rate table as well as in the Task Group 43 dose-calculation formalism. In this study, excellent agreement was found between the point source theoretical model using fitted polynomial values and Monte Carlo calculations of the dose rate distribution for the pellet source. A comparison study was also made between the dose rate distribution obtained from a complete Monte Carlo simulation (Cs-137 pellet sources + remote afterloading system plastic guide tube + gynaecological applicator) and that calculated by applying the superposition principle to Monte Carlo data of the individual pellet sources. The data were obtained for a portion of uterine tandem of typical train source configurations. Significant differences with a strong dependence on polar angle have been found that must be kept in mind for clinical dosimetry. PMID:15070245
NASA Technical Reports Server (NTRS)
Kylling, Arve; Stamnes, Knut
1992-01-01
The present solutions to the linear transport equation pertain to monoenergetic particles' interaction with a multiple scattering/absorbing layered medium with a general anisotropic internal source term. Attention is given to a novel exponential-linear approximation to the internal source, as a function of scattering depth, which furnishes an at-once efficient and accurate solution to the linear transport equation through its reduction of the spatial mesh size. The great superiority of the proposed method is demonstrated by the numerical results obtained in the illustrative cases of (1) an embedded thermal source and (2) a rapidly varying beam pseudosource.
Assessment tools are being developed to predict diffuse NPS effects from watershed development and distinguish these from local impacts (point sources, contaminated sediments). Using EMAP data from the New England Wadeable Stream Survey and two state datasets (CT, ME), we are de...
Technology Transfer Automated Retrieval System (TEKTRAN)
Thermal conductivity and thermal diffusivity are two important physical properties essential for designing any food engineering processes. Recently a new transient plane-source method was developed to measure a variety of materials, but its application in foods has not been documented. Therefore, ...
Hα line measurements from ten diffuse galactic sources using the DEFPOS facility
NASA Astrophysics Data System (ADS)
Sahan, M.; Oflaz, F. M.; Yegingil, I.; Tel, E.
2015-08-01
The hydrogen Balmer-α emission line spectrum of ten diffuse ionization sources in the Milk Way - NGC 40 (WC8), NGC 2022, NGC 6210, NGC 6618 (M17, Sh2-45), NGC 6720 (M57), NGC 6781, NGC 6888 (Sh2-105), NGC 6992 (Sh2-103), NGC 7635 (Sh2-162,) and IC 1848 (Sh2-199) - has been investigated using a dual etalon Fabry-Pérot optical spectrometer (DEFPOS) aatached to the 150 cm RTT150 telescope at TUBITAK National Observatory (TUG, Antalya, Turkey: 36° 51' N; 30° 20' E; elevation: 2547 m). All of our galactic Hα observations discussed in this paper were carried out during the nights of 2013 June 21-24 with exposure time of 3600 s. As main results the intensity, the full width at half maximum, and the radial velocity with respect to the LSR have been determined for each data set. The intensities, the radial velocities, and the line widths of the Hα emission line vary from 59.15 to 8923.44 R, -46.72 to +54.07 km s-1, and 31.4 to 48.01 km s-1, respectively. The radial velocities and the half-widths of the H II regions and planetary nebulae determined from our measurements are found to be consistent with values given in literature, especially with those in Schneider et al. (1983) and Fich et al. (1990).
Water quality along a river continuum subject to point and diffuse sources
NASA Astrophysics Data System (ADS)
Neal, Colin; Jarvie, Helen P.; Love, Alison; Neal, Margaret; Wickham, Heather; Harman, Sarah
2008-02-01
SummaryThe water quality along the River Kennet, in the Thames basin of southern England, was examined in terms of the influence of point- and diffuse-nutrient inputs. The river is supplied mainly from a Cretaceous Chalk aquifer and hence the waters are of a calcium bicarbonate type. The nitrate largely comes from agricultural sources, with concentrations decreasing downstream due to plant uptake and probable denitrification. In contrast, soluble reactive phosphorus (SRP) is largely associated with sewage inputs and concentrations increase downstream in line with effluents from major towns such as Newbury and Reading. Adjacent to the river in the lower half of the catchment is the Kennet and Avon Canal and the two are in places hydrologically connected. The canal inputs may influence calcium carbonate (calcite) precipitation and increase suspended sediment and particulate phosphorus concentrations in the river. Monitoring upstream and downstream of Marlborough sewage treatment works (STW) showed that SRP concentrations in the effluent were highly variable due to variable efficiency of P stripping and still sufficiently concentrated to dominate downstream river SRP with potential impacts on stream ecology. Biological recovery in this river following P stripping at STWs is complex and controlling those spikes in SRP that are above a threshold of 100 μg l -1 may be a critical requirement. More stringent effluent targets than are currently recommended may be needed (less than 800 μg RP l -1) to achieve good ecological status in this river depending on SRP concentrations upstream.
NASA Astrophysics Data System (ADS)
Katzarov, Ivaylo H.; Pashov, Dimitar L.; Paxton, Anthony T.
2013-08-01
We present calculations of free energy barriers and diffusivities as functions of temperature for the diffusion of hydrogen in α-Fe. This is a fully quantum mechanical approach since the total energy landscape is computed using a self-consistent, transferable tight binding model for interstitial impurities in magnetic iron. Also the hydrogen nucleus is treated quantum mechanically and we compare here two approaches in the literature both based in the Feynman path integral formulation of statistical mechanics. We find that the quantum transition state theory which admits greater freedom for the proton to explore phase space gives result in better agreement with experiment than the alternative which is based on fixed centroid calculations of the free energy barrier. This will have an impact on future modeling and the simulation of hydrogen trapping and diffusion.
Chen, Hao; Zheng, Bing-Hui; Zhang, Lei
2013-02-01
Dissolved organic matter collected in Daning River (China) in July 2009 was investigated with parallel factor analysis (PARAFAC) and fluorescence spectroscopy with the aim of identifying the origin of dissolved humic substance (HS) components. Two HS-like fluorescence components (peak M and C) with excitation/emission (ex/em) maxima at 305/406 nm and 360/464 nm showed relatively uniform distribution in the vertical direction for each sampling site but a trend of accumulation down the river, independent of the highly heterogeneous water environment as implicated by water quality parameters (i.e., water temperature, algae density, chlorophyll a, dissolved oxygen, dissolved organic carbon, pH, conductivity and turbidity), while an amino acid/protein-like component (peak T; ex/em = 280/334 nm) was quite variable in its spatial distribution, implying strong influence from point sources (e.g. sewage discharge) and local microbial activities. The fluorescence intensity (F max in Raman units) at these ex/em wavelength pairs fell in the range of 0.031-0.358, 0.051-0.224 and 0.026-0.115 for peak T, M and C, respectively. In addition, the F max values of peak C covaried with M (i.e. C = 0.503 ×M, p < 0.01, R (2) = 0.973). Taken together, these results indicate that peak M and C originated primarily and directly from the same soil sources that were diffusive in the catchment, but peak T was more influenced by local point sources (e.g. wastewater discharge) and in situ microbial activities. This study presents new insights into the currently controversial origin of some HS components (e.g."peak M", as commonly referred to in the literature). This study highlights that natural water samples should be collected at various depths in addition to along a river/stream flow path so as to better evaluate the origin of HS fluorescence components. PMID:25208714
Soils as sinks or sources for diffuse pollution of the water cycle
NASA Astrophysics Data System (ADS)
Grathwohl, Peter
2010-05-01
Numerous chemical compounds have been released into the environment by human activities and can nowadays be found everywhere, i.e. in the compartments water, soil, and air, at the poles and in high mountains. Examples for a global distribution of toxic compounds are the persistent organic pollutants (PCB, dioxins, PAH, fluorinated surfactants and flame retardants, etc.: "the Stockholm dirty dozen") but also mercury and other metals. Many of these compounds reached a global distribution via the atmo¬sphere; others have been and are still directly applied to top soils at the large scale by agriculture or are released into groundwater at landfill sites or by discharge of treated or untreated waste waters. Sooner or later such compounds end up in the water cycle - often via an intermediate storage in soils. Pollutants in soils are leached by seepage waters, transferred to ground¬water, and transported to rivers via groundwater flow. Adsorbed compounds may be transported from soils into surface waters by erosion processes and will end up in the sediments. Diffuse pollution of the subsurface environment not only reflects the history of the economic development of the modern society but it is still ongoing - e.g. the number of organic pollutants released into the environment is increasing even though the con¬centrations may decrease compared to the past. Evidence shows that many compounds are persistent in the subsurface environment at large time scales (up to centuries). Thus polluted soils already are or may become a future source for pollution of adjacent compartments such as the atmosphere and groundwater. A profound understanding on how diffuse pollutants are stored and processed in the subsurface environment is crucial to assess their long term fate and transport at large scales. Thus integrated studies e.g. at the catchment scale and models are needed which couple not only the relevant compartments (soil - atmosphere - groundwater/surface waters) but also flow
FRETBursts: An Open Source Toolkit for Analysis of Freely-Diffusing Single-Molecule FRET
Lerner, Eitan; Chung, SangYoon; Weiss, Shimon; Michalet, Xavier
2016-01-01
Single-molecule Förster Resonance Energy Transfer (smFRET) allows probing intermolecular interactions and conformational changes in biomacromolecules, and represents an invaluable tool for studying cellular processes at the molecular scale. smFRET experiments can detect the distance between two fluorescent labels (donor and acceptor) in the 3-10 nm range. In the commonly employed confocal geometry, molecules are free to diffuse in solution. When a molecule traverses the excitation volume, it emits a burst of photons, which can be detected by single-photon avalanche diode (SPAD) detectors. The intensities of donor and acceptor fluorescence can then be related to the distance between the two fluorophores. While recent years have seen a growing number of contributions proposing improvements or new techniques in smFRET data analysis, rarely have those publications been accompanied by software implementation. In particular, despite the widespread application of smFRET, no complete software package for smFRET burst analysis is freely available to date. In this paper, we introduce FRETBursts, an open source software for analysis of freely-diffusing smFRET data. FRETBursts allows executing all the fundamental steps of smFRET bursts analysis using state-of-the-art as well as novel techniques, while providing an open, robust and well-documented implementation. Therefore, FRETBursts represents an ideal platform for comparison and development of new methods in burst analysis. We employ modern software engineering principles in order to minimize bugs and facilitate long-term maintainability. Furthermore, we place a strong focus on reproducibility by relying on Jupyter notebooks for FRETBursts execution. Notebooks are executable documents capturing all the steps of the analysis (including data files, input parameters, and results) and can be easily shared to replicate complete smFRET analyzes. Notebooks allow beginners to execute complex workflows and advanced users to
FRETBursts: An Open Source Toolkit for Analysis of Freely-Diffusing Single-Molecule FRET.
Ingargiola, Antonino; Lerner, Eitan; Chung, SangYoon; Weiss, Shimon; Michalet, Xavier
2016-01-01
Single-molecule Förster Resonance Energy Transfer (smFRET) allows probing intermolecular interactions and conformational changes in biomacromolecules, and represents an invaluable tool for studying cellular processes at the molecular scale. smFRET experiments can detect the distance between two fluorescent labels (donor and acceptor) in the 3-10 nm range. In the commonly employed confocal geometry, molecules are free to diffuse in solution. When a molecule traverses the excitation volume, it emits a burst of photons, which can be detected by single-photon avalanche diode (SPAD) detectors. The intensities of donor and acceptor fluorescence can then be related to the distance between the two fluorophores. While recent years have seen a growing number of contributions proposing improvements or new techniques in smFRET data analysis, rarely have those publications been accompanied by software implementation. In particular, despite the widespread application of smFRET, no complete software package for smFRET burst analysis is freely available to date. In this paper, we introduce FRETBursts, an open source software for analysis of freely-diffusing smFRET data. FRETBursts allows executing all the fundamental steps of smFRET bursts analysis using state-of-the-art as well as novel techniques, while providing an open, robust and well-documented implementation. Therefore, FRETBursts represents an ideal platform for comparison and development of new methods in burst analysis. We employ modern software engineering principles in order to minimize bugs and facilitate long-term maintainability. Furthermore, we place a strong focus on reproducibility by relying on Jupyter notebooks for FRETBursts execution. Notebooks are executable documents capturing all the steps of the analysis (including data files, input parameters, and results) and can be easily shared to replicate complete smFRET analyzes. Notebooks allow beginners to execute complex workflows and advanced users to
NASA Astrophysics Data System (ADS)
Ge, J.; Everett, M. E.; Weiss, C. J.
2010-12-01
An interpretation based on the Continuous Time Random Walk theory (CTRW) to the diffusion of electromagnetic fields generated by a transient straight current source over a porous geological media is presented here. The CTRW theory is demonstrated to be a powerful tool to concisely and more accurately model a transport process in a fractal medium with complex structures, comparing to the classical transport theory. In the controlled-source electromagnetic (EM) induction setting, the time dependent evolution of the EM field of some sources over a rough medium are governed by the fractional diffusion EM equation in a CTRW sense. The master equation can be solved for a uniform conducting half-space in the Laplace domain semi-analytically. We use 2D finite difference method to calculate the solution numerically for the assigned space and transform to time domain with Gaver-Stehfest algorithm. Here we adopt a spatially uniform roughness parameter β in the solution to characterize the complexity of the geoelectrical structure of the geological medium. To introduce the heterogeneity to our model, we set up the space as several 2 D blocks with different conductivities and βs. Then we compare our results with the synthetic data we got from the high resolution numerical simulations. We are able to show that by introducing the heterogeneity to the fractional diffusion perspective, our approach is competent for tracing the diffusion process with less model parameters.
Kolpin, Dana W.; Schenzel, Judith; Meyer, Michael T.; Phillips, Patrick J.; Hubbard, Laura E.; Scott, Tia-Marie; Bucheli, Thomas D.
2014-01-01
To determine the prevalence of mycotoxins in streams, 116 water samples from 32 streams and three wastewater treatment plant effluents were collected in 2010 providing the broadest investigation on the spatial and temporal occurrence of mycotoxins in streams conducted in the United States to date. Out of the 33 target mycotoxins measured, nine were detected at least once during this study. The detections of mycotoxins were nearly ubiquitous during this study even though the basin size spanned four orders of magnitude. At least one mycotoxin was detected in 94% of the 116 samples collected. Deoxynivalenol was the most frequently detected mycotoxin (77%), followed by nivalenol (59%), beauvericin (43%), zearalenone (26%), β-zearalenol (20%), 3-acetyl-deoxynivalenol (16%), α-zearalenol (10%), diacetoxyscirpenol (5%), and verrucarin A (1%). In addition, one or more of the three known estrogenic compounds (i.e. zearalenone, α-zearalenol, and β-zearalenol) were detected in 43% of the samples, with maximum concentrations substantially higher than observed in previous research. While concentrations were generally low (i.e. < 50 ng/L) during this study, concentrations exceeding 1000 ng/L were measured during spring snowmelt conditions in agricultural settings and in wastewater treatment plant effluent. Results of this study suggest that both diffuse (e.g. release from infected plants and manure applications from exposed livestock) and point (e.g. wastewater treatment plants and food processing plants) sources are important environmental pathways for mycotoxin transport to streams. The ecotoxicological impacts from the long-term, low-level exposures to mycotoxins alone or in combination with complex chemical mixtures are unknown
Tang Jie; Jiang Weiman; Zhao Wei; Wang Yishan; Li Shibo; Wang Haojing; Duan Yixiang
2013-01-21
A stable diffuse large-volume air plasma source was developed by using argon-induced dielectric-barrier discharges at atmospheric pressure. This plasma source can be operated in a filamentary discharge with the average areal power density of 0.27 W/cm{sup 2} and the gas temperature of 315{+-}3 K. Spatial measurement of emission spectrum and temperature indicates that this plasma is uniform in the central region along the transverse direction. It is also found that the formation of diffuse air plasma mainly lies in the creation of sufficient seed electrons by the Penning effect through collisions between two argon or nitrogen metastables at low electric fields.
NASA Technical Reports Server (NTRS)
Maoz, Eyal; Grindlay, Jonathan E.
1995-01-01
The incompatibility of the properties of the X-ray background (XRB) with active galactic nuclei (AGNs) contributing approximately greater than 60% at energies of a few keV has often been interpreted as being due to a substantial contribution of a new population of yet unrecognized X-ray sources. The existence of such population has been recently suggested also by an analysis of very deep ROSAT observations which revealed a considerable excess of faint X-ray sources over that expected from QSO evolution models, and that the average spectrum of the resolved sources becomes harder with decreasing flux limit. These sources could be extragalactic in origin, but if they make a substantial contribution to the XRB then they must exhibit much weaker clustering than galaxies or QSOs in order to be consistent with the stringent constraints on source clustering imposed by autocorrelation analyses of the unresolved XRB. We investigate the possibility that the indicated new population of X-ray sources is Galactic in origin. Examining spherical halo and thick disk distributions, we derive the allowed properties of such populations which would resolve the discrepancy found in the number counts of faint sources and be consistent with observational constraints on the total background intensity, the XRB anisotropy, the number of unidentified bright sources, the Galaxy's total X-ray luminosity, and with the results of fluctuation analyses of the unresolved XRB. We find that a flattened Galactic halo (or a thick disk) distribution with a scale height of a few kpc is consistent with all the above requirements. The typical X-ray luminosity of the sources is approximately equal to 10(exp 30-31)ergs/s in the 0.5-2 keV band, the number density of sources in the solar vicinity is approximately 10(exp -4.5)pc(exp -3), their total number in the Galaxy is approximately 10(exp 8.5), and their total contribution to the Galaxy's X-ray luminosity is approximately 10(exp 39) ergs/s. We discuss the
NASA Astrophysics Data System (ADS)
Sujatha, N.; Nivetha, K. Bala; Singhal, Akshay
2014-05-01
Localization and determination of blood region parameters in skin tissue can serve as a valuable supplement to standard non invasive techniques, especially in accessing the degree of depth of burns on skin and for the classification of vascular malformations. Quantitative optical examination of skin local blood region requires the use of depth sensitive techniques and preferential probing for assessment of data from specific layers of skin tissue. This work incorporates the depth sensitivity of diffuse reflectance spectroscopy and optimal source to detector fiber separation for maximum reflectance collection efficiency from local blood region in skin. Monte Carlo simulation for diffuse reflectance was performed on a multi layered skin tissue model consisting of epidermis, perfused dermis and localized blood region. It was found that the slope of the spatially resolved reflectance curve plotted with respect to the source to detector separation distance in semi log scale varies with the depth of the local blood region at specific wavelengths corresponding to the absorption wavelengths of hemoglobin. From the depth information obtained from the spatially resolved reflectance data, the optimum source to detector separation (SDS) is determined for maximum collection efficiency from the chromophore layer. The results obtained from simulation suggest the design of a linearly variable source to detector separation probe for preferential analysis of the depth of a specific tissue layer and subsequent determination of optimal source to detector separation for extracting the layer information.
NASA Astrophysics Data System (ADS)
Park, Kyoungsu; Lee, Minseok; Lee, Seung-ha; Cerussi, Albert E.; Chung, Phil-sang; Kim, Sehwan
2015-03-01
Many biomedical applications require an efficient combination and localization of multiple discrete light sources. In this paper, we present a compact six-channel combiner of optical sub-assembly type that couples the output of independent solid-state light sources into a single 400 μm diameter optical fiber. It is equipped with six discrete laser diodes, 658, 690, 705, 785, 830, and 850 nm for the measurement of the tissue optical properties from optical spectroscopy and imaging. We demonstrate coupling efficiencies ≥ 77% and output optical power ≥ 20 mW for each of the 6 laser diodes installed into the prototype. The design supports the use of continuous wave and intensity modulated laser diodes (with bandwidth ≥ 3 GHz). The developed light source could be used to construct custom multi-wavelength sources for tissue oximeters, diffuse optical imaging, and molecular imaging technologies.
Zushi, Yasuyuki; Masunaga, Shigeki
2011-11-01
To efficiently reduce perfluorinated compound (PFC) pollution, it is important to have an understanding of PFC sources and their contribution to the pollution. In this study, source identification of diffuse water pollution by PFCs was conducted using a GIS-based approach. Major components of the source identification were collection of the monitoring data and preparation of the corresponding geographic information that was extracted from a constructed GIS database. The spatially distributed pollution factors were then explored by multiple linear regression analysis, after which they were visually expressed using GIS. Among the 35 PFC homologues measured in a survey of the Tokyo Bay basin, 18 homologues were analyzed. Pollution by perfluorooctane sulfonate (PFOS) was explained well by the percentage of arterial traffic area in the basin, and the 84% variance of the measured PFOS concentration was explained by two geographic variables, arterial traffic area and population. Source apportionment between point and nonpoint sources was conducted based on the results of the analysis. The contribution of PFOS from nonpoint sources was comparable to that from point sources in several major rivers flowing into Tokyo Bay. Source identification and apportionment using the GIS-based approach was shown to be effective, especially for ubiquitous types of pollution, such as PFC pollution. PMID:21885084
A diffusion source for sodium and potassium in the atmospheres of Mercury and the moon
NASA Astrophysics Data System (ADS)
Sprague, A. L.
1990-03-01
Deep grain-boundary diffusion and regolith diffusion through a fractured crust and regolith can account not only for the Na/K ratios observed in the Mercurian and lunar atmospheres, but the large Na abundance enhancement of Mercury over lunar levels. A hot component of Na and K at Mercury is noted to be smaller in proportion to the total abundances of these two constituents than at the moon; this hot component is consistent with a population of meteoritic substances similar to lunar ones, as well as with a surface composition which has undergone no greater K depletion than that of the moon.
A diffusion source for sodium and potassium in the atmospheres of Mercury and the moon
NASA Technical Reports Server (NTRS)
Sprague, Ann L.
1990-01-01
Deep grain-boundary diffusion and regolith diffusion through a fractured crust and regolith can account not only for the Na/K ratios observed in the Mercurian and lunar atmospheres, but the large Na abundance enhancement of Mercury over lunar levels. A hot component of Na and K at Mercury is noted to be smaller in proportion to the total abundances of these two constituents than at the moon; this hot component is consistent with a population of meteoritic substances similar to lunar ones, as well as with a surface composition which has undergone no greater K depletion than that of the moon.
NASA Astrophysics Data System (ADS)
Sivakoff, Gregory R.; Sarazin, Craig L.; Carlin, Jeffrey L.
2004-12-01
We observed the X-ray-bright E3 galaxy NGC 1600 and nearby members of the NGC 1600 group with the Chandra X-Ray Observatory ACIS-S3 to study their X-ray properties. Unresolved emission dominates the observation; however, we resolved some of the emission into 71 sources, most of which are low-mass X-ray binaries associated with NGC 1600. Twenty-one of the sources have LX>2×1039 ergs s-1 (0.3-10.0 keV; assuming they are at the distance of NGC 1600), marking them as ultraluminous X-ray point source (ULX) candidates; we expect that only 11+/-2 are unrelated foreground/background sources. NGC 1600 may have the largest number of ULX candidates in an early-type galaxy to date; however, cosmic variance in the number of background active galactic nuclei cannot be ruled out. The spectrum and luminosity function (LF) of the resolved sources are more consistent with sources found in other early-type galaxies than with sources found in star-forming regions of galaxies. The source LF and the spectrum of the unresolved emission both indicate that there are a large number of unresolved point sources. We propose that these sources are associated with globular clusters (GCs) and that NGC 1600 has a large GC specific frequency. Observations of the GC population in NGC 1600 would be very useful for testing this prediction. Approximately 50%-75% of the unresolved flux comes from diffuse gaseous emission. The spectral fits, hardness ratios, and X-ray surface brightness profile all point to two gas components. We interpret the soft inner component (a<~25'', kT~0.85 keV) as the interstellar medium of NGC 1600 and the hotter outer component (a>~25'', kT~1.5 keV) as the intragroup medium of the NGC 1600 group. The X-ray image shows several interesting structures. First, there is a central region of excess emission that is roughly cospatial with Hα and dust filaments immediately west of the center of NGC 1600. There appear to be holes in the X-ray emission to the north and south of the
NASA Astrophysics Data System (ADS)
Healy, R. W.
2015-12-01
Water movement through soils is often dominated by preferential flow processes such as fingering and macropore flow. Traditional models of flow in the unsaturated zone are based on the diffusion or Richards equation and assume that diffusive (surface-tension viscous) flow is the only flow process. These models are incapable of accurately simulating preferential flow. Several alternative approaches, including kinematic wave, transfer function, and water-content wave models, have been suggested for simulating water movement through preferential flow paths. The source-responsive model proposed by Nimmo (2010) and Nimmo and Mitchell (2013) is unique among such models in that water transfer from land surface to depth is controlled by the water-application rate at land surface. The source-responsive model has been coupled with a one-dimensional version of the Richards-equation based model of variably saturated flow, VS2DT. The new model, can simulate flow within the preferential (S) domain alone, within the diffuse (D) domain alone, or within both the S and D domains simultaneously. Water exchange between the two domains is treated as a first-order diffusive process, with the exchange coefficient being a function of soil-water content. The new model was used to simulate field and laboratory infiltration experiments described in the literature. Simulations were calibrated against measured soil water contents with the PEST parameter estimation package; values for hydraulic conductivity and 3 van Genuchten and 3 source-responsive parameters were optimized. Although exact matches between measured and simulated water contents were not obtained, the simulation results captured the salient characteristics of the published data sets, including features typical of preferential as well as diffusive flow. Results obtained from simulating flow simultaneously in both the S and D domain provided better matches to measured data than results obtained from simulating flow independently
Diffusion of radon through concrete block walls: A significant source of indoor radon
Lively, R.S.; Goldberg, L.F.
1999-01-01
Basement modules located in southern Minnesota have been the site of continuous radon and environmental measurements during heating seasons since 1993. Concentrations of radon within the basement modules ranged from 70 Bq.m-3 to over 4000 Bq.m-3 between November to April during the three measurement periods. In the soil gas for the same times, concentrations of radon ranged between 25,000 and 70,000 Bq.m-3. Levels of radon within the basement modules changed by factors of five or more within 24 h, in concert with pressure gradients of 4 to 20 Pa that developed between the basement modules and their surroundings. Diffusion is identified as the principal method by which radon is transferred into and out of the basement modules, and appears to be relatively independent of insulating materials and vapour retarders. The variability of radon and correlations with differential pressure gradients may be related to air currents in the block walls and soil that interrupt radon diffusing inward. This yields a net decrease of radon in the basement modules by decay and outward diffusion. Levels of radon within the basement modules increase when the pressure differential is zero and air flow ceases, allowing diffusion gradients to be re-established. Radon levels in both the soil and the basement modules then increase until an equilibrium is achieved.
NASA Astrophysics Data System (ADS)
Wadle, Stephen; Wuest, Daniel; Cantalupo, John; Lakes, Roderic S.
1994-01-01
Holographic diffusers are prepared using silver halide (Agfa 8E75 and Kodak 649F) and photopolymer (Polaroid DMP 128 and DuPont 600, 705, and 150 series) media. It is possible to control the diffusion angle in three ways: by selection of the properties of the source diffuser, by control of its subtended angle, and by selection of the holographic medium. Several conventional diffusers based on refraction or scattering of light are examined for comparison.
NASA Astrophysics Data System (ADS)
Chang, Xiao-Chuan; Liu, Ruo-Yu; Wang, Xiang-Yu
2016-07-01
The nearly isotropic distribution of teraelectronvolt to petaelectronvolt neutrinos recently detected by the IceCube Collaboration suggests that they come from sources at a distance beyond our Galaxy, but how far away they are is largely unknown because of a lack of any associations with known sources. In this paper, we propose that the cumulative TeV gamma-ray emission accompanying the production of neutrinos can be used to constrain the distance of these neutrino sources, since the opacity of TeV gamma rays due to absorption by the extragalactic background light depends on the distance these TeV gamma rays have traveled. As the diffuse extragalactic TeV background measured by Fermi is much weaker than the expected cumulative flux associated with IceCube neutrinos, the majority of IceCube neutrinos, if their sources are transparent to TeV gamma rays, must come from distances larger than the horizon of TeV gamma rays. We find that above 80% of the IceCube neutrinos should come from sources at redshift z > 0.5. Thus, the chance of finding nearby sources correlated with IceCube neutrinos would be small. We also find that, to explain the flux of neutrinos under the TeV gamma-ray emission constraint, the redshift evolution of neutrino source density must be at least as fast as the cosmic star formation rate.
NASA Astrophysics Data System (ADS)
Di Mauro, Mattia;
2016-05-01
The isotropic diffuse γ-ray background (IGRB) has been detected by various experiments and recently the Fermi-LAT Collaboration has precisely measured its spectrum in a wide energy range. The origin of the IGRB is still unclear and we show in this paper the significative improvements that have been done, thanks to the new Fermi-LAT catalogs, to solve this mystery. We demonstrate that the γ-ray intensity and spectrum of the IGRB is fully consistent with the unresolved emission from extragalactic point sources, namely Active Galactic Nuclei and Star Forming Galaxies. We show also that the IGRB can be employed to derive sever constraints for the γ-ray emission from diffuse processes such as annihilation of Dark Matter (DM) particles. Our method is able to provide low bounds for the thermal annihilation cross section for a wide range of DM masses.
Atoms in carbon cages as a source of interstellar diffuse lines
NASA Technical Reports Server (NTRS)
Ballester, J. L.; Antoniewicz, P. R.; Smoluchowski, R.
1990-01-01
A model to describe the resonance absorption lines of various atoms trapped in closed carbon cages is presented. These systems may be responsible for some of the as yet unexplained diffuse interstellar bands. Model potentials for possible atom-C60 systems are obtained and used to calculate the resonance lines. The trapped atoms considered are O, N, Si, Mg, Al, Na, and S, and in all cases the resonance lines are shifted toward the red as compared to the isolated atoms. The calculated wavelengths are compared to the range of wavelengths observed for the diffuse interstellar bands, and good agreement is found for Mg and Si resonance lines. Other lines may be caused by other than resonance transitions or by trapped molecules. The oscillator strengths and the abundances are evaluated and compared with observation. Mechanisms to explain the observed band width of the lines and the existence of certain correlated pairs of lines are discussed.
Shi, Zhenzhi; Fan, Ying; Zhao, Huijuan; Xu, Kexin
2012-06-01
Accurate determination of the optical properties (the absorption coefficient μ(a) and the reduced scattering coefficient μ(s) (')) of tissues is very important in a variety of diagnostic and therapeutic procedures. Optical diffusion theory is frequently used as the forward model for describing the photon transfer in media with large reduced albedos (a(')) and in large source-detector separations (SDS). Several other methods (PN approximation, hybrid diffusion-P3 approximation) have also been published that describe photon transfer in media with low a(') or small SDSs. We studied the theoretical models for the steady-state spatially resolved diffuse reflectance measurement to accurately determine μ(a) and μ(s) (') at large a(') range but small SDSs. Instead of using a single model, a joint derivation method is proposed. The developed method uses one of the best aforementioned theoretical methods separately in five ranges of a(') determined from several forward models. In the region of small SDSs (the range between 0.4 and 8 mm) and large a(') range (between 0.5 and 0.99), the best theoretical derivation model was determined. The results indicate that the joint derivation method can improve the derivation accuracy and that a(') range can be determined by the steady-state spatially resolved diffuse reflectance measurement. PMID:22734782
NASA Astrophysics Data System (ADS)
Suard, Maxime
We studied the near field dispersion of natural gas plumes leaking from transmission lines and diffusing in a convective Atmospheric Boundary Layer (ABL), with the intent of providing an aerial system of leak detection and pinpointing, as well as quantitative leak rate estimation. We used high frequency measurements of methane and ethane concentrations on a fixed wing aircraft using high rate spectroscopic gas concentration measurements. We looked for characteristics of the effluent concentration field which can be related to the distance from the effluent source, and developed an empirical approach to effluent source position estimation from airborne effluent concentration measurements. From a mass-balance approach we developed a practical method of effluent leak rate estimation based on airborne effluent concentration measurements. Since gathering experimental data was costly and time-expensive, Large Eddy Simulation (LES) results were also investigated. Results showed that analysis of effluent concentration variability is likely to provide information about the position of the effluent source. The developed leak rate estimation method provided encouraging results showing that such an approach is able to yield relatively accurate leak rate estimates. LES results proved to be very useful as they helped to provide guidelines for experiments as well as to deepen our understanding of the diffusion dynamics of turbulent effluent plumes.
Greiner, R; Herr, A; Brodie, J; Haynes, D
2005-01-01
This paper presents a multi-criteria based tool for assessing the relative impact of diffuse-source pollution to the Great Barrier Reef (GBR) from the river basins draining into the GBR lagoon. The assessment integrates biophysical and ecological data of water quality and pollutant concentrations with socio-economic information pertaining to non-point source pollution and (potential) pollutant impact. The tool generates scores for each river basin against four criteria, thus profiling the basins and enabling prioritization of management alternatives between and within basins. The results support policy development for pollution control through community participation, scientific data integration and expert knowledge contributed by people from across the catchment. The results specifically provided support for the Reef Water Quality Protection Plan, released in October 2003. The aim of the plan is to provide a framework for reducing discharge of sediment, nutrient and other diffuse-source loads and (potential) impact of that discharge and for prioritising management actions both between and within river basins. PMID:15757715
NASA Technical Reports Server (NTRS)
Cohen, Martin; Jones, B. F.
1987-01-01
Spectroscopy of the starlike optical counterpart to IRAS 21282+5050, a source with the hydrocarbon infrared emission band spectrum, shows an 07(f)-(WC11) planetary nebula nucleus suffering an extinction of 5.7 mag. Emission line widths in the WC spectrum are only approx. 100 km/s, indicating a very slow stellar wind. Optical diffuse interstellar bands (DIBs) are prominent. Five DIBs are strongly enhanced, namely lamda lamda 5797, 6196, 6203, 6283, and 6613. The presence of circumstellar hydrocarbon molecules may explain both the infrared emission bands and the enhanced DIBs.
NASA Astrophysics Data System (ADS)
Cernak, M.; Kovacik, D.; Zahoranova, A.; Rahel, J.
2008-07-01
Collaborating Czech and Slovakian university teams have recently developed an innovative plasma source, the so-called Diffuse Coplanar Surface Barrier Discharge (DCSBD), which has the potential to move a step closer to the industry requirement for in-line treatment of low-added-value materials using a highly-nonequlibrium ambient air plasma (Simor et al. 2002, The idea is to generate a thin (on the order of 0.1 mm) layer of highly-nonequlibrium plasma with a high power density (up to 100 W/cm^3) in the immediate vicinity of the treated surface and bring it into a close contact with the treated surface. Comparing to atmospheric-pressure glow discharge, volume dielectric barrier discharge, and plasma jet plasmas, such a diffuse plasma layer is believed to provide substantial advantages in energy consumption, exposure time, and technical simplicity. A brief outline of physical mechanism and basic properties of DCSBD will given using the results of emission spectroscopy, high-speed camera, and spatially resolved cross-correlation spectroscopy studies. The presentation will review also a current state of the art in in-line plasma treatment of low-cost materials and opportunities for the use of the so-called Diffuse Coplanar Surface Dielectric Barrier Discharge (DCSBD). The results obtained on the ambient air plasma treatments of textile, paper, wood, and glass illustrate that DCSBD offers outstanding performance with extremely low energy consumption for large area, uniform surface modifications of materials under continuous process conditions.
NASA Astrophysics Data System (ADS)
Baksht, E. Kh; Burachenko, A. G.; Lomaev, M. I.; Panchenko, A. N.; Tarasenko, V. F.
2015-04-01
An extended repetitively pulsed source of spontaneous UV radiation is fabricated, which may also be used for producing laser radiation. Voltage pulses with an incident wave amplitude of up to 30 kV, a half-amplitude duration of ~4 ns and a rise time of ~2.5 ns are applied to a gap with a nonuniform electric field. For an excitation region length of 35 cm and a nitrogen pressure of 30 - 760 Torr, a diffusive discharge up to a pulse repetition rate of 2 kHz is produced without using an additional system for gap preionisation. An investigation is made of the plasma of the run-away electron preionised diffuse discharge. Using a CCD camera it is found that the dense diffused plasma fills the gap in a time shorter than 1 ns. X-ray radiation is recorded from behind the foil anode throughout the pressure range under study; a supershort avalanche electron beam is recorded by the collector electrode at pressures below 100 Torr.
Analysis of XMM-Newton Data from Extended Sources and the Diffuse X-Ray Background
NASA Technical Reports Server (NTRS)
Snowden, Steven
2011-01-01
Reduction of X-ray data from extended objects and the diffuse background is a complicated process that requires attention to the details of the instrumental response as well as an understanding of the multiple background components. We present methods and software that we have developed to reduce data from XMM-Newton EPIC imaging observations for both the MOS and PN instruments. The software has now been included in the Science Analysis System (SAS) package available through the XMM-Newton Science Operations Center (SOC).
NASA Astrophysics Data System (ADS)
Petlenko, A. V.; Kopytenko, Yu. A.
2015-01-01
The coherence conditions for local ionospheric sources of Pi2 high-latitude geomagnetic pulsations in the 6-8 mHz frequency band were studied with a method for representing their component spatial distributions by vertical magnetic dipole field integrals over random walk trajectories in the ionospheric horizontal plane. High-quality representations are reached when the field intensity is proportional to a self-intersection density that has spatial and time walk trajectory discontinuities. According to anomalous diffusion theory, the trajectory of the subordinate Levy process is characterized by the stabilization of turning points grouped near pulsation local ionospheric sources. This makes it possible to simulate disturbances of the Pi2 high-latitude pulsation field by a uniform motion along the Hall ionospheric source Levy trajectories with short-term pulsed current switching on-off. A coherent redistribution of the Pi2 local ionospheric source intensity and several observed specific features in the source field distribution dynamics are explained by the independent processes of local field-aligned current formation and intensification-weakening of pulsations caused by the current carrier shift relative to wandering Levy turning points and a stationary observer.
Diffuse radio emission around FR II sources as exemplified by 3C452
NASA Astrophysics Data System (ADS)
Wiita, Paul J.; Sirothia, S. K.; Gopal-Krishna, ..
2014-01-01
We have discovered a pair of megaparsec size radio lobes of extremely steep spectrum straddling the well-known classical double radio source 3C452. For the past several decades 3C452 has been regarded as a textbook example of an edge-brightened double radio source of Fanaroff-Riley type II (FR II) but we show it to be a bonafide "double-double" radio galaxy (DDRG). The inner double fed by the jets has evolved into a perfectly normal FR II radio source. Thus, 3C452 presents a uniquely robust example of recurrent nuclear activity in which the restarted jets are expanding non-relativistically within the relic synchrotron plasma from an earlier active phase. This situation contrasts markedly with the strikingly narrow inner doubles observed in a few other DDRGs that have been interpreted in terms of compression of the synchrotron plasma of the relic outer lobes at the relativistic bow-shocks driven by the near ballistic propagation of the two inner jets through the relic plasma. We also present additional examples of the occurrence of faded outer lobes around well defined FRII sources, using our deep GMRT images at meter wavelengths processed with AIPS++ software. We also examine the statistics of the occurrence of such sources using a flux density limited sample. A key ramification of our findings are that they caution against the use of FR II classical double radio sources for testing cosmological models and unification schemes for active galactic nuclei.
NASA Astrophysics Data System (ADS)
Ge, J.; Everett, M. E.; Weiss, C. J.
2011-12-01
A 2.5D finite difference (FD) frequency-domain modeling based on the Continuous Time Random Walk theory (CTRW) to the diffusion of electromagnetic fields generated by a transient loop source over a porous geological media is presented. The CTRW theory is demonstrated to be a powerful tool to concisely describe a transport process in fractal mediums with complex structures (Scher and Lax 1973, Metzler and Klafter 2000). In the controlled-source electromagnetic (CSEM) induction setting, the diffusion of Eddy currents in a rough medium is governed by the fractional differential EM equations in a CTRW sense (Everett 2009). We use staggered finite difference (FD) approach to discretize the equations in frequency domain. The biconjugate gradient iterative solver is implemented to solve sparse matrix equations. We compare the FD solutions with analytical solutions to a half space model to verify the accuracy of our code. Then FD solutions are calculated for a model with a 2D block fractured structures included in a homogeneous space. Via displaying the high sensitivity of the spatial variation of Ey fields to the geoelectrical properties and roughness of the block, we are able to show the potential of combination of CTRW theory with CSEM method in detecting unknown fractured targets in subsurface.
Interacting biochemical and diffusive controls on trace gas sources in unsaturated soils
NASA Astrophysics Data System (ADS)
Rubol, S.; Manzoni, S.; Bellin, A.; Porporato, A. M.
2011-12-01
Microbes react to environmental conditions on different timescales. When conditions improve (e.g., rewetting, substrate amendment), the residing population exits the dormant state, becomes active and starts synthesizing extra-cellular enzymes. If substrate availability, and hence energy, is sufficient, microbes may start to reproduce and increase the size of their population. These dynamics make it complicated to interpret measured relationships between microbial activity (e.g., respiration, denitrification, N mineralization) and environmental conditions. In particular, the relationship between bacterial activity and soil moisture, which is derived by incubating soil samples at constant soil moisture levels, seems to vary under dynamic hydrological conditions. This may be related to both soil physical properties and the resilience of bacteria to adapt to rapid changes in soil moisture. We present a process-based model that includes both the above effects and test the hypothesis that the ratio of the time scale of biological versus physical factors determines the shape describing the relationship between microbial activity and soil moisture. In particular, we focus on the role of oxygen dynamics, which regulate the prevalence of aerobic versus anaerobic conditions and thus the prevalence of nitrification versus denitrification. We identify and compare the time scale of the biological oxygen consumption with the time scale of physical diffusion. Starting from well-aerated conditions, as bacteria consume O2 in solution, more oxygen dissolves from the atmosphere - depending on gas-filled porosity. If water dynamics or tillage limits re-aeration, this can affect the equilibrium between the aqueous and the gaseous phase and thus alter the time scale of the reactions. This balance between consumption and re-aeration by diffusion ultimately controls the water quality as well the production of trace gases.
CONTRIBUTION OF UNRESOLVED POINT SOURCES TO THE DIFFUSE X-RAY BACKGROUND BELOW 1 keV
Gupta, A.; Galeazzi, M.
2009-09-01
We present here the analysis of X-ray point sources detected in several observations available in the XMM-Newton public archive. We focused, in particular, on energies below 1 keV, which are of particular relevance to the understanding of the diffuse X-ray background (DXB). The average field of all the exposures is 0.09 deg{sup -2}. We reached an average flux sensitivity of 5.8 x 10{sup -16}ergs{sup -1}cm{sup -2} in the soft band (0.5-2.0 keV) and 2.5 x 10{sup -16}ergs{sup -1}cm{sup -2} in the very soft band (0.4-0.6 keV). In this paper, we discuss the log N-log S results, the contribution to the integrated X-ray sky flux, and the properties of the cumulative spectrum from all sources. In particular, we found an excess flux at around 0.5 keV in the composite spectrum of faint sources. The excess seems to be a general property of all the fields observed suggesting an additional class of weak sources is contributing to the X-ray emission at these energies. Combining our results with previous investigations, we have also quantified the contribution of the individual components of the DXB in the 3/4 keV band.
An Exploratory Investigation of Diffused Point Arrival Time and Source Credibility.
ERIC Educational Resources Information Center
Baxter, Leslie A.; Ward, Jean M.
Because much conjecture (but limited empirical research) exists about the nonverbal variable to time, this exploratory study investigated the effects of differential arrival times on four dimensions of source credibility: sociability, dynamism, competence, and composure. Subjects were 84 educational secretaries with a mean age of 49 years. Each…
Source term evaluation for UF{sub 6} release event in feed facility at gaseous diffusion plants
Kim, S.H.; Taleyarkhan, R.P.
1997-01-30
An assessment of UF{sub 6} release accidents was conducted for the feed facility of a gaseous diffusion plant (GDP). Release rates from pig-tail connections were estimated from CYLIND code predictions, whereas, MELCOR was utilized for simulating reactions of UF{sub 6} with moisture and consequent transport of UO{sub 2}F{sub 2} aerosols and HF vapor through the building and to the environment. Two wind speeds were utilized. At the high end (Case 1) a wind speed of {approximately} 1 m/s (200 fpm) was assumed to flow parallel to the building length. At the low end (Case 2) to represent stagnant conditions a corresponding wind speed of 1 cm/s (2 fpm) was utilized. A further conservative assumption was made to specify no closure of crane and train doors at either end of the building. Relaxation of this assumption should provide for additional margins. Results indicated that, for the high (200 fpm) wind speed, close to 66% of the UO{sub 2}F{sub 2} aerosols and 100% of the HF gas get released to the environment over a 10-minute period. However, for the low (2 fpm) wind speed, negligible amount ({approximately} 1% UO{sub 2}F{sub 2}) of aerosols get released even over a 2 hour period.
Source stabilization for high quality time-domain diffuse optical tomography
NASA Astrophysics Data System (ADS)
Mo, Weirong; Chen, Nanguang
2009-02-01
We report a new close-loop feedback control method to keep a Mach-Zehnder electro-optic modulator (MZ-EOM) biased at the quadrature point and simultaneously correct the bias drift caused by the temperature changes as well as the inherent photorefractive effect. The modulator is a key part of our high speed time-domain diffuse optical tomography system. It modulates the dual-wavelength near-infrared light with the high speed pseudorandom bit sequence (PRBS) signal for the temporal point spread function (TPSF) measurements. Our method applies a periodical low frequency square wave with 50% duty cycle as the pilot tone upon the MZ-EOM together with the PRBS and sweep the bias voltage of the MZ-EOM in a self-adaptive step. A constant fraction of the modulated output power is measured by a photodiode via a tap coupler. After demodulation, the modulation depth versus the bias voltage can be measured from which the peak value corresponding to the quadrature point can be located quickly by curve fitting. Our stabilization technique is simple, fast and cost effective and is effective to correct the bias drift caused by the photorefractive and the change of ambient conditions. The experiment results show the TPSFs measurements can be stabilized to within +/-2% in an hour duration, which helps improved the image quality.
Delayed plastic relaxation limit in SiGe islands grown by Ge diffusion from a local source
Vanacore, G. M.; Zani, M.; Tagliaferri, A.; Nicotra, G.; Bollani, M.; Bonera, E.; Montalenti, F.; Picco, A.; Boioli, F.; Capellini, G.; Isella, G.; Osmond, J.
2015-03-14
The hetero-epitaxial strain relaxation in nano-scale systems plays a fundamental role in shaping their properties. Here, the elastic and plastic relaxation of self-assembled SiGe islands grown by surface-thermal-diffusion from a local Ge solid source on Si(100) are studied by atomic force and transmission electron microscopies, enabling the simultaneous investigation of the strain relaxation in different dynamical regimes. Islands grown by this technique remain dislocation-free and preserve a structural coherence with the substrate for a base width as large as 350 nm. The results indicate that a delay of the plastic relaxation is promoted by an enhanced Si-Ge intermixing, induced by the surface-thermal-diffusion, which takes place already in the SiGe overlayer before the formation of a critical nucleus. The local entropy of mixing dominates, leading the system toward a thermodynamic equilibrium, where non-dislocated, shallow islands with a low residual stress are energetically stable. These findings elucidate the role of the interface dynamics in modulating the lattice distortion at the nano-scale, and highlight the potential use of our growth strategy to create composition and strain-controlled nano-structures for new-generation devices.
Single-layer model to predict the source/sink behavior of diffusion-controlled building materials.
Kumar, Deept; Little, John C
2003-09-01
Building materials may act as both sources of and sinks forvolatile organic compounds (VOCs) in indoor air. A strategy to characterize the rate of absorption and desorption of VOCs by diffusion-controlled building materials is validated. A previously developed model that predicts mass transfer between a flat slab of material and the well-mixed air within a chamber or room is extended. The generalized model allows a nonuniform initial material-phase concentration and a transient influent gas-phase concentration to be simultaneously considered. An analytical solution to the more general model is developed. Experimental data are obtained by placing samples of vinyl flooring inside a small stainless steel chamber and exposing them to absorption/desorption cycles of n-dodecane and phenol. Measured values for the material-air partition coefficient and the material-phase diffusion coefficient were obtained previously in a series of completely independent experiments. The a priori model predictions are in close agreement with the observed experimental data. PMID:12967101
Delayed plastic relaxation limit in SiGe islands grown by Ge diffusion from a local source
NASA Astrophysics Data System (ADS)
Vanacore, G. M.; Nicotra, G.; Zani, M.; Bollani, M.; Bonera, E.; Montalenti, F.; Capellini, G.; Isella, G.; Osmond, J.; Picco, A.; Boioli, F.; Tagliaferri, A.
2015-03-01
The hetero-epitaxial strain relaxation in nano-scale systems plays a fundamental role in shaping their properties. Here, the elastic and plastic relaxation of self-assembled SiGe islands grown by surface-thermal-diffusion from a local Ge solid source on Si(100) are studied by atomic force and transmission electron microscopies, enabling the simultaneous investigation of the strain relaxation in different dynamical regimes. Islands grown by this technique remain dislocation-free and preserve a structural coherence with the substrate for a base width as large as 350 nm. The results indicate that a delay of the plastic relaxation is promoted by an enhanced Si-Ge intermixing, induced by the surface-thermal-diffusion, which takes place already in the SiGe overlayer before the formation of a critical nucleus. The local entropy of mixing dominates, leading the system toward a thermodynamic equilibrium, where non-dislocated, shallow islands with a low residual stress are energetically stable. These findings elucidate the role of the interface dynamics in modulating the lattice distortion at the nano-scale, and highlight the potential use of our growth strategy to create composition and strain-controlled nano-structures for new-generation devices.
Diffuse-source pesticide inputs in surface waters: online risk assessment at field scale
NASA Astrophysics Data System (ADS)
Reichenberger, S.; Röpke, B.; Bach, M.; Frede, H.-G.
2003-04-01
Diffuse pesticide input from agricultural fields into surface waters depends, apart from substance properties, also on soil properties, site hydrology, tillage and application practices, weather conditions, and distance to water bodies. A product which can be used safely at one site may pose unacceptable risk to aquatic life at another site. Thus, there is a great need for site-specific risk assessment approaches. Our objective is therefore to develop a tool for assessing pesticide inputs via runoff/erosion, drainflow, and spraydrift for single fields. The tool will be part of the web-based Information System for Integrated Plant Production (ISIP, www.isip.de), but also able to run independently. ISIP is an interactive portal for both farmers and advisors and will become the leading decision support system for plant production and plant protection in Germany over the next years. The aim is an approach which i) is applicable for the majority of European countries, ii) needs only input data readily available for large areas (such as provided by hydrological and soil maps), iii) predicts, with daily resolution, realistic (not worst-case) loads and resulting aquatic concentrations (PECsw), and iv) provides a substance-specific risk assessment on the basis of model results for consulting and management purposes. Special emphasis is laid on scientifically up-to-date model approaches and robust, but site-specific parameterization. For instance, drainage inputs will be calculated using the preferential flow model MACRO 5 (Nick Jarvis, SLU, Sweden). In the later stages of the project, the tool will provide a fully GIS-based risk assessment in ISIP for single fields in Germany.
NASA Astrophysics Data System (ADS)
Hoffmann, T. L.; Lieb, S.; Pauldrach, A. W. A.; Lesch, H.; Hultzsch, P. J. N.; Birk, G. T.
2012-08-01
Aims: The aim of this work is to verify whether turbulent magnetic reconnection can provide the additional energy input required to explain the up to now only poorly understood ionization mechanism of the diffuse ionized gas (DIG) in galaxies and its observed emission line spectra. Methods: We use a detailed non-LTE radiative transfer code that does not make use of the usual restrictive gaseous nebula approximations to compute synthetic spectra for gas at low densities. Excitation of the gas is via an additional heating term in the energy balance as well as by photoionization. Numerical values for this heating term are derived from three-dimensional resistive magnetohydrodynamic two-fluid plasma-neutral-gas simulations to compute energy dissipation rates for the DIG under typical conditions. Results: Our simulations show that magnetic reconnection can liberate enough energy to by itself fully or partially ionize the gas. However, synthetic spectra from purely thermally excited gas are incompatible with the observed spectra; a photoionization source must additionally be present to establish the correct (observed) ionization balance in the gas.
Moyes, Andrew B; Gaines, Sarah J; Siegwolf, Rolf T W; Bowling, David R
2010-11-01
Carbon isotope ratios (δ¹³C) of heterotrophic and rhizospheric sources of soil respiration under deciduous trees were evaluated over two growing seasons. Fluxes and δ¹³C of soil respiratory CO₂ on trenched and untrenched plots were calculated from closed chambers, profiles of soil CO₂ mole fraction and δ¹³C and continuous open chambers. δ¹³C of respired CO₂ and bulk carbon were measured from excised leaves and roots and sieved soil cores. Large diel variations (>5‰) in δ¹³C of soil respiration were observed when diel flux variability was large relative to average daily fluxes, independent of trenching. Soil gas transport modelling supported the conclusion that diel surface flux δ¹³C variation was driven by non-steady state gas transport effects. Active roots were associated with high summertime soil respiration rates and around 1‰ enrichment in the daily average δ¹³C of the soil surface CO₂ flux. Seasonal δ¹³C variability of about 4‰ (most enriched in summer) was observed on all plots and attributed to the heterotrophic CO₂ source. PMID:20545887
NASA Astrophysics Data System (ADS)
Tǎnasǎ, C.; Bosioc, A. I.; Susan-Resiga, R. F.; Muntean, S.
2012-11-01
The previous experimental and numerical investigations of decelerated swirling flows in conical diffusers have demonstrated that water injection along to the axis mitigates the pressure fluctuations associated to the precessing vortex rope [1]. However, for swirling flows similar to Francis turbines operated at partial discharge, the water jet becomes effective when the jet discharge is larger than 10% from the turbine discharge, leading to large volumetric losses when the jet is supplied from upstream the runner. As a result, it was introduced a new approach for supplying the jet by using a fraction of the discharge collected downstream the conical diffuser [2]. This is called flow-feedback control technique (FFCT) and it was investigated experimentally in order to assess its capability [3]. The FFCT approach not requires additional energy to supply the jet. Consequently, the turbine efficiency is not diminished due to the volumetric losses injected even if around 10% of the main flow is used. However, the equivalent amplitude of the pressure pulsations associated to the vortex rope decreases with 30% if 10% jet discharge is applied [3]. Using 12% water jet discharge from upstream then the equivalent amplitude of the pressure pulsations is mitigated with 70% according to Bosioc et al. [4]. In our case, an extra 2% jet discharge is required in order to obtain similar results with FFCT. This extra discharge is provided using an additional energy source. Therefore, the paper presents experimental investigation performed with FFCT with additional energy source. The experimental results obtained with this technique are compared against FFCT and the swirling flow with vortex rope, respectively.
Beach Sands Along the California Coast are Diffuse Sources of Fecal Bacteria to Coastal Waters
NASA Astrophysics Data System (ADS)
Boehm, A. B.; Yamahara, K.; Layton, B.
2007-05-01
Fecal indicator bacteria (FIB) are nearly ubiquitous in California (CA) beach sands. Sands were collected from 55 beaches along the CA coast. Ninety-one percent of the beaches had detectable enterococci (ENT) while 62% had detectable E. coli (EC) in their sands. The presence of a putative bacterial source (such as a river), the degree of wave shelter, and surrounding land use explained a significant (p<0.05) fraction of the variation in both ENT and EC densities between beaches. Sand characteristics including moisture content, organic carbon, and percent fines, significantly (p<0.05) influenced only EC densities in beach sand. We assayed 34 of 163 sand samples for salmonellae, but did not detect this bacterial pathogen. The potential for FIB to be transported from the sand to sea was investigated at a single wave-sheltered beach with high densities of ENT in beach sand: Lovers Point, CA (LP). We collected samples of exposed and submerged sands as well as water over a 24 h period in order to compare the disappearance or appearance of ENT in sand and the water column. Exposed sands had significantly higher densities of ENT than submerged sands with the highest densities located near the high tide line. Water column ENT densities began low, increased sharply during the first flood tide and slowly decreased over the remainder of the study. During the first flood tide, the number of ENT that entered the water column was nearly equivalent to the number of ENT lost from exposed sands when they were submerged by seawater. The decrease in nearshore ENT concentrations after the initial influx can be explained by ENT die-off and dilution with clean ocean water. A source tracking study at LP indicated that ENT were likely of human origin because they were positive for the esp gene.
Schreiber, H; Behrendt, H; Constantinescu, L T; Cvitanic, I; Drumea, D; Jabucar, D; Juran, S; Pataki, B; Snishko, S; Zessner, M
2005-01-01
Nutrient emissions by point and diffuse sources were estimated for 388 sub-catchments of the Danube river basin for the period 1998-2000 by means of the Model MONERIS. For nitrogen total emissions of 684 kt/a N were estimated for the Danube basin. 80% of these emissions were caused by diffuse sources (mainly groundwater, urban areas and tile drainage). For phosphorus the emission was 57 kt/a P, with a contribution of diffuse sources to this sum of 58%. The comparison of calculated and observed loads shows that the mean deviation for the investigated sub-catchments of the Danube river basin is 20% for dissolved inorganic nitrogen and 34% for phosphorus. The spatial resolution of the emission calculations allows the identification of regional hot spots and the derivation of specific regional measures to reduce the emissions into the Danube and consequently into the Western Black Sea. PMID:15850201
NASA Technical Reports Server (NTRS)
Boggs, S. E.; Lin, R. P.; Coburn, W.; Feffer, P.; Pelling, R. M.; Schroeder, P.; Slassi-Sennou, S.
1997-01-01
The balloon-borne high resolution gamma ray and X-ray germanium spectrometer (HIREGS) was used to observe the Galactic center and two positions along the Galactic plane from Antarctica in January 1995. For its flight, the collimators were configured to measure the Galactic diffuse hard X-ray continuum between 20 and 200 keV by directly measuring the point source contributions to the wide field of view flux for subtraction. The hard X-ray spectra of GX 1+4 and GRO J1655-40 were measured with the diffuse continuum subtracted off. The analysis technique for source separation is discussed and the preliminary separated spectra for these point sources and the Galactic diffuse emission are presented.
Benes, V.; Pekny, V.; Skorepa, J.; Vrba, J. )
1989-11-01
In several regions of Czechoslovakia with intensive agricultural production, the correlation between the amount of nitrogen fertilizer applied and the nitrate content in groundwater has been recognized. Nitrate pollution of groundwater is considered to be the most serious source of nonpoint pollution in Czechoslovakia. A program of research into the effects of farming activities on groundwater quality in Czechoslovakia is under way on experimental fields (20 to 30 hectares) and, simultaneously, in regions in which shallow, vulnerable aquifers occur. The importance of the soil organic matter's stability for maintaining the groundwater quality is emphasized. Regional investigation of the impact of farming on shallow aquifers in the fluvial deposits of the Elbe River in Bohemia has proved the hydrochemical instability and vertical hydrochemical heterogeneity of these aquifers. The WASTEN deterministic model was used for modeling the transport and transformation of various types of inorganic fertilizers. The input data is based on laboratory and field measurements. Special topics are the verification of model calculations and the time and spatial variability of input data with respect to the unsaturated zone. The research results are being used for making regional and national agro-groundwater managerial schemes more precise, as well as for decision-making.
Benes, V; Pĕkný, V; Skorepa, J; Vrba, J
1989-01-01
In several regions of Czechoslovakia with intensive agricultural production, the correlation between the amount of nitrogen fertilizer applied and the nitrate content in groundwater has been recognized. Nitrate pollution of groundwater is considered to be the most serious source of nonpoint pollution in Czechoslovakia. A program of research into the effects of farming activities on groundwater quality in Czechoslovakia is under way on experimental fields (20 to 30 hectares) and, simultaneously, in regions in which shallow, vulnerable aquifers occur. The importance of the soil organic matter's stability for maintaining the groundwater quality is emphasized. Research based on nitrogen and organic carbon balance has shown that the restoration of a soil-groundwater system is a complicated process that usually requires changes in the extent and intensity of agricultural activities and consistent attention to the effects produced by natural conditions. Regional investigation of the impact of farming on shallow aquifers in the fluvial deposits of the Elbe River in Bohemia has proved the hydrochemical instability and vertical hydrochemical heterogeneity of these aquifers. The WASTEN deterministic model was used for modeling the transport and transformation of various types of inorganic fertilizers. The input data is based on laboratory and field measurements. Special topics are the verification of model calculations and the time and spatial variability of input data with respect to the unsaturated zone. The research results are being used for making regional and national agro-groundwater managerial schemes more precise, as well as for decision-making. PMID:2559844
NASA Astrophysics Data System (ADS)
Santos, Frederico P.; Filho, Hermes Alves; Barros, Ricardo C.
2013-10-01
The scattering source iterative (SI) scheme is traditionally applied to converge fine-mesh numerical solutions to fixed-source discrete ordinates (SN) neutron transport problems. The SI scheme is very simple to implement under a computational viewpoint. However, the SI scheme may show very slow convergence rate, mainly for diffusive media (low absorption) with several mean free paths in extent (low leakage). In this work we describe an acceleration technique based on an improved initial guess for the scattering source distribution within the slab. In other words, we use as initial guess for the fine-mesh scattering source, the coarse-mesh solution of the neutron diffusion equation with special boundary conditions to account for the classical SN prescribed boundary conditions, including vacuum boundary conditions. Therefore, we first implement a spectral nodal method that generates coarse-mesh diffusion solution that is completely free from spatial truncation errors, then we reconstruct this coarse-mesh solution within each spatial cell of the discretization grid, to further yield the initial guess for the fine-mesh scattering source in the first SN transport sweep (forward and backward) across the spatial grid. We consider a number of numerical experiments to illustrate the efficiency of the offered diffusion synthetic acceleration (DSA) technique.
NASA Astrophysics Data System (ADS)
Thomas, Ian; Jordan, Phil; Mellander, Per-Erik; Fenton, Owen; Shine, Oliver; hUallacháin, Daire Ó.; Creamer, Rachel; McDonald, Noeleen; Dunlop, Paul; Murphy, Paul
2016-04-01
Identifying critical source areas (CSAs) of diffuse pollution in agricultural catchments requires the accurate identification of hydrologically sensitive areas (HSAs) at highest propensity for generating surface runoff and transporting pollutants such as phosphorus (P). A new GIS-based HSA Index is presented that identifies HSAs at the sub-field scale. It uses a soil topographic index (STI) and accounts for the hydrological disconnection of overland flow via topographic impediment from flow sinks such as hedgerows and depressions. High resolution (0.25-2 m) LiDAR Digital Elevation Models (DEMs) are utilised to capture these microtopographic controls on flow pathways and hydrological connectivity. The HSA Index was applied to four agricultural catchments (~7.5-12 km2) with contrasting topography and soil types. Catchment HSA sizes were estimated using high resolution rainfall-quickflow measurements during saturated winter storm events in 2009-2014, and mapped using the HSA Index. HSA sizes ranged from 1.6-3.4% of the catchment area during median storm events and 2.9-8.5% during upper quartile events depending on whether well or poorly drained soils dominated, which validated HSA Index value distributions. Total flow sink volume capacities ranged from 8,298-59,584 m3 and caused 8.5-24.2% of overland-flow-generating-areas and 16.8-33.4% of catchment areas to become hydrologically disconnected from the open drainage channel network. HSA maps identified 'delivery points' along surface runoff pathways where transported pollutants such as P are delivered to the open drainage network. Using these as proposed locations for targeting mitigation measures such as riparian buffer strips (RBS) reduced costs compared to blanket implementation within an example agri-environment scheme by 66% and 91% over 1 and 5 years respectively, which included LiDAR DEM acquisition costs. Considering that HSAs are often the dominant P CSA factor in agricultural catchments and can override
Diffusion of monochromatic classical waves.
Gerritsen, Sijmen; Bauer, Gerrit E W
2006-01-01
We study the diffusion of monochromatic classical waves in a disordered acoustic medium by scattering theory. In order to avoid artifacts associated with mathematical point scatterers, we model the randomness by small but finite insertions. We derive expressions for the configuration-averaged energy flux, energy density, and intensity for one-, two-, and three-dimensional (3D) systems with an embedded monochromatic source using the ladder approximation to the Bethe-Salpeter equation. We study the transition from ballistic to diffusive wave propagation and obtain results for the frequency dependence of the medium properties such as mean free path and diffusion coefficient as a function of the scattering parameters. We discover characteristic differences of the diffusion in 2D as compared to the conventional 3D case, such as an explicit dependence of the energy flux on the mean free path and quite different expressions for the effective transport velocity. PMID:16486306
NASA Astrophysics Data System (ADS)
Davis, J. A.; Guo, H.; Lai, B.; Kemner, K. M.; Ercius, P.; Fox, P. M.; Singer, D. M.; Minor, A.; Waychunas, G.
2012-12-01
) tomography, it can be seen that there are large numbers of pore throat sizes less than 10 nm within the coatings. We hypothesize that diffusion through these pores, which likely have electrically charged surfaces, controls the observed macroscopic rates of U(VI) sorption in batch experiments with sand grains. Evidence to support this hypothesis was observed by studying U and Fe fluorescence spatial variation within FIB samples (1 micron thick) at 200 nm spatial resolution. With this greater spatial resolution, it is possible to see U concentration variations within the coatings that are dependent on the time of sorption reaction, and illustrates how the coating environment constitutes a diffusion constraint to achieve adsorptive equilibrium between an aqueous phase and the mineral surfaces. Including this diffusion constraint within conceptual models for reactive contaminant transport may be significant at the field scale, because secondary mineral coatings are potentially both sinks and sources of contaminants depending on the history of a contaminated site. This is important in resolving long-term transport predictions at DOE sites, such as Hanford and Savannah River, where equilibrium versus kinetic reactive transport models are being evaluated.
NASA Astrophysics Data System (ADS)
Chen, Yung-Sheng; Wang, Jeng-Yau
2015-09-01
Light source plays a significant role to acquire a qualified image from objects for facilitating the image processing and pattern recognition. For objects possessing specular surface, the phenomena of reflection and halo appearing in the acquired image will increase the difficulty of information processing. Such a situation may be improved by the assistance of valuable diffuse light source. Consider reading resistor via computer vision, due to the resistor's specular reflective surface it will face with a severe non-uniform luminous intensity on image yielding a higher error rate in recognition without a well-controlled light source. A measurement system including mainly a digital microscope embedded in a replaceable diffuse cover, a ring-type LED embedded onto a small pad carrying a resistor for evaluation, and Arduino microcontrollers connected with PC, is presented in this paper. Several replaceable cost-effective diffuse covers made by paper bowl, cup and box inside pasted with white paper are presented for reducing specular reflection and halo effects and compared with a commercial diffuse some. The ring-type LED can be flexibly configured to be a full or partial lighting based on the application. For each self-made diffuse cover, a set of resistors with 4 or 5 color bands are captured via digital microscope for experiments. The signal-to-noise ratio from the segmented resistor-image is used for performance evaluation. The detected principal axis of resistor body is used for the partial LED configuration to further improve the lighting condition. Experimental results confirm that the proposed mechanism can not only evaluate the cost-effective diffuse light source but also be extended as an automatic recognition system for resistor reading.
Diffuse sources of human fecal pollution allow for the direct discharge of waste into receiving waters with minimal or no treatment. Traditional culture-based methods are commonly used to characterize fecal pollution in ambient waters, however these methods do not discern between...
Yang, Hao; Zhang, Tao; Zhou, Junli; Carney, Paul R.; Jiang, Huabei
2014-01-01
SUMMARY Objective The goal of this work is to establish a new dual-modal brain mapping technique based on diffuse optical tomography (DOT) and electroencephalographic source localization (ESL) that can chronically/intracranially record optical/EEG data to precisely map seizures and localize the seizure onset zone and associated epileptic brain network. Methods The dual-modal imaging system was employed to image seizures in an experimental acute bicuculline methiodide rat model of focal epilepsy. Depth information derived from DOT was used as constraint in ESL to enhance the image reconstruction. Groups of animals were compared based on localization of seizure foci, either at different positions or at different depths. Results This novel imaging technique successfully localized the seizure onset zone in rat induced by bicuculline methiodide injected at a depth of 1mm, 2mm and 3mm, respectively. The results demonstrated that the incorporation of the depth information from DOT into the ESL image reconstruction resulted in more accurate and reliable ESL images. Although the ESL images showed a horizontal shift of the source localization, the DOT identified the seizure focus accurately. In one case, when the BMI was injected at a site outside the field of view (FOV) of the DOT/ESL interface, ESL gives false positive detection of the focus while DOT shows negative detection. Significance This study represents the first to identify seizure onset zone using implantable DOT. In addition, the combination of DOT/ESL has never been documented in neuroscience and epilepsy imaging. This technology will enable us to precisely measure the neural activity and hemodynamic response at exactly the same tissue site and at both cortical and sub cortical levels. PMID:25524046
NASA Astrophysics Data System (ADS)
Karawacki, Ernest; Suleiman, Bashir M.; ul-Haq, Izhar; Nhi, Bui-Thi
1992-10-01
The recently developed dynamic plane source (DPS) technique for simultaneous determination of the thermal properties of fast thermally conducting materials with thermal conductivities between 200 and 2 W/mK has now been extended for studying relatively slow conducting materials with thermal conductivities equal or below 2 W/mK. The method is self-checking since the thermal conductivity, thermal diffusivity specific heat, and effusivity of the material are obtained independently from each other. The theory of the technique and the experimental arrangement are given in detail. The data evaluation procedure is simple and makes it possible to reveal the distortions due to the nonideal experimental conditions. The extension to the DPS technique has been implemented at room temperature to study the samples of cordierite-based ceramic Cecorite 130P (thermal conductivity equal to 1.48 W/mK), rubber (0.403 W/mK), and polycarbonate (0.245 W/mK). The accuracy of the method is within ±5%.
NASA Technical Reports Server (NTRS)
Nicolet, M. A.
1983-01-01
The choice of the metallic film for the contact to a semiconductor device is discussed. One way to try to stabilize a contact is by interposing a thin film of a material that has low diffusivity for the atoms in question. This thin film application is known as a diffusion barrier. Three types of barriers can be distinguished. The stuffed barrier derives its low atomic diffusivity to impurities that concentrate along the extended defects of a polycrystalline layer. Sacrificial barriers exploit the fact that some (elemental) thin films react in a laterally uniform and reproducible fashion. Sacrificial barriers have the advantage that the point of their failure is predictable. Passive barriers are those most closely approximating an ideal barrier. The most-studied case is that of sputtered TiN films. Stuffed barriers may be viewed as passive barriers whose low diffusivity material extends along the defects of the polycrystalline host.
Cuthbert, M.O.; Mackay, R.; Nimmo, J.R.
2012-01-01
Results are presented of a detailed study into the vadose zone and shallow water table hydrodynamics of a field site in Shropshire, UK. A conceptual model is developed and tested using a range of numerical models, including a modified soil moisture balance model (SMBM) for estimating groundwater recharge in the presence of both diffuse and preferential flow components. Tensiometry reveals that the loamy sand topsoil wets up via macropore flow and subsequent redistribution of moisture into the soil matrix. Recharge does not occur until near-positive pressures are achieved at the top of the sandy glaciofluvial outwash material that underlies the topsoil, about 1 m above the water table. Once this occurs, very rapid water table rises follow. This threshold behaviour is attributed to the vertical discontinuity in the macropore system due to seasonal ploughing of the topsoil, and a lower permeability plough/iron pan restricting matrix flow between the topsoil and the lower outwash deposits. Although the wetting process in the topsoil is complex, a SMBM is shown to be effective in predicting the initiation of preferential flow from the base of the topsoil into the lower outwash horizon. The rapidity of the response at the water table and a water table rise during the summer period while flow gradients in the unsaturated profile were upward suggest that preferential flow is also occurring within the outwash deposits below the topsoil. A variation of the source-responsive model proposed by Nimmo (2010) is shown to reproduce the observed water table dynamics well in the lower outwash horizon when linked to a SMBM that quantifies the potential recharge from the topsoil. The results reveal new insights into preferential flow processes in cultivated soils and provide a useful and practical approach to accounting for preferential flow in studies of groundwater recharge estimation.
NASA Astrophysics Data System (ADS)
Kobulnicky, Henry A.; Babler, Brian L.; Alexander, Michael J.; Meade, Marilyn R.; Whitney, Barbara A.; Churchwell, Edward B.
2013-07-01
We characterize the completeness of point source lists from Spitzer Space Telescope surveys in the four Infrared Array Camera (IRAC) bandpasses, emphasizing the Galactic Legacy Infrared Mid-Plane Survey Extraordinaire (GLIMPSE) programs (GLIMPSE I, II, 3D, 360; Deep GLIMPSE) and their resulting point source Catalogs and Archives. The analysis separately addresses effects of incompleteness resulting from high diffuse background emission and incompleteness resulting from point source confusion (i.e., crowding). An artificial star addition and extraction analysis demonstrates that completeness is strongly dependent on local background brightness and structure, with high-surface-brightness regions suffering up to five magnitudes of reduced sensitivity to point sources. This effect is most pronounced at the IRAC 5.8 and 8.0 μm bands where UV-excited polycyclic aromatic hydrocarbon emission produces bright, complex structures (photodissociation regions). With regard to diffuse background effects, we provide the completeness as a function of stellar magnitude and diffuse background level in graphical and tabular formats. These data are suitable for estimating completeness in the low-source-density limit in any of the four IRAC bands in GLIMPSE Catalogs and Archives and some other Spitzer IRAC programs that employ similar observational strategies and are processed by the GLIMPSE pipeline. By performing the same analysis on smoothed images we show that the point source incompleteness is primarily a consequence of structure in the diffuse background emission rather than photon noise. With regard to source confusion in the high-source-density regions of the Galactic Plane, we provide figures illustrating the 90% completeness levels as a function of point source density at each band. We caution that completeness of the GLIMPSE 360/Deep GLIMPSE Catalogs is suppressed relative to the corresponding Archives as a consequence of rejecting stars that lie in the point
Kobulnicky, Henry A.; Alexander, Michael J.; Babler, Brian L.; Meade, Marilyn R.; Whitney, Barbara A.; Churchwell, Edward B. E-mail: malexan9@uwyo.edu E-mail: meade@astro.wisc.edu E-mail: ebc@astro.wisc.edu
2013-07-01
We characterize the completeness of point source lists from Spitzer Space Telescope surveys in the four Infrared Array Camera (IRAC) bandpasses, emphasizing the Galactic Legacy Infrared Mid-Plane Survey Extraordinaire (GLIMPSE) programs (GLIMPSE I, II, 3D, 360; Deep GLIMPSE) and their resulting point source Catalogs and Archives. The analysis separately addresses effects of incompleteness resulting from high diffuse background emission and incompleteness resulting from point source confusion (i.e., crowding). An artificial star addition and extraction analysis demonstrates that completeness is strongly dependent on local background brightness and structure, with high-surface-brightness regions suffering up to five magnitudes of reduced sensitivity to point sources. This effect is most pronounced at the IRAC 5.8 and 8.0 {mu}m bands where UV-excited polycyclic aromatic hydrocarbon emission produces bright, complex structures (photodissociation regions). With regard to diffuse background effects, we provide the completeness as a function of stellar magnitude and diffuse background level in graphical and tabular formats. These data are suitable for estimating completeness in the low-source-density limit in any of the four IRAC bands in GLIMPSE Catalogs and Archives and some other Spitzer IRAC programs that employ similar observational strategies and are processed by the GLIMPSE pipeline. By performing the same analysis on smoothed images we show that the point source incompleteness is primarily a consequence of structure in the diffuse background emission rather than photon noise. With regard to source confusion in the high-source-density regions of the Galactic Plane, we provide figures illustrating the 90% completeness levels as a function of point source density at each band. We caution that completeness of the GLIMPSE 360/Deep GLIMPSE Catalogs is suppressed relative to the corresponding Archives as a consequence of rejecting stars that lie in the point
Vanacore, G M; Zani, M; Bollani, M; Bonera, E; Nicotra, G; Osmond, J; Capellini, Giovanni; Isella, G; Tagliaferri, A
2014-04-01
In this paper we experimentally study the growth of self-assembled SiGe islands formed on Si(001) by exploiting the thermally activated surface diffusion of Ge atoms from a local Ge source stripe in the temperature range 600-700 °C. This new growth strategy allows us to vary continuously the Ge coverage from 8 to 0 monolayers as the distance from the source increases, and thus enables the investigation of the island growth over a wide range of dynamical regimes at the same time, providing a unique birds eye view of the factors governing the growth process and the dominant mechanism for the mass collection by a critical nucleus. Our results give experimental evidence that the nucleation process evolves within a diffusion limited regime. At a given annealing temperature, we find that the nucleation density depends only on the kinetics of the Ge surface diffusion resulting in a universal scaling distribution depending only on the Ge coverage. An analytical model is able to reproduce quantitatively the trend of the island density. Following the nucleation, the growth process appears to be driven mainly by short-range interactions between an island and the atoms diffusing within its vicinities. The islands volume distribution is, in fact, well described in the whole range of parameters by the Mulheran's capture zone model. The complex growth mechanism leads to a strong intermixing of Si and Ge within the island volume. Our growth strategy allows us to directly investigate the correlation between the Si incorporation and the Ge coverage in the same experimental conditions: higher intermixing is found for lower Ge coverage. This confirms that, besides the Ge gathering from the surface, also the Si incorporation from the substrate is driven by the diffusion kinetics, thus imposing a strict constraint on the initial Ge coverage, its diffusion properties and the final island volume. PMID:24594569
NASA Technical Reports Server (NTRS)
1981-01-01
A diffuse celestial radiation which is isotropic at least on a course scale were measured from the soft X-ray region to about 150 MeV, at which energy the intensity falls below that of the galactic emission for most galactic latitudes. The spectral shape, the intensity, and the established degree of isotropy of this diffuse radiation already place severe constraints on the possible explanations for this radiation. Among the extragalactic theories, the more promising explanations of the isotropic diffuse emission appear to be radiation from exceptional galaxies from matter antimatter annihilation at the boundaries of superclusters of galaxies of matter and antimatter in baryon symmetric big bang models. Other possible sources for extragalactic diffuse gamma radiation are discussed and include normal galaxies, clusters of galaxies, primordial cosmic rays interacting with intergalactic matter, primordial black holes, and cosmic ray leakage from galaxies.
NASA Astrophysics Data System (ADS)
Plucinsky, P. P.; Snowden, S. L.; Briel, U. G.; Hasinger, G.; Pfeffermann, E.
1993-11-01
In order to permit quantitative studies of the cosmic diffuse X-ray background (DXRB) and of extended X-ray sources, we present updated calibrations of the particle-induced background of the Position Sensitive Proportional Counters (PSPCs) on board the Röntgen Satellite (ROSAT). We present new parameterizations of the temporal, spectral, and spatial distributions of the particle-induced events following closely the analysis discussed in Snowden et al. (1992). The ROSAT Guest Observer (GO) may find a step-by-step method for applying these parameterizations to a GO observation in § 3.4. Except for a variable contamination which is present in channels ≤ 18 and a change in our understanding of the externally produced components, the current parameterizations are quite similar to the previous results. We have used the spectral information available on the variable contamination to formulate a method for determining the level of this contamination in a given observation. The PSPC rejection efficiency for particle background events in the pulse-height range 18 ≤ CH ≤ 249 is 99.90%, with a typical count rate of 4 × 10-6 counts s-1 arcmin-2 keV-1. During typical conditions, the count rate of residual events is well correlated with the Master Veto count rate. The spectrum in the pulse-height range 18 ≤ CH ≤ 249 is well described by a power law, a flat component, and an Al Kα line at 1.5 keV. The spatial distribution of counts with pulse heights ≥ 18 is uniform over the field of view except for a small radial gradient and shadowing of the Al Kα line and part of the flat continuum by the window support structure. During an astronomical observation in low-gain mode (after 1991 October 11), the particle background can also be monitored by the count rate in channels 260 ≤ CH ≤ 370, since in most cases all these events are produced by particles. We have used a 54 ks observation of the Ursa Major region to verify the accuracy of our model. We have also
ERIC Educational Resources Information Center
Michaelides, Michalis P.; Haertel, Edward H.
2014-01-01
The standard error of equating quantifies the variability in the estimation of an equating function. Because common items for deriving equated scores are treated as fixed, the only source of variability typically considered arises from the estimation of common-item parameters from responses of samples of examinees. Use of alternative, equally…
NASA Astrophysics Data System (ADS)
Putze, A.; Maurin, D.; Donato, F.
2011-02-01
Context. The source spectrum of cosmic rays is not well determined by diffusive shock acceleration models. The propagated fluxes of proton, helium, and heavier primary cosmic-ray species (up to Fe) are a means to indirectly access it. But how robust are the constraints, and how degenerate are the source and transport parameters? Aims: We check the compatibility of the primary fluxes with the transport parameters derived from the B/C analysis, but also ask whether they add further constraints. We study whether the spectral shapes of these fluxes and their ratios are mostly driven by source or propagation effects. We then derive the source parameters (slope, abundance, and low-energy shape). Methods: Simple analytical formulae are used to address the issue of degeneracies between source/transport parameters, and to understand the shape of the p/He and C/O to Fe/O data. The full analysis relies on the USINE propagation package, the MINUIT minimisation routines (χ2 analysis) and a Markov Chain Monte Carlo (MCMC) technique. Results: Proton data are well described in the simplest model defined by a power-law source spectrum and plain diffusion. They can also be accommodated by models with, e.g., convection and/or reacceleration. There is no need for breaks in the source spectral indices below ~1 TeV/n. Fits to the primary fluxes alone do not provide physical constraints on the transport parameters. If we leave the source spectrum free, parametrised by the form dQ/dE = q βη_S R-α, and fix the diffusion coefficient K(R) = K_0βη_T Rδ so as to reproduce the B/C ratio, the MCMC analysis constrains the source spectral index α to be in the range 2.2-2.5 for all primary species up to Fe, regardless of the value of the diffusion slope δ. The values of the parameter ηS describing the low-energy shape of the source spectrum are degenerate with the parameter ηT describing the low-energy shape of the diffusion coefficient: we find ηS - ηT ≈ 0 for p and He data, but
Technology Transfer Automated Retrieval System (TEKTRAN)
Legacy phosphorus (P), the accumulation of P in soils and sediments due to past agricultural management activities, represents an emerging challenge to ongoing efforts to mitigate diffuse P pollution from agriculture. Nutrient management programs, already tasked with minimizing the effects of today...
The paper proposes three alternative, diffusion-limited mathematical models to account for volatile organic compound (VOC) interactions with indoor sinks, using the linear isotherm model as a reference point. (NOTE: Recent reports by both the U.S. EPA and a study committee of the...
Chen, Yu-Wen; Chen, Chien-Chih; Huang, Po-Jung; Tseng, Sheng-Hao
2016-01-01
Diffuse reflectance spectroscopy (DRS) based on the frequency-domain (FD) technique has been employed to investigate the optical properties of deep tissues such as breast and brain using source to detector separation up to 40 mm. Due to the modeling and system limitations, efficient and precise determination of turbid sample optical properties from the FD diffuse reflectance acquired at a source-detector separation (SDS) of around 1 mm has not been demonstrated. In this study, we revealed that at SDS of 1 mm, acquiring FD diffuse reflectance at multiple frequencies is necessary for alleviating the influence of inevitable measurement uncertainty on the optical property recovery accuracy. Furthermore, we developed artificial neural networks (ANNs) trained by Monte Carlo simulation generated databases that were capable of efficiently determining FD reflectance at multiple frequencies. The ANNs could work in conjunction with a least-square optimization algorithm to rapidly (within 1 second), accurately (within 10%) quantify the sample optical properties from FD reflectance measured at SDS of 1 mm. In addition, we demonstrated that incorporating the steady-state apparatus into the FD DRS system with 1 mm SDS would enable obtaining broadband absorption and reduced scattering spectra of turbid samples in the wavelength range from 650 to 1000 nm. PMID:27446671
Zeinali-Rafsanjani, B.; Mosleh-Shirazi, M. A.; Faghihi, R.; Karbasi, S.; Mosalaei, A.
2015-01-01
To accurately recompute dose distributions in chest-wall radiotherapy with 120 kVp kilovoltage X-rays, an MCNP4C Monte Carlo model is presented using a fast method that obviates the need to fully model the tube components. To validate the model, half-value layer (HVL), percentage depth doses (PDDs) and beam profiles were measured. Dose measurements were performed for a more complex situation using thermoluminescence dosimeters (TLDs) placed within a Rando phantom. The measured and computed first and second HVLs were 3.8, 10.3 mm Al and 3.8, 10.6 mm Al, respectively. The differences between measured and calculated PDDs and beam profiles in water were within 2 mm/2% for all data points. In the Rando phantom, differences for majority of data points were within 2%. The proposed model offered an approximately 9500-fold reduced run time compared to the conventional full simulation. The acceptable agreement, based on international criteria, between the simulations and the measurements validates the accuracy of the model for its use in treatment planning and radiobiological modeling studies of superficial therapies including chest-wall irradiation using kilovoltage beam. PMID:26170553
Zeinali-Rafsanjani, B; Mosleh-Shirazi, M A; Faghihi, R; Karbasi, S; Mosalaei, A
2015-01-01
To accurately recompute dose distributions in chest-wall radiotherapy with 120 kVp kilovoltage X-rays, an MCNP4C Monte Carlo model is presented using a fast method that obviates the need to fully model the tube components. To validate the model, half-value layer (HVL), percentage depth doses (PDDs) and beam profiles were measured. Dose measurements were performed for a more complex situation using thermoluminescence dosimeters (TLDs) placed within a Rando phantom. The measured and computed first and second HVLs were 3.8, 10.3 mm Al and 3.8, 10.6 mm Al, respectively. The differences between measured and calculated PDDs and beam profiles in water were within 2 mm/2% for all data points. In the Rando phantom, differences for majority of data points were within 2%. The proposed model offered an approximately 9500-fold reduced run time compared to the conventional full simulation. The acceptable agreement, based on international criteria, between the simulations and the measurements validates the accuracy of the model for its use in treatment planning and radiobiological modeling studies of superficial therapies including chest-wall irradiation using kilovoltage beam. PMID:26170553
Suzaku observations of two diffuse hard X-ray source regions, G22.0+0.0 and G23.5+0.1
NASA Astrophysics Data System (ADS)
Yamauchi, Shigeo; Sumita, Mayu; Bamba, Aya
2016-06-01
G22.0+0.0 and G23.5+0.1 are diffuse hard X-ray sources discovered in the ASCA Galactic Plane Survey. We present Suzaku results of spectral analysis for these sources. G22.0+0.0 is confirmed to be a largely extended emission. Its spectra were represented by a highly absorbed power-law model with a photon index of 1.7 ± 0.3 and a moderately absorbed thermal emission with a temperature of 0.34^{+0.11}_{-0.08}keV. The difference in the NH values between the two components suggests that the thermal component is unrelated to the power-law component and is a foreground emission located in the same line of sight. G23.5+0.1 is an extended source with a size of ˜3{^'.}5. Its spectra were fitted with an absorbed power-law model with a photon index of 2.4^{+0.5}_{-0.4}. The spatial and spectral properties show that both sources are candidates for old pulsar wind nebulae (PWNe). In addition to the extended sources, we analyzed spectra of three point sources found in the observed fields. Based on the spectral features, we discuss the origin of the sources.
Handbook on atmospheric diffusion
Hanna, S.R.; Briggs, G.A.; Hosker, R.P. Jr.
1982-01-01
Basic meteorological concepts are covered as well as plume rise, source effects, and diffusion models. Chapters are included on cooling tower plumes and urban diffusion. Suggestions are given for calculating diffusion in special situations, such as for instantaneous releases over complex terrain, over long distances, and during times when chemical reactions or dry or wet deposition are important. (PSB)
Doping of ZnO nanowires using phosphorus diffusion from a spin-on doped glass source
Bocheux, A.; Robin, I. C.; Bonaimé, J.; Hyot, B.; Feuillet, G.; Kolobov, A. V.; Fons, P.; Mitrofanov, K. V.; Tominaga, J.; Tamenori, Y.
2014-05-21
In this article, we report on ZnO nanowires that were phosphorus doped using a spin on dopant glass deposition and diffusion method. Photoluminescence measurements suggest that this process yields p-doped ZnO. The spatial location of P atoms was studied using x-ray near-edge absorption structure spectroscopy and it is concluded that the doping is amphoteric with P atoms located on both Zn and O sites.
Doping of ZnO nanowires using phosphorus diffusion from a spin-on doped glass source
NASA Astrophysics Data System (ADS)
Bocheux, A.; Robin, I. C.; Bonaimé, J.; Hyot, B.; Kolobov, A. V.; Mitrofanov, K. V.; Fons, P.; Tominaga, J.; Tamenori, Y.; Feuillet, G.
2014-05-01
In this article, we report on ZnO nanowires that were phosphorus doped using a spin on dopant glass deposition and diffusion method. Photoluminescence measurements suggest that this process yields p-doped ZnO. The spatial location of P atoms was studied using x-ray near-edge absorption structure spectroscopy and it is concluded that the doping is amphoteric with P atoms located on both Zn and O sites.