High-order solution methods for grey discrete ordinates thermal radiative transfer

DOE PAGES

Maginot, Peter G.; Ragusa, Jean C.; Morel, Jim E.

2016-09-29

This paper presents a solution methodology for solving the grey radiative transfer equations that is both spatially and temporally more accurate than the canonical radiative transfer solution technique of linear discontinuous finite element discretization in space with implicit Euler integration in time. We solve the grey radiative transfer equations by fully converging the nonlinear temperature dependence of the material specific heat, material opacities, and Planck function. The grey radiative transfer equations are discretized in space using arbitrary-order self-lumping discontinuous finite elements and integrated in time with arbitrary-order diagonally implicit Runge–Kutta time integration techniques. Iterative convergence of the radiation equation ismore » accelerated using a modified interior penalty diffusion operator to precondition the full discrete ordinates transport operator.« less

Modeling boundary measurements of scattered light using the corrected diffusion approximation

PubMed Central

Lehtikangas, Ossi; Tarvainen, Tanja; Kim, Arnold D.

2012-01-01

We study the modeling and simulation of steady-state measurements of light scattered by a turbid medium taken at the boundary. In particular, we implement the recently introduced corrected diffusion approximation in two spatial dimensions to model these boundary measurements. This implementation uses expansions in plane wave solutions to compute boundary conditions and the additive boundary layer correction, and a finite element method to solve the diffusion equation. We show that this corrected diffusion approximation models boundary measurements substantially better than the standard diffusion approximation in comparison to numerical solutions of the radiative transport equation. PMID:22435102

Entropy-based artificial viscosity stabilization for non-equilibrium Grey Radiation-Hydrodynamics

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

Delchini, Marc O., E-mail: delchinm@email.tamu.edu; Ragusa, Jean C., E-mail: jean.ragusa@tamu.edu; Morel, Jim, E-mail: jim.morel@tamu.edu

2015-09-01

The entropy viscosity method is extended to the non-equilibrium Grey Radiation-Hydrodynamic equations. The method employs a viscous regularization to stabilize the numerical solution. The artificial viscosity coefficient is modulated by the entropy production and peaks at shock locations. The added dissipative terms are consistent with the entropy minimum principle. A new functional form of the entropy residual, suitable for the Radiation-Hydrodynamic equations, is derived. We demonstrate that the viscous regularization preserves the equilibrium diffusion limit. The equations are discretized with a standard Continuous Galerkin Finite Element Method and a fully implicit temporal integrator within the MOOSE multiphysics framework. The methodmore » of manufactured solutions is employed to demonstrate second-order accuracy in both the equilibrium diffusion and streaming limits. Several typical 1-D radiation-hydrodynamic test cases with shocks (from Mach 1.05 to Mach 50) are presented to establish the ability of the technique to capture and resolve shocks.« less

Accurate radiative transfer calculations for layered media.

PubMed

Selden, Adrian C

2016-07-01

Simple yet accurate results for radiative transfer in layered media with discontinuous refractive index are obtained by the method of K-integrals. These are certain weighted integrals applied to the angular intensity distribution at the refracting boundaries. The radiative intensity is expressed as the sum of the asymptotic angular intensity distribution valid in the depth of the scattering medium and a transient term valid near the boundary. Integrated boundary equations are obtained, yielding simple linear equations for the intensity coefficients, enabling the angular emission intensity and the diffuse reflectance (albedo) and transmittance of the scattering layer to be calculated without solving the radiative transfer equation directly. Examples are given of half-space, slab, interface, and double-layer calculations, and extensions to multilayer systems are indicated. The K-integral method is orders of magnitude more accurate than diffusion theory and can be applied to layered scattering media with a wide range of scattering albedos, with potential applications to biomedical and ocean optics.

Quantifying equation-of-state and opacity errors using integrated supersonic diffusive radiation flow experiments on the National Ignition Facility

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

Guymer, T. M., E-mail: Thomas.Guymer@awe.co.uk; Moore, A. S.; Morton, J.

A well diagnosed campaign of supersonic, diffusive radiation flow experiments has been fielded on the National Ignition Facility. These experiments have used the accurate measurements of delivered laser energy and foam density to enable an investigation into SESAME's tabulated equation-of-state values and CASSANDRA's predicted opacity values for the low-density C{sub 8}H{sub 7}Cl foam used throughout the campaign. We report that the results from initial simulations under-predicted the arrival time of the radiation wave through the foam by ≈22%. A simulation study was conducted that artificially scaled the equation-of-state and opacity with the intended aim of quantifying the systematic offsets inmore » both CASSANDRA and SESAME. Two separate hypotheses which describe these errors have been tested using the entire ensemble of data, with one being supported by these data.« less

On the renormalisation of the diffusion asymptotics in the problem of reflection of a narrow optical beam from a biological medium

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

Appanov, A Yu; Barabanenkov, Yu N

2005-12-31

An analytic hybrid method is considered for solving the stationary radiation transfer equation in the problem on reflection of a narrow laser beam from biological media such as the 2% aqueous solution of intralipid and erythrocyte suspension with the volume concentration (hematocrit) H=0.41. The method is based on the reciprocity of the Green function in the radiation transfer theory and on the iteration solution of the integral equation for this function. As a result, the ray intensity is represented as a sum of two terms. The first of them describes the contribution of finite-order scattering to the intensity of amore » beam diffusely reflected from the medium. The second term contains the explicit analytic expression for a spatially distributed effective source of diffuse radiation emerging from the deep layers of the medium to the surface. This approach substantially improves the diffusion approximation for the problem under study and allows one to obtain the uniform asymptotics of the reflection coefficient at the specified interval of distances between the radiation source and detector on the medium surface with the relative error within {+-}6% for the 2% intralipid emulsion and erythrocyte suspension (H=0.41). (radiation scattering)« less

Numerical applications of the advective-diffusive codes for the inner magnetosphere

NASA Astrophysics Data System (ADS)

Aseev, N. A.; Shprits, Y. Y.; Drozdov, A. Y.; Kellerman, A. C.

2016-11-01

In this study we present analytical solutions for convection and diffusion equations. We gather here the analytical solutions for the one-dimensional convection equation, the two-dimensional convection problem, and the one- and two-dimensional diffusion equations. Using obtained analytical solutions, we test the four-dimensional Versatile Electron Radiation Belt code (the VERB-4D code), which solves the modified Fokker-Planck equation with additional convection terms. The ninth-order upwind numerical scheme for the one-dimensional convection equation shows much more accurate results than the results obtained with the third-order scheme. The universal limiter eliminates unphysical oscillations generated by high-order linear upwind schemes. Decrease in the space step leads to convergence of a numerical solution of the two-dimensional diffusion equation with mixed terms to the analytical solution. We compare the results of the third- and ninth-order schemes applied to magnetospheric convection modeling. The results show significant differences in electron fluxes near geostationary orbit when different numerical schemes are used.

An asymptotic-preserving stochastic Galerkin method for the radiative heat transfer equations with random inputs and diffusive scalings

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

Jin, Shi, E-mail: sjin@wisc.edu; Institute of Natural Sciences, Department of Mathematics, MOE-LSEC and SHL-MAC, Shanghai Jiao Tong University, Shanghai 200240; Lu, Hanqing, E-mail: hanqing@math.wisc.edu

2017-04-01

In this paper, we develop an Asymptotic-Preserving (AP) stochastic Galerkin scheme for the radiative heat transfer equations with random inputs and diffusive scalings. In this problem the random inputs arise due to uncertainties in cross section, initial data or boundary data. We use the generalized polynomial chaos based stochastic Galerkin (gPC-SG) method, which is combined with the micro–macro decomposition based deterministic AP framework in order to handle efficiently the diffusive regime. For linearized problem we prove the regularity of the solution in the random space and consequently the spectral accuracy of the gPC-SG method. We also prove the uniform (inmore » the mean free path) linear stability for the space-time discretizations. Several numerical tests are presented to show the efficiency and accuracy of proposed scheme, especially in the diffusive regime.« less

A study of the radiative transfer equation using a spherical harmonics-nodal collocation method

NASA Astrophysics Data System (ADS)

Capilla, M. T.; Talavera, C. F.; Ginestar, D.; Verdú, G.

2017-03-01

Optical tomography has found many medical applications that need to know how the photons interact with the different tissues. The majority of the photon transport simulations are done using the diffusion approximation, but this approximation has a limited validity when optical properties of the different tissues present large gradients, when structures near the photons source are studied or when anisotropic scattering has to be taken into account. As an alternative to the diffusion model, the PL equations for the radiative transfer problem are studied. These equations are discretized in a rectangular mesh using a nodal collocation method. The performance of this model is studied by solving different 1D and 2D benchmark problems of light propagation in tissue having media with isotropic and anisotropic scattering.

Cross diffusion effect on MHD mixed convection flow of nonlinear radiative heat and mass transfer of Casson fluid over a vertical plate

NASA Astrophysics Data System (ADS)

Ganesh Kumar, K.; Archana, M.; Gireesha, B. J.; Krishanamurthy, M. R.; Rudraswamy, N. G.

2018-03-01

A study on magnetohydrodynamic mixed convection flow of Casson fluid over a vertical plate has been modelled in the presence of Cross diffusion effect and nonlinear thermal radiation. The governing partial differential equations are remodelled into ordinary differential equations by using similarity transformation. The accompanied differential equations are resolved numerically by using Runge-Kutta-Fehlberg forth-fifth order along with shooting method (RKF45 Method). The results of various physical parameters on velocity and temperature profiles are given diagrammatically. The numerical values of the local skin friction coefficient, local Nusselt number and local Sherwood number also are shown in a tabular form. It is found that, effect of Dufour and Soret parameter increases the temperature and concentration component correspondingly.

Theory of transformation thermal convection for creeping flow in porous media: Cloaking, concentrating, and camouflage

NASA Astrophysics Data System (ADS)

Dai, Gaole; Shang, Jin; Huang, Jiping

2018-02-01

Heat can transfer via thermal conduction, thermal radiation, and thermal convection. All the existing theories of transformation thermotics and optics can treat thermal conduction and thermal radiation, respectively. Unfortunately, thermal convection has seldom been touched in transformation theories due to the lack of a suitable theory, thus limiting applications associated with heat transfer through fluids (liquid or gas). Here, we develop a theory of transformation thermal convection by considering the convection-diffusion equation, the equation of continuity, and the Darcy law. By introducing porous media, we get a set of equations keeping their forms under coordinate transformation. As model applications, the theory helps to show the effects of cloaking, concentrating, and camouflage. Our finite-element simulations confirm the theoretical findings. This work offers a transformation theory for thermal convection, thus revealing novel behaviors associated with potential applications; it not only provides different hints on how to control heat transfer by combining thermal conduction, thermal convection, and thermal radiation, but also benefits mass diffusion and other related fields that contain a set of equations and need to transform velocities at the same time.

Effects of Drift-Shell Splitting by Chorus Waves on Radiation Belt Electrons

NASA Astrophysics Data System (ADS)

Chan, A. A.; Zheng, L.; O'Brien, T. P., III; Tu, W.; Cunningham, G.; Elkington, S. R.; Albert, J.

2015-12-01

Drift shell splitting in the radiation belts breaks all three adiabatic invariants of charged particle motion via pitch angle scattering, and produces new diffusion terms that fully populate the diffusion tensor in the Fokker-Planck equation. Based on the stochastic differential equation method, the Radbelt Electron Model (REM) simulation code allows us to solve such a fully three-dimensional Fokker-Planck equation, and to elucidate the sources and transport mechanisms behind the phase space density variations. REM has been used to perform simulations with an empirical initial phase space density followed by a seed electron injection, with a Tsyganenko 1989 magnetic field model, and with chorus wave and ULF wave diffusion models. Our simulation results show that adding drift shell splitting changes the phase space location of the source to smaller L shells, which typically reduces local electron energization (compared to neglecting drift-shell splitting effects). Simulation results with and without drift-shell splitting effects are compared with Van Allen Probe measurements.

An iterative phase-space explicit discontinuous Galerkin method for stellar radiative transfer in extended atmospheres

NASA Astrophysics Data System (ADS)

de Almeida, Valmor F.

2017-07-01

A phase-space discontinuous Galerkin (PSDG) method is presented for the solution of stellar radiative transfer problems. It allows for greater adaptivity than competing methods without sacrificing generality. The method is extensively tested on a spherically symmetric, static, inverse-power-law scattering atmosphere. Results for different sizes of atmospheres and intensities of scattering agreed with asymptotic values. The exponentially decaying behavior of the radiative field in the diffusive-transparent transition region, and the forward peaking behavior at the surface of extended atmospheres were accurately captured. The integrodifferential equation of radiation transfer is solved iteratively by alternating between the radiative pressure equation and the original equation with the integral term treated as an energy density source term. In each iteration, the equations are solved via an explicit, flux-conserving, discontinuous Galerkin method. Finite elements are ordered in wave fronts perpendicular to the characteristic curves so that elemental linear algebraic systems are solved quickly by sweeping the phase space element by element. Two implementations of a diffusive boundary condition at the origin are demonstrated wherein the finite discontinuity in the radiation intensity is accurately captured by the proposed method. This allows for a consistent mechanism to preserve photon luminosity. The method was proved to be robust and fast, and a case is made for the adequacy of parallel processing. In addition to classical two-dimensional plots, results of normalized radiation intensity were mapped onto a log-polar surface exhibiting all distinguishing features of the problem studied.

Stability of stationary inverse transport equation in diffusion scaling

NASA Astrophysics Data System (ADS)

Chen, Ke; Li, Qin; Wang, Li

2018-02-01

We consider the inverse problem of reconstructing the optical parameters for the stationary radiative transfer equation (RTE) from velocity-averaged measurement. The RTE often contains multiple scales, characterized by the magnitude of a dimensionless parameter—the Knudsen number ( \

Polarization radiation in the planetary atmosphere delimited by a heterogeneous diffusely reflecting surface

NASA Technical Reports Server (NTRS)

Strelkov, S. A.; Sushkevich, T. A.

1983-01-01

Spatial frequency characteristics (SFC) and the scattering functions were studied in the two cases of a uniform horizontal layer with absolutely black bottom, and an isolated layer. The mathematical model for these examples describes the horizontal heterogeneities in a light field with regard to radiation polarization in a three dimensional planar atmosphere, delimited by a heterogeneous surface with diffuse reflection. The perturbation method was used to obtain vector transfer equations which correspond to the linear and nonlinear systems of polarization radiation transfer. The boundary value tasks for the vector transfer equation that is a parametric set and one dimensional are satisfied by the SFC of the nonlinear system, and are expressed through the SFC of linear approximation. As a consequence of the developed theory, formulas were obtained for analytical calculation of albedo in solving the task of dissemination of polarization radiation in the planetary atmosphere with uniform Lambert bottom.

Discrete ordinates solutions of nongray radiative transfer with diffusely reflecting walls

NASA Technical Reports Server (NTRS)

Menart, J. A.; Lee, Haeok S.; Kim, Tae-Kuk

1993-01-01

Nongray gas radiation in a plane parallel slab bounded by gray, diffusely reflecting walls is studied using the discrete ordinates method. The spectral equation of transfer is averaged over a narrow wavenumber interval preserving the spectral correlation effect. The governing equations are derived by considering the history of multiple reflections between two reflecting wails. A closure approximation is applied so that only a finite number of reflections have to be explicitly included. The closure solutions express the physics of the problem to a very high degree and show relatively little error. Numerical solutions are obtained by applying a statistical narrow-band model for gas properties and a discrete ordinates code. The net radiative wail heat fluxes and the radiative source distributions are obtained for different temperature profiles. A zeroth-degree formulation, where no wall reflection is handled explicitly, is sufficient to predict the radiative transfer accurately for most cases considered, when compared with increasingly accurate solutions based on explicitly tracing a larger number of wail reflections without any closure approximation applied.

Using hybrid implicit Monte Carlo diffusion to simulate gray radiation hydrodynamics

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

Cleveland, Mathew A., E-mail: cleveland7@llnl.gov; Gentile, Nick

This work describes how to couple a hybrid Implicit Monte Carlo Diffusion (HIMCD) method with a Lagrangian hydrodynamics code to evaluate the coupled radiation hydrodynamics equations. This HIMCD method dynamically applies Implicit Monte Carlo Diffusion (IMD) [1] to regions of a problem that are opaque and diffusive while applying standard Implicit Monte Carlo (IMC) [2] to regions where the diffusion approximation is invalid. We show that this method significantly improves the computational efficiency as compared to a standard IMC/Hydrodynamics solver, when optically thick diffusive material is present, while maintaining accuracy. Two test cases are used to demonstrate the accuracy andmore » performance of HIMCD as compared to IMC and IMD. The first is the Lowrie semi-analytic diffusive shock [3]. The second is a simple test case where the source radiation streams through optically thin material and heats a thick diffusive region of material causing it to rapidly expand. We found that HIMCD proves to be accurate, robust, and computationally efficient for these test problems.« less

Diffusion approximation of the radiative-conductive heat transfer model with Fresnel matching conditions

NASA Astrophysics Data System (ADS)

Chebotarev, Alexander Yu.; Grenkin, Gleb V.; Kovtanyuk, Andrey E.; Botkin, Nikolai D.; Hoffmann, Karl-Heinz

2018-04-01

The paper is concerned with a problem of diffraction type. The study starts with equations of complex (radiative and conductive) heat transfer in a multicomponent domain with Fresnel matching conditions at the interfaces. Applying the diffusion, P1, approximation yields a pair of coupled nonlinear PDEs describing the radiation intensity and temperature for each component of the domain. Matching conditions for these PDEs, imposed at the interfaces between the domain components, are derived. The unique solvability of the obtained problem is proven, and numerical experiments are conducted.

An AMR capable finite element diffusion solver for ALE hydrocodes [An AMR capable diffusion solver for ALE-AMR

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

Fisher, A. C.; Bailey, D. S.; Kaiser, T. B.

2015-02-01

Here, we present a novel method for the solution of the diffusion equation on a composite AMR mesh. This approach is suitable for including diffusion based physics modules to hydrocodes that support ALE and AMR capabilities. To illustrate, we proffer our implementations of diffusion based radiation transport and heat conduction in a hydrocode called ALE-AMR. Numerical experiments conducted with the diffusion solver and associated physics packages yield 2nd order convergence in the L 2 norm.

A Radiation Chemistry Code Based on the Greens Functions of the Diffusion Equation

NASA Technical Reports Server (NTRS)

Plante, Ianik; Wu, Honglu

2014-01-01

Ionizing radiation produces several radiolytic species such as.OH, e-aq, and H. when interacting with biological matter. Following their creation, radiolytic species diffuse and chemically react with biological molecules such as DNA. Despite years of research, many questions on the DNA damage by ionizing radiation remains, notably on the indirect effect, i.e. the damage resulting from the reactions of the radiolytic species with DNA. To simulate DNA damage by ionizing radiation, we are developing a step-by-step radiation chemistry code that is based on the Green's functions of the diffusion equation (GFDE), which is able to follow the trajectories of all particles and their reactions with time. In the recent years, simulations based on the GFDE have been used extensively in biochemistry, notably to simulate biochemical networks in time and space and are often used as the "gold standard" to validate diffusion-reaction theories. The exact GFDE for partially diffusion-controlled reactions is difficult to use because of its complex form. Therefore, the radial Green's function, which is much simpler, is often used. Hence, much effort has been devoted to the sampling of the radial Green's functions, for which we have developed a sampling algorithm This algorithm only yields the inter-particle distance vector length after a time step; the sampling of the deviation angle of the inter-particle vector is not taken into consideration. In this work, we show that the radial distribution is predicted by the exact radial Green's function. We also use a technique developed by Clifford et al. to generate the inter-particle vector deviation angles, knowing the inter-particle vector length before and after a time step. The results are compared with those predicted by the exact GFDE and by the analytical angular functions for free diffusion. This first step in the creation of the radiation chemistry code should help the understanding of the contribution of the indirect effect in the formation of DNA damage and double-strand breaks.

A Fast Vector Radiative Transfer Model for Atmospheric and Oceanic Remote Sensing

NASA Astrophysics Data System (ADS)

Ding, J.; Yang, P.; King, M. D.; Platnick, S. E.; Meyer, K.

2017-12-01

A fast vector radiative transfer model is developed in support of atmospheric and oceanic remote sensing. This model is capable of simulating the Stokes vector observed at the top of the atmosphere (TOA) and the terrestrial surface by considering absorption, scattering, and emission. The gas absorption is parameterized in terms of atmospheric gas concentrations, temperature, and pressure. The parameterization scheme combines a regression method and the correlated-K distribution method, and can easily integrate with multiple scattering computations. The approach is more than four orders of magnitude faster than a line-by-line radiative transfer model with errors less than 0.5% in terms of transmissivity. A two-component approach is utilized to solve the vector radiative transfer equation (VRTE). The VRTE solver separates the phase matrices of aerosol and cloud into forward and diffuse parts and thus the solution is also separated. The forward solution can be expressed by a semi-analytical equation based on the small-angle approximation, and serves as the source of the diffuse part. The diffuse part is solved by the adding-doubling method. The adding-doubling implementation is computationally efficient because the diffuse component needs much fewer spherical function expansion terms. The simulated Stokes vector at both the TOA and the surface have comparable accuracy compared with the counterparts based on numerically rigorous methods.

Modelling the radiotherapy effect in the reaction-diffusion equation.

PubMed

Borasi, Giovanni; Nahum, Alan

2016-09-01

In recent years, the reaction-diffusion (Fisher-Kolmogorov) equation has received much attention from the oncology research community due to its ability to describe the infiltrating nature of glioblastoma multiforme and its extraordinary resistance to any type of therapy. However, in a number of previous papers in the literature on applications of this equation, the term (R) expressing the 'External Radiotherapy effect' was incorrectly derived. In this note we derive an analytical expression for this term in the correct form to be included in the reaction-diffusion equation. The R term has been derived starting from the Linear-Quadratic theory of cell killing by ionizing radiation. The correct definition of R was adopted and the basic principles of differential calculus applied. The compatibility of the R term derived here with the reaction-diffusion equation was demonstrated. Referring to a typical glioblastoma tumour, we have compared the results obtained using our expression for the R term with the 'incorrect' expression proposed by other authors. Copyright © 2016 Associazione Italiana di Fisica Medica. Published by Elsevier Ltd. All rights reserved.

Relationships between diffuse reflectance and vegetation canopy variables based on the radiative transfer theory

NASA Technical Reports Server (NTRS)

Park, J. K.; Deering, D. W.

1981-01-01

Out of the lengthy original expression of the diffuse reflectance formula, simple working equations were derived by employing characteristic parameters, which are independent of the canopy coverage and identifiable by field observations. The typical asymptotic nature of reflectance data that is usually observed in biomass studies was clearly explained. The usefulness of the simplified equations was demonstrated by the exceptionally close fit of the theoretical curves to two separately acquired data sets for alfalfa and shortgrass prairie canopies.

Thermal radiation and mass transfer effects on unsteady MHD free convection flow past a vertical oscillating plate

NASA Astrophysics Data System (ADS)

Rana, B. M. Jewel; Ahmed, Rubel; Ahmmed, S. F.

2017-06-01

Unsteady MHD free convection flow past a vertical porous plate in porous medium with radiation, diffusion thermo, thermal diffusion and heat source are analyzed. The governing non-linear, partial differential equations are transformed into dimensionless by using non-dimensional quantities. Then the resultant dimensionless equations are solved numerically by applying an efficient, accurate and conditionally stable finite difference scheme of explicit type with the help of a computer programming language Compaq Visual Fortran. The stability and convergence analysis has been carried out to establish the effect of velocity, temperature, concentration, skin friction, Nusselt number, Sherwood number, stream lines and isotherms line. Finally, the effects of various parameters are presented graphically and discussed qualitatively.

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

Liemert, André, E-mail: andre.liemert@ilm.uni-ulm.de; Kienle, Alwin

Purpose: Explicit solutions of the monoenergetic radiative transport equation in the P{sub 3} approximation have been derived which can be evaluated with nearly the same computational effort as needed for solving the standard diffusion equation (DE). In detail, the authors considered the important case of a semi-infinite medium which is illuminated by a collimated beam of light. Methods: A combination of the classic spherical harmonics method and the recently developed method of rotated reference frames is used for solving the P{sub 3} equations in closed form. Results: The derived solutions are illustrated and compared to exact solutions of the radiativemore » transport equation obtained via the Monte Carlo (MC) method as well as with other approximated analytical solutions. It is shown that for the considered cases which are relevant for biomedical optics applications, the P{sub 3} approximation is close to the exact solution of the radiative transport equation. Conclusions: The authors derived exact analytical solutions of the P{sub 3} equations under consideration of boundary conditions for defining a semi-infinite medium. The good agreement to Monte Carlo simulations in the investigated domains, for example, in the steady-state and time domains, as well as the short evaluation time needed suggests that the derived equations can replace the often applied solutions of the diffusion equation for the homogeneous semi-infinite medium.« less

Refractive Index Effects on Radiation in an Absorbing, Emitting, and Scattering Laminated Layer

NASA Technical Reports Server (NTRS)

Siegel, R.; Spuckler, C. M.

1993-01-01

A simple set of equations is derived for predicting temperature radiative energy flow in a two-region semitransparent laminated layer in the limit of zero heat conduction. The composite is heated on its two sides by unequal amounts of incident radiation. The two layers of the composite have different refractive indices, and each material absorbs, emits, and isotropically scatters radiation. The interfaces are diffuse, and all interface reflections are included. To illustrate the thermal behavior that is readily calculated from the equations, typical results an given for various optical thicknesses and refractive indices of the layers. Internal reflections have a substantial effect on the temperature distribution and radiative heat flow.

Two-Flux Green's Function Analysis for Transient Spectral Radiation in a Composite

NASA Technical Reports Server (NTRS)

Siegel, Robert

1996-01-01

An analysis is developed for obtaining transient temperatures in a two-layer semitransparent composite with spectrally dependent properties. Each external boundary of the composite is subjected to radiation and convection. The two-flux radiative transfer equations are solved by deriving a Green's function. This yields the local radiative heat source needed to numerically solve the transient energy equation. An advantage of the two-flux method is that isotropic scattering is included without added complexity. The layer refractive indices are larger than one. This produces internal reflections at the boundaries and the internal interface; the reflections are assumed diffuse. Spectral results using the Green's function method are verified by comparing with numerical solutions using the exact radiative transfer equations. Transient temperature distributions are given to illustrate the effect of radiative heating on one side of a composite with external convective cooling. The protection of a material from incident radiation is illustrated by adding scattering to the layer adjacent to the radiative source.

A second order radiative transfer equation and its solution by meshless method with application to strongly inhomogeneous media

NASA Astrophysics Data System (ADS)

Zhao, J. M.; Tan, J. Y.; Liu, L. H.

2013-01-01

A new second order form of radiative transfer equation (named MSORTE) is proposed, which overcomes the singularity problem of a previously proposed second order radiative transfer equation [J.E. Morel, B.T. Adams, T. Noh, J.M. McGhee, T.M. Evans, T.J. Urbatsch, Spatial discretizations for self-adjoint forms of the radiative transfer equations, J. Comput. Phys. 214 (1) (2006) 12-40 (where it was termed SAAI), J.M. Zhao, L.H. Liu, Second order radiative transfer equation and its properties of numerical solution using finite element method, Numer. Heat Transfer B 51 (2007) 391-409] in dealing with inhomogeneous media where some locations have very small/zero extinction coefficient. The MSORTE contains a naturally introduced diffusion (or second order) term which provides better numerical property than the classic first order radiative transfer equation (RTE). The stability and convergence characteristics of the MSORTE discretized by central difference scheme is analyzed theoretically, and the better numerical stability of the second order form radiative transfer equations than the RTE when discretized by the central difference type method is proved. A collocation meshless method is developed based on the MSORTE to solve radiative transfer in inhomogeneous media. Several critical test cases are taken to verify the performance of the presented method. The collocation meshless method based on the MSORTE is demonstrated to be capable of stably and accurately solve radiative transfer in strongly inhomogeneous media, media with void region and even with discontinuous extinction coefficient.

Application of linear multifrequency-grey acceleration to preconditioned Krylov iterations for thermal radiation transport

DOE PAGES

Till, Andrew T.; Warsa, James S.; Morel, Jim E.

2018-06-15

The thermal radiative transfer (TRT) equations comprise a radiation equation coupled to the material internal energy equation. Linearization of these equations produces effective, thermally-redistributed scattering through absorption-reemission. In this paper, we investigate the effectiveness and efficiency of Linear-Multi-Frequency-Grey (LMFG) acceleration that has been reformulated for use as a preconditioner to Krylov iterative solution methods. We introduce two general frameworks, the scalar flux formulation (SFF) and the absorption rate formulation (ARF), and investigate their iterative properties in the absence and presence of true scattering. SFF has a group-dependent state size but may be formulated without inner iterations in the presence ofmore » scattering, while ARF has a group-independent state size but requires inner iterations when scattering is present. We compare and evaluate the computational cost and efficiency of LMFG applied to these two formulations using a direct solver for the preconditioners. Finally, this work is novel because the use of LMFG for the radiation transport equation, in conjunction with Krylov methods, involves special considerations not required for radiation diffusion.« less

Heuristic Green's function of the time dependent radiative transfer equation for a semi-infinite medium.

PubMed

Martelli, Fabrizio; Sassaroli, Angelo; Pifferi, Antonio; Torricelli, Alessandro; Spinelli, Lorenzo; Zaccanti, Giovanni

2007-12-24

The Green's function of the time dependent radiative transfer equation for the semi-infinite medium is derived for the first time by a heuristic approach based on the extrapolated boundary condition and on an almost exact solution for the infinite medium. Monte Carlo simulations performed both in the simple case of isotropic scattering and of an isotropic point-like source, and in the more realistic case of anisotropic scattering and pencil beam source, are used to validate the heuristic Green's function. Except for the very early times, the proposed solution has an excellent accuracy (> 98 % for the isotropic case, and > 97 % for the anisotropic case) significantly better than the diffusion equation. The use of this solution could be extremely useful in the biomedical optics field where it can be directly employed in conditions where the use of the diffusion equation is limited, e.g. small volume samples, high absorption and/or low scattering media, short source-receiver distances and early times. Also it represents a first step to derive tools for other geometries (e.g. slab and slab with inhomogeneities inside) of practical interest for noninvasive spectroscopy and diffuse optical imaging. Moreover the proposed solution can be useful to several research fields where the study of a transport process is fundamental.

Steady state model for the thermal regimes of shells of airships and hot air balloons

NASA Astrophysics Data System (ADS)

Luchev, Oleg A.

1992-10-01

A steady state model of the temperature regime of airships and hot air balloons shells is developed. The model includes three governing equations: the equation of the temperature field of airships or balloons shell, the integral equation for the radiative fluxes on the internal surface of the shell, and the integral equation for the natural convective heat exchange between the shell and the internal gas. In the model the following radiative fluxes on the shell external surface are considered: the direct and the earth reflected solar radiation, the diffuse solar radiation, the infrared radiation of the earth surface and that of the atmosphere. For the calculations of the infrared external radiation the model of the plane layer of the atmosphere is used. The convective heat transfer on the external surface of the shell is considered for the cases of the forced and the natural convection. To solve the mentioned set of the equations the numerical iterative procedure is developed. The model and the numerical procedure are used for the simulation study of the temperature fields of an airship shell under the forced and the natural convective heat transfer.

Particle Acceleration and Radiative Losses at Relativistic Shocks

NASA Astrophysics Data System (ADS)

Dempsey, P.; Duffy, P.

A semi-analytic approach to the relativistic transport equation with isotropic diffusion and consistent radiative losses is presented. It is based on the eigenvalue method first introduced in Kirk & Schneider [5]and Heavens & Drury [3]. We demonstrate the pitch-angle dependence of the cut-off in relativistic shocks.

Approximate Solution Methods for Spectral Radiative Transfer in High Refractive Index Layers

NASA Technical Reports Server (NTRS)

Siegel, R.; Spuckler, C. M.

1994-01-01

Some ceramic materials for high temperature applications are partially transparent for radiative transfer. The refractive indices of these materials can be substantially greater than one which influences internal radiative emission and reflections. Heat transfer behavior of single and laminated layers has been obtained in the literature by numerical solutions of the radiative transfer equations coupled with heat conduction and heating at the boundaries by convection and radiation. Two-flux and diffusion methods are investigated here to obtain approximate solutions using a simpler formulation than required for exact numerical solutions. Isotropic scattering is included. The two-flux method for a single layer yields excellent results for gray and two band spectral calculations. The diffusion method yields a good approximation for spectral behavior in laminated multiple layers if the overall optical thickness is larger than about ten. A hybrid spectral model is developed using the two-flux method in the optically thin bands, and radiative diffusion in bands that are optically thick.

SAS-2 observations of the diffuse gamma radiation in the galactic latitude interval 10 deg absolute b or equal to 90 deg

NASA Technical Reports Server (NTRS)

Fichtel, C. E.; Hartman, R. C.; Kniffen, D. A.; Thompson, D. J.; Oegelman, H. B.; Oezel, M. E.; Tuemer, T.

1977-01-01

An analysis of all of the second small astronomy satellite gamma-ray data for galactic latitudes with the absolute value of b 10 deg has shown that the intensity varies with galactic latitude, being larger near 10 deg than 90 deg. For energies above 100 MeV the gamma-ray data are consistent with a latitude distribution of the form I(b) = C sub 1 + C sub 2/sin b, with the second term being dominant. This result suggests that the radiation above 100 MeV is coming largely from local regions of the galactic disk. Between 35 and 100 MeV, a similar equation is also a good representation of the data, but here the two terms are comparable. These results indicate that the diffuse radiation above 35 MeV consists of two parts, one with a relatively hard galactic component and the other an isotropic, steep spectral component which extrapolates back well to the low energy diffuse radiation. The steepness of the diffuse isotropic component places significant constraints on possible theoretical models of this radiation.

Photon migration through a turbid slab described by a model based on diffusion approximation. I. Theory.

PubMed

Contini, D; Martelli, F; Zaccanti, G

1997-07-01

The diffusion approximation of the radiative transfer equation is a model used widely to describe photon migration in highly diffusing media and is an important matter in biological tissue optics. An analysis of the time-dependent diffusion equation together with its solutions for the slab geometry and for a semi-infinite diffusing medium are reported. These solutions, presented for both the time-dependent and the continuous wave source, account for the refractive index mismatch between the turbid medium and the surrounding medium. The results have been compared with those obtained when different boundary conditions were assumed. The comparison has shown that the effect of the refractive index mismatch cannot be disregarded. This effect is particularly important for the transmittance. The discussion of results also provides an analysis of the role of the absorption coefficient in the expression of the diffusion coefficient.

Cattaneo-Christov heat flux effect on hydromagnetic radiative Oldroyd-B liquid flow across a cone/wedge in the presence of cross-diffusion

NASA Astrophysics Data System (ADS)

Gnaneswara Reddy, M.

2018-01-01

The present article scrutinizes the prominent characteristics of the Cattaneo-Christov heat flux on magnetohydrodynamic Oldroyd-B radiative liquid flow over two different geometries. The effects of cross-diffusion are considered in the modeling of species and energy equations. Similarity transformations are employed to transmute the governing flow, species and energy equations into a set of nonlinear ordinary differential equations (ODEs) with the appropriate boundary conditions. The final system of dimensionless equations is resolved numerically by utilizing the R-K-Fehlberg numerical approach. The behaviors of all physical pertinent flow controlling variables on the three flow distributions are analyzed through plots. The obtained numerical results have been compared with earlier published work and reveal good agreement. The Deborah numbers γ1 and γ2 have quite opposite effects on velocity and energy fields. The increase in thermal relaxation parameter β corresponds to a decrease in the fluid temperature. This study has salient applications in heat and mass transfer manufacturing system processing for energy conversion.

Three-temperature plasma shock solutions with gray radiation diffusion

DOE PAGES

Johnson, Bryan M.; Klein, Richard I.

2016-04-19

Here we discuss the effects of radiation on the structure of shocks in a fully ionized plasma are investigated by solving the steady-state fluid equations for ions, electrons, and radiation. The electrons and ions are assumed to have the same bulk velocity but separate temperatures, and the radiation is modeled with the gray diffusion approximation. Both electron and ion conduction are included, as well as ion viscosity. When the material is optically thin, three-temperature behavior occurs. When the diffusive flux of radiation is important but radiation pressure is not, two-temperature behavior occurs, with the electrons strongly coupled to the radiation.more » Since the radiation heats the electrons on length scales that are much longer than the electron–ion Coulomb coupling length scale, these solutions resemble radiative shock solutions rather than plasma shock solutions that neglect radiation. When radiation pressure is important, all three components are strongly coupled. Results with constant values for the transport and coupling coefficients are compared to a full numerical simulation with a good match between the two, demonstrating that steady shock solutions constitute a straightforward and comprehensive verification test methodology for multi-physics numerical algorithms.« less

Three-temperature plasma shock solutions with gray radiation diffusion

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

Johnson, Bryan M.; Klein, Richard I.

Here we discuss the effects of radiation on the structure of shocks in a fully ionized plasma are investigated by solving the steady-state fluid equations for ions, electrons, and radiation. The electrons and ions are assumed to have the same bulk velocity but separate temperatures, and the radiation is modeled with the gray diffusion approximation. Both electron and ion conduction are included, as well as ion viscosity. When the material is optically thin, three-temperature behavior occurs. When the diffusive flux of radiation is important but radiation pressure is not, two-temperature behavior occurs, with the electrons strongly coupled to the radiation.more » Since the radiation heats the electrons on length scales that are much longer than the electron–ion Coulomb coupling length scale, these solutions resemble radiative shock solutions rather than plasma shock solutions that neglect radiation. When radiation pressure is important, all three components are strongly coupled. Results with constant values for the transport and coupling coefficients are compared to a full numerical simulation with a good match between the two, demonstrating that steady shock solutions constitute a straightforward and comprehensive verification test methodology for multi-physics numerical algorithms.« less

Quantitative computational infrared imaging of buoyant diffusion flames

NASA Astrophysics Data System (ADS)

Newale, Ashish S.

Studies of infrared radiation from turbulent buoyant diffusion flames impinging on structural elements have applications to the development of fire models. A numerical and experimental study of radiation from buoyant diffusion flames with and without impingement on a flat plate is reported. Quantitative images of the radiation intensity from the flames are acquired using a high speed infrared camera. Large eddy simulations are performed using fire dynamics simulator (FDS version 6). The species concentrations and temperature from the simulations are used in conjunction with a narrow-band radiation model (RADCAL) to solve the radiative transfer equation. The computed infrared radiation intensities rendered in the form of images and compared with the measurements. The measured and computed radiation intensities reveal necking and bulging with a characteristic frequency of 7.1 Hz which is in agreement with previous empirical correlations. The results demonstrate the effects of stagnation point boundary layer on the upstream buoyant shear layer. The coupling between these two shear layers presents a model problem for sub-grid scale modeling necessary for future large eddy simulations.

An iterative phase-space explicit discontinuous Galerkin method for stellar radiative transfer in extended atmospheres

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

de Almeida, Valmor F.

In this work, a phase-space discontinuous Galerkin (PSDG) method is presented for the solution of stellar radiative transfer problems. It allows for greater adaptivity than competing methods without sacrificing generality. The method is extensively tested on a spherically symmetric, static, inverse-power-law scattering atmosphere. Results for different sizes of atmospheres and intensities of scattering agreed with asymptotic values. The exponentially decaying behavior of the radiative field in the diffusive-transparent transition region, and the forward peaking behavior at the surface of extended atmospheres were accurately captured. The integrodifferential equation of radiation transfer is solved iteratively by alternating between the radiative pressure equationmore » and the original equation with the integral term treated as an energy density source term. In each iteration, the equations are solved via an explicit, flux-conserving, discontinuous Galerkin method. Finite elements are ordered in wave fronts perpendicular to the characteristic curves so that elemental linear algebraic systems are solved quickly by sweeping the phase space element by element. Two implementations of a diffusive boundary condition at the origin are demonstrated wherein the finite discontinuity in the radiation intensity is accurately captured by the proposed method. This allows for a consistent mechanism to preserve photon luminosity. The method was proved to be robust and fast, and a case is made for the adequacy of parallel processing. In addition to classical two-dimensional plots, results of normalized radiation intensity were mapped onto a log-polar surface exhibiting all distinguishing features of the problem studied.« less

An iterative phase-space explicit discontinuous Galerkin method for stellar radiative transfer in extended atmospheres

DOE PAGES

de Almeida, Valmor F.

2017-04-19

In this work, a phase-space discontinuous Galerkin (PSDG) method is presented for the solution of stellar radiative transfer problems. It allows for greater adaptivity than competing methods without sacrificing generality. The method is extensively tested on a spherically symmetric, static, inverse-power-law scattering atmosphere. Results for different sizes of atmospheres and intensities of scattering agreed with asymptotic values. The exponentially decaying behavior of the radiative field in the diffusive-transparent transition region, and the forward peaking behavior at the surface of extended atmospheres were accurately captured. The integrodifferential equation of radiation transfer is solved iteratively by alternating between the radiative pressure equationmore » and the original equation with the integral term treated as an energy density source term. In each iteration, the equations are solved via an explicit, flux-conserving, discontinuous Galerkin method. Finite elements are ordered in wave fronts perpendicular to the characteristic curves so that elemental linear algebraic systems are solved quickly by sweeping the phase space element by element. Two implementations of a diffusive boundary condition at the origin are demonstrated wherein the finite discontinuity in the radiation intensity is accurately captured by the proposed method. This allows for a consistent mechanism to preserve photon luminosity. The method was proved to be robust and fast, and a case is made for the adequacy of parallel processing. In addition to classical two-dimensional plots, results of normalized radiation intensity were mapped onto a log-polar surface exhibiting all distinguishing features of the problem studied.« less

Modelling thermal radiation in buoyant turbulent diffusion flames

NASA Astrophysics Data System (ADS)

Consalvi, J. L.; Demarco, R.; Fuentes, A.

2012-10-01

This work focuses on the numerical modelling of radiative heat transfer in laboratory-scale buoyant turbulent diffusion flames. Spectral gas and soot radiation is modelled by using the Full-Spectrum Correlated-k (FSCK) method. Turbulence-Radiation Interactions (TRI) are taken into account by considering the Optically-Thin Fluctuation Approximation (OTFA), the resulting time-averaged Radiative Transfer Equation (RTE) being solved by the Finite Volume Method (FVM). Emission TRIs and the mean absorption coefficient are then closed by using a presumed probability density function (pdf) of the mixture fraction. The mean gas flow field is modelled by the Favre-averaged Navier-Stokes (FANS) equation set closed by a buoyancy-modified k-ɛ model with algebraic stress/flux models (ASM/AFM), the Steady Laminar Flamelet (SLF) model coupled with a presumed pdf approach to account for Turbulence-Chemistry Interactions, and an acetylene-based semi-empirical two-equation soot model. Two sets of experimental pool fire data are used for validation: propane pool fires 0.3 m in diameter with Heat Release Rates (HRR) of 15, 22 and 37 kW and methane pool fires 0.38 m in diameter with HRRs of 34 and 176 kW. Predicted flame structures, radiant fractions, and radiative heat fluxes on surrounding surfaces are found in satisfactory agreement with available experimental data across all the flames. In addition further computations indicate that, for the present flames, the gray approximation can be applied for soot with a minor influence on the results, resulting in a substantial gain in Computer Processing Unit (CPU) time when the FSCK is used to treat gas radiation.

Characterization of supersonic radiation diffusion waves

DOE PAGES

Moore, Alastair S.; Guymer, Thomas M.; Morton, John; ...

2015-02-27

Supersonic and diffusive radiation flow is an important test problem for the radiative transfer models used in radiation-hydrodynamics computer codes owing to solutions being accessible via analytic and numeric methods. We present experimental results with which we compare these solutions by studying supersonic and diffusive flow in the laboratory. Here, we present results of higher-accuracy experiments than previously possible studying radiation flow through up to 7 high-temperature mean free paths of low-density, chlorine-doped polystyrene foam and silicon dioxide aerogel contained by an Au tube. Measurements of the heat front position and absolute measurements of the x-ray emission arrival at themore » end of the tube are used to test numerical and analytical models. We find excellent absolute agreement with simulations provided that the opacity and the equation of state are adjusted within expected uncertainties; analytical models provide a good phenomenological match to measurements but are not in quantitative agreement due to their limited scope.« less

Characterization of supersonic radiation diffusion waves

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

Moore, Alastair S.; Guymer, Thomas M.; Morton, John

Supersonic and diffusive radiation flow is an important test problem for the radiative transfer models used in radiation-hydrodynamics computer codes owing to solutions being accessible via analytic and numeric methods. We present experimental results with which we compare these solutions by studying supersonic and diffusive flow in the laboratory. Here, we present results of higher-accuracy experiments than previously possible studying radiation flow through up to 7 high-temperature mean free paths of low-density, chlorine-doped polystyrene foam and silicon dioxide aerogel contained by an Au tube. Measurements of the heat front position and absolute measurements of the x-ray emission arrival at themore » end of the tube are used to test numerical and analytical models. We find excellent absolute agreement with simulations provided that the opacity and the equation of state are adjusted within expected uncertainties; analytical models provide a good phenomenological match to measurements but are not in quantitative agreement due to their limited scope.« less

Evaluation of the telegrapher's equation and multiple-flux theories for calculating the transmittance and reflectance of a diffuse absorbing slab.

PubMed

Kong, Steven H; Shore, Joel D

2007-03-01

We study the propagation of light through a medium containing isotropic scattering and absorption centers. With a Monte Carlo simulation serving as the benchmark solution to the radiative transfer problem of light propagating through a turbid slab, we compare the transmission and reflection density computed from the telegrapher's equation, the diffusion equation, and multiple-flux theories such as the Kubelka-Munk and four-flux theories. Results are presented for both normally incident light and diffusely incident light. We find that we can always obtain very good results from the telegrapher's equation provided that two parameters that appear in the solution are set appropriately. We also find an interesting connection between certain solutions of the telegrapher's equation and solutions of the Kubelka-Munk and four-flux theories with a small modification to how the phenomenological parameters in those theories are traditionally related to the optical scattering and absorption coefficients of the slab. Finally, we briefly explore how well the theories can be extended to the case of anisotropic scattering by multiplying the scattering coefficient by a simple correction factor.

A multidimensional unified gas-kinetic scheme for radiative transfer equations on unstructured mesh

NASA Astrophysics Data System (ADS)

Sun, Wenjun; Jiang, Song; Xu, Kun

2017-12-01

In order to extend the unified gas kinetic scheme (UGKS) to solve radiative transfer equations in a complex geometry, a multidimensional asymptotic preserving implicit method on unstructured mesh is constructed in this paper. With an implicit formulation, the CFL condition for the determination of the time step in UGKS can be much relaxed, and a large time step is used in simulations. Differently from previous direction-by-direction UGKS on orthogonal structured mesh, on unstructured mesh the interface flux transport takes into account multi-dimensional effect, where gradients of radiation intensity and material temperature in both normal and tangential directions of a cell interface are included in the flux evaluation. The multiple scale nature makes the UGKS be able to capture the solutions in both optically thin and thick regions seamlessly. In the optically thick region the condition of cell size being less than photon's mean free path is fully removed, and the UGKS recovers a solver for diffusion equation in such a limit on unstructured mesh. For a distorted quadrilateral mesh, the UGKS goes to a nine-point scheme for the diffusion equation, and it naturally reduces to the standard five-point scheme for a orthogonal quadrilateral mesh. Numerical computations covering a wide range of transport regimes on unstructured and distorted quadrilateral meshes will be presented to validate the current approach.

SAS-2 observations of the diffuse gamma radiation in the galactic latitude interval from 10 to 90 deg in both hemispheres

NASA Technical Reports Server (NTRS)

Fichtel, C. E.; Hartman, R. C.; Kniffen, D. A.; Thompson, D. J.; Ogelman, H. B.; Ozel, M. E.; Tumer, T.

1977-01-01

An analysis of all the second Small Astronomy Satellite (SAS-2) gamma-ray data for galactic latitudes higher than 10 deg in both hemispheres has shown that the intensity varies with galactic latitude, being larger near 10 deg than 90 deg. For energies above 100 MeV the gamma-ray data are consistent with a latitude distribution of the form I(b) = C1 + C2/sin b, with the second term being dominant. This result suggests that the radiation above 100 MeV is coming largely from local regions of the galactic disk. Between 35 and 100 MeV, a similar equation is also a good representation of the data, but here the two terms are comparable. These results indicate that the diffuse radiation above 35 MeV consists of two parts, one with a relatively hard galactic component and the other an isotropic steep spectral component which extrapolates back well to the low-energy (less than 10 MeV) diffuse radiation. The steepness of the diffuse isotropic component places significant constraints on possible theoretical models of this radiation.

Recursive recovery of Markov transition probabilities from boundary value data

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

Patch, Sarah Kathyrn

1994-04-01

In an effort to mathematically describe the anisotropic diffusion of infrared radiation in biological tissue Gruenbaum posed an anisotropic diffusion boundary value problem in 1989. In order to accommodate anisotropy, he discretized the temporal as well as the spatial domain. The probabilistic interpretation of the diffusion equation is retained; radiation is assumed to travel according to a random walk (of sorts). In this random walk the probabilities with which photons change direction depend upon their previous as well as present location. The forward problem gives boundary value data as a function of the Markov transition probabilities. The inverse problem requiresmore » finding the transition probabilities from boundary value data. Problems in the plane are studied carefully in this thesis. Consistency conditions amongst the data are derived. These conditions have two effects: they prohibit inversion of the forward map but permit smoothing of noisy data. Next, a recursive algorithm which yields a family of solutions to the inverse problem is detailed. This algorithm takes advantage of all independent data and generates a system of highly nonlinear algebraic equations. Pluecker-Grassmann relations are instrumental in simplifying the equations. The algorithm is used to solve the 4 x 4 problem. Finally, the smallest nontrivial problem in three dimensions, the 2 x 2 x 2 problem, is solved.« less

Simulation of a fast diffuse optical tomography system based on radiative transfer equation

NASA Astrophysics Data System (ADS)

Motevalli, S. M.; Payani, A.

2016-12-01

Studies show that near-infrared (NIR) light (light with wavelength between 700nm and 1300nm) undergoes two interactions, absorption and scattering, when it penetrates a tissue. Since scattering is the predominant interaction, the calculation of light distribution in the tissue and the image reconstruction of absorption and scattering coefficients are very complicated. Some analytical and numerical methods, such as radiative transport equation and Monte Carlo method, have been used for the simulation of light penetration in tissue. Recently, some investigators in the world have tried to develop a diffuse optical tomography system. In these systems, NIR light penetrates the tissue and passes through the tissue. Then, light exiting the tissue is measured by NIR detectors placed around the tissue. These data are collected from all the detectors and transferred to the computational parts (including hardware and software), which make a cross-sectional image of the tissue after performing some computational processes. In this paper, the results of the simulation of an optical diffuse tomography system are presented. This simulation involves two stages: a) Simulation of the forward problem (or light penetration in the tissue), which is performed by solving the diffusion approximation equation in the stationary state using FEM. b) Simulation of the inverse problem (or image reconstruction), which is performed by the optimization algorithm called Broyden quasi-Newton. This method of image reconstruction is faster compared to the other Newton-based optimization algorithms, such as the Levenberg-Marquardt one.

Synchronous flowering of the rubber tree (Hevea brasiliensis) induced by high solar radiation intensity.

PubMed

Yeang, Hoong-Yeet

2007-01-01

How tropical trees flower synchronously near the equator in the absence of significant day length variation or other meteorological cues has long been a puzzle. The rubber tree (Hevea brasiliensis) is used as a model to investigate this phenomenon. The annual cycle of solar radiation intensity is shown to correspond closely with the flowering of the rubber tree planted near the equator and in the subtropics. Unlike in temperate regions, where incoming solar radiation (insolation) is dependent on both day length and radiation intensity, insolation at the equator is due entirely to the latter. Insolation at the upper atmosphere peaks twice a year during the spring and autumn equinoxes, but the actual solar radiation that reaches the ground is attenuated to varying extents in different localities. The rubber tree shows one or two flowering seasons a year (with major and minor seasons in the latter) in accordance with the solar radiation intensity received. High solar radiation intensity, and in particular bright sunshine (as distinct from prolonged diffuse radiation), induces synchronous anthesis and blooming in Hevea around the time of the equinoxes. The same mechanism may be operational in other tropical tree species.

Four decades of implicit Monte Carlo

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

Wollaber, Allan B.

In 1971, Fleck and Cummings derived a system of equations to enable robust Monte Carlo simulations of time-dependent, thermal radiative transfer problems. Denoted the “Implicit Monte Carlo” (IMC) equations, their solution remains the de facto standard of high-fidelity radiative transfer simulations. Over the course of 44 years, their numerical properties have become better understood, and accuracy enhancements, novel acceleration methods, and variance reduction techniques have been suggested. In this review, we rederive the IMC equations—explicitly highlighting assumptions as they are made—and outfit the equations with a Monte Carlo interpretation. We put the IMC equations in context with other approximate formsmore » of the radiative transfer equations and present a new demonstration of their equivalence to another well-used linearization solved with deterministic transport methods for frequency-independent problems. We discuss physical and numerical limitations of the IMC equations for asymptotically small time steps, stability characteristics and the potential of maximum principle violations for large time steps, and solution behaviors in an asymptotically thick diffusive limit. We provide a new stability analysis for opacities with general monomial dependence on temperature. Here, we consider spatial accuracy limitations of the IMC equations and discussion acceleration and variance reduction techniques.« less

Four decades of implicit Monte Carlo

DOE PAGES

Wollaber, Allan B.

2016-02-23

In 1971, Fleck and Cummings derived a system of equations to enable robust Monte Carlo simulations of time-dependent, thermal radiative transfer problems. Denoted the “Implicit Monte Carlo” (IMC) equations, their solution remains the de facto standard of high-fidelity radiative transfer simulations. Over the course of 44 years, their numerical properties have become better understood, and accuracy enhancements, novel acceleration methods, and variance reduction techniques have been suggested. In this review, we rederive the IMC equations—explicitly highlighting assumptions as they are made—and outfit the equations with a Monte Carlo interpretation. We put the IMC equations in context with other approximate formsmore » of the radiative transfer equations and present a new demonstration of their equivalence to another well-used linearization solved with deterministic transport methods for frequency-independent problems. We discuss physical and numerical limitations of the IMC equations for asymptotically small time steps, stability characteristics and the potential of maximum principle violations for large time steps, and solution behaviors in an asymptotically thick diffusive limit. We provide a new stability analysis for opacities with general monomial dependence on temperature. Here, we consider spatial accuracy limitations of the IMC equations and discussion acceleration and variance reduction techniques.« less

Effect of the scattering delay on time-dependent photon migration in turbid media.

PubMed

Yaroslavsky, I V; Yaroslavsky, A N; Tuchin, V V; Schwarzmaier, H J

1997-09-01

We modified the diffusion approximation of the time-dependent radiative transfer equation to account for a finite scattering delay time. Under the usual assumptions of the diffusion approximation, the effect of the scattering delay leads to a simple renormalization of the light velocity that appears in the diffusion equation. Accuracy of the model was evaluated by comparison with Monte Carlo simulations in the frequency domain for a semi-infinite geometry. A good agreement is demonstrated for both matched and mismatched boundary conditions when the distance from the source is sufficiently large. The modified diffusion model predicts that the neglect of the scattering delay when the optical properties of the turbid material are derived from normalized frequency- or time-domain measurements should result in an underestimation of the absorption coefficient and an overestimation of the transport coefficient. These observations are consistent with the published experimental data.

Computation of diffuse sky irradiance from multidirectional radiance measurements

NASA Technical Reports Server (NTRS)

Ahmad, Suraiya P.; Middleton, Elizabeth M.; Deering, Donald W.

1987-01-01

Accurate determination of the diffuse solar spectral irradiance directly above the land surface is important in characterizing the reflectance properties of these surfaces, especially vegetation canopies. This determination is also needed to infer the net radiation budget of the earth-atmosphere system above these surfaces. An algorithm is developed here for the computation of hemispheric diffuse irradiance using the measurements from an instrument called PARABOLA, which rapidly measures upwelling and downwelling radiances in three selected wavelength bands. The validity of the algorithm is established from simulations. The standard reference data set of diffuse radiances of Dave (1978), obtained by solving the radiative transfer equation numerically for realistic atmospheric models, is used to simulate PARABOLA radiances. Hemispheric diffuse irradiance is estimated from a subset of simulated radiances by using the algorithm described. The algorithm is validated by comparing the estimated diffuse irradiance with the true diffuse irradiance of the standard data set. The validations include sensitivity studies for two wavelength bands (visible, 0.65-0.67 micron; near infrared, 0.81-0.84 micron), different atmospheric conditions, solar elevations, and surface reflectances. In most cases the hemispheric diffuse irradiance computed from simulated PARABOLA radiances and the true irradiance obtained from radiative transfer calculations agree within 1-2 percent. This technique can be applied to other sampling instruments designed to estimate hemispheric diffuse sky irradiance.

Dynamic Theory of Relativistic Electrons Stochastic Heating by Whistler Mode Waves with Application to the Earth Magnetosphere

NASA Technical Reports Server (NTRS)

Khazanov, G. V.; Tel'nikhin, A. A.; Kronberg, T. K.

2007-01-01

In the Hamiltonian approach an electron motion in a coherent packet of the whistler mode waves propagating along the direction of an ambient magnetic field is studied. The physical processes by which these particles are accelerated to high energy are established. Equations governing a particle motion were transformed in to a closed pair of nonlinear difference equations. The solutions of these equations have shown there exists the energetic threshold below that the electron motion is regular, and when the initial energy is above the threshold an electron moves stochastically. Particle energy spectra and pitch angle electron scattering are described by the Fokker-Planck-Kolmogorov equations. Calculating the stochastic diffusion of electrons due to a spectrum of whistler modes is presented. The parametric dependence of the diffusion coefficients on the plasma particle density, magnitude of wave field, and the strength of magnetic field is studies. It is shown that significant pitch angle diffusion occurs for the Earth radiation belt electrons with energies from a few keV up to a few MeV.

THE EFFECT OF DIFFUSION ON THE PARTICLE SPECTRA IN PULSAR WIND NEBULAE

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

Vorster, M. J.; Moraal, H., E-mail: 12792322@nwu.ac.za

2013-03-01

A possible way to calculate particle spectra as a function of position in pulsar wind nebulae is to solve a Fokker-Planck transport equation. This paper presents numerical solutions to the transport equation with the processes of convection, diffusion, adiabatic losses, and synchrotron radiation included. In the first part of the paper, the steady-state version of the transport equation is solved as a function of position and energy. This is done to distinguish the various effects of the aforementioned processes on the solutions to the transport equation. The second part of the paper deals with a time-dependent solution to the transportmore » equation, specifically taking into account the effect of a moving outer boundary. The paper highlights the fact that diffusion can play a significant role in reducing the amount of synchrotron losses, leading to a modification in the expected particle spectra. These modified spectra can explain the change in the photon index of the synchrotron emission as a function of position. The solutions presented in this paper are not limited to pulsar wind nebulae, but can be applied to any similar central source system, e.g., globular clusters.« less

Diffuse reflectance relations based on diffusion dipole theory for large absorption and reduced scattering

NASA Astrophysics Data System (ADS)

Bremmer, Rolf H.; van Gemert, Martin J. C.; Faber, Dirk J.; van Leeuwen, Ton G.; Aalders, Maurice C. G.

2013-08-01

Diffuse reflectance spectra are used to determine the optical properties of biological samples. In medicine and forensic science, the turbid objects under study often possess large absorption and/or scattering properties. However, data analysis is frequently based on the diffusion approximation to the radiative transfer equation, implying that it is limited to tissues where the reduced scattering coefficient dominates over the absorption coefficient. Nevertheless, up to absorption coefficients of 20 m at reduced scattering coefficients of 1 and 11.5 mm-1, we observed excellent agreement (r2=0.994) between reflectance measurements of phantoms and the diffuse reflectance equation proposed by Zonios et al. [Appl. Opt. 38, 6628-6637 (1999)], derived as an approximation to one of the diffusion dipole equations of Farrell et al. [Med. Phys. 19, 879-888 (1992)]. However, two parameters were fitted to all phantom experiments, including strongly absorbing samples, implying that the reflectance equation differs from diffusion theory. Yet, the exact diffusion dipole approximation at high reduced scattering and absorption also showed agreement with the phantom measurements. The mathematical structure of the diffuse reflectance relation used, derived by Zonios et al. [Appl. Opt. 38, 6628-6637 (1999)], explains this observation. In conclusion, diffuse reflectance relations derived as an approximation to the diffusion dipole theory of Farrell et al. can analyze reflectance ratios accurately, even for much larger absorption than reduced scattering coefficients. This allows calibration of fiber-probe set-ups so that the object's diffuse reflectance can be related to its absorption even when large. These findings will greatly expand the application of diffuse reflection spectroscopy. In medicine, it may allow the use of blue/green wavelengths and measurements on whole blood, and in forensic science, it may allow inclusion of objects such as blood stains and cloth at crime scenes.

Inverse atmospheric radiative transfer problems - A nonlinear minimization search method of solution. [aerosol pollution monitoring

NASA Technical Reports Server (NTRS)

Fymat, A. L.

1976-01-01

The paper studies the inversion of the radiative transfer equation describing the interaction of electromagnetic radiation with atmospheric aerosols. The interaction can be considered as the propagation in the aerosol medium of two light beams: the direct beam in the line-of-sight attenuated by absorption and scattering, and the diffuse beam arising from scattering into the viewing direction, which propagates more or less in random fashion. The latter beam has single scattering and multiple scattering contributions. In the former case and for single scattering, the problem is reducible to first-kind Fredholm equations, while for multiple scattering it is necessary to invert partial integrodifferential equations. A nonlinear minimization search method, applicable to the solution of both types of problems has been developed, and is applied here to the problem of monitoring aerosol pollution, namely the complex refractive index and size distribution of aerosol particles.

Radial Diffusion study of the 1 June 2013 CME event using MHD simulations.

NASA Astrophysics Data System (ADS)

Patel, M.; Hudson, M.; Wiltberger, M. J.; Li, Z.; Boyd, A. J.

2016-12-01

The June 1, 2013 storm was a CME-shock driven geomagnetic storm (Dst = -119 nT) that caused a dropout affecting all radiation belt electron energies measured by the Energetic Particle, Composition and Thermal Plasma Suite (ECT) instrument on Van Allen Probes at higher L-shells following dynamic pressure enhancement in the solar wind. Lower energies (up to about 700 keV) were enhanced by the storm while MeV electrons were depleted throughout the belt. We focus on depletion through radial diffusion caused by the enhanced ULF wave activity due to the CME-shock. This study utilities the Lyon-Fedder-Mobarry (LFM) model, a 3D global magnetospheric simulation code based on the ideal MHD equations, coupled with the Magnetosphere Ionosphere Coupler (MIX) and Rice Convection Model (RCM). The MHD electric and magnetic fields with equations described by Fei et al. [JGR, 2006] are used to calculate radial diffusion coefficients (DLL). These DLL values are input into a radial diffusion code to recreate the dropouts observed by the Van Allen Probes. The importance of understanding the complex role that ULF waves play in radial transport and the effects of CME-driven storms on the relativistic energy electrons in the radiation belts can be accomplished using MHD simulations to obtain diffusion coefficients, initial phase space density and the outer boundary condition from the ECT instrument suite and a radial diffusion model to reproduce observed fluxes which compare favorably with Van Allen Probes ECT measurements.

A Two-moment Radiation Hydrodynamics Module in ATHENA Using a Godunov Method

NASA Astrophysics Data System (ADS)

Skinner, M. A.; Ostriker, E. C.

2013-04-01

We describe a module for the Athena code that solves the grey equations of radiation hydrodynamics (RHD) using a local variable Eddington tensor (VET) based on the M1 closure of the two-moment hierarchy of the transfer equation. The variables are updated via a combination of explicit Godunov methods to advance the gas and radiation variables including the non-stiff source terms, and a local implicit method to integrate the stiff source terms. We employ the reduced speed of light approximation (RSLA) with subcycling of the radiation variables in order to reduce computational costs. The streaming and diffusion limits are well-described by the M1 closure model, and our implementation shows excellent behavior for problems containing both regimes simultaneously. Our operator-split method is ideally suited for problems with a slowly-varying radiation field and dynamical gas flows, in which the effect of the RSLA is minimal.

Development of a robust modeling tool for radiation-induced segregation in austenitic stainless steels

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

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

2015-09-01

Irradiation-assisted stress corrosion cracking (IASCC) of austenitic stainless steels in Light Water Reactor (LWR) components has been linked to changes in grain boundary composition due to irradiation induced segregation (RIS). This work developed a robust RIS modeling tool to account for thermodynamics and kinetics of the atom and defect transportation under combined thermal and radiation conditions. The diffusion flux equations were based on the Perks model formulated through the linear theory of the thermodynamics of irreversible processes. Both cross and non-cross phenomenological diffusion coefficients in the flux equations were considered and correlated to tracer diffusion coefficients through Manning’s relation. Themore » preferential atomvacancy coupling was described by the mobility model, whereas the preferential atom-interstitial coupling was described by the interstitial binding model. The composition dependence of the thermodynamic factor was modeled using the CALPHAD approach. Detailed analysis on the diffusion fluxes near and at grain boundaries of irradiated austenitic stainless steels suggested the dominant diffusion mechanism for chromium and iron is via vacancy, while that for nickel can swing from the vacancy to the interstitial dominant mechanism. The diffusion flux in the vicinity of a grain boundary was found to be greatly influenced by the composition gradient formed from the transient state, leading to the oscillatory behavior of alloy compositions in this region. This work confirms that both vacancy and interstitial diffusion, and segregation itself, have important roles in determining the microchemistry of Fe, Cr, and Ni at irradiated grain boundaries in austenitic stainless steels.« less

Modeling charge collection efficiency degradation in partially depleted GaAs photodiodes using the 1- and 2-carrier Hecht equations

DOE PAGES

Auden, E. C.; Vizkelethy, G.; Serkland, D. K.; ...

2017-03-24

Here, the Hecht equation can be used to model the nonlinear degradation of charge collection efficiency (CCE) in response to radiation-induced displacement damage in both fully and partially depleted GaAs photodiodes. CCE degradation is measured for laser-generated photocurrent as a function of fluence and bias in Al 0.3Ga 0.7As/GaAs/Al 0.25Ga 0.75As p-i-n photodiodes which have been irradiated with 12 MeV C and 7.5 MeV Si ions. CCE is observed to degrade more rapidly with fluence in partially depleted photodiodes than in fully depleted photodiodes. When the intrinsic GaAs layer is fully depleted, the 2-carrier Hecht equation describes CCE degradation asmore » photogenerated electrons and holes recombine at defect sites created by radiation damage in the depletion region. If the GaAs layer is partially depleted, CCE degradation is more appropriately modeled as the sum of the 2-carrier Hecht equation applied to electrons and holes generated within the depletion region and the 1-carrier Hecht equation applied to minority carriers that diffuse from the field-free (non-depleted) region into the depletion region. Enhanced CCE degradation is attributed to holes that recombine within the field-free region of the partially depleted intrinsic GaAs layer before they can diffuse into the depletion region.« less

Modeling charge collection efficiency degradation in partially depleted GaAs photodiodes using the 1- and 2-carrier Hecht equations

NASA Astrophysics Data System (ADS)

Auden, E. C.; Vizkelethy, G.; Serkland, D. K.; Bossert, D. J.; Doyle, B. L.

2017-05-01

The Hecht equation can be used to model the nonlinear degradation of charge collection efficiency (CCE) in response to radiation-induced displacement damage in both fully and partially depleted GaAs photodiodes. CCE degradation is measured for laser-generated photocurrent as a function of fluence and bias in Al0.3Ga0.7As/GaAs/Al0.25Ga0.75As p-i-n photodiodes which have been irradiated with 12 MeV C and 7.5 MeV Si ions. CCE is observed to degrade more rapidly with fluence in partially depleted photodiodes than in fully depleted photodiodes. When the intrinsic GaAs layer is fully depleted, the 2-carrier Hecht equation describes CCE degradation as photogenerated electrons and holes recombine at defect sites created by radiation damage in the depletion region. If the GaAs layer is partially depleted, CCE degradation is more appropriately modeled as the sum of the 2-carrier Hecht equation applied to electrons and holes generated within the depletion region and the 1-carrier Hecht equation applied to minority carriers that diffuse from the field-free (non-depleted) region into the depletion region. Enhanced CCE degradation is attributed to holes that recombine within the field-free region of the partially depleted intrinsic GaAs layer before they can diffuse into the depletion region.

Modeling charge collection efficiency degradation in partially depleted GaAs photodiodes using the 1- and 2-carrier Hecht equations

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

Auden, E. C.; Vizkelethy, G.; Serkland, D. K.

Here, the Hecht equation can be used to model the nonlinear degradation of charge collection efficiency (CCE) in response to radiation-induced displacement damage in both fully and partially depleted GaAs photodiodes. CCE degradation is measured for laser-generated photocurrent as a function of fluence and bias in Al 0.3Ga 0.7As/GaAs/Al 0.25Ga 0.75As p-i-n photodiodes which have been irradiated with 12 MeV C and 7.5 MeV Si ions. CCE is observed to degrade more rapidly with fluence in partially depleted photodiodes than in fully depleted photodiodes. When the intrinsic GaAs layer is fully depleted, the 2-carrier Hecht equation describes CCE degradation asmore » photogenerated electrons and holes recombine at defect sites created by radiation damage in the depletion region. If the GaAs layer is partially depleted, CCE degradation is more appropriately modeled as the sum of the 2-carrier Hecht equation applied to electrons and holes generated within the depletion region and the 1-carrier Hecht equation applied to minority carriers that diffuse from the field-free (non-depleted) region into the depletion region. Enhanced CCE degradation is attributed to holes that recombine within the field-free region of the partially depleted intrinsic GaAs layer before they can diffuse into the depletion region.« less

Diffuse reflection from a stochastically bounded, semi-infinite medium

NASA Technical Reports Server (NTRS)

Lumme, K.; Peltoniemi, J. I.; Irvine, W. M.

1990-01-01

In order to determine the diffuse reflection from a medium bounded by a rough surface, the problem of radiative transfer in a boundary layer characterized by a statistical distribution of heights is considered. For the case that the surface is defined by a multivariate normal probability density, the propagation probability for rays traversing the boundary layer is derived and, from that probability, a corresponding radiative transfer equation. A solution of the Eddington (two stream) type is found explicitly, and examples are given. The results should be applicable to reflection from the regoliths of solar system bodies, as well as from a rough ocean surface.

Modulated wave formation in myocardial cells under electromagnetic radiation

NASA Astrophysics Data System (ADS)

Takembo, Clovis N.; Mvogo, A.; Ekobena Fouda, H. P.; Kofané, T. C.

2018-06-01

We exclusively analyze the onset and condition of formation of modulated waves in a diffusive FitzHugh-Nagumo model for myocardial cell excitations. The cells are connected through gap junction coupling. An additive magnetic flux variable is used to describe the effect of electromagnetic induction, while electromagnetic radiation is imposed on the magnetic flux variable as a periodic forcing. We used the discrete multiple scale expansion and obtained, from the model equations, a single differential-difference amplitude nonlinear equation. We performed the linear stability analysis of this equation and found that instability features are importantly influenced by the induced electromagnetic gain. We present the unstable and stable regions of modulational instability (MI). The resulting analytic predictions are confirmed by numerical experiments of the generic equations. The results reveal that due to MI, an initial steady state that consisted of a plane wave with low amplitude evolves into a modulated localized wave patterns, soliton-like in shape, with features of synchronization. Furthermore, the formation of periodic pulse train with breathing motion presents a disappearing pattern in the presence of electromagnetic radiation. This could provide guidance and better understanding of sudden heart failure exposed to heavily electromagnetic radiation.

Fermi-Compton scattering due to magnetopause surface fluctuations in Jupiter's magnetospheric cavity

NASA Technical Reports Server (NTRS)

Barbosa, D. D.

1981-01-01

The effects of boundary surface fluctuations on a spectrum of electromagnetic radiation trapped in a high Q (quality) cavity are considered. Undulating walls introduce small frequency shifts at reflection to the radiation, and it is argued that the process is entirely analogous to both Fermi (particle) acceleration and inverse Compton scattering. A Fokker-Planck formalism is pursued; it yields a diffusion equation in frequency for which the Green's function and steady-state solutions are found. Applying this analysis to the Jovian continuum radiation discovered by Voyager spacecraft, it is suggested that characteristic diffusion times are greater than 1 year, and that in order to account for the steep frequency spectra observed, an unidentified loss mechanism must operate in the cavity with a decay time constant approximately equal to the characteristic diffusion time divided by 28. A radiator-reactor model of the cavity is investigated to provide an estimate for the intrinsic luminosity of the low frequency (approximately 100 Hz) continuum source whose power is approximately 7 x 10 to the 6th W.

Estimation of optimal hologram recording modes on photothermal materials

NASA Astrophysics Data System (ADS)

Dzhamankyzov, Nasipbek Kurmanalievich; Ismanov, Yusupzhan Khakimzhanovich; Zhumaliev, Kubanychbek Myrzabekovich; Alymkulov, Samsaly Amanovich

2018-01-01

A theoretical analysis of the hologram recording process on photothermal media to estimate the required laser radiation power for the information recording as the function of the spatial frequency and radiation exposure duration is considered. Results of the analysis showed that materials with a low thermal diffusivity are necessary to increase the recording density in these media and the recording should be performed with short pulses to minimize the thermal diffusion length. A solution for the heat conduction equation for photothermal materials heated by an interference laser field was found. The solution obtained allows one to determine the required value of the recording temperature for given spatial frequencies, depending on the thermal physical parameters of the medium and on the power and duration of the heating radiation.

Simulation of the radiolysis of water using Green's functions of the diffusion equation.

PubMed

Plante, I; Cucinotta, F A

2015-09-01

Radiation chemistry is of fundamental importance in the understanding of the effects of ionising radiation, notably with regard to DNA damage by indirect effect (e.g. damage by ·OH radicals created by the radiolysis of water). In the recent years, Green's functions of the diffusion equation (GFDEs) have been used extensively in biochemistry, notably to simulate biochemical networks in time and space. In the present work, an approach based on the GFDE will be used to refine existing models on the indirect effect of ionising radiation on DNA. As a starting point, the code RITRACKS (relativistic ion tracks) will be used to simulate the radiation track structure and calculate the position of all radiolytic species formed during irradiation. The chemical reactions between these radiolytic species and with DNA will be done by using an efficient Monte Carlo sampling algorithm for the GFDE of reversible reactions with an intermediate state that has been developed recently. These simulations should help the understanding of the contribution of the indirect effect in the formation of DNA damage, particularly with regards to the formation of double-strand breaks. © The Author 2015. Published by Oxford University Press. All rights reserved. For Permissions, please email: journals.permissions@oup.com.

Relativistic theory of particles in a scattering flow III: photon transport.

NASA Astrophysics Data System (ADS)

Achterberg, A.; Norman, C. A.

2018-06-01

We use the theory developed in Achterberg & Norman (2018a) and Achterberg & Norman (2018b) to calculate the stress due to photons that are scattered elastically by a relativistic flow. We show that the energy-momentum tensor of the radiation takes the form proposed by Eckart (1940). In particular we show that no terms associated with a bulk viscosity appear if one makes the diffusion approximation for radiation transport and treats the radiation as a separate fluid. We find only shear (dynamic) viscosity terms and heat flow terms in our expression for the energy-momentum tensor. This conclusion holds quite generally for different forms of scattering: Krook-type integral scattering, diffusive (Fokker-Planck) scattering and Thomson scattering. We also derive the transport equation in the diffusion approximation that shows the effects of the flow on the photon gas in the form of a combination of adiabatic heating and an irreversible heating term. We find no diffusive changes to the comoving number density and energy density of the scattered photons, in contrast with some published results in Radiation Hydrodynamics. It is demonstrated that these diffusive corrections to the number- and energy density of the photons are in fact higher-order terms that can (and should) be neglected in the diffusion approximation. Our approach eliminates these terms at the root of the expansion that yields the anisotropic terms in the phase-space density of particles and photons, the terms responsible for the photon viscosity.

Specular reflection treatment for the 3D radiative transfer equation solved with the discrete ordinates method

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

Le Hardy, D.; Favennec, Y., E-mail: yann.favennec@univ-nantes.fr; Rousseau, B.

The contribution of this paper relies in the development of numerical algorithms for the mathematical treatment of specular reflection on borders when dealing with the numerical solution of radiative transfer problems. The radiative transfer equation being integro-differential, the discrete ordinates method allows to write down a set of semi-discrete equations in which weights are to be calculated. The calculation of these weights is well known to be based on either a quadrature or on angular discretization, making the use of such method straightforward for the state equation. Also, the diffuse contribution of reflection on borders is usually well taken intomore » account. However, the calculation of accurate partition ratio coefficients is much more tricky for the specular condition applied on arbitrary geometrical borders. This paper presents algorithms that calculate analytically partition ratio coefficients needed in numerical treatments. The developed algorithms, combined with a decentered finite element scheme, are validated with the help of comparisons with analytical solutions before being applied on complex geometries.« less

Variable Refractive Index Effects on Radiation in Semitransparent Scattering Multilayered Regions

NASA Technical Reports Server (NTRS)

Siegel, R.; Spuckler, C. M.

1993-01-01

A simple set of equations is derived for predicting the temperature distribution and radiative energy flow in a semitransparent layer consisting of an arbitrary number of laminated sublayers that absorb, emit, and scatter radiation. Each sublayer can have a different refractive index and optical thickness. The plane composite region is heated on each exterior side by a different amount of incident radiation. The results are for the limiting case where heat conduction within the layers is very small relative to radiative transfer, and is neglected. The interfaces are assumed diffuse, and all interface reflections are included in the analysis. The thermal behavior is readily calculated from the analytical expressions that are obtained. By using many sublayers, expressions provide the temperature distribution and heat flow for a diffusing medium with a continually varying refractive index, including internal reflection effects caused by refractive index gradients. Temperature and heat flux results are given to show the effect of variations in refractive index and optical thickness through the multilayer laminate.

Variable Refractive Index Effects on Radiation in Semitransparent Scattering Multilayered Regions

NASA Technical Reports Server (NTRS)

Siegel, R.; Spuckler, C. M.

1993-01-01

A simple set of equations is derived for predicting the temperature distribution and radiative energy flow in a semitransparent layer consisting of an arbitrary number of laminated sublayers that absorb, emit, and scatter radiation. Each sublayer can have a different refractive index and optical thickness. The plane composite region is heated on each exterior side by a different amount of incident radiation. The results are for the limiting case where heat conduction within the layers is very small relative to radiative transfer, and is neglected. The interfaces are assumed diffuse, and all interface reflections are included in the analysis. The thermal behavior is readily calculated from the analytical expressions that are obtained. By using many sublayers, the analytical expressions provide the temperature distribution and heat flow for a diffusing medium with a continuously varying refractive index, including internal reflection effects caused by refractive index gradients. Temperature and heat flux results are given to show the effect of variations in refractive index and optical thickness through the multilayer laminate.

Radiative transfer in dusty nebulae. III - The effects of dust albedo

NASA Technical Reports Server (NTRS)

Petrosian, V.; Dana, R. A.

1980-01-01

The effects of an albedo of internal dust, such as ionization structure and temperature of dust grain, were studied by the quasi-diffusion method with an iterative technique for solving the radiative heat transfer equations. It was found that the generalized on-the-spot approximation solution is adequate for most astrophysical applications for a zero albedo; for a nonzero albedo, the Eddington approximation is more accurate. The albedo increases the average energy of the diffuse photons, increasing the ionization level of hydrogen and heavy elements if the Eddington approximation is applied; the dust thermal gradient is reduced so that the infrared spectrum approaches blackbody spectrum with an increasing albedo.

Numerical study of photon migration in the presence of a void region using the radiative transfer and diffusion equations

NASA Astrophysics Data System (ADS)

Miyakawa, Erina; Fujii, Hiroyuki; Hattori, Kiyohito; Tatekura, Yuki; Kobayashi, Kazumichi; Watanabe, Masao

2016-12-01

Diffuse optical tomography (DOT), which is still under development, has a potential to enable non-invasive diagnoses of thyroid cancers in the human neck using the near-infrared light. This modality needs a photon migration model because scattered light is used. There are two types of photon migration models: the radiative transport equation (RTE) and diffusion equation (DE). The RTE can describe photon migration in the human neck with accuracy, while the DE enables an efficient calculation. For developing the accurate and efficient model of photon migration, it is crucial to investigate a condition where the DE holds in a scattering medium including a void region under the refractive-index mismatch at the void boundary because the human neck has a trachea (void region) and the refractive indices are different between the human neck and trachea. Hence, in this paper, we compare photon migration using the RTE with that using the DE in the medium. The numerical results show that the DE is valid under the refractive-index match at the void boundary even though the void region is near the source and detector positions. Under the refractive-index mismatch at the boundary, the numerical results using the DE disagree with those using the RTE when the void region is near the source and detector positions. This is probably because the anisotropy of the light scattering remains around the void boundary.

Measurements and Modeling of Soot Formation and Radiation in Microgravity Jet Diffusion Flames. Volume 4

NASA Technical Reports Server (NTRS)

Ku, Jerry C.; Tong, Li; Greenberg, Paul S.

1996-01-01

This is a computational and experimental study for soot formation and radiative heat transfer in jet diffusion flames under normal gravity (1-g) and microgravity (0-g) conditions. Instantaneous soot volume fraction maps are measured using a full-field imaging absorption technique developed by the authors. A compact, self-contained drop rig is used for microgravity experiments in the 2.2-second drop tower facility at NASA Lewis Research Center. On modeling, we have coupled flame structure and soot formation models with detailed radiation transfer calculations. Favre-averaged boundary layer equations with a k-e-g turbulence model are used to predict the flow field, and a conserved scalar approach with an assumed Beta-pdf are used to predict gaseous species mole fraction. Scalar transport equations are used to describe soot volume fraction and number density distributions, with formation and oxidation terms modeled by one-step rate equations and thermophoretic effects included. An energy equation is included to couple flame structure and radiation analyses through iterations, neglecting turbulence-radiation interactions. The YIX solution for a finite cylindrical enclosure is used for radiative heat transfer calculations. The spectral absorption coefficient for soot aggregates is calculated from the Rayleigh solution using complex refractive index data from a Drude- Lorentz model. The exponential-wide-band model is used to calculate the spectral absorption coefficient for H20 and C02. It is shown that when compared to results from true spectral integration, the Rosseland mean absorption coefficient can provide reasonably accurate predictions for the type of flames studied. The soot formation model proposed by Moss, Syed, and Stewart seems to produce better fits to experimental data and more physically sound than the simpler model by Khan et al. Predicted soot volume fraction and temperature results agree well with published data for a normal gravity co-flow laminar flames and turbulent jet flames. Predicted soot volume fraction results also agree with our data for 1-g and 0-g laminar jet names as well as 1-g turbulent jet flames.

Steady-state solutions of a diffusive energy-balance climate model and their stability

NASA Technical Reports Server (NTRS)

Ghil, M.

1975-01-01

A diffusive energy-balance climate model, governed by a nonlinear parabolic partial differential equation, was studied. Three positive steady-state solutions of this equation are found; they correspond to three possible climates of our planet: an interglacial (nearly identical to the present climate), a glacial, and a completely ice-covered earth. Models similar to the main one are considered, and the number of their steady states was determined. All the models have albedo continuously varying with latitude and temperature, and entirely diffusive horizontal heat transfer. The stability under small perturbations of the main model's climates was investigated. A stability criterion is derived, and its application shows that the present climate and the deep freeze are stable, whereas the model's glacial is unstable. The dependence was examined of the number of steady states and of their stability on the average solar radiation.

Sweat Rate Prediction Equations for Outdoor Exercise with Transient Solar Radiation

DTIC Science & Technology

2012-01-01

AD] 15 Interchangeable variables gSL W/m2 Global solar load Direct weather station data; pyranometer values 25 Direct measurement from weather station ...Fanger equations 2, 4, 13, Direct or weather station values Rdif W Diffuse irradiance Rref W Reflected irradiance AD m2 Body surface area (BSA) from DuBois...assuming the given weather station uses standard meteorological measuring instru- ments. In the heat flow form expressed by Matthew et al. (25

Thermal radiation from large bolides and impact plumes

NASA Astrophysics Data System (ADS)

Svetsov, V.; Shuvalov, V.

2017-09-01

Numerical simulations of the impacts of asteroids and comets from 20 m to 3 km in diameter have been carried out and thermal radiation fluxes on the ground and luminous efficiencies of the impacts have been calculated. It was assumed that the cosmic objects have no strength, deform, fragment, and vaporize in the atmosphere. After the impact on the ground, formation of craters and plumes was simulated taking into account internal friction of destroyed rocks and a wake formed in the atmosphere. The equations of radiative transfer, added to the equations of gas dynamics, were used in the approximation of radiative heat diffusion or, if the Rosseland optical depth of a radiating volume of gas and vapor was less than unity, in the approximation of volume emission. Radiation fluxes on the Earth's surface were calculated by integrating the equation of radiative transfer along rays passing through a luminous area. Direct thermal radiation from fireballs and impact plumes produced by asteroids and comets larger than 50 m in diameter is dangerous for people, animals, plants, economic objects. Forest fires can be ignited on the ground within a radius of roughly 1000 times the body's diameter (for diameters of the order or smaller than 1 km), 50-m-diameter bodies can ignite forest fires within a radius of up to 40 km and 3-km asteroids - within 1700 km.

3D unstructured-mesh radiation transport codes

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

Morel, J.

1997-12-31

Three unstructured-mesh radiation transport codes are currently being developed at Los Alamos National Laboratory. The first code is ATTILA, which uses an unstructured tetrahedral mesh in conjunction with standard Sn (discrete-ordinates) angular discretization, standard multigroup energy discretization, and linear-discontinuous spatial differencing. ATTILA solves the standard first-order form of the transport equation using source iteration in conjunction with diffusion-synthetic acceleration of the within-group source iterations. DANTE is designed to run primarily on workstations. The second code is DANTE, which uses a hybrid finite-element mesh consisting of arbitrary combinations of hexahedra, wedges, pyramids, and tetrahedra. DANTE solves several second-order self-adjoint forms of the transport equation including the even-parity equation, the odd-parity equation, and a new equation called the self-adjoint angular flux equation. DANTE also offers three angular discretization options:more » $$S{_}n$$ (discrete-ordinates), $$P{_}n$$ (spherical harmonics), and $$SP{_}n$$ (simplified spherical harmonics). DANTE is designed to run primarily on massively parallel message-passing machines, such as the ASCI-Blue machines at LANL and LLNL. The third code is PERICLES, which uses the same hybrid finite-element mesh as DANTE, but solves the standard first-order form of the transport equation rather than a second-order self-adjoint form. DANTE uses a standard $$S{_}n$$ discretization in angle in conjunction with trilinear-discontinuous spatial differencing, and diffusion-synthetic acceleration of the within-group source iterations. PERICLES was initially designed to run on workstations, but a version for massively parallel message-passing machines will be built. The three codes will be described in detail and computational results will be presented.« less

Development of response models for the Earth Radiation Budget Experiment (ERBE) sensors. Part 2: Analysis of the ERBE integrating sphere ground calibration

NASA Technical Reports Server (NTRS)

Halyo, Nesim; Taylor, Deborah B.

1987-01-01

An explicit solution of the spectral radiance leaving an arbitrary point on the wall of a spherical cavity with diffuse reflectivity is obtained. The solution is applicable to spheres with an arbitrary number of openings of any size and shape, an arbitrary number of light sources with possible non-diffuse characteristics, a non-uniform sphere wall temperature distribution, non-uniform and non-diffuse sphere wall emissivity and non-uniform but diffuse sphere wall spectral reflectivity. A general measurement equation describing the output of a sensor with a given field of view, angular and spectral response measuring the sphere output is obtained. The results are applied to the Earth Radiation Budget Experiment (ERBE) integrating sphere. The sphere wall radiance uniformity, loading effects and non-uniform wall temperature effects are investigated. It is shown that using appropriate interpretation and processing, a high-accuracy short-wave calibration of the ERBE sensors can be achieved.

Photon diffusion coefficient in scattering and absorbing media.

PubMed

Pierrat, Romain; Greffet, Jean-Jacques; Carminati, Rémi

2006-05-01

We present a unified derivation of the photon diffusion coefficient for both steady-state and time-dependent transport in disordered absorbing media. The derivation is based on a modal analysis of the time-dependent radiative transfer equation. This approach confirms that the dynamic diffusion coefficient is given by the random-walk result D = cl(*)/3, where l(*) is the transport mean free path and c is the energy velocity, independent of the level of absorption. It also shows that the diffusion coefficient for steady-state transport, often used in biomedical optics, depends on absorption, in agreement with recent theoretical and experimental works. These two results resolve a recurrent controversy in light propagation and imaging in scattering media.

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

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

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

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

On the Green's function of the partially diffusion-controlled reversible ABCD reaction for radiation chemistry codes

NASA Astrophysics Data System (ADS)

Plante, Ianik; Devroye, Luc

2015-09-01

Several computer codes simulating chemical reactions in particles systems are based on the Green's functions of the diffusion equation (GFDE). Indeed, many types of chemical systems have been simulated using the exact GFDE, which has also become the gold standard for validating other theoretical models. In this work, a simulation algorithm is presented to sample the interparticle distance for partially diffusion-controlled reversible ABCD reaction. This algorithm is considered exact for 2-particles systems, is faster than conventional look-up tables and uses only a few kilobytes of memory. The simulation results obtained with this method are compared with those obtained with the independent reaction times (IRT) method. This work is part of our effort in developing models to understand the role of chemical reactions in the radiation effects on cells and tissues and may eventually be included in event-based models of space radiation risks. However, as many reactions are of this type in biological systems, this algorithm might play a pivotal role in future simulation programs not only in radiation chemistry, but also in the simulation of biochemical networks in time and space as well.

ULF Waves and Diffusive Radial Transport of Charged Particles

NASA Astrophysics Data System (ADS)

Ali, Ashar Fawad

The Van Allen radiation belts contain highly energetic particles which interact with a variety of plasma and magnetohydrodynamic (MHD) waves. Waves in the ultra low-frequency (ULF) range play an important role in the loss and acceleration of energetic particles. Considering the geometry of the geomagnetic field, charged particles trapped in the inner magnetosphere undergo three distinct types of periodic motions; an adiabatic invariant is associated with each type of motion. The evolution of the phase space density of charged particles in the magnetosphere in the coordinate space of the three adiabatic invariants is modeled by the Fokker-Planck equation. If we assume that the first two adiabatic invariants are conserved while the third invariant is violated, then the general Fokker-Planck equation reduces to a radial diffusion equation with the radial diffusion coefficient quantifying the rate of the radial diffusion of charged particles, including contributions from perturbations in both the magnetic and the electric fields. This thesis investigates two unanswered questions about ULF wave-driven radial transport of charged particles. First, how important are the ULF fluctuations in the magnetic field compared with the ULF fluctuations in the electric field in driving the radial diffusion of charged particles in the Earth's inner magnetosphere? It has generally been accepted that magnetic field perturbations dominate over electric field perturbations, but several recently published studies suggest otherwise. Second, what is the distribution of ULF wave power in azimuth, and how does ULF wave power depend upon radial distance and the level of geomagnetic activity? Analytic treatments of the diffusion coefficients generally assume uniform distribution of power in azimuth, but in situ measurements suggest that this may not be the case. We used the magnetic field data from the Combined Release and Radiation Effects Satellite (CRRES) and the electric and the magnetic field data from the Radiation Belt Storm Probes (RBSP) to compute the electric and the magnetic component of the radial diffusion coefficient using the Fei et al. [2006] formulation. We conclude that contrary to prior notions, the electric component is dominant in driving radial diffusion of charged particles in the Earth's inner magnetosphere instead of the magnetic component. The electric component can be up to two orders of magnitude larger than the magnetic component. In addition, we see that ULF wave power in both the electric and the magnetic fields has a clear dependence on Kp with wave power decreasing as radial distance decreases. For both fields, the noon sectors generally contain more ULF wave power than the dawn, dusk, and the midnight magnetic local time (MLT) sectors. There is no significant difference between ULF wave power in the dawn, dusk, and the midnight sectors.

Brownian diffusion and thermophoresis mechanisms in Casson fluid over a moving wedge

NASA Astrophysics Data System (ADS)

Ullah, Imran; Shafie, Sharidan; Khan, Ilyas; Hsiao, Kai Long

2018-06-01

The effect of Brownian diffusion and thermophoresis on electrically conducting mixed convection flow of Casson fluid induced by moving wedge is investigated in this paper. It is assumed that the wedge is saturated in a porous medium and experiences the thermal radiation and chemical reaction effects. The transformed nonlinear governing equations are solved numerically by Keller box scheme. Findings reveal that increase in Casson and magnetic parameters reduced the boundary layer thickness. The effect of Brownian motion and thermophoresis parameters are more pronounced on temperature profile as compared to nanoparticles concentration. The presence of thermal radiation assisted the heat transfer rate significantly. The influence of magnetic parameter is observed less significant on temperature and nanoparticles concentration.

Variable mass diffusion effects on free convection flow past an impulsively started infinite vertical plate

NASA Astrophysics Data System (ADS)

Rushi Kumar, B.; Jayakar, R.; Vijay Kumar, A. G.

2017-11-01

An exact analysis of the problem of free convection flow of a viscous incompressible chemically reacting fluid past an infinite vertical plate with the flow due to impulsive motion of the plate with Newtonian heating in the presence of thermal radiation and variable mass diffusion is performed. The resulting governing equations were tackled by Laplace transform technique. Finally the effects of pertinent flow parameters such as the radiation parameter, chemical reaction parameter, buoyancy ratio parameter, thermal Grashof number, Schmidt number, Prandtl number and time on the velocity, temperature, concentration and skin friction for both aiding and opposing flows were examined in detail when Pr=0.71(conducting air) and Pr=7.0(water).

Comparison of the Radiative Two-Flux and Diffusion Approximations

NASA Technical Reports Server (NTRS)

Spuckler, Charles M.

2006-01-01

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

Compton Scattering by Static and Moving Media. Part 1; The Transfer Equation and its Moments

NASA Technical Reports Server (NTRS)

Psaltis, Dimitrios; Lamb, Frederick K.

1997-01-01

Compton scattering of photons by nonrelativistic particles is thought to play an important role in forming the radiation spectrum of many astrophysical systems. Here we derive the time-dependent photon kinetic equation that describes spontaneous and induced Compton scattering, as well as absorption and emission by static and moving media, the corresponding radiative transfer equation, and their zeroth and first angular moments, both in the system frame and in the frame comoving with the medium. We show that it is necessary to use the correct relativistic differential scattering cross section in order to obtain a photon kinetic equation that is correct to first order in Epsilon/m(sub e), T(sub e)/m(sub e), and V, where Epsilon is the photon energy, T(sub e) and m(sub e) are the electron temperature and rest mass, and V is the electron bulk velocity in units of the speed of light. We also demonstrate that the terms in the radiative transfer equation that are second order in V should usually be retained, because if the radiation energy density is sufficiently large, compared to the radiation flux, the effects of bulk Comptonization described by the terms that are second order in V can be as important as the effects described by the terms that are first order in V, even when V is small. The system- and fluid-frame equations that we derive are correct to first order in Epsilon/m(sub e). Our system-frame equations, which are correct to second order in V, may be used when V is not too large. Our fluid-frame equations, which are exact in V, may be used when V approaches 1. Both sets of equations are valid for systems of arbitrary optical depth and can therefore be used in both the free-streaming and diffusion regimes. We demonstrate that Comptonization by the electron bulk motion occurs whether or not the radiation field is isotropic or the bulk flow converges and that it is more important than thermal Comptonization if V(sup 2) is greater than 3T(sub e)/m(sub e).

A stochastic method for Brownian-like optical transport calculations in anisotropic biosuspensions and blood

NASA Astrophysics Data System (ADS)

Miller, Steven

1998-03-01

A generic stochastic method is presented that rapidly evaluates numerical bulk flux solutions to the one-dimensional integrodifferential radiative transport equation, for coherent irradiance of optically anisotropic suspensions of nonspheroidal bioparticles, such as blood. As Fermat rays or geodesics enter the suspension, they evolve into a bundle of random paths or trajectories due to scattering by the suspended bioparticles. Overall, this can be interpreted as a bundle of Markov trajectories traced out by a "gas" of Brownian-like point photons being scattered and absorbed by the homogeneous distribution of uncorrelated cells in suspension. By considering the cumulative vectorial intersections of a statistical bundle of random trajectories through sets of interior data planes in the space containing the medium, the effective equivalent information content and behavior of the (generally unknown) analytical flux solutions of the radiative transfer equation rapidly emerges. The fluxes match the analytical diffuse flux solutions in the diffusion limit, which verifies the accuracy of the algorithm. The method is not constrained by the diffusion limit and gives correct solutions for conditions where diffuse solutions are not viable. Unlike conventional Monte Carlo and numerical techniques adapted from neutron transport or nuclear reactor problems that compute scalar quantities, this vectorial technique is fast, easily implemented, adaptable, and viable for a wide class of biophotonic scenarios. By comparison, other analytical or numerical techniques generally become unwieldy, lack viability, or are more difficult to utilize and adapt. Illustrative calculations are presented for blood medias at monochromatic wavelengths in the visible spectrum.

Functional imaging of small tissue volumes with diffuse optical tomography

NASA Astrophysics Data System (ADS)

Klose, Alexander D.; Hielscher, Andreas H.

2006-03-01

Imaging of dynamic changes in blood parameters, functional brain imaging, and tumor imaging are the most advanced application areas of diffuse optical tomography (DOT). When dealing with the image reconstruction problem one is faced with the fact that near-infrared photons, unlike X-rays, are highly scattered when they traverse biological tissue. Image reconstruction schemes are required that model the light propagation inside biological tissue and predict measurements on the tissue surface. By iteratively changing the tissue-parameters until the predictions agree with the real measurements, a spatial distribution of optical properties inside the tissue is found. The optical properties can be related to the tissue oxygenation, inflammation, or to the fluorophore concentration of a biochemical marker. If the model of light propagation is inaccurate, the reconstruction process will lead to an inaccurate result as well. Here, we focus on difficulties that are encountered when DOT is employed for functional imaging of small tissue volumes, for example, in cancer studies involving small animals, or human finger joints for early diagnosis of rheumatoid arthritis. Most of the currently employed image reconstruction methods rely on the diffusion theory that is an approximation to the equation of radiative transfer. But, in the cases of small tissue volumes and tissues that contain low scattering regions diffusion theory has been shown to be of limited applicability Therefore, we employ a light propagation model that is based on the equation of radiative transfer, which promises to overcome the limitations.

Fractional Diffusion Equations and Anomalous Diffusion

NASA Astrophysics Data System (ADS)

Evangelista, Luiz Roberto; Kaminski Lenzi, Ervin

2018-01-01

Preface; 1. Mathematical preliminaries; 2. A survey of the fractional calculus; 3. From normal to anomalous diffusion; 4. Fractional diffusion equations: elementary applications; 5. Fractional diffusion equations: surface effects; 6. Fractional nonlinear diffusion equation; 7. Anomalous diffusion: anisotropic case; 8. Fractional Schrödinger equations; 9. Anomalous diffusion and impedance spectroscopy; 10. The Poisson–Nernst–Planck anomalous (PNPA) models; References; Index.

Modeling and Simulation for Particle Radiation Damage to Electronic and Opto-Electronic Devices

DTIC Science & Technology

2018-01-25

AFRL-RV-PS- AFRL-RV-PS- TR-2018-0001 TR-2018-0001 MODELING & SIMULATION FOR PARTICLE RADIATION DAMAGE TO ELECTRONIC AND OPTO- ELECTRONIC DEVICES... Electronic and Opto- Electronic Devices 5a. CONTRACT NUMBER FA9453-14-1-0248 5b. GRANT NUMBER 5c. PROGRAM ELEMENT NUMBER 62601F 6. AUTHOR(S) Sanjay...nBp Diode Assuming the light is incident on the n-side of the photodiode, the drift-diffusion equation for the minority electron of the p-type

Similarity solutions of reaction–diffusion equation with space- and time-dependent diffusion and reaction terms

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

Ho, C.-L.; Lee, C.-C., E-mail: chieh.no27@gmail.com

2016-01-15

We consider solvability of the generalized reaction–diffusion equation with both space- and time-dependent diffusion and reaction terms by means of the similarity method. By introducing the similarity variable, the reaction–diffusion equation is reduced to an ordinary differential equation. Matching the resulting ordinary differential equation with known exactly solvable equations, one can obtain corresponding exactly solvable reaction–diffusion systems. Several representative examples of exactly solvable reaction–diffusion equations are presented.

An Infrared Survey of the Diffuse Emission within 5 deg of the Galactic Plane.

DTIC Science & Technology

1980-06-05

t O ±60. Over the region of 100 to 3 0 oi longitude along the galactic equator, this emission can be fit by 500( K black -body emission with a dilution...from the AFGL catalog, which they classify as stars. The assumed background is, therefore, composed of black -body radiators with a characteristic...SUPPLEMENTARY NOTES 19 KEY WORDS (c-nIIl, ,l IY ,I. AIIId-1, hI MI’< A III-15SI, Infra red Diffuse emission Galactic structure 1111 regions yI 40

FOKKER-PLANCK ANALYSIS OF TRANSVERSE COLLECTIVE INSTABILITIES IN ELECTRON STORAGE RINGS

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

Lindberg, R. R.

We analyze single bunch transverse instabilities due to wakefields using a Fokker-Planck model. We expand on the work of Suzuki [1], writing out the linear matrix equation including chromaticity, both dipolar and quadrupolar transverse wakefields, and the effects of damping and diffusion due to the synchrotron radiation. The eigenvalues and eigenvectors determine the collective stability of the beam, and we show that the predicted threshold current for transverse instability and the profile of the unstable agree well with tracking simulations. In particular, we find that predicting collective stability for high energy electron beams at moderate to large values of chromaticitymore » requires the full Fokker-Planck analysis to properly account for the effects of damping and diffusion due to synchrotron radiation.« less

A Theoretical Understanding of Circular Polarization Memory in Random Media

NASA Astrophysics Data System (ADS)

Dark, Julia

Radiative transport theory describes the propagation of light in random media that absorb, scatter, and emit radiation. To describe the propagation of light, the full polarization state is quantified using the Stokes parameters. For the sake of mathematical convenience, the polarization state of light is often neglected leading to the scalar radiative transport equation for the intensity only. For scalar transport theory, there is a well-established body of literature on numerical and analytic approximations to the radiative transport equation. We extend the scalar theory to the vector radiative transport equation (vRTE). In particular, we are interested in the theoretical basis for a phenomena called circular polarization memory. Circular polarization memory is the physical phenomena whereby circular polarization retains its ellipticity and handedness when propagating in random media. This is in contrast to the propagation of linear polarization in random media, which depolarizes at a faster rate, and specular reflection of circular polarization, whereby the circular polarization handedness flips. We investigate two limits that are of known interest in the phenomena of circular polarization memory. The first limit we investigate is that of forward-peaked scattering, i.e. the limit where most scattering events occur in the forward or near-forward directions. The second limit we consider is that of strong scattering and weak absorption. In the forward-peaked scattering limit we approximate the vRTE by a system of partial differential equations motivated by the scalar Fokker-Planck approximation. We call the leading order approximation the vector Fokker-Planck approximation. The vector Fokker Planck approximation predicts that strongly forward-peaked media exhibit circular polarization memory where the strength of the effect can be calculated from the expansion of the scattering matrix in special functions. In addition, we find in this limit that total intensity, linear polarization, and circular polarization decouple. From this result we conclude, that in the Fokker-Planck limit the scalar approximation is an appropriate leading order approximation. In the strong scattering and weak absorbing limit the vector radiative transport equation can be analyzed using boundary layer theory. In this case, the problem of light scattering in an optically thick medium is reduced to a 1D vRTE near the boundary and a 3D diffusion equation in the interior. We develop and implement a numerical solver for the boundary layer problem by using a discrete ordinate solver in the boundary layer and a spectral method to solve the diffusion approximation in the interior. We implement the method in Fortran 95 with external dependencies on BLAS, LAPACK, and FFTW. By analyzing the spectrum of the discretized vRTE in the boundary layer, we are able to predict the presence of circular polarization memory in a given medium.

Marangoni Convection during Free Electron Laser Nitriding of Titanium

NASA Astrophysics Data System (ADS)

Höche, Daniel; Müller, Sven; Rapin, Gerd; Shinn, Michelle; Remdt, Elvira; Gubisch, Maik; Schaaf, Peter

2009-08-01

Pure titanium was treated by free electron laser (FEL) radiation in a nitrogen atmosphere. As a result, nitrogen diffusion occurs and a TiN coating was synthesized. Local gradients of interfacial tension due to the local heating lead to a Marangoni convection, which determines the track properties. Because of the experimental inaccessibility of time-dependent occurrences, finite element calculations were performed, to determine the physical processes such as heat transfer, melt flow, and mass transport. In order to calculate the surface deformation of the gas-liquid interface, the level set approach was used. The equations were modified and coupled with heat-transfer and diffusion equations. The process was characterized by dimensionless numbers such as the Reynolds, Peclet, and capillary numbers, to obtain more information about the acting forces and the coating development. Moreover, the nitrogen distribution was calculated using the corresponding transport equation. The simulations were compared with cross-sectional micrographs of the treated titanium sheets and checked for their validity. Finally, the process presented is discussed and compared with similar laser treatments.

An analysis of the flow field near the fuel injection location in a gas core reactor.

NASA Technical Reports Server (NTRS)

Weinstein, H.; Murty, B. G. K.; Porter, R. W.

1971-01-01

An analytical study is presented which shows the effects of large energy release and the concurrent high acceleration of inner stream fluid on the coaxial flow field in a gas core reactor. The governing equations include the assumptions of only radial radiative transport of energy represented as an energy diffusion term in the Euler equations. The method of integral relations is used to obtain the numerical solution. Results show that the rapidly accelerating, heat generating inner stream actually shrinks in radius as it expands axially.

Relationship between the Kubelka-Munk scattering and radiative transfer coefficients.

PubMed

Thennadil, Suresh N

2008-07-01

The relationship between the Kubelka-Munk (K-M) and the transport scattering coefficient is obtained through a semi-empirical approach. This approach gives the same result as that given by Gate [Appl. Opt.13, 236 (1974)] when the incident beam is diffuse. This result and those given by Star et al. [Phys. Med. Biol.33, 437 (1988)] and Brinkworth [Appl. Opt.11, 1434 (1972)] are compared with the exact solution of the radiative transfer equation over a large range of optical properties. It is found that the latter expressions, which include an absorption component, do not give accurate results over the range considered. Using the semi-empirical approach, the relationship between the K-M and the transport scattering coefficient is derived for the case where the incident light is collimated. It is shown that although the K-M equation is derived based on diffuse incident light, it can also represent very well the reflectance from a slab of infinite thickness when the incident light is collimated. However, in this case the relationship between the coefficients has to include a function that is dependent on the anisotropy factor. Analysis indicates that the K-M transform achieves the objective of obtaining a measure that gives the ratio of absorption to scattering effects for both diffuse and collimated incident beams over a large range of optical properties.

An analytic solution of the radiative transfer equation for a gray scattering atmosphere in motion

NASA Astrophysics Data System (ADS)

Pistinner, Shlomi; Shaviv, Giora

1994-12-01

We provide a formal analytic solution of the radiative transfer equation for a gray moving atmosphere in a plane parallel geometry. A formal solution in the diffusion and the free-streaming limit is also provided in the case of a spherically extended atmosphere. The formal solutions are written explicitly for scattering atmospheres in which the density and the velocity fields are given by a power law. A self-consistent temperature profile accurate to O(Beta = v/c) is provided for the case in which the absorption or the scattering are temperature independent. The gray extinction temperature profile is considerably simplified in the case of a scattering atmosphere. Steady state flow and homologous expansion are special cases that are considered in detail.

An analytic solution of the radiative transfer equation for a gray scattering atmosphere in motion

NASA Technical Reports Server (NTRS)

Pistinner, Shlomi; Shaviv, Giora

1994-01-01

We provide a formal analytic solution of the radiative transfer equation for a gray moving atmosphere in a plane parallel geometry. A formal solution in the diffusion and the free-streaming limit is also provided in the case of a spherically extended atmosphere. The formal solutions are written explicitly for scattering atmospheres in which the density and the velocity fields are given by a power law. A self-consistent temperature profile accurate to O(Beta = v/c) is provided for the case in which the absorption or the scattering are temperature independent. The gray extinction temperature profile is considerably simplified in the case of a scattering atmosphere. Steady state flow and homologous expansion are special cases that are considered in detail.

Two-Flux and Green's Function Method for Transient Radiative Transfer in a Semi-Transparent Layer

NASA Technical Reports Server (NTRS)

Siegel, Robert

1995-01-01

A method using a Green's function is developed for computing transient temperatures in a semitransparent layer by using the two-flux method coupled with the transient energy equation. Each boundary of the layer is exposed to a hot or cold radiative environment, and is heated or cooled by convection. The layer refractive index is larger than one, and the effect of internal reflections is included with the boundaries assumed diffuse. The analysis accounts for internal emission, absorption, heat conduction, and isotropic scattering. Spectrally dependent radiative properties are included, and transient results are given to illustrate two-band spectral behavior with optically thin and thick bands. Transient results using the present Green's function method are verified for a gray layer by comparison with a finite difference solution of the exact radiative transfer equations; excellent agreement is obtained. The present method requires only moderate computing times and incorporates isotropic scattering without additional complexity. Typical temperature distributions are given to illustrate application of the method by examining the effect of strong radiative heating on one side of a layer with convective cooling on the other side, and the interaction of strong convective heating with radiative cooling from the layer interior.

Some New Lidar Equations for Laser Pulses Scattered Back from Optically Thick Media Such as Clouds, Dense Aerosol Plumes, Sea Ice, Snow, and Turbid Coastal Waters

NASA Technical Reports Server (NTRS)

Davis, Anthony B.

2013-01-01

I survey the theoretical foundations of the slowly-but-surely emerging field of multiple scattering lidar, which has already found applications in atmospheric and cryospheric optics that I also discuss. In multiple scattering lidar, returned pulses are stretched far beyond recognition, and there is no longer a one-to-one connection between range and return-trip timing. Moreover, one can exploit the radial profile of the diffuse radiance field excited by the laser source that, by its very nature, is highly concentrated in space and collimated in direction. One needs, however, a new class of lidar equations to explore this new phenomenology. A very useful set is derived from radiative diffusion theory, which is found at the opposite asymptotic limit of radiative transfer theory than the conventional (single-scattering) limit used to derive the standard lidar equation. In particular, one can use it to show that, even if the simple time-of-flight-to-range connection is irretrievably lost, multiply-scattered lidar light can be used to restore a unique profiling capability with coarser resolution but much deeper penetration into a wide variety of optical thick media in nature. Several new applications are proposed, including a laser bathymetry technique that should work for highly turbid coastal waters.

Diffusion of chemically reactive species in MHD oscillatory flow with thermal radiation in the presence of constant suction and injection

NASA Astrophysics Data System (ADS)

Sasikumar, J.; Bhuvaneshwari, S.; Govindarajan, A.

2018-04-01

In this project, it is proposed to investigate the effect of suction/injection on the unsteady oscillatory flow of an incompressible viscous electrically conducting fluid through a channel filled with porous medium and non-uniform wall temperature. The fluid is subjected to a uniform magnetic field normal to the channel and the velocity slip at the cold plate is taken into consideration. With the assumption of magnetic Reynolds number to be very small, the induced magnetic field is neglected. Assuming pressure gradient to be oscillatory across the ends of the channel, resulting flow as unsteady oscillatory flow. Under the usual Bousinessq approximation, a mathematical model representing this fluid flow consisting of governing equations with boundary conditions will be developed. Closed form solutions of the dimensionless governing equations of the fluid flow, namely momentum equation, energy equation and species concentration can be obtained . The effects of heat radiation and chemical reaction with suction and injection on temperature, velocity and species concentration profiles will be analysed with tables and graphs.

Radiation Heat Transfer Between Diffuse-Gray Surfaces Using Higher Order Finite Elements

NASA Technical Reports Server (NTRS)

Gould, Dana C.

2000-01-01

This paper presents recent work on developing methods for analyzing radiation heat transfer between diffuse-gray surfaces using p-version finite elements. The work was motivated by a thermal analysis of a High Speed Civil Transport (HSCT) wing structure which showed the importance of radiation heat transfer throughout the structure. The analysis also showed that refining the finite element mesh to accurately capture the temperature distribution on the internal structure led to very large meshes with unacceptably long execution times. Traditional methods for calculating surface-to-surface radiation are based on assumptions that are not appropriate for p-version finite elements. Two methods for determining internal radiation heat transfer are developed for one and two-dimensional p-version finite elements. In the first method, higher-order elements are divided into a number of sub-elements. Traditional methods are used to determine radiation heat flux along each sub-element and then mapped back to the parent element. In the second method, the radiation heat transfer equations are numerically integrated over the higher-order element. Comparisons with analytical solutions show that the integration scheme is generally more accurate than the sub-element method. Comparison to results from traditional finite elements shows that significant reduction in the number of elements in the mesh is possible using higher-order (p-version) finite elements.

A Fast and Accurate Method of Radiation Hydrodynamics Calculation in Spherical Symmetry

NASA Astrophysics Data System (ADS)

Stamer, Torsten; Inutsuka, Shu-ichiro

2018-06-01

We develop a new numerical scheme for solving the radiative transfer equation in a spherically symmetric system. This scheme does not rely on any kind of diffusion approximation, and it is accurate for optically thin, thick, and intermediate systems. In the limit of a homogeneously distributed extinction coefficient, our method is very accurate and exceptionally fast. We combine this fast method with a slower but more generally applicable method to describe realistic problems. We perform various test calculations, including a simplified protostellar collapse simulation. We also discuss possible future improvements.

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

NASA Astrophysics Data System (ADS)

Löwe, H.; Helbig, N.

2012-04-01

We provide a new quasi-analytical method to compute the topographic influence on the effective albedo of complex topography as required for meteorological, land-surface or climate models. We investigate radiative transfer in complex terrain via the radiosity equation on isotropic Gaussian random fields. Under controlled approximations we derive expressions for domain averages of direct, diffuse and terrain radiation and the sky view factor. Domain averaged quantities are related to a type of level-crossing probability of the random field which is approximated by longstanding results developed for acoustic scattering at ocean boundaries. This allows us to express all non-local horizon effects in terms of a local terrain parameter, namely the mean squared slope. Emerging integrals are computed numerically and fit formulas are given for practical purposes. As an implication of our approach we provide an expression for the effective albedo of complex terrain in terms of the sun elevation angle, mean squared slope, the area averaged surface albedo, and the direct-to-diffuse ratio of solar radiation. As an application, we compute the effective albedo for the Swiss Alps and discuss possible generalizations of the method.

Numerical modelling of soot formation and oxidation in laminar coflow non-smoking and smoking ethylene diffusion flames

NASA Astrophysics Data System (ADS)

Liu, Fengshan; Guo, Hongsheng; Smallwood, Gregory J.; Gülder, Ömer L.

2003-06-01

A numerical study of soot formation and oxidation in axisymmetric laminar coflow non-smoking and smoking ethylene diffusion flames was conducted using detailed gas-phase chemistry and complex thermal and transport properties. A modified two-equation soot model was employed to describe soot nucleation, growth and oxidation. Interaction between the gas-phase chemistry and soot chemistry was taken into account. Radiation heat transfer by both soot and radiating gases was calculated using the discrete-ordinates method coupled with a statistical narrow-band correlated-k based band model, and was used to evaluate the simple optically thin approximation. The governing equations in fully elliptic form were solved. The current models in the literature describing soot oxidation by O2 and OH have to be modified in order to predict the smoking flame. The modified soot oxidation model has only moderate effects on the calculation of the non-smoking flame, but dramatically affects the soot oxidation near the flame tip in the smoking flame. Numerical results of temperature, soot volume fraction and primary soot particle size and number density were compared with experimental data in the literature. Relatively good agreement was found between the prediction and the experimental data. The optically thin approximation radiation model significantly underpredicts temperatures in the upper portion of both flames, seriously affecting the soot prediction.

Solution Methods for Certain Evolution Equations

NASA Astrophysics Data System (ADS)

Vega-Guzman, Jose Manuel

Solution methods for certain linear and nonlinear evolution equations are presented in this dissertation. Emphasis is placed mainly on the analytical treatment of nonautonomous differential equations, which are challenging to solve despite the existent numerical and symbolic computational software programs available. Ideas from the transformation theory are adopted allowing one to solve the problems under consideration from a non-traditional perspective. First, the Cauchy initial value problem is considered for a class of nonautonomous and inhomogeneous linear diffusion-type equation on the entire real line. Explicit transformations are used to reduce the equations under study to their corresponding standard forms emphasizing on natural relations with certain Riccati(and/or Ermakov)-type systems. These relations give solvability results for the Cauchy problem of the parabolic equation considered. The superposition principle allows to solve formally this problem from an unconventional point of view. An eigenfunction expansion approach is also considered for this general evolution equation. Examples considered to corroborate the efficacy of the proposed solution methods include the Fokker-Planck equation, the Black-Scholes model and the one-factor Gaussian Hull-White model. The results obtained in the first part are used to solve the Cauchy initial value problem for certain inhomogeneous Burgers-type equation. The connection between linear (the Diffusion-type) and nonlinear (Burgers-type) parabolic equations is stress in order to establish a strong commutative relation. Traveling wave solutions of a nonautonomous Burgers equation are also investigated. Finally, it is constructed explicitly the minimum-uncertainty squeezed states for quantum harmonic oscillators. They are derived by the action of corresponding maximal kinematical invariance group on the standard ground state solution. It is shown that the product of the variances attains the required minimum value only at the instances that one variance is a minimum and the other is a maximum, when the squeezing of one of the variances occurs. Such explicit construction is possible due to the relation between the diffusion-type equation studied in the first part and the time-dependent Schrodinger equation. A modication of the radiation field operators for squeezed photons in a perfect cavity is also suggested with the help of a nonstandard solution of Heisenberg's equation of motion.

An asymptotic preserving unified gas kinetic scheme for frequency-dependent radiative transfer equations

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

Sun, Wenjun, E-mail: sun_wenjun@iapcm.ac.cn; Jiang, Song, E-mail: jiang@iapcm.ac.cn; Xu, Kun, E-mail: makxu@ust.hk

This paper presents an extension of previous work (Sun et al., 2015 [22]) of the unified gas kinetic scheme (UGKS) for the gray radiative transfer equations to the frequency-dependent (multi-group) radiative transfer system. Different from the gray radiative transfer equations, where the optical opacity is only a function of local material temperature, the simulation of frequency-dependent radiative transfer is associated with additional difficulties from the frequency-dependent opacity. For the multiple frequency radiation, the opacity depends on both the spatial location and the frequency. For example, the opacity is typically a decreasing function of frequency. At the same spatial region themore » transport physics can be optically thick for the low frequency photons, and optically thin for high frequency ones. Therefore, the optical thickness is not a simple function of space location. In this paper, the UGKS for frequency-dependent radiative system is developed. The UGKS is a finite volume method and the transport physics is modeled according to the ratio of the cell size to the photon's frequency-dependent mean free path. When the cell size is much larger than the photon's mean free path, a diffusion solution for such a frequency radiation will be obtained. On the other hand, when the cell size is much smaller than the photon's mean free path, a free transport mechanism will be recovered. In the regime between the above two limits, with the variation of the ratio between the local cell size and photon's mean free path, the UGKS provides a smooth transition in the physical and frequency space to capture the corresponding transport physics accurately. The seemingly straightforward extension of the UGKS from the gray to multiple frequency radiation system is due to its intrinsic consistent multiple scale transport modeling, but it still involves lots of work to properly discretize the multiple groups in order to design an asymptotic preserving (AP) scheme in all regimes. The current scheme is tested in a few frequency-dependent radiation problems, and the results are compared with the solutions from the well-defined implicit Monte Carlo (IMC) method. The UGKS is much more efficient than IMC, and the computational times of both schemes for all test cases are listed. The UGKS seems to be the first discrete ordinate method (DOM) for the accurate capturing of multiple frequency radiative transport physics from ballistic particle motion to the diffusive wave propagation.« less

Data fitting and image fine-tuning approach to solve the inverse problem in fluorescence molecular imaging

NASA Astrophysics Data System (ADS)

Gorpas, Dimitris; Politopoulos, Kostas; Yova, Dido; Andersson-Engels, Stefan

2008-02-01

One of the most challenging problems in medical imaging is to "see" a tumour embedded into tissue, which is a turbid medium, by using fluorescent probes for tumour labeling. This problem, despite the efforts made during the last years, has not been fully encountered yet, due to the non-linear nature of the inverse problem and the convergence failures of many optimization techniques. This paper describes a robust solution of the inverse problem, based on data fitting and image fine-tuning techniques. As a forward solver the coupled radiative transfer equation and diffusion approximation model is proposed and compromised via a finite element method, enhanced with adaptive multi-grids for faster and more accurate convergence. A database is constructed by application of the forward model on virtual tumours with known geometry, and thus fluorophore distribution, embedded into simulated tissues. The fitting procedure produces the best matching between the real and virtual data, and thus provides the initial estimation of the fluorophore distribution. Using this information, the coupled radiative transfer equation and diffusion approximation model has the required initial values for a computational reasonable and successful convergence during the image fine-tuning application.

Dynamic implicit 3D adaptive mesh refinement for non-equilibrium radiation diffusion

NASA Astrophysics Data System (ADS)

Philip, B.; Wang, Z.; Berrill, M. A.; Birke, M.; Pernice, M.

2014-04-01

The time dependent non-equilibrium radiation diffusion equations are important for solving the transport of energy through radiation in optically thick regimes and find applications in several fields including astrophysics and inertial confinement fusion. The associated initial boundary value problems that are encountered often exhibit a wide range of scales in space and time and are extremely challenging to solve. To efficiently and accurately simulate these systems we describe our research on combining techniques that will also find use more broadly for long term time integration of nonlinear multi-physics systems: implicit time integration for efficient long term time integration of stiff multi-physics systems, local control theory based step size control to minimize the required global number of time steps while controlling accuracy, dynamic 3D adaptive mesh refinement (AMR) to minimize memory and computational costs, Jacobian Free Newton-Krylov methods on AMR grids for efficient nonlinear solution, and optimal multilevel preconditioner components that provide level independent solver convergence.

Modelling of soot formation in laminar diffusion flames using a comprehensive CFD-PBE model with detailed gas-phase chemistry

NASA Astrophysics Data System (ADS)

Akridis, Petros; Rigopoulos, Stelios

2017-01-01

A discretised population balance equation (PBE) is coupled with an in-house computational fluid dynamics (CFD) code in order to model soot formation in laminar diffusion flames. The unsteady Navier-Stokes, species and enthalpy transport equations and the spatially-distributed discretised PBE for the soot particles are solved in a coupled manner, together with comprehensive gas-phase chemistry and an optically thin radiation model, thus yielding the complete particle size distribution of the soot particles. Nucleation, surface growth and oxidation are incorporated into the PBE using an acetylene-based soot model. The potential of the proposed methodology is investigated by comparing with experimental results from the Santoro jet burner [Santoro, Semerjian and Dobbins, Soot particle measurements in diffusion flames, Combustion and Flame, Vol. 51 (1983), pp. 203-218; Santoro, Yeh, Horvath and Semerjian, The transport and growth of soot particles in laminar diffusion flames, Combustion Science and Technology, Vol. 53 (1987), pp. 89-115] for three laminar axisymmetric non-premixed ethylene flames: a non-smoking, an incipient smoking and a smoking flame. Overall, good agreement is observed between the numerical and the experimental results.

Hybrid diffusion-P3 equation in N-layered turbid media: steady-state domain.

PubMed

Shi, Zhenzhi; Zhao, Huijuan; Xu, Kexin

2011-10-01

This paper discusses light propagation in N-layered turbid media. The hybrid diffusion-P3 equation is solved for an N-layered finite or infinite turbid medium in the steady-state domain for one point source using the extrapolated boundary condition. The Fourier transform formalism is applied to derive the analytical solutions of the fluence rate in Fourier space. Two inverse Fourier transform methods are developed to calculate the fluence rate in real space. In addition, the solutions of the hybrid diffusion-P3 equation are compared to the solutions of the diffusion equation and the Monte Carlo simulation. For the case of small absorption coefficients, the solutions of the N-layered diffusion equation and hybrid diffusion-P3 equation are almost equivalent and are in agreement with the Monte Carlo simulation. For the case of large absorption coefficients, the model of the hybrid diffusion-P3 equation is more precise than that of the diffusion equation. In conclusion, the model of the hybrid diffusion-P3 equation can replace the diffusion equation for modeling light propagation in the N-layered turbid media for a wide range of absorption coefficients.

Large-scale magnetic field in the accretion discs of young stars: the influence of magnetic diffusion, buoyancy and Hall effect

NASA Astrophysics Data System (ADS)

Khaibrakhmanov, S. A.; Dudorov, A. E.; Parfenov, S. Yu.; Sobolev, A. M.

2017-01-01

We investigate the fossil magnetic field in the accretion and protoplanetary discs using the Shakura and Sunyaev approach. The distinguishing feature of this study is the accurate solution of the ionization balance equations and the induction equation with Ohmic diffusion, magnetic ambipolar diffusion, buoyancy and the Hall effect. We consider the ionization by cosmic rays, X-rays and radionuclides, radiative recombinations, recombinations on dust grains and also thermal ionization. The buoyancy appears as the additional mechanism of magnetic flux escape in the steady-state solution of the induction equation. Calculations show that Ohmic diffusion and magnetic ambipolar diffusion constraint the generation of the magnetic field inside the `dead' zones. The magnetic field in these regions is quasi-vertical. The buoyancy constraints the toroidal magnetic field strength close to the disc inner edge. As a result, the toroidal and vertical magnetic fields become comparable. The Hall effect is important in the regions close to the borders of the `dead' zones because electrons are magnetized there. The magnetic field in these regions is quasi-radial. We calculate the magnetic field strength and geometry for the discs with accretion rates (10^{-8}-10^{-6}) {M}_{⊙} {yr}^{-1}. The fossil magnetic field geometry does not change significantly during the disc evolution while the accretion rate decreases. We construct the synthetic maps of dust emission polarized due to the dust grain alignment by the magnetic field. In the polarization maps, the `dead' zones appear as the regions with the reduced values of polarization degree in comparison to those in the adjacent regions.

The drag of airplane radiators with special reference to air heating : comparison of theory and experiment

NASA Technical Reports Server (NTRS)

Gothert, B

1939-01-01

This report contains a survey of past radiator research. This report also is intended as a systematic comparison of theoretical and experimental radiator drag, with the object of ascertaining the most important loss sources and their interaction in different cases of installation, and to separate the radiator systems which are amenable to calculation, both as regards axial flow and drag. The sources of loss due to the diffuser are to be looked into closely as in many cases they can be of preeminent magnitude and their customary appraisal, according to Fliegner's formula, does not meet actual conditions. Besides, generally applicable equations and charts are developed for the rapid determination of the heating effect of radiators as regards flow and drag, and then checked by routine tests on hot radiators.

Considerations for the independent reaction times and step-by-step methods for radiation chemistry simulations

NASA Astrophysics Data System (ADS)

Plante, Ianik; Devroye, Luc

2017-10-01

Ionizing radiation interacts with the water molecules of the tissues mostly by ionizations and excitations, which result in the formation of the radiation track structure and the creation of radiolytic species such as H.,.OH, H2, H2O2, and e-aq. After their creation, these species diffuse and may chemically react with the neighboring species and with the molecules of the medium. Therefore radiation chemistry is of great importance in radiation biology. As the chemical species are not distributed homogeneously, the use of conventional models of homogeneous reactions cannot completely describe the reaction kinetics of the particles. Actually, many simulations of radiation chemistry are done using the Independent Reaction Time (IRT) method, which is a very fast technique to calculate radiochemical yields but which do not calculate the positions of the radiolytic species as a function of time. Step-by-step (SBS) methods, which are able to provide such information, have been used only sparsely because these are time-consuming in terms of calculation. Recent improvements in computer performance now allow the regular use of the SBS method in radiation chemistry. The SBS and IRT methods are both based on the Green's functions of the diffusion equation (GFDE). In this paper, several sampling algorithms of the GFDE and for the IRT method are presented. We show that the IRT and SBS methods are exactly equivalent for 2-particles systems for diffusion and partially diffusion-controlled reactions between non-interacting particles. We also show that the results obtained with the SBS simulation method with periodic boundary conditions are in agreement with the predictions by classical reaction kinetics theory, which is an important step towards using this method for modelling of biochemical networks and metabolic pathways involved in oxidative stress. Finally, the first simulation results obtained with the code RITRACKS (Relativistic Ion Tracks) are presented.

Turbulent Flow past High Temperature Surfaces

NASA Astrophysics Data System (ADS)

Mehmedagic, Igbal; Thangam, Siva; Carlucci, Pasquale; Buckley, Liam; Carlucci, Donald

2014-11-01

Flow over high-temperature surfaces subject to wall heating is analyzed with applications to projectile design. In this study, computations are performed using an anisotropic Reynolds-stress model to study flow past surfaces that are subject to radiative flux. The model utilizes a phenomenological treatment of the energy spectrum and diffusivities of momentum and heat to include the effects of wall heat transfer and radiative exchange. The radiative transport is modeled using Eddington approximation including the weighted effect of nongrayness of the fluid. The time-averaged equations of motion and energy are solved using the modeled form of transport equations for the turbulence kinetic energy and the scalar form of turbulence dissipation with an efficient finite-volume algorithm. The model is applied for available test cases to validate its predictive capabilities for capturing the effects of wall heat transfer. Computational results are compared with experimental data available in the literature. Applications involving the design of projectiles are summarized. Funded in part by U.S. Army, ARDEC.

Radiation effects studies for the high-resolution spectrograph

NASA Technical Reports Server (NTRS)

Smith, L. C.; Becher, J.

1982-01-01

The generation and collection of charge carriers created during the passage of energetic protons through a silicon photodiode array are modeled. Pulse height distributions of noise charge collected during exposure of a digicon type diode array to 21 and 75 MeV protons were obtained. The magnitude of charge collected by a diode from each proton event is determined not only by diffusion, but by statistical considerations involving the ionization process itself. Utilizing analytical solutions to the diffusion equation for transport of minority carriers, together with the Vavilov theory of energy loss fluctuations in thin absorbers, simulations of the pulse height spectra which follow the experimental distributions fairly well are presented and an estimate for the minority carrier diffusion length L sub d is provided.

A Rate-Theory-Phase-Field Model of Irradiation-Induced Recrystallization in UMo Nuclear Fuels

NASA Astrophysics Data System (ADS)

Hu, Shenyang; Joshi, Vineet; Lavender, Curt A.

2017-12-01

In this work, we developed a recrystallization model to study the effect of microstructures and radiation conditions on recrystallization kinetics in UMo fuels. The model integrates the rate theory of intragranular gas bubble and interstitial loop evolutions and a phase-field model of recrystallization zone evolution. A first passage method is employed to describe one-dimensional diffusion of interstitials with a diffusivity value several orders of magnitude larger than that of fission gas xenons. With the model, the effect of grain sizes on recrystallization kinetics is simulated. The results show that (1) recrystallization in large grains starts earlier than that in small grains, (2) the recrystallization kinetics (recrystallization volume fraction) decrease as the grain size increases, (3) the predicted recrystallization kinetics are consistent with the experimental results, and (4) the recrystallization kinetics can be described by the modified Avrami equation, but the parameters of the Avrami equation strongly depend on the grain size.

Dynamics of the diffusive DM-DE interaction – Dynamical system approach

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

Haba, Zbigniew; Stachowski, Aleksander; Szydłowski, Marek, E-mail: zhab@ift.uni.wroc.pl, E-mail: aleksander.stachowski@uj.edu.pl, E-mail: marek.szydlowski@uj.edu.pl

We discuss dynamics of a model of an energy transfer between dark energy (DE) and dark matter (DM) . The energy transfer is determined by a non-conservation law resulting from a diffusion of dark matter in an environment of dark energy. The relativistic invariance defines the diffusion in a unique way. The system can contain baryonic matter and radiation which do not interact with the dark sector. We treat the Friedman equation and the conservation laws as a closed dynamical system. The dynamics of the model is examined using the dynamical systems methods for demonstration how solutions depend on initialmore » conditions. We also fit the model parameters using astronomical observation: SNIa, H ( z ), BAO and Alcock-Paczynski test. We show that the model with diffuse DM-DE is consistent with the data.« less

Effects of whistler mode hiss waves on the radiation belts structure during quiet times

NASA Astrophysics Data System (ADS)

Ripoll, J. F.; Santolik, O.; Reeves, G. D.; Kurth, W. S.; Denton, M.; Loridan, V.; Thaller, S. A.; Cunningham, G.; Kletzing, C.; Turner, D. L.; Henderson, M. G.; Ukhorskiy, S.; Drozdov, A.; Cervantes Villa, J. S.; Shprits, Y.

2017-12-01

We present dynamic Fokker-Planck simulations of the electron radiation belts and slot formation during the quiet days that can follow a storm. Simulations are made for all energies and L-shells between 2 and 6 in the view of recovering the observations of two particular events. Pitch angle diffusion is essential to energy structure of the belts and slot region. Pitch angle diffusion is computed from data-driven spatially and temporally-resolved whistler mode hiss wave and ambient plasma observations from the Van Allen Probes satellites. The simulations are performed either with a 3D formulation that uses pitch angle diffusion coefficients or with a simpler 1D Fokker-Planck equation based on losses computed from a lifetime. Validation is carried out globally against Magnetic Electron and Ion Spectrometer observations of the belts at all energy. Results are complemented with a sensitivity study involving different radial diffusion coefficients, electron lifetimes, and pitch angle diffusion coefficients. We discuss which models allow to recover the observed "S-shaped" energy-dependent inner boundary to the outer zone that results from the competition between diffusive radial transport and losses. Periods when the plasmasphere extends beyond L 5 favor long-lasting hiss losses from the outer belt. Through these simulations, we explain the full structure in energy and L-shell of the belts and the slot formation by hiss scattering during quiet storm recovery.

Diffuse Optical Tomography for Brain Imaging: Continuous Wave Instrumentation and Linear Analysis Methods

NASA Astrophysics Data System (ADS)

Giacometti, Paolo; Diamond, Solomon G.

Diffuse optical tomography (DOT) is a functional brain imaging technique that measures cerebral blood oxygenation and blood volume changes. This technique is particularly useful in human neuroimaging measurements because of the coupling between neural and hemodynamic activity in the brain. DOT is a multichannel imaging extension of near-infrared spectroscopy (NIRS). NIRS uses laser sources and light detectors on the scalp to obtain noninvasive hemodynamic measurements from spectroscopic analysis of the remitted light. This review explains how NIRS data analysis is performed using a combination of the modified Beer-Lambert law (MBLL) and the diffusion approximation to the radiative transport equation (RTE). Laser diodes, photodiode detectors, and optical terminals that contact the scalp are the main components in most NIRS systems. Placing multiple sources and detectors over the surface of the scalp allows for tomographic reconstructions that extend the individual measurements of NIRS into DOT. Mathematically arranging the DOT measurements into a linear system of equations that can be inverted provides a way to obtain tomographic reconstructions of hemodynamics in the brain.

A Radiation Transfer Solver for Athena Using Short Characteristics

NASA Astrophysics Data System (ADS)

Davis, Shane W.; Stone, James M.; Jiang, Yan-Fei

2012-03-01

We describe the implementation of a module for the Athena magnetohydrodynamics (MHD) code that solves the time-independent, multi-frequency radiative transfer (RT) equation on multidimensional Cartesian simulation domains, including scattering and non-local thermodynamic equilibrium (LTE) effects. The module is based on well known and well tested algorithms developed for modeling stellar atmospheres, including the method of short characteristics to solve the RT equation, accelerated Lambda iteration to handle scattering and non-LTE effects, and parallelization via domain decomposition. The module serves several purposes: it can be used to generate spectra and images, to compute a variable Eddington tensor (VET) for full radiation MHD simulations, and to calculate the heating and cooling source terms in the MHD equations in flows where radiation pressure is small compared with gas pressure. For the latter case, the module is combined with the standard MHD integrators using operator splitting: we describe this approach in detail, including a new constraint on the time step for stability due to radiation diffusion modes. Implementation of the VET method for radiation pressure dominated flows is described in a companion paper. We present results from a suite of test problems for both the RT solver itself and for dynamical problems that include radiative heating and cooling. These tests demonstrate that the radiative transfer solution is accurate and confirm that the operator split method is stable, convergent, and efficient for problems of interest. We demonstrate there is no need to adopt ad hoc assumptions of questionable accuracy to solve RT problems in concert with MHD: the computational cost for our general-purpose module for simple (e.g., LTE gray) problems can be comparable to or less than a single time step of Athena's MHD integrators, and only few times more expensive than that for more general (non-LTE) problems.

A unified radiative magnetohydrodynamics code for lightning-like discharge simulations

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

Chen, Qiang, E-mail: cq0405@126.com; Chen, Bin, E-mail: emcchen@163.com; Xiong, Run

2014-03-15

A two-dimensional Eulerian finite difference code is developed for solving the non-ideal magnetohydrodynamic (MHD) equations including the effects of self-consistent magnetic field, thermal conduction, resistivity, gravity, and radiation transfer, which when combined with specified pulse current models and plasma equations of state, can be used as a unified lightning return stroke solver. The differential equations are written in the covariant form in the cylindrical geometry and kept in the conservative form which enables some high-accuracy shock capturing schemes to be equipped in the lightning channel configuration naturally. In this code, the 5-order weighted essentially non-oscillatory scheme combined with Lax-Friedrichs fluxmore » splitting method is introduced for computing the convection terms of the MHD equations. The 3-order total variation diminishing Runge-Kutta integral operator is also equipped to keep the time-space accuracy of consistency. The numerical algorithms for non-ideal terms, e.g., artificial viscosity, resistivity, and thermal conduction, are introduced in the code via operator splitting method. This code assumes the radiation is in local thermodynamic equilibrium with plasma components and the flux limited diffusion algorithm with grey opacities is implemented for computing the radiation transfer. The transport coefficients and equation of state in this code are obtained from detailed particle population distribution calculation, which makes the numerical model is self-consistent. This code is systematically validated via the Sedov blast solutions and then used for lightning return stroke simulations with the peak current being 20 kA, 30 kA, and 40 kA, respectively. The results show that this numerical model consistent with observations and previous numerical results. The population distribution evolution and energy conservation problems are also discussed.« less

The Earth’s Radiation Belts.

DTIC Science & Technology

1983-09-20

on -- 1d. It --- . d id-er, c S, blck -.o~b.1) ’Trapped radiation Steady-state miodels Adiabatic invariants Empirical flux models Diffusion equations...Shell -splitting, Transport theory Nuclear detonations Wave-oarticle interactions Effects on microelectronics 20 ABSTRACT ( C -0- n OR e -~ d . It -~e-lay...olo -i t i os5 at 500 ke\\% live lrtI’m i s pt, eOI iS .1: litv, ,Ie It if 5)ht* stIweo f iul-’t, wi te thle hie av itk, i il - il v t’il 1 Ltt sI c a

Multiscale solutions of radiative heat transfer by the discrete unified gas kinetic scheme

NASA Astrophysics Data System (ADS)

Luo, Xiao-Ping; Wang, Cun-Hai; Zhang, Yong; Yi, Hong-Liang; Tan, He-Ping

2018-06-01

The radiative transfer equation (RTE) has two asymptotic regimes characterized by the optical thickness, namely, optically thin and optically thick regimes. In the optically thin regime, a ballistic or kinetic transport is dominant. In the optically thick regime, energy transport is totally dominated by multiple collisions between photons; that is, the photons propagate by means of diffusion. To obtain convergent solutions to the RTE, conventional numerical schemes have a strong dependence on the number of spatial grids, which leads to a serious computational inefficiency in the regime where the diffusion is predominant. In this work, a discrete unified gas kinetic scheme (DUGKS) is developed to predict radiative heat transfer in participating media. Numerical performances of the DUGKS are compared in detail with conventional methods through three cases including one-dimensional transient radiative heat transfer, two-dimensional steady radiative heat transfer, and three-dimensional multiscale radiative heat transfer. Due to the asymptotic preserving property, the present method with relatively coarse grids gives accurate and reliable numerical solutions for large, small, and in-between values of optical thickness, and, especially in the optically thick regime, the DUGKS demonstrates a pronounced computational efficiency advantage over the conventional numerical models. In addition, the DUGKS has a promising potential in the study of multiscale radiative heat transfer inside the participating medium with a transition from optically thin to optically thick regimes.

Fokker-Planck analysis of transverse collective instabilities in electron storage rings

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

Lindberg, Ryan R.

We analyze single bunch transverse instabilities due to wakefields using a Fokker-Planck model. We first expand on the work of T. Suzuki, Part. Accel. 12, 237 (1982) to derive the theoretical model including chromaticity, both dipolar and quadrupolar transverse wakefields, and the effects of damping and diffusion due to the synchrotron radiation. We reduce the problem to a linear matrix equation, whose eigenvalues and eigenvectors determine the collective stability of the beam. We then show that various predictions of the theory agree quite well with results from particle tracking simulations, including the threshold current for transverse instability and the profilemore » of the unstable mode. In particular, we find that predicting collective stability for high energy electron beams at moderate to large values of chromaticity requires the full Fokker-Planck analysis to properly account for the effects of damping and diffusion due to synchrotron radiation.« less

Cellular Automata for Spatiotemporal Pattern Formation from Reaction-Diffusion Partial Differential Equations

NASA Astrophysics Data System (ADS)

Ohmori, Shousuke; Yamazaki, Yoshihiro

2016-01-01

Ultradiscrete equations are derived from a set of reaction-diffusion partial differential equations, and cellular automaton rules are obtained on the basis of the ultradiscrete equations. Some rules reproduce the dynamical properties of the original reaction-diffusion equations, namely, bistability and pulse annihilation. Furthermore, other rules bring about soliton-like preservation and periodic pulse generation with a pacemaker, which are not obtained from the original reaction-diffusion equations.

RITRACKS: A Software for Simulation of Stochastic Radiation Track Structure, Micro and Nanodosimetry, Radiation Chemistry and DNA Damage for Heavy Ions

NASA Technical Reports Server (NTRS)

Plante, I; Wu, H

2014-01-01

The code RITRACKS (Relativistic Ion Tracks) has been developed over the last few years at the NASA Johnson Space Center to simulate the effects of ionizing radiations at the microscopic scale, to understand the effects of space radiation at the biological level. The fundamental part of this code is the stochastic simulation of radiation track structure of heavy ions, an important component of space radiations. The code can calculate many relevant quantities such as the radial dose, voxel dose, and may also be used to calculate the dose in spherical and cylindrical targets of various sizes. Recently, we have incorporated DNA structure and damage simulations at the molecular scale in RITRACKS. The direct effect of radiations is simulated by introducing a slight modification of the existing particle transport algorithms, using the Binary-Encounter-Bethe model of ionization cross sections for each molecular orbitals of DNA. The simulation of radiation chemistry is done by a step-by-step diffusion-reaction program based on the Green's functions of the diffusion equation]. This approach is also used to simulate the indirect effect of ionizing radiation on DNA. The software can be installed independently on PC and tablets using the Windows operating system and does not require any coding from the user. It includes a Graphic User Interface (GUI) and a 3D OpenGL visualization interface. The calculations are executed simultaneously (in parallel) on multiple CPUs. The main features of the software will be presented.

Hourly global and diffuse radiation of Lagos, Nigeria-correlation with some atmospheric parameters

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

Chendo, M.A.C.; Maduekwe, A.A.L.

1994-03-01

The influence of four climatic parameters on the hourly diffuse fraction in Lagos, Nigeria, has been studied. Using data for two years, new correlations were established. The standard error of the Liu and Jordan-type equation was reduced by 12.83% when solar elevation, ambient temperature, and relative humidity were used together as predictor variables for the entire data set. Ambient temperature and relative humidity proved to be very important variables for predicting the diffuse fraction of the solar radiation passing through the humid atmosphere of the coastal and tropic city of Lagos. Seasonal analysis carried out with the data showed improvementsmore » on the standard errors for the new seasonal correlations. In the case of the dry season, the improvement was 18.37%, whole for the wet season, this was 12.37%. Comparison with existing correlations showed that the performance of the one parameter model (namely K[sub t]), of Orgill and Hollands and Reindl, Beckman, and Duffie were very different from the Liu and Jordan-type model obtained for Lagos.« less

Photon migration through a turbid slab described by a model based on diffusion approximation. II. Comparison with Monte Carlo results.

PubMed

Martelli, F; Contini, D; Taddeucci, A; Zaccanti, G

1997-07-01

In our companion paper we presented a model to describe photon migration through a diffusing slab. The model, developed for a homogeneous slab, is based on the diffusion approximation and is able to take into account reflection at the boundaries resulting from the refractive index mismatch. In this paper the predictions of the model are compared with solutions of the radiative transfer equation obtained by Monte Carlo simulations in order to determine the applicability limits of the approximated theory in different physical conditions. A fitting procedure, carried out with the optical properties as fitting parameters, is used to check the application of the model to the inverse problem. The results show that significant errors can be made if the effect of the refractive index mismatch is not properly taken into account. Errors are more important when measurements of transmittance are used. The effects of using a receiver with a limited angular field of view and the angular distribution of the radiation that emerges from the slab have also been investigated.

Long-term evolution of electron distribution function due to nonlinear resonant interaction with whistler mode waves

NASA Astrophysics Data System (ADS)

Artemyev, Anton V.; Neishtadt, Anatoly I.; Vasiliev, Alexei A.

2018-04-01

Accurately modelling and forecasting of the dynamics of the Earth's radiation belts with the available computer resources represents an important challenge that still requires significant advances in the theoretical plasma physics field of wave-particle resonant interaction. Energetic electron acceleration or scattering into the Earth's atmosphere are essentially controlled by their resonances with electromagnetic whistler mode waves. The quasi-linear diffusion equation describes well this resonant interaction for low intensity waves. During the last decade, however, spacecraft observations in the radiation belts have revealed a large number of whistler mode waves with sufficiently high intensity to interact with electrons in the nonlinear regime. A kinetic equation including such nonlinear wave-particle interactions and describing the long-term evolution of the electron distribution is the focus of the present paper. Using the Hamiltonian theory of resonant phenomena, we describe individual electron resonance with an intense coherent whistler mode wave. The derived characteristics of such a resonance are incorporated into a generalized kinetic equation which includes non-local transport in energy space. This transport is produced by resonant electron trapping and nonlinear acceleration. We describe the methods allowing the construction of nonlinear resonant terms in the kinetic equation and discuss possible applications of this equation.

A New Numerical Scheme for Cosmic-Ray Transport

NASA Astrophysics Data System (ADS)

Jiang, Yan-Fei; Oh, S. Peng

2018-02-01

Numerical solutions of the cosmic-ray (CR) magnetohydrodynamic equations are dogged by a powerful numerical instability, which arises from the constraint that CRs can only stream down their gradient. The standard cure is to regularize by adding artificial diffusion. Besides introducing ad hoc smoothing, this has a significant negative impact on either computational cost or complexity and parallel scalings. We describe a new numerical algorithm for CR transport, with close parallels to two-moment methods for radiative transfer under the reduced speed of light approximation. It stably and robustly handles CR streaming without any artificial diffusion. It allows for both isotropic and field-aligned CR streaming and diffusion, with arbitrary streaming and diffusion coefficients. CR transport is handled explicitly, while source terms are handled implicitly. The overall time step scales linearly with resolution (even when computing CR diffusion) and has a perfect parallel scaling. It is given by the standard Courant condition with respect to a constant maximum velocity over the entire simulation domain. The computational cost is comparable to that of solving the ideal MHD equation. We demonstrate the accuracy and stability of this new scheme with a wide variety of tests, including anisotropic streaming and diffusion tests, CR-modified shocks, CR-driven blast waves, and CR transport in multiphase media. The new algorithm opens doors to much more ambitious and hitherto intractable calculations of CR physics in galaxies and galaxy clusters. It can also be applied to other physical processes with similar mathematical structure, such as saturated, anisotropic heat conduction.

Method for calculation of light field characteristics in optical diagnosis problems and personalized laser treatment of biological tissues

NASA Astrophysics Data System (ADS)

Lisenko, S. A.; Kugeiko, M. M.

2013-05-01

We have developed a simple method for solving the radiation transport equation, permitting us to rapidly calculate (with accuracy acceptable in practice) the diffuse reflection coeffi cient for a broad class of biological tissues in the spectral region of strong and weak absorption of light, and also the light flux distribution over the depth of the tissue. We show that it is feasible to use the proposed method for quantitative estimates of tissue parameters from its diffuse reflectance spectrum and also for selecting the irradiation dose which is optimal for a specifi c patient in laser therapy for various diseases.

Fluorescence molecular imaging based on the adjoint radiative transport equation

NASA Astrophysics Data System (ADS)

Asllanaj, Fatmir; Addoum, Ahmad; Rodolphe Roche, Jean

2018-07-01

A new reconstruction algorithm for fluorescence diffuse optical tomography of biological tissues is proposed. The radiative transport equation in the frequency domain is used to model light propagation. The adjoint method studied in this work provides an efficient way for solving the inverse problem. The methodology is applied to a 2D tissue-like phantom subjected to a collimated laser beam. Indocyanine Green is used as fluorophore. Reconstructed images of the spatial fluorophore absorption distribution is assessed taking into account the residual fluorescence in the medium. We show that illuminating the tissue surface from a collimated centered direction near the inclusion gaves a better reconstruction quality. Two closely positioned inclusions can be accurately localized. Additionally, their fluorophore absorption coefficients can be quantified. However, the algorithm failes to reconstruct smaller or deeper inclusions. This is due to light attenuation in the medium. Reconstructions with noisy data are also achieved with a reasonable accuracy.

3D Radiative Hydrodynamics Simulations of Protoplanetary Disks: A Comparison Between Two Radiative Cooling Algorithms

NASA Astrophysics Data System (ADS)

Lord, Jesse W.; Boley, A. C.; Durisen, R. H.

2006-12-01

We present a comparison between two three-dimensional radiative hydrodynamics simulations of a gravitationally unstable 0.07 Msun protoplanetary disk around a 0.5 Msun star. The first simulation is the radiatively cooled disk described in Boley et al. (2006, ApJ, 651). This simulation employed an algorithm that uses 3D flux-limited diffusion wherever the vertical Rosseland optical depth is greater than 2/3, which defines the optically thick region. The optically thin atmosphere of the disk, which cools according to its emissivity, is coupled to the optically thick region through an Eddington-like boundary condition. The second simulation employed an algorithm that uses a combination of solving the radiative transfer equation along rays in the z direction and flux limited diffusion in the r and phi directions on a cylindrical grid. We compare the following characteristics of the disk simulations: the mass transport and torques induced by gravitational instabilities, the effective temperature profiles of the disks, the gravitational and Reynolds stresses measured in the disk and those expected in an alpha-disk, and the amplitudes of the Fourier modes. This work has been supported by the National Science Foundation through grant AST-0452975 (astronomy REU program to Indiana University).

Radiative and Convective Heat Transfer over Ablating Composite Flat Surface in Hypersonic Flow Regime.

DTIC Science & Technology

1987-04-22

absorptivity in the presence of scatteringsc B Defined in equation (40) B wBE Diffuse surface radiosity C Mass fraction of injected species D. jiCoefficient of...Then 20 A eb)x 8 eb- (49) where B and B., are the surface radiosities . It follows invnediately that wX 0 T to d 2e (50) ~ f ~ b W 2 L 3 ( ) 2 1 - 1

Transformed Fourier and Fick equations for the control of heat and mass diffusion

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

Guenneau, S.; Petiteau, D.; Zerrad, M.

We review recent advances in the control of diffusion processes in thermodynamics and life sciences through geometric transforms in the Fourier and Fick equations, which govern heat and mass diffusion, respectively. We propose to further encompass transport properties in the transformed equations, whereby the temperature is governed by a three-dimensional, time-dependent, anisotropic heterogeneous convection-diffusion equation, which is a parabolic partial differential equation combining the diffusion equation and the advection equation. We perform two dimensional finite element computations for cloaks, concentrators and rotators of a complex shape in the transient regime. We precise that in contrast to invisibility cloaks for waves,more » the temperature (or mass concentration) inside a diffusion cloak crucially depends upon time, its distance from the source, and the diffusivity of the invisibility region. However, heat (or mass) diffusion outside cloaks, concentrators and rotators is unaffected by their presence, whatever their shape or position. Finally, we propose simplified designs of layered cylindrical and spherical diffusion cloaks that might foster experimental efforts in thermal and biochemical metamaterials.« less

Simulations of the stratocumulus-topped boundary layer with a third-order closure model

NASA Technical Reports Server (NTRS)

Moeng, C. H.; Randall, D. A.

1984-01-01

A third order closure model is proposed by Andre et al. (1982), in which the time rate of change terms, the relaxation and rapid effects for the pressure related terms, and the clipping approximation are included along with the quasi-normal closure, to study turbulence in a cloudy layer which is cooled radiatively from above. A spurious oscillation which is strongest near the inversion occurs. An analysis of the problem shows that the oscillation arises from the mean gradient and buoyancy terms of the triple moment equations; these terms are largest near the cloud top. The oscillation is physical, rather than computational. In nature the oscillation is effectively damped, by a mechanism which apparently is not included in our model. In the stably stratified layer just above the mixed layer top, turbulence can excite gravity waves, whose energy is radiated away. Because the closure assumption for the pressure terms does not take into account the transport of wave energy, the model generates spurious oscillations. Damping of the oscillations is possible by introducing diffusion terms into the triple moment equations. With a large enough choice for the diffusion coefficient, the oscillation is effectively eliminated. The results are quite sensitive to the ad hoc eddy coefficient.

Fractional Diffusion Processes: Probability Distributions and Continuous Time Random Walk

NASA Astrophysics Data System (ADS)

Gorenflo, R.; Mainardi, F.

A physical-mathematical approach to anomalous diffusion may be based on generalized diffusion equations (containing derivatives of fractional order in space or/and time) and related random walk models. By the space-time fractional diffusion equation we mean an evolution equation obtained from the standard linear diffusion equation by replacing the second-order space derivative with a Riesz-Feller derivative of order alpha in (0,2] and skewness theta (\\verttheta\\vertlemin \\{alpha ,2-alpha \\}), and the first-order time derivative with a Caputo derivative of order beta in (0,1] . The fundamental solution (for the Cauchy problem) of the fractional diffusion equation can be interpreted as a probability density evolving in time of a peculiar self-similar stochastic process. We view it as a generalized diffusion process that we call fractional diffusion process, and present an integral representation of the fundamental solution. A more general approach to anomalous diffusion is however known to be provided by the master equation for a continuous time random walk (CTRW). We show how this equation reduces to our fractional diffusion equation by a properly scaled passage to the limit of compressed waiting times and jump widths. Finally, we describe a method of simulation and display (via graphics) results of a few numerical case studies.

Multiphysics modelling of the separation of suspended particles via frequency ramping of ultrasonic standing waves.

PubMed

Trujillo, Francisco J; Eberhardt, Sebastian; Möller, Dirk; Dual, Jurg; Knoerzer, Kai

2013-03-01

A model was developed to determine the local changes of concentration of particles and the formations of bands induced by a standing acoustic wave field subjected to a sawtooth frequency ramping pattern. The mass transport equation was modified to incorporate the effect of acoustic forces on the concentration of particles. This was achieved by balancing the forces acting on particles. The frequency ramping was implemented as a parametric sweep for the time harmonic frequency response in time steps of 0.1s. The physics phenomena of piezoelectricity, acoustic fields and diffusion of particles were coupled and solved in COMSOL Multiphysics™ (COMSOL AB, Stockholm, Sweden) following a three step approach. The first step solves the governing partial differential equations describing the acoustic field by assuming that the pressure field achieves a pseudo steady state. In the second step, the acoustic radiation force is calculated from the pressure field. The final step allows calculating the locally changing concentration of particles as a function of time by solving the modified equation of particle transport. The diffusivity was calculated as function of concentration following the Garg and Ruthven equation which describes the steep increase of diffusivity when the concentration approaches saturation. However, it was found that this steep increase creates numerical instabilities at high voltages (in the piezoelectricity equations) and high initial particle concentration. The model was simplified to a pseudo one-dimensional case due to computation power limitations. The predicted particle distribution calculated with the model is in good agreement with the experimental data as it follows accurately the movement of the bands in the centre of the chamber. Crown Copyright © 2012. Published by Elsevier B.V. All rights reserved.

Anisotropic diffusion in mesh-free numerical magnetohydrodynamics

NASA Astrophysics Data System (ADS)

Hopkins, Philip F.

2017-04-01

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

Some Properties of the Fractional Equation of Continuity and the Fractional Diffusion Equation

NASA Astrophysics Data System (ADS)

Fukunaga, Masataka

2006-05-01

The fractional equation of continuity (FEC) and the fractional diffusion equation (FDE) show peculiar behaviors that are in the opposite sense to those expected from the equation of continuity and the diffusion equation, respectively. The behaviors are interpreted in terms of the memory effect of the fractional time derivatives included in the equations. Some examples are given by solutions of the FDE.

Solution of a modified fractional diffusion equation

NASA Astrophysics Data System (ADS)

Langlands, T. A. M.

2006-07-01

Recently, a modified fractional diffusion equation has been proposed [I. Sokolov, J. Klafter, From diffusion to anomalous diffusion: a century after Einstein's brownian motion, Chaos 15 (2005) 026103; A.V. Chechkin, R. Gorenflo, I.M. Sokolov, V.Yu. Gonchar, Distributed order time fractional diffusion equation, Frac. Calc. Appl. Anal. 6 (3) (2003) 259279; I.M. Sokolov, A.V. Checkin, J. Klafter, Distributed-order fractional kinetics, Acta. Phys. Pol. B 35 (2004) 1323.] for describing processes that become less anomalous as time progresses by the inclusion of a second fractional time derivative acting on the diffusion term. In this letter we give the solution of the modified equation on an infinite domain. In contrast to the solution of the traditional fractional diffusion equation, the solution of the modified equation requires an infinite series of Fox functions instead of a single Fox function.

Re-ionization and decaying dark matter

NASA Technical Reports Server (NTRS)

Dodelson, Scott; Jubas, Jay M.

1991-01-01

Gunn-Peterson tests suggest that the Universe was reionized after the standard recombination epoch. A systematic treatment is presented of the ionization process by deriving the Boltzmann equations appropriate to this regime. A compact solution for the photon spectrum is found in terms of the ionization ratio. These equations are then solved numerically for the Decaying Dark Matter scenario, wherein neutrinos with mass of order 30 eV radiatively decay producing photons which ionize the intergalactic medium. It was found that the neutrino mass and lifetime are severely constrained by Gunn-Peterson tests, observations of the diffuse photon spectrum in the ultraviolet regime, and the Hubble parameter.

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

Prinja, A. K.

The Karhunen-Loeve stochastic spectral expansion of a random binary mixture of immiscible fluids in planar geometry is used to explore asymptotic limits of radiation transport in such mixtures. Under appropriate scalings of mixing parameters - correlation length, volume fraction, and material cross sections - and employing multiple- scale expansion of the angular flux, previously established atomic mix and diffusion limits are reproduced. When applied to highly contrasting material properties in the small cor- relation length limit, the methodology yields a nonstandard reflective medium transport equation that merits further investigation. Finally, a hybrid closure is proposed that produces both small andmore » large correlation length limits of the closure condition for the material averaged equations.« less

Fully implicit moving mesh adaptive algorithm

NASA Astrophysics Data System (ADS)

Chacon, Luis

2005-10-01

In many problems of interest, the numerical modeler is faced with the challenge of dealing with multiple time and length scales. The former is best dealt with with fully implicit methods, which are able to step over fast frequencies to resolve the dynamical time scale of interest. The latter requires grid adaptivity for efficiency. Moving-mesh grid adaptive methods are attractive because they can be designed to minimize the numerical error for a given resolution. However, the required grid governing equations are typically very nonlinear and stiff, and of considerably difficult numerical treatment. Not surprisingly, fully coupled, implicit approaches where the grid and the physics equations are solved simultaneously are rare in the literature, and circumscribed to 1D geometries. In this study, we present a fully implicit algorithm for moving mesh methods that is feasible for multidimensional geometries. A crucial element is the development of an effective multilevel treatment of the grid equation.ootnotetextL. Chac'on, G. Lapenta, A fully implicit, nonlinear adaptive grid strategy, J. Comput. Phys., accepted (2005) We will show that such an approach is competitive vs. uniform grids both from the accuracy (due to adaptivity) and the efficiency standpoints. Results for a variety of models 1D and 2D geometries, including nonlinear diffusion, radiation-diffusion, Burgers equation, and gas dynamics will be presented.

Groups in the radiative transfer theory

NASA Astrophysics Data System (ADS)

Nikoghossian, Arthur

2016-11-01

The paper presents a group-theoretical description of radiation transfer in inhomogeneous and multi-component atmospheres with the plane-parallel geometry. It summarizes and generalizes the results obtained recently by the author for some standard transfer problems of astrophysical interest with allowance of the angle and frequency distributions of the radiation field. We introduce the concept of composition groups for media with different optical and physical properties. Group representations are derived for two possible cases of illumination of a composite finite atmosphere. An algorithm for determining the reflectance and transmittance of inhomogeneous and multi-component atmospheres is described. The group theory is applied also to determining the field of radiation inside an inhomogeneous atmosphere. The concept of a group of optical depth translations is introduced. The developed theory is illustrated with the problem of radiation diffusion with partial frequency distribution assuming that the inhomogeneity is due to depth-variation of the scattering coefficient. It is shown that once reflectance and transmittance of a medium are determined, the internal field of radiation in the source-free atmosphere is found without solving any new equations. The transfer problems for a semi-infinite atmosphere and an atmosphere with internal sources of energy are discussed. The developed theory allows to derive summation laws for the mean number of scattering events underwent by the photons in the course of diffusion in the atmosphere.

Preliminary study: Moisture-polymer interaction. Stuby objectives

NASA Technical Reports Server (NTRS)

Wen, L. C.

1985-01-01

The problems associated with mathematically modeling water-module interaction phenomena, including sorption and desorption, diffusion, and permeation are discussed. With reliable analytical models, an extensive materials data base, and solar radiation surface meteorological observations (SOLMET) weather data, predicting module lifetimes in realistic environments can become a practical reality. The status of the present techniques of simulating the various transport mechanisms was reported. The Dent model (a modified Brunauer-Emmet-Teller) approach represented polyvinyl butyral (PVB) sorption data. A 100-layer material model and Fick's diffusion model gave diffusivity values exhibiting adequate agreement with those measured for PVB. Diffusivity of PVB is concentration dependent, decreasing as the water content in PVB increases. The temperature dependence of diffusion in PVB is well modeled by the Arrhenius rate equation. Equilibrium conductivity and leakage current data are well represented by Hearle's model for bulk ionic conductivity. A nodal network analysis using the Systems Improved Numerical Differencing Analyzer (SINDA) Thermal Analyzer gave reasonable correlation with measurable data. It is concluded that realistic lifetime predictions seem to be feasible.

Two-dimensional Radiative Magnetohydrodynamic Simulations of Partial Ionization in the Chromosphere. II. Dynamics and Energetics of the Low Solar Atmosphere

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

Martínez-Sykora, Juan; Pontieu, Bart De; Hansteen, Viggo H.

2017-09-20

We investigate the effects of interactions between ions and neutrals on the chromosphere and overlying corona using 2.5D radiative MHD simulations with the Bifrost code. We have extended the code capabilities implementing ion–neutral interaction effects using the generalized Ohm’s law, i.e., we include the Hall term and the ambipolar diffusion (Pedersen dissipation) in the induction equation. Our models span from the upper convection zone to the corona, with the photosphere, chromosphere, and transition region partially ionized. Our simulations reveal that the interactions between ionized particles and neutral particles have important consequences for the magnetothermodynamics of these modeled layers: (1) ambipolarmore » diffusion increases the temperature in the chromosphere; (2) sporadically the horizontal magnetic field in the photosphere is diffused into the chromosphere, due to the large ambipolar diffusion; (3) ambipolar diffusion concentrates electrical currents, leading to more violent jets and reconnection processes, resulting in (3a) the formation of longer and faster spicules, (3b) heating of plasma during the spicule evolution, and (3c) decoupling of the plasma and magnetic field in spicules. Our results indicate that ambipolar diffusion is a critical ingredient for understanding the magnetothermodynamic properties in the chromosphere and transition region. The numerical simulations have been made publicly available, similar to previous Bifrost simulations. This will allow the community to study realistic numerical simulations with a wider range of magnetic field configurations and physics modules than previously possible.« less

Two-dimensional Radiative Magnetohydrodynamic Simulations of Partial Ionization in the Chromosphere. II. Dynamics and Energetics of the Low Solar Atmosphere

NASA Astrophysics Data System (ADS)

Martínez-Sykora, Juan; De Pontieu, Bart; Carlsson, Mats; Hansteen, Viggo H.; Nóbrega-Siverio, Daniel; Gudiksen, Boris V.

2017-09-01

We investigate the effects of interactions between ions and neutrals on the chromosphere and overlying corona using 2.5D radiative MHD simulations with the Bifrost code. We have extended the code capabilities implementing ion-neutral interaction effects using the generalized Ohm’s law, I.e., we include the Hall term and the ambipolar diffusion (Pedersen dissipation) in the induction equation. Our models span from the upper convection zone to the corona, with the photosphere, chromosphere, and transition region partially ionized. Our simulations reveal that the interactions between ionized particles and neutral particles have important consequences for the magnetothermodynamics of these modeled layers: (1) ambipolar diffusion increases the temperature in the chromosphere; (2) sporadically the horizontal magnetic field in the photosphere is diffused into the chromosphere, due to the large ambipolar diffusion; (3) ambipolar diffusion concentrates electrical currents, leading to more violent jets and reconnection processes, resulting in (3a) the formation of longer and faster spicules, (3b) heating of plasma during the spicule evolution, and (3c) decoupling of the plasma and magnetic field in spicules. Our results indicate that ambipolar diffusion is a critical ingredient for understanding the magnetothermodynamic properties in the chromosphere and transition region. The numerical simulations have been made publicly available, similar to previous Bifrost simulations. This will allow the community to study realistic numerical simulations with a wider range of magnetic field configurations and physics modules than previously possible.

Radiative entropy generation in a gray absorbing, emitting, and scattering planar medium at radiative equilibrium

NASA Astrophysics Data System (ADS)

Sadeghi, Pegah; Safavinejad, Ali

2017-11-01

Radiative entropy generation through a gray absorbing, emitting, and scattering planar medium at radiative equilibrium with diffuse-gray walls is investigated. The radiative transfer equation and radiative entropy generation equations are solved using discrete ordinates method. Components of the radiative entropy generation are considered for two different boundary conditions: two walls are at a prescribed temperature and mixed boundary conditions, which one wall is at a prescribed temperature and the other is at a prescribed heat flux. The effect of wall emissivities, optical thickness, single scattering albedo, and anisotropic-scattering factor on the entropy generation is attentively investigated. The results reveal that entropy generation in the system mainly arises from irreversible radiative transfer at wall with lower temperature. Total entropy generation rate for the system with prescribed temperature at walls remarkably increases as wall emissivity increases; conversely, for system with mixed boundary conditions, total entropy generation rate slightly decreases. Furthermore, as the optical thickness increases, total entropy generation rate remarkably decreases for the system with prescribed temperature at walls; nevertheless, for the system with mixed boundary conditions, total entropy generation rate increases. The variation of single scattering albedo does not considerably affect total entropy generation rate. This parametric analysis demonstrates that the optical thickness and wall emissivities have a significant effect on the entropy generation in the system at radiative equilibrium. Considering the parameters affecting radiative entropy generation significantly, provides an opportunity to optimally design or increase overall performance and efficiency by applying entropy minimization techniques for the systems at radiative equilibrium.

DREAM3D simulations of inner-belt dynamics

NASA Astrophysics Data System (ADS)

Cunningham, G.

2015-12-01

A 1973 paper by Lyons and Thorne explains the two-belt structure for electrons in the inner magnetosphere as a balance between inward radial diffusion and loss to the atmosphere due to pitch-angle scattering from Coulomb and VLF wave-particle interactions. In this paper, equilibrium solutions to a set of 1D radial diffusion equations, one for each value of the first invariant of motion, μ, were computed to produce the equilibrium structure. Each diffusion equation incorporated an L- and μ-dependent `lifetime' due to the Coulomb and wave-particle interactions. This model is appropriate under the assumption that radial diffusion is slow in comparison to pitch-angle scattering, and that there is no acceleration caused by the VLF wave-particle interactions. We have revisited this model using our DREAM3D 3D diffusion code, which allows the user to explicitly model the diffusion in pitch-angle and momentum rather than using a lifetime. We find that a) replacing the lifetimes with an explicit model of pitch-angle diffusion, thus allowing for coupling between radial and pitch-angle diffusion, affects the equilibrium structure, and b) over the long time scales needed to reach equilibrium, significant acceleration due to VLF wave particle interactions takes place due to the 'cross-terms' in pitch-angle and momentum and the sharp gradient in the equilibrium pitch-angle distributions. We also find that the equilibrium solutions are quite sensitive to various aspects of the physics model employed in the 1973 paper that can be improved, suggesting that additional work needs to be done to fully understand the equilibirum nature of the trapped electron radiation belts.

A reconsideration of the noise equivalent power and the data analysis procedure for the infrared imaging video bolometers

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

Pandya, Shwetang N., E-mail: pandya.shwetang@LHD.nifs.ac.jp; Sano, Ryuichi; Peterson, Byron J.

The infrared imaging video bolometer (IRVB) used for measurement of the two-dimensional (2D) radiation profiles from the Large Helical Device has been significantly upgraded recently to improve its signal to noise ratio, sensitivity, and calibration, which ultimately provides quantitative measurements of the radiation from the plasma. The reliability of the quantified data needs to be established by various checks. The noise estimates also need to be revised and more realistic values need to be established. It is shown that the 2D heat diffusion equation can be used for estimating the power falling on the IRVB foil, even with a significantmore » amount of spatial variation in the thermal diffusivity across the area of the platinum foil found experimentally during foil calibration. The equation for the noise equivalent power density (NEPD) is re-derived to include the errors in the measurement of the thermophysical and the optical properties of the IRVB foil. The theoretical value estimated using this newly derived equation matches closely, within 5.5%, with the mean experimental value. The change in the contribution of each error term of the NEPD equation with rising foil temperature is also studied and the blackbody term is found to dominate the other terms at elevated operating temperatures. The IRVB foil is also sensitive to the charge exchange (CX) neutrals escaping from the plasma. The CX neutral contribution is estimated to be marginally higher than the noise equivalent power (NEP) of the IRVB. It is also established that the radiation measured by the IRVB originates from the impurity line radiation from the plasma and not from the heated divertor tiles. The change in the power density due to noise reduction measures such as data smoothing and averaging is found to be comparable to the IRVB NEPD. The precautions that need to be considered during background subtraction are also discussed with experimental illustrations. Finally, the analysis algorithm with all the improvements is validated and found to reproduce the input power well within 10% accuracy. This article answers many fundamental questions relevant to the IRVB and illustrates the care to be exercised while processing the IRVB data.« less

Influence of a Simple Heat Loss Profile on a Pure Diffusion Flame

NASA Technical Reports Server (NTRS)

Ray, Anjan; Wichman, Indrek S.

1996-01-01

The presence of soot on the fuel side of a diffusion flame results in significant radiative heat losses. The influence of a fuel side heat loss zone on a pure diffusion flame established between a fuel and an oxidizer wall is investigated by assuming a hypothetical sech(sup 2) heat loss profile. The intensity and width of the loss zone are parametrically varied. The loss zone is placed at different distances from the Burke-Schumann flame location. The migration of the temperature and reactivity peaks are examined for a variety of situations. For certain cases the reaction zone breaks through the loss zone and relocates itself on the fuel side of the loss zone. In all cases the temperature and reactivity peaks move toward the fuel side with increased heat losses. The flame structure reveals that the primary balance for the energy equation is between the reaction term and the diffusion term. Extinction plots are generated for a variety of situations. The heat transfer from the flame to the walls and the radiative fraction is also investigated, and an analytical correlation formula, derived in a previous study, is shown to produce excellent predictions of our numerical results when an O(l) numerical multiplicative constant is employed.

Diffusion Coefficients from Molecular Dynamics Simulations in Binary and Ternary Mixtures

NASA Astrophysics Data System (ADS)

Liu, Xin; Schnell, Sondre K.; Simon, Jean-Marc; Krüger, Peter; Bedeaux, Dick; Kjelstrup, Signe; Bardow, André; Vlugt, Thijs J. H.

2013-07-01

Multicomponent diffusion in liquids is ubiquitous in (bio)chemical processes. It has gained considerable and increasing interest as it is often the rate limiting step in a process. In this paper, we review methods for calculating diffusion coefficients from molecular simulation and predictive engineering models. The main achievements of our research during the past years can be summarized as follows: (1) we introduced a consistent method for computing Fick diffusion coefficients using equilibrium molecular dynamics simulations; (2) we developed a multicomponent Darken equation for the description of the concentration dependence of Maxwell-Stefan diffusivities. In the case of infinite dilution, the multicomponent Darken equation provides an expression for [InlineEquation not available: see fulltext.] which can be used to parametrize the generalized Vignes equation; and (3) a predictive model for self-diffusivities was proposed for the parametrization of the multicomponent Darken equation. This equation accurately describes the concentration dependence of self-diffusivities in weakly associating systems. With these methods, a sound framework for the prediction of mutual diffusion in liquids is achieved.

Systematic derivation of reaction-diffusion equations with distributed delays and relations to fractional reaction-diffusion equations and hyperbolic transport equations: application to the theory of Neolithic transition.

PubMed

Vlad, Marcel Ovidiu; Ross, John

2002-12-01

We introduce a general method for the systematic derivation of nonlinear reaction-diffusion equations with distributed delays. We study the interactions among different types of moving individuals (atoms, molecules, quasiparticles, biological organisms, etc). The motion of each species is described by the continuous time random walk theory, analyzed in the literature for transport problems, whereas the interactions among the species are described by a set of transformation rates, which are nonlinear functions of the local concentrations of the different types of individuals. We use the time interval between two jumps (the transition time) as an additional state variable and obtain a set of evolution equations, which are local in time. In order to make a connection with the transport models used in the literature, we make transformations which eliminate the transition time and derive a set of nonlocal equations which are nonlinear generalizations of the so-called generalized master equations. The method leads under different specified conditions to various types of nonlocal transport equations including a nonlinear generalization of fractional diffusion equations, hyperbolic reaction-diffusion equations, and delay-differential reaction-diffusion equations. Thus in the analysis of a given problem we can fit to the data the type of reaction-diffusion equation and the corresponding physical and kinetic parameters. The method is illustrated, as a test case, by the study of the neolithic transition. We introduce a set of assumptions which makes it possible to describe the transition from hunting and gathering to agriculture economics by a differential delay reaction-diffusion equation for the population density. We derive a delay evolution equation for the rate of advance of agriculture, which illustrates an application of our analysis.

Diffuse correlation tomography in the transport regime: A theoretical study of the sensitivity to Brownian motion.

PubMed

Tricoli, Ugo; Macdonald, Callum M; Durduran, Turgut; Da Silva, Anabela; Markel, Vadim A

2018-02-01

Diffuse correlation tomography (DCT) uses the electric-field temporal autocorrelation function to measure the mean-square displacement of light-scattering particles in a turbid medium over a given exposure time. The movement of blood particles is here estimated through a Brownian-motion-like model in contrast to ordered motion as in blood flow. The sensitivity kernel relating the measurable field correlation function to the mean-square displacement of the particles can be derived by applying a perturbative analysis to the correlation transport equation (CTE). We derive an analytical expression for the CTE sensitivity kernel in terms of the Green's function of the radiative transport equation, which describes the propagation of the intensity. We then evaluate the kernel numerically. The simulations demonstrate that, in the transport regime, the sensitivity kernel provides sharper spatial information about the medium as compared with the correlation diffusion approximation. Also, the use of the CTE allows one to explore some additional degrees of freedom in the data such as the collimation direction of sources and detectors. Our results can be used to improve the spatial resolution of DCT, in particular, with applications to blood flow imaging in regions where the Brownian motion is dominant.

Diffuse correlation tomography in the transport regime: A theoretical study of the sensitivity to Brownian motion

NASA Astrophysics Data System (ADS)

Tricoli, Ugo; Macdonald, Callum M.; Durduran, Turgut; Da Silva, Anabela; Markel, Vadim A.

2018-02-01

Diffuse correlation tomography (DCT) uses the electric-field temporal autocorrelation function to measure the mean-square displacement of light-scattering particles in a turbid medium over a given exposure time. The movement of blood particles is here estimated through a Brownian-motion-like model in contrast to ordered motion as in blood flow. The sensitivity kernel relating the measurable field correlation function to the mean-square displacement of the particles can be derived by applying a perturbative analysis to the correlation transport equation (CTE). We derive an analytical expression for the CTE sensitivity kernel in terms of the Green's function of the radiative transport equation, which describes the propagation of the intensity. We then evaluate the kernel numerically. The simulations demonstrate that, in the transport regime, the sensitivity kernel provides sharper spatial information about the medium as compared with the correlation diffusion approximation. Also, the use of the CTE allows one to explore some additional degrees of freedom in the data such as the collimation direction of sources and detectors. Our results can be used to improve the spatial resolution of DCT, in particular, with applications to blood flow imaging in regions where the Brownian motion is dominant.

Light radiation pressure upon an optically orthotropic surface

NASA Astrophysics Data System (ADS)

Nerovny, Nikolay A.; Lapina, Irina E.; Grigorjev, Anton S.

2017-11-01

In this paper, we discuss the problem of determination of light radiation pressure force upon an anisotropic surface. The optical parameters of such a surface are considered to have major and minor axes, so the model is called an orthotropic model. We derive the equations for force components from emission, absorption, and reflection, utilizing a modified Maxwell's specular-diffuse model. The proposed model can be used to model a flat solar sail with wrinkles. By performing Bayesian analysis for example of a wrinkled surface, we show that there are cases in which an orthotropic model of the optical parameters of a surface may be more accurate than an isotropic model.

The role of fractional time-derivative operators on anomalous diffusion

NASA Astrophysics Data System (ADS)

Tateishi, Angel A.; Ribeiro, Haroldo V.; Lenzi, Ervin K.

2017-10-01

The generalized diffusion equations with fractional order derivatives have shown be quite efficient to describe the diffusion in complex systems, with the advantage of producing exact expressions for the underlying diffusive properties. Recently, researchers have proposed different fractional-time operators (namely: the Caputo-Fabrizio and Atangana-Baleanu) which, differently from the well-known Riemann-Liouville operator, are defined by non-singular memory kernels. Here we proposed to use these new operators to generalize the usual diffusion equation. By analyzing the corresponding fractional diffusion equations within the continuous time random walk framework, we obtained waiting time distributions characterized by exponential, stretched exponential, and power-law functions, as well as a crossover between two behaviors. For the mean square displacement, we found crossovers between usual and confined diffusion, and between usual and sub-diffusion. We obtained the exact expressions for the probability distributions, where non-Gaussian and stationary distributions emerged. This former feature is remarkable because the fractional diffusion equation is solved without external forces and subjected to the free diffusion boundary conditions. We have further shown that these new fractional diffusion equations are related to diffusive processes with stochastic resetting, and to fractional diffusion equations with derivatives of distributed order. Thus, our results suggest that these new operators may be a simple and efficient way for incorporating different structural aspects into the system, opening new possibilities for modeling and investigating anomalous diffusive processes.

Instability of turing patterns in reaction-diffusion-ODE systems.

PubMed

Marciniak-Czochra, Anna; Karch, Grzegorz; Suzuki, Kanako

2017-02-01

The aim of this paper is to contribute to the understanding of the pattern formation phenomenon in reaction-diffusion equations coupled with ordinary differential equations. Such systems of equations arise, for example, from modeling of interactions between cellular processes such as cell growth, differentiation or transformation and diffusing signaling factors. We focus on stability analysis of solutions of a prototype model consisting of a single reaction-diffusion equation coupled to an ordinary differential equation. We show that such systems are very different from classical reaction-diffusion models. They exhibit diffusion-driven instability (turing instability) under a condition of autocatalysis of non-diffusing component. However, the same mechanism which destabilizes constant solutions of such models, destabilizes also all continuous spatially heterogeneous stationary solutions, and consequently, there exist no stable Turing patterns in such reaction-diffusion-ODE systems. We provide a rigorous result on the nonlinear instability, which involves the analysis of a continuous spectrum of a linear operator induced by the lack of diffusion in the destabilizing equation. These results are extended to discontinuous patterns for a class of nonlinearities.

A tutorial solution to scattering of radiation in a thin atmosphere bounded below by a diffusely reflecting, absorbing surface

NASA Technical Reports Server (NTRS)

Buglia, J. J.

1982-01-01

A simple tutorial method, based on a photon tracking procedure, is described to determine the spherical albedo for a thin atmosphere overlying a reflecting surface. This procedure is used to provide a physical structure with which to interpret the more detailed but highly mathematical analyses presented. The final equations are shown to be in good numerical agreement with more exact solutions for thin atmospheres.

Effect of localized surface-plasmon mode on exciton transport and radiation emission in carbon nanotubes.

PubMed

Roslyak, Oleksiy; Cherqui, Charles; Dunlap, David H; Piryatinski, Andrei

2014-07-17

We report on a general theoretical approach to study exciton transport and emission in a single-walled carbon nanotube (SWNT) in the presence of a localized surface-plasmon (SP) mode within a metal nanoparticle interacting via near-field coupling. We derive a set of quantum mechanical equations of motion and approximate rate equations that account for the exciton, SP, and the environmental degrees of freedom. The material equations are complemented by an expression for the radiated power that depends on the exciton and SP populations and coherences, allowing for an examination of the angular distribution of the emitted radiation that would be measured in experiment. Numerical simulations for a (6,5) SWNT and cone-shaped Ag metal tip (MT) have been performed using this methodology. Comparison with physical parameters shows that the near-field interaction between the exciton-SP occurs in a weak coupling regime, with the diffusion processes being much faster than the exciton-SP population exchange. In such a case, the effect of the exciton population transfer to the MT with its subsequent dissipation (i.e., the Förster energy transfer) is to modify the exciton steady state distribution while reducing the equilibration time for excitons to reach a steady sate distribution. We find that the radiation distribution is dominated by SP emission for a SWNT-MT separation of a few tens of nanometers due to the fast SP emission rate, whereas the exciton-SP coherences can cause its rotation.

Influence of thermal radiation on soot production in Laminar axisymmetric diffusion flames

NASA Astrophysics Data System (ADS)

Demarco, R.; Nmira, F.; Consalvi, J. L.

2013-05-01

The aim of this paper is to study the effect of radiative heat transfer on soot production in laminar axisymmetric diffusion flames. Twenty-four C1-C3 hydrocarbon-air flames, consisting of normal (NDF) and inverse (IDF) diffusion flames at both normal gravity (1 g) and microgravity (0 g), and covering a wide range of conditions affecting radiative heat transfer, were simulated. The numerical model is based on the Steady Laminar Flamelet (SLF) model, a semi-empirical two-equation acetylene/benzene based soot model and the Statistical Narrow Band Correlated K (SNBCK) model coupled to the Finite Volume Method (FVM) to compute thermal radiation. Predictions relative to velocity, temperature, soot volume fraction and radiative losses are on the whole in good agreement with the available experimental data. Model results show that, for all the flames considered, thermal radiation is a crucial process with a view to providing accurate predictions for temperatures and soot concentrations. It becomes increasingly significant from IDFs to NDFs and its influence is much greater as gravity is reduced. The radiative contribution of gas prevails in the weakly-sooting IDFs and in the methane and ethane NDFs, whereas soot radiation dominates in the other flames. However, both contributions are significant in all cases, with the exception of the 1 g IDFs investigated where soot radiation can be ignored. The optically-thin approximation (OTA) was also tested and found to be applicable as long as the optical thickness, based on flame radius and Planck mean absorption coefficient, is less than 0.05. The OTA is reasonable for the IDFs and for most of the 1 g NDFs, but it fails to predict the radiative heat transfer for the 0 g NDFs. The accuracy of radiative-property models was then assessed in the latter cases. Simulations show that the gray approximation can be applied to soot but not to combustion gases. Both the non-gray and gray soot versions of the Full Spectrum Correlated k (FSCK) model can be then substituted for the SNBCK with a reduction in CPU time by a factor of about 20 in the latter case.

Background-Error Correlation Model Based on the Implicit Solution of a Diffusion Equation

DTIC Science & Technology

2010-01-01

1 Background- Error Correlation Model Based on the Implicit Solution of a Diffusion Equation Matthew J. Carrier* and Hans Ngodock...4. TITLE AND SUBTITLE Background- Error Correlation Model Based on the Implicit Solution of a Diffusion Equation 5a. CONTRACT NUMBER 5b. GRANT...2001), which sought to model error correlations based on the explicit solution of a generalized diffusion equation. The implicit solution is

Improved understanding of the acoustophoretic focusing of dense suspensions in a microchannel

NASA Astrophysics Data System (ADS)

Karthick, S.; Sen, A. K.

2017-11-01

We provide improved understanding of acoustophoretic focusing of a dense suspension (volume fraction φ >10 % ) in a microchannel subjected to an acoustic standing wave using a proposed theoretical model and experiments. The model is based on the theory of interacting continua and utilizes a momentum transport equation for the mixture, continuity equation, and transport equation for the solid phase. The model demonstrates the interplay between acoustic radiation and shear-induced diffusion (SID) forces that is critical in the focusing of dense suspensions. The shear-induced particle migration model of Leighton and Acrivos, coupled with the acoustic radiation force, is employed to simulate the continuum behavior of particles. In the literature, various closures for the diffusion coefficient Dφ* are available for rigid spheres at high concentrations and nonspherical deformable particles [e.g., red blood cells (RBCs)] at low concentrations. Here we propose a closure for Dφ* for dense suspension of RBCs and validate the proposed model with experimental data. While the available closures for Dφ* fail to predict the acoustic focusing of a dense suspension of nonspherical deformable particles like RBCs, the predictions of the proposed model match experimental data within 15%. Both the model and experiments reveal a competition between acoustic radiation and SID forces that gives rise to an equilibrium width w* of a focused stream of particles at some distance Leq* along the flow direction. Using different shear rates, acoustic energy densities, and particle concentrations, we show that the equilibrium width is governed by Péclet number Pe and Strouhal number Stas w*=1.4(PeSt) -0.5 while the length required to obtain the equilibrium-focused width depends on St as Leq*=3.8 /(St)0.6 . The proposed model and correlations would find significance in the design of microchannels for acoustic focusing of dense suspensions such as undiluted blood.

Diffusion Forecasting Model with Basis Functions from QR-Decomposition

NASA Astrophysics Data System (ADS)

Harlim, John; Yang, Haizhao

2018-06-01

The diffusion forecasting is a nonparametric approach that provably solves the Fokker-Planck PDE corresponding to Itô diffusion without knowing the underlying equation. The key idea of this method is to approximate the solution of the Fokker-Planck equation with a discrete representation of the shift (Koopman) operator on a set of basis functions generated via the diffusion maps algorithm. While the choice of these basis functions is provably optimal under appropriate conditions, computing these basis functions is quite expensive since it requires the eigendecomposition of an N× N diffusion matrix, where N denotes the data size and could be very large. For large-scale forecasting problems, only a few leading eigenvectors are computationally achievable. To overcome this computational bottleneck, a new set of basis functions constructed by orthonormalizing selected columns of the diffusion matrix and its leading eigenvectors is proposed. This computation can be carried out efficiently via the unpivoted Householder QR factorization. The efficiency and effectiveness of the proposed algorithm will be shown in both deterministically chaotic and stochastic dynamical systems; in the former case, the superiority of the proposed basis functions over purely eigenvectors is significant, while in the latter case forecasting accuracy is improved relative to using a purely small number of eigenvectors. Supporting arguments will be provided on three- and six-dimensional chaotic ODEs, a three-dimensional SDE that mimics turbulent systems, and also on the two spatial modes associated with the boreal winter Madden-Julian Oscillation obtained from applying the Nonlinear Laplacian Spectral Analysis on the measured Outgoing Longwave Radiation.

Diffusion Forecasting Model with Basis Functions from QR-Decomposition

NASA Astrophysics Data System (ADS)

Harlim, John; Yang, Haizhao

2017-12-01

The diffusion forecasting is a nonparametric approach that provably solves the Fokker-Planck PDE corresponding to Itô diffusion without knowing the underlying equation. The key idea of this method is to approximate the solution of the Fokker-Planck equation with a discrete representation of the shift (Koopman) operator on a set of basis functions generated via the diffusion maps algorithm. While the choice of these basis functions is provably optimal under appropriate conditions, computing these basis functions is quite expensive since it requires the eigendecomposition of an N× N diffusion matrix, where N denotes the data size and could be very large. For large-scale forecasting problems, only a few leading eigenvectors are computationally achievable. To overcome this computational bottleneck, a new set of basis functions constructed by orthonormalizing selected columns of the diffusion matrix and its leading eigenvectors is proposed. This computation can be carried out efficiently via the unpivoted Householder QR factorization. The efficiency and effectiveness of the proposed algorithm will be shown in both deterministically chaotic and stochastic dynamical systems; in the former case, the superiority of the proposed basis functions over purely eigenvectors is significant, while in the latter case forecasting accuracy is improved relative to using a purely small number of eigenvectors. Supporting arguments will be provided on three- and six-dimensional chaotic ODEs, a three-dimensional SDE that mimics turbulent systems, and also on the two spatial modes associated with the boreal winter Madden-Julian Oscillation obtained from applying the Nonlinear Laplacian Spectral Analysis on the measured Outgoing Longwave Radiation.

Radiation Transport for Explosive Outflows: Opacity Regrouping

NASA Astrophysics Data System (ADS)

Wollaeger, Ryan T.; van Rossum, Daniel R.

2014-10-01

Implicit Monte Carlo (IMC) and Discrete Diffusion Monte Carlo (DDMC) are methods used to stochastically solve the radiative transport and diffusion equations, respectively. These methods combine into a hybrid transport-diffusion method we refer to as IMC-DDMC. We explore a multigroup IMC-DDMC scheme that in DDMC, combines frequency groups with sufficient optical thickness. We term this procedure "opacity regrouping." Opacity regrouping has previously been applied to IMC-DDMC calculations for problems in which the dependence of the opacity on frequency is monotonic. We generalize opacity regrouping to non-contiguous groups and implement this in SuperNu, a code designed to do radiation transport in high-velocity outflows with non-monotonic opacities. We find that regrouping of non-contiguous opacity groups generally improves the speed of IMC-DDMC radiation transport. We present an asymptotic analysis that informs the nature of the Doppler shift in DDMC groups and summarize the derivation of the Gentile-Fleck factor for modified IMC-DDMC. We test SuperNu using numerical experiments including a quasi-manufactured analytic solution, a simple 10 group problem, and the W7 problem for Type Ia supernovae. We find that opacity regrouping is necessary to make our IMC-DDMC implementation feasible for the W7 problem and possibly Type Ia supernova simulations in general. We compare the bolometric light curves and spectra produced by the SuperNu and PHOENIX radiation transport codes for the W7 problem. The overall shape of the bolometric light curves are in good agreement, as are the spectra and their evolution with time. However, for the numerical specifications we considered, we find that the peak luminosity of the light curve calculated using SuperNu is ~10% less than that calculated using PHOENIX.

Direct Collapse to Supermassive Black Hole Seeds with Radiative Transfer: Isolated Halos

NASA Astrophysics Data System (ADS)

Luo, Yang; Ardaneh, Kazem; Shlosman, Isaac; Nagamine, Kentaro; Wise, John H.; Begelman, Mitchell C.

2018-05-01

Direct collapse within dark matter haloes is a promising path to form supermassive black hole seeds at high redshifts. The outer part of this collapse remains optically thin. However, the innermost region of the collapse is expected to become optically thick and requires to follow the radiation field in order to understand its evolution. So far, the adiabatic approximation has been used exclusively for this purpose. We apply radiative transfer in the flux-limited diffusion (FLD) approximation to solve the evolution of coupled gas and radiation for isolated haloes. We find that (1) the photosphere forms at 10-6 pc and rapidly expands outwards. (2) A central core forms, with a mass of 1 M⊙, supported by gas pressure gradients and rotation. (3) Growing gas and radiation pressure gradients dissolve it. (4) This process is associated with a strong anisotropic outflow; another core forms nearby and grows rapidly. (5) Typical radiation luminosity emerging from the photosphere is 5 × 1037-5 × 1038 erg s-1, of the order the Eddington luminosity. (6) Two variability time-scales are associated with this process: a long one, which is related to the accretion flow within the central 10-4-10-3 pc, and 0.1 yr, related to radiation diffusion. (7) Adiabatic models evolution differs profoundly from that of the FLD models, by forming a geometrically thick disc. Overall, an adiabatic equation of state is not a good approximation to the advanced stage of direct collapse, because the radiation is capable of escaping due to anisotropy in the optical depth and associated gradients.

THERMOHALINE INSTABILITIES INSIDE STARS: A SYNTHETIC STUDY INCLUDING EXTERNAL TURBULENCE AND RADIATIVE LEVITATION

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

Vauclair, Sylvie; Theado, Sylvie, E-mail: sylvie.vauclair@irap.omp.eu

2012-07-01

We have derived a new expression for the thermohaline mixing coefficient in stars, including the effects of radiative levitation and external turbulence, by solving Boussinesq equations in a nearly incompressible stratified fluid with a linear approximation. It is well known that radiative levitation of individual elements can lead to their accumulation in specific stellar layers. In some cases, it can induce important effects on the stellar structure. Here we confirm that this accumulation is moderated by thermohaline convection due to the resulting inverse {mu}-gradient. The new coefficient that we have derived shows that the effect of radiative accelerations on themore » thermohaline instability itself is small. This effect must however be checked in all computations. We also confirm that the presence of large horizontal turbulence can reduce or even suppress the thermohaline convection. These results are important as they concern all the cases of heavy element accumulation in stars. Computations of radiative diffusion must be revisited to include thermohaline convection and its consequences. It may be one of the basic reasons for the fact that the observed abundances are always smaller than those predicted by pure atomic diffusion. In any case, these processes have to compete with rotation-induced mixing, but this competition is more complex than previously thought due to their mutual interaction.« less

Efficient Charge Collection in Coplanar-Grid Radiation Detectors

NASA Astrophysics Data System (ADS)

Kunc, J.; Praus, P.; Belas, E.; Dědič, V.; Pekárek, J.; Grill, R.

2018-05-01

We model laser-induced transient-current waveforms in radiation coplanar-grid detectors. Poisson's equation is solved by the finite-element method and currents induced by a photogenerated charge are obtained using the Shockley-Ramo theorem. The spectral response on a radiation flux is modeled by Monte Carlo simulations. We show a 10 × improved spectral resolution of the coplanar-grid detector using differential signal sensing. We model the current waveform dependence on the doping, depletion width, diffusion, and detector shielding, and their mutual dependence is discussed in terms of detector optimization. The numerical simulations are successfully compared to experimental data, and further model simplifications are proposed. The space charge below electrodes and a nonhomogeneous electric field on a coplanar-grid anode are found to be the dominant contributions to laser-induced transient-current waveforms.

Diffusion of Charged Species in Liquids

NASA Astrophysics Data System (ADS)

Del Río, J. A.; Whitaker, S.

2016-11-01

In this study the laws of mechanics for multi-component systems are used to develop a theory for the diffusion of ions in the presence of an electrostatic field. The analysis begins with the governing equation for the species velocity and it leads to the governing equation for the species diffusion velocity. Simplification of this latter result provides a momentum equation containing three dominant forces: (a) the gradient of the partial pressure, (b) the electrostatic force, and (c) the diffusive drag force that is a central feature of the Maxwell-Stefan equations. For ideal gas mixtures we derive the classic Nernst-Planck equation. For liquid-phase diffusion we encounter a situation in which the Nernst-Planck contribution to diffusion differs by several orders of magnitude from that obtained for ideal gases.

Diffusion of Charged Species in Liquids.

PubMed

Del Río, J A; Whitaker, S

2016-11-04

In this study the laws of mechanics for multi-component systems are used to develop a theory for the diffusion of ions in the presence of an electrostatic field. The analysis begins with the governing equation for the species velocity and it leads to the governing equation for the species diffusion velocity. Simplification of this latter result provides a momentum equation containing three dominant forces: (a) the gradient of the partial pressure, (b) the electrostatic force, and (c) the diffusive drag force that is a central feature of the Maxwell-Stefan equations. For ideal gas mixtures we derive the classic Nernst-Planck equation. For liquid-phase diffusion we encounter a situation in which the Nernst-Planck contribution to diffusion differs by several orders of magnitude from that obtained for ideal gases.

Diffusion of Charged Species in Liquids

PubMed Central

del Río, J. A.; Whitaker, S.

2016-01-01

In this study the laws of mechanics for multi-component systems are used to develop a theory for the diffusion of ions in the presence of an electrostatic field. The analysis begins with the governing equation for the species velocity and it leads to the governing equation for the species diffusion velocity. Simplification of this latter result provides a momentum equation containing three dominant forces: (a) the gradient of the partial pressure, (b) the electrostatic force, and (c) the diffusive drag force that is a central feature of the Maxwell-Stefan equations. For ideal gas mixtures we derive the classic Nernst-Planck equation. For liquid-phase diffusion we encounter a situation in which the Nernst-Planck contribution to diffusion differs by several orders of magnitude from that obtained for ideal gases. PMID:27811959

Effects of Gravity on Soot Formation in a Coflow Laminar Methane/Air Diffusion Flame

NASA Astrophysics Data System (ADS)

Kong, Wenjun; Liu, Fengshan

2010-04-01

Simulations of a laminar coflow methane/air diffusion flame at atmospheric pressure are conducted to gain better understanding of the effects of gravity on soot formation by using detailed gas-phase chemistry, complex thermal and transport properties coupled with a semiempirical two-equation soot model and a nongray radiation model. Soot oxidation by O2, OH and O was considered. Thermal radiation was calculated using the discrete ordinate method coupled with a statistical narrow-band correlated-K model. The spectral absorption coefficient of soot was obtained by Rayleigh's theory for small particles. The results show that the peak temperature decreases with the decrease of the gravity level. The peak soot volume fraction in microgravity is about twice of that in normal gravity under the present conditions. The numerical results agree very well with available experimental results. The predicted results also show that gravity affects the location and intensity for soot nucleation and surface growth.

Direct Radiative Effect and Heating Rate of black carbon aerosol: high time resolution measurements and source-identified forcing effects

NASA Astrophysics Data System (ADS)

Ferrero, Luca; Mocnik, Grisa; Cogliati, Sergio; Comi, Alberto; Degni, Francesca; Di Mauro, Biagio; Colombo, Roberto; Bolzacchini, Ezio

2016-04-01

Black carbon (BC) absorbs sunlight in the atmosphere heating it. However, up to now, heating rate (HR) calculations from the divergence of the net radiative flux with altitude or from the modelling activity are too sparse. This work fills the aforementioned gap presenting a new methodology based on a full set of physical equations to experimentally determine both the radiative power density absorbed into a ground-based atmospheric layer (ADRE), and the consequent HR induced by the absorptive component of aerosol. In urban context, it is essentially related to the BC. The methodology is also applicable to natural components (i.e. dust) and is obtained solving the first derivative of the main radiative transfer equations. The ADRE and the consequent HR can be determined coupling spectral aerosol absorption measurements with the spectrally resolved measurements of the direct, diffuse downward radiation and the surface reflected radiance components. Moreover, the spectral absorption of BC aerosol allows its source apportionment (traffic and biomass burning (BB)) allowing the same apportionment on HR. This work reports one year of high-time resolution measurements (5 min) of sunlight absorption and HR induced by BC aerosol over Milan. A unique sampling site was set up from March 2015 with: 1) Aethalometer (AE-31, Magee Scientific, 7-λ), 2) the Multiplexer-Radiometer-Irradiometer which detects downward and reflected radiance (350-1000 nm in 3648 spectral bands) coupled with a rotating shadow-band to measure spectrally-resolved global and diffuse radiation (thus direct), 3) a meteorological station (LSI-Lastem) equipped with 3 pyranometers (global, diffuse and refrected radiation; 300-3000 nm), a thermohygrometer, a barometer, an anemometer, 4) condensation and optical particle counters (TSI 3775 and Grimm 1.107), 5) low volume sampler (FAI Hydra dual sampler, PM2.5 and PM10) for sample collection and chemistry determination. Results concerning the radiative power density absorbed by BC and the consequent HR allowed to determine: 1) the mean monthly values along one year (i.e. October: 14.5±0.2 mW/m3, 1.04±0.01 K/day of HR, 3.0±0.1 μg/m3 of BC); 2) the importance of the direct, diffuse and reflected radiation and thus of sky conditions (clear/cloudy) on the HR (i.e. October HR: 0.42±0.10 K/day for direct, 0.44±0.10 K/day for diffuse, 0.18±0.10 K/day for reflected); 3) the daily cycle influence of BC and radiation on HR; 4) the influence of anthropogenic activity studying the daily cycle along working and non-working days (i.e. September HR: 1.00±0.06 K/day for working days, 0.35±0.02 K/day for non-working days); 5) determine the radiative effect of traffic and BB sources of BC in function of the domestic heating operation. In particular, the domestic heating is allowed in Italy starting from 15th October: traffic BC for 1-15 and 15-31 October was 1.3±0.1 μg/m3 and 2.3±0.1 μg/m3 respectively, while BB BC was 0.7±0.1 μg/m3 and 1.5±0.1 μg/m3. In terms of HR, traffic BC for 1-15 and 15-31 October contributed with 0.46±0.01 K/day and 0.74±0.02 K/day while BB BC was 0.28±0.01 K/day and 0.61±0.02 K/day. All of the aforementioned results will be detailed in the presentation using the full set of data collected.

Calibration of a thin metal foil for infrared imaging video bolometer to estimate the spatial variation of thermal diffusivity using a photo-thermal technique.

PubMed

Pandya, Shwetang N; Peterson, Byron J; Sano, Ryuichi; Mukai, Kiyofumi; Drapiko, Evgeny A; Alekseyev, Andrey G; Akiyama, Tsuyoshi; Itomi, Muneji; Watanabe, Takashi

2014-05-01

A thin metal foil is used as a broad band radiation absorber for the InfraRed imaging Video Bolometer (IRVB), which is a vital diagnostic for studying three-dimensional radiation structures from high temperature plasmas in the Large Helical Device. The two-dimensional (2D) heat diffusion equation of the foil needs to be solved numerically to estimate the radiation falling on the foil through a pinhole geometry. The thermal, physical, and optical properties of the metal foil are among the inputs to the code besides the spatiotemporal variation of temperature, for reliable estimation of the exhaust power from the plasma illuminating the foil. The foil being very thin and of considerable size, non-uniformities in these properties need to be determined by suitable calibration procedures. The graphite spray used for increasing the surface emissivity also contributes to a change in the thermal properties. This paper discusses the application of the thermographic technique for determining the spatial variation of the effective in-plane thermal diffusivity of the thin metal foil and graphite composite. The paper also discusses the advantages of this technique in the light of limitations and drawbacks presented by other calibration techniques being practiced currently. The technique is initially applied to a material of known thickness and thermal properties for validation and finally to thin foils of gold and platinum both with two different thicknesses. It is observed that the effect of the graphite layer on the estimation of the thermal diffusivity becomes more pronounced for thinner foils and the measured values are approximately 2.5-3 times lower than the literature values. It is also observed that the percentage reduction in thermal diffusivity due to the coating is lower for high thermal diffusivity materials such as gold. This fact may also explain, albeit partially, the higher sensitivity of the platinum foil as compared to gold.

Calibration of a thin metal foil for infrared imaging video bolometer to estimate the spatial variation of thermal diffusivity using a photo-thermal technique

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

Pandya, Shwetang N., E-mail: pandya.shwetang@LHD.nifs.ac.jp; Sano, Ryuichi; Peterson, Byron J.

A thin metal foil is used as a broad band radiation absorber for the InfraRed imaging Video Bolometer (IRVB), which is a vital diagnostic for studying three-dimensional radiation structures from high temperature plasmas in the Large Helical Device. The two-dimensional (2D) heat diffusion equation of the foil needs to be solved numerically to estimate the radiation falling on the foil through a pinhole geometry. The thermal, physical, and optical properties of the metal foil are among the inputs to the code besides the spatiotemporal variation of temperature, for reliable estimation of the exhaust power from the plasma illuminating the foil.more » The foil being very thin and of considerable size, non-uniformities in these properties need to be determined by suitable calibration procedures. The graphite spray used for increasing the surface emissivity also contributes to a change in the thermal properties. This paper discusses the application of the thermographic technique for determining the spatial variation of the effective in-plane thermal diffusivity of the thin metal foil and graphite composite. The paper also discusses the advantages of this technique in the light of limitations and drawbacks presented by other calibration techniques being practiced currently. The technique is initially applied to a material of known thickness and thermal properties for validation and finally to thin foils of gold and platinum both with two different thicknesses. It is observed that the effect of the graphite layer on the estimation of the thermal diffusivity becomes more pronounced for thinner foils and the measured values are approximately 2.5–3 times lower than the literature values. It is also observed that the percentage reduction in thermal diffusivity due to the coating is lower for high thermal diffusivity materials such as gold. This fact may also explain, albeit partially, the higher sensitivity of the platinum foil as compared to gold.« less

Microscopic Interpretation and Generalization of the Bloch-Torrey Equation for Diffusion Magnetic Resonance

PubMed Central

Seroussi, Inbar; Grebenkov, Denis S.; Pasternak, Ofer; Sochen, Nir

2017-01-01

In order to bridge microscopic molecular motion with macroscopic diffusion MR signal in complex structures, we propose a general stochastic model for molecular motion in a magnetic field. The Fokker-Planck equation of this model governs the probability density function describing the diffusion-magnetization propagator. From the propagator we derive a generalized version of the Bloch-Torrey equation and the relation to the random phase approach. This derivation does not require assumptions such as a spatially constant diffusion coefficient, or ad-hoc selection of a propagator. In particular, the boundary conditions that implicitly incorporate the microstructure into the diffusion MR signal can now be included explicitly through a spatially varying diffusion coefficient. While our generalization is reduced to the conventional Bloch-Torrey equation for piecewise constant diffusion coefficients, it also predicts scenarios in which an additional term to the equation is required to fully describe the MR signal. PMID:28242566

Study on monostable and bistable reaction-diffusion equations by iteration of travelling wave maps

NASA Astrophysics Data System (ADS)

Yi, Taishan; Chen, Yuming

2017-12-01

In this paper, based on the iterative properties of travelling wave maps, we develop a new method to obtain spreading speeds and asymptotic propagation for monostable and bistable reaction-diffusion equations. Precisely, for Dirichlet problems of monostable reaction-diffusion equations on the half line, by making links between travelling wave maps and integral operators associated with the Dirichlet diffusion kernel (the latter is NOT invariant under translation), we obtain some iteration properties of the Dirichlet diffusion and some a priori estimates on nontrivial solutions of Dirichlet problems under travelling wave transformation. We then provide the asymptotic behavior of nontrivial solutions in the space-time region for Dirichlet problems. These enable us to develop a unified method to obtain results on heterogeneous steady states, travelling waves, spreading speeds, and asymptotic spreading behavior for Dirichlet problem of monostable reaction-diffusion equations on R+ as well as of monostable/bistable reaction-diffusion equations on R.

A NUMERICAL ALGORITHM FOR MODELING MULTIGROUP NEUTRINO-RADIATION HYDRODYNAMICS IN TWO SPATIAL DIMENSIONS

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

Swesty, F. Douglas; Myra, Eric S.

It is now generally agreed that multidimensional, multigroup, neutrino-radiation hydrodynamics (RHD) is an indispensable element of any realistic model of stellar-core collapse, core-collapse supernovae, and proto-neutron star instabilities. We have developed a new, two-dimensional, multigroup algorithm that can model neutrino-RHD flows in core-collapse supernovae. Our algorithm uses an approach similar to the ZEUS family of algorithms, originally developed by Stone and Norman. However, this completely new implementation extends that previous work in three significant ways: first, we incorporate multispecies, multigroup RHD in a flux-limited-diffusion approximation. Our approach is capable of modeling pair-coupled neutrino-RHD, and includes effects of Pauli blocking inmore » the collision integrals. Blocking gives rise to nonlinearities in the discretized radiation-transport equations, which we evolve implicitly in time. We employ parallelized Newton-Krylov methods to obtain a solution of these nonlinear, implicit equations. Our second major extension to the ZEUS algorithm is the inclusion of an electron conservation equation that describes the evolution of electron-number density in the hydrodynamic flow. This permits calculating deleptonization of a stellar core. Our third extension modifies the hydrodynamics algorithm to accommodate realistic, complex equations of state, including those having nonconvex behavior. In this paper, we present a description of our complete algorithm, giving sufficient details to allow others to implement, reproduce, and extend our work. Finite-differencing details are presented in appendices. We also discuss implementation of this algorithm on state-of-the-art, parallel-computing architectures. Finally, we present results of verification tests that demonstrate the numerical accuracy of this algorithm on diverse hydrodynamic, gravitational, radiation-transport, and RHD sample problems. We believe our methods to be of general use in a variety of model settings where radiation transport or RHD is important. Extension of this work to three spatial dimensions is straightforward.« less

Prediction of stream volatilization coefficients

USGS Publications Warehouse

Rathbun, Ronald E.

1990-01-01

Equations are developed for predicting the liquid-film and gas-film reference-substance parameters for quantifying volatilization of organic solutes from streams. Molecular weight and molecular-diffusion coefficients of the solute are used as correlating parameters. Equations for predicting molecular-diffusion coefficients of organic solutes in water and air are developed, with molecular weight and molal volume as parameters. Mean absolute errors of prediction for diffusion coefficients in water are 9.97% for the molecular-weight equation, 6.45% for the molal-volume equation. The mean absolute error for the diffusion coefficient in air is 5.79% for the molal-volume equation. Molecular weight is not a satisfactory correlating parameter for diffusion in air because two equations are necessary to describe the values in the data set. The best predictive equation for the liquid-film reference-substance parameter has a mean absolute error of 5.74%, with molal volume as the correlating parameter. The best equation for the gas-film parameter has a mean absolute error of 7.80%, with molecular weight as the correlating parameter.

Causal Diffusion and the Survival of Charge Fluctuations

NASA Astrophysics Data System (ADS)

Abdel-Aziz, Mohamed; Gavin, Sean

2004-10-01

Diffusion may obliterate fluctuation signals of the QCD phase transition in nuclear collisions at SPS and RHIC energies. We propose a hyperbolic diffusion equation to study the dissipation of net charge fluctuations [1]. This equation is needed in a relativistic context, because the classic parabolic diffusion equation violates causality. We find that causality substantially limits the extent to which diffusion can dissipate these fluctuations. [1] M. Abdel-Aziz and S. Gavin, nucl-th/0404058

Theoretical Relationships between Luminescence and Hillslope Soil Vertical Diffusivity: a Numerical Modeling Approach

NASA Astrophysics Data System (ADS)

Gray, H. J.; Tucker, G. E.; Mahan, S.

2017-12-01

Luminescence is a property of matter that can be used to obtain depositional ages from fine sand. Luminescence generates due to exposure to background ionizing radiation and is removed by sunlight exposure in a process known as bleaching. There is evidence to suggest that luminescence can also serve as a sediment tracer in fluvial and hillslope environments. For hillslope environments, it has been suggested that the magnitude of luminescence as a function of soil depth is related to the strength of soil mixing. Hillslope soils with a greater extent of mixing will have previously surficial sand grains moved to greater depths in a soil column. These previously surface-exposed grains will contain a lower luminescence than those which have never seen the surface. To attempt to connect luminescence profiles with soil mixing rate, here defined as the soil vertical diffusivity, I conduct numerical modelling of particles in hillslope soils coupled with equations describing the physics of luminescence. I use recently published equations describing the trajectories of particles under both exponential and uniform soil velocity soils profiles and modify them to include soil diffusivity. Results from the model demonstrates a strong connection between soil diffusivity and luminescence. Both the depth profiles of luminescence and the total percent of surface exposed grains will change drastically based on the magnitude of the diffusivity. This suggests that luminescence could potentially be used to infer the magnitude of soil diffusivity. However, I test other variables such as the soil production rate, e-folding length of soil velocity, background dose rate, and soil thickness, and I find these other variables can also affect the relationship between luminescence and diffusivity. This suggests that these other variables may need to be constrained prior to any inferences of soil diffusivity from luminescence measurements. Further field testing of the model in areas where the soil vertical diffusivity and other parameters are independently known will provide a test of this potential new method.

Global dynamics of a nonlocal delayed reaction-diffusion equation on a half plane

NASA Astrophysics Data System (ADS)

Hu, Wenjie; Duan, Yueliang

2018-04-01

We consider a delayed reaction-diffusion equation with spatial nonlocality on a half plane that describes population dynamics of a two-stage species living in a semi-infinite environment. A Neumann boundary condition is imposed accounting for an isolated domain. To describe the global dynamics, we first establish some a priori estimate for nontrivial solutions after investigating asymptotic properties of the nonlocal delayed effect and the diffusion operator, which enables us to show the permanence of the equation with respect to the compact open topology. We then employ standard dynamical system arguments to establish the global attractivity of the nontrivial equilibrium. The main results are illustrated by the diffusive Nicholson's blowfly equation and the diffusive Mackey-Glass equation.

FRACTIONAL PEARSON DIFFUSIONS.

PubMed

Leonenko, Nikolai N; Meerschaert, Mark M; Sikorskii, Alla

2013-07-15

Pearson diffusions are governed by diffusion equations with polynomial coefficients. Fractional Pearson diffusions are governed by the corresponding time-fractional diffusion equation. They are useful for modeling sub-diffusive phenomena, caused by particle sticking and trapping. This paper provides explicit strong solutions for fractional Pearson diffusions, using spectral methods. It also presents stochastic solutions, using a non-Markovian inverse stable time change.

A nonlinear equation for ionic diffusion in a strong binary electrolyte

PubMed Central

Ghosal, Sandip; Chen, Zhen

2010-01-01

The problem of the one-dimensional electro-diffusion of ions in a strong binary electrolyte is considered. The mathematical description, known as the Poisson–Nernst–Planck (PNP) system, consists of a diffusion equation for each species augmented by transport owing to a self-consistent electrostatic field determined by the Poisson equation. This description is also relevant to other important problems in physics, such as electron and hole diffusion across semiconductor junctions and the diffusion of ions in plasmas. If concentrations do not vary appreciably over distances of the order of the Debye length, the Poisson equation can be replaced by the condition of local charge neutrality first introduced by Planck. It can then be shown that both species diffuse at the same rate with a common diffusivity that is intermediate between that of the slow and fast species (ambipolar diffusion). Here, we derive a more general theory by exploiting the ratio of the Debye length to a characteristic length scale as a small asymptotic parameter. It is shown that the concentration of either species may be described by a nonlinear partial differential equation that provides a better approximation than the classical linear equation for ambipolar diffusion, but reduces to it in the appropriate limit. PMID:21818176

An Ab Initio and Kinetic Monte Carlo Simulation Study of Lithium Ion Diffusion on Graphene

PubMed Central

Zhong, Kehua; Yang, Yanmin; Xu, Guigui; Zhang, Jian-Min; Huang, Zhigao

2017-01-01

The Li+ diffusion coefficients in Li+-adsorbed graphene systems were determined by combining first-principle calculations based on density functional theory with Kinetic Monte Carlo simulations. The calculated results indicate that the interactions between Li ions have a very important influence on lithium diffusion. Based on energy barriers directly obtained from first-principle calculations for single-Li+ and two-Li+ adsorbed systems, a new equation predicting energy barriers with more than two Li ions was deduced. Furthermore, it is found that the temperature dependence of Li+ diffusion coefficients fits well to the Arrhenius equation, rather than meeting the equation from electrochemical impedance spectroscopy applied to estimate experimental diffusion coefficients. Moreover, the calculated results also reveal that Li+ concentration dependence of diffusion coefficients roughly fits to the equation from electrochemical impedance spectroscopy in a low concentration region; however, it seriously deviates from the equation in a high concentration region. So, the equation from electrochemical impedance spectroscopy technique could not be simply used to estimate the Li+ diffusion coefficient for all Li+-adsorbed graphene systems with various Li+ concentrations. Our work suggests that interactions between Li ions, and among Li ion and host atoms will influence the Li+ diffusion, which determines that the Li+ intercalation dependence of Li+ diffusion coefficient should be changed and complex. PMID:28773122

Chemical vapor deposition fluid flow simulation modelling tool

NASA Technical Reports Server (NTRS)

Bullister, Edward T.

1992-01-01

Accurate numerical simulation of chemical vapor deposition (CVD) processes requires a general purpose computational fluid dynamics package combined with specialized capabilities for high temperature chemistry. In this report, we describe the implementation of these specialized capabilities in the spectral element code NEKTON. The thermal expansion of the gases involved is shown to be accurately approximated by the low Mach number perturbation expansion of the incompressible Navier-Stokes equations. The radiative heat transfer between multiple interacting radiating surfaces is shown to be tractable using the method of Gebhart. The disparate rates of reaction and diffusion in CVD processes are calculated via a point-implicit time integration scheme. We demonstrate the use above capabilities on prototypical CVD applications.

A RADIATION TRANSFER SOLVER FOR ATHENA USING SHORT CHARACTERISTICS

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

Davis, Shane W.; Stone, James M.; Jiang Yanfei

2012-03-01

We describe the implementation of a module for the Athena magnetohydrodynamics (MHD) code that solves the time-independent, multi-frequency radiative transfer (RT) equation on multidimensional Cartesian simulation domains, including scattering and non-local thermodynamic equilibrium (LTE) effects. The module is based on well known and well tested algorithms developed for modeling stellar atmospheres, including the method of short characteristics to solve the RT equation, accelerated Lambda iteration to handle scattering and non-LTE effects, and parallelization via domain decomposition. The module serves several purposes: it can be used to generate spectra and images, to compute a variable Eddington tensor (VET) for full radiationmore » MHD simulations, and to calculate the heating and cooling source terms in the MHD equations in flows where radiation pressure is small compared with gas pressure. For the latter case, the module is combined with the standard MHD integrators using operator splitting: we describe this approach in detail, including a new constraint on the time step for stability due to radiation diffusion modes. Implementation of the VET method for radiation pressure dominated flows is described in a companion paper. We present results from a suite of test problems for both the RT solver itself and for dynamical problems that include radiative heating and cooling. These tests demonstrate that the radiative transfer solution is accurate and confirm that the operator split method is stable, convergent, and efficient for problems of interest. We demonstrate there is no need to adopt ad hoc assumptions of questionable accuracy to solve RT problems in concert with MHD: the computational cost for our general-purpose module for simple (e.g., LTE gray) problems can be comparable to or less than a single time step of Athena's MHD integrators, and only few times more expensive than that for more general (non-LTE) problems.« less

Inverse Compton Scattering in Mildly Relativistic Plasma

NASA Technical Reports Server (NTRS)

Molnar, S. M.; Birkinshaw, M.

1998-01-01

We investigated the effect of inverse Compton scattering in mildly relativistic static and moving plasmas with low optical depth using Monte Carlo simulations, and calculated the Sunyaev-Zel'dovich effect in the cosmic background radiation. Our semi-analytic method is based on a separation of photon diffusion in frequency and real space. We use Monte Carlo simulation to derive the intensity and frequency of the scattered photons for a monochromatic incoming radiation. The outgoing spectrum is determined by integrating over the spectrum of the incoming radiation using the intensity to determine the correct weight. This method makes it possible to study the emerging radiation as a function of frequency and direction. As a first application we have studied the effects of finite optical depth and gas infall on the Sunyaev-Zel'dovich effect (not possible with the extended Kompaneets equation) and discuss the parameter range in which the Boltzmann equation and its expansions can be used. For high temperature clusters (k(sub B)T(sub e) greater than or approximately equal to 15 keV) relativistic corrections based on a fifth order expansion of the extended Kompaneets equation seriously underestimate the Sunyaev-Zel'dovich effect at high frequencies. The contribution from plasma infall is less important for reasonable velocities. We give a convenient analytical expression for the dependence of the cross-over frequency on temperature, optical depth, and gas infall speed. Optical depth effects are often more important than relativistic corrections, and should be taken into account for high-precision work, but are smaller than the typical kinematic effect from cluster radial velocities.

Thermal diffusion effect on MHD mixed convective flow along a vertically inclined plate: A casson fluid flow

NASA Astrophysics Data System (ADS)

Prasad, D. V. V. Krishna; Chaitanya, G. S. Krishna; Raju, R. Srinivasa

2018-05-01

The nature of Casson fluid on MHD free convective flow of over an impulsively started infinite vertically inclined plate in presence of thermal diffusion (Soret), thermal radiation, heat and mass transfer effects is studied. The basic governing nonlinear coupled partial differential equations are solved numerically using finite element method. The relevant physical parameters appearing in velocity, temperature and concentration profiles are analyzed and discussed through graphs. Finally, the results for velocity profiles and the reduced Nusselt and Sherwood numbers are obtained and compared with previous results in the literature and are found to be in excellent agreement. Applications of the present study would be useful in magnetic material processing and chemical engineering systems.

Diffusion equations and the time evolution of foreign exchange rates

NASA Astrophysics Data System (ADS)

Figueiredo, Annibal; de Castro, Marcio T.; da Fonseca, Regina C. B.; Gleria, Iram

2013-10-01

We investigate which type of diffusion equation is most appropriate to describe the time evolution of foreign exchange rates. We modify the geometric diffusion model assuming a non-exponential time evolution and the stochastic term is the sum of a Wiener noise and a jump process. We find the resulting diffusion equation to obey the Kramers-Moyal equation. Analytical solutions are obtained using the characteristic function formalism and compared with empirical data. The analysis focus on the first four central moments considering the returns of foreign exchange rate. It is shown that the proposed model offers a good improvement over the classical geometric diffusion model.

A numerical study of three-dimensional flame propagation over thin solids in purely forced concurrent flow including gas-phase radiation

NASA Astrophysics Data System (ADS)

Feier, Ioan I., Jr.

The effect of flame radiation on concurrent-flow flame spread over a thin solid sample of finite width in a low-speed wind tunnel is modeled using three-dimensional full Navier-Stokes equations and three-dimensional flame radiation transfer equations. The formulation includes the conservation of mass, momentum, energy, and species: fuel vapor, oxygen, carbon dioxide and water vapor. The SN discrete ordinates method is used to solve the radiation transfer equation with a mean absorption coefficient kappa = Ckappa p, where kappap is the Planck mean absorption coefficient of the gas mixture. The varying parameter C has a value between 0 and 1; C represents the strength of flame radiation. In addition, the solid fuel absorptivity alpha is varied to ascertain the effect of flame radiation heat feedback to the solid. The flow tunnel modeled has a dimension of 10x10x30 cm, the solid fuel has a width of 6-cm with two 1-cm inert strips as edges. Incoming forced flow velocity (5 cm/s) of 21% oxygen is assumed. For comparison with the three-dimensional results, corresponding two-dimensional computations are also performed. Detailed spatial flame profiles, solid surface profiles, and heat fluxes are presented. Increasing the flame radiation strength decreases the flame length. Although flame radiation provides an additional heat transfer mechanism to preheat the solid, it is insufficient to offset the decreased convective heating due to the shorter flame; the net effect is a slower spread rate. The percentage of unreacted fuel vapor that escapes from the flame is under 2%. It is theorized that some of the pyrolyzed fuel vapor diffuses sideway and reacts at the flame edges. A radiative energy balance is analyzed also. Flame radiative feedback to the solid plays a more important role in two-dimensional flames. With high solid fuel absorptivity, a peak in the flame spread rate occurs at an intermediate value of flame radiation strength---due to the competition between two mechanisms: gas-radiation heat loss weakening the flame and the radiative feedback boosting the solid pyrolysis. Two-dimensional calculations suggest that a larger percentage of unreacted fuel vapor can escape from the flame when the flame radiation strength is high.

The exit-time problem for a Markov jump process

NASA Astrophysics Data System (ADS)

Burch, N.; D'Elia, M.; Lehoucq, R. B.

2014-12-01

The purpose of this paper is to consider the exit-time problem for a finite-range Markov jump process, i.e, the distance the particle can jump is bounded independent of its location. Such jump diffusions are expedient models for anomalous transport exhibiting super-diffusion or nonstandard normal diffusion. We refer to the associated deterministic equation as a volume-constrained nonlocal diffusion equation. The volume constraint is the nonlocal analogue of a boundary condition necessary to demonstrate that the nonlocal diffusion equation is well-posed and is consistent with the jump process. A critical aspect of the analysis is a variational formulation and a recently developed nonlocal vector calculus. This calculus allows us to pose nonlocal backward and forward Kolmogorov equations, the former equation granting the various moments of the exit-time distribution.

Effects of curved midline and varying width on the description of the effective diffusivity of Brownian particles

NASA Astrophysics Data System (ADS)

Chávez, Yoshua; Chacón-Acosta, Guillermo; Dagdug, Leonardo

2018-05-01

Axial diffusion in channels and tubes of smoothly-varying geometry can be approximately described as one-dimensional diffusion in the entropy potential with a position-dependent effective diffusion coefficient, by means of the modified Fick–Jacobs equation. In this work, we derive analytical expressions for the position-dependent effective diffusivity for two-dimensional asymmetric varying-width channels, and for three-dimensional curved midline tubes, formed by straight walls. To this end, we use a recently developed theoretical framework using the Frenet–Serret moving frame as the coordinate system (2016 J. Chem. Phys. 145 074105). For narrow tubes and channels, an effective one-dimensional description reducing the diffusion equation to a Fick–Jacobs-like equation in general coordinates is used. From this last equation, one can calculate the effective diffusion coefficient applying Neumann boundary conditions.

General PFG signal attenuation expressions for anisotropic anomalous diffusion by modified-Bloch equations

NASA Astrophysics Data System (ADS)

Lin, Guoxing

2018-05-01

Anomalous diffusion exists widely in polymer and biological systems. Pulsed-field gradient (PFG) anomalous diffusion is complicated, especially in the anisotropic case where limited research has been reported. A general PFG signal attenuation expression, including the finite gradient pulse (FGPW) effect for free general anisotropic fractional diffusion { 0 < α , β ≤ 2 } based on the fractional derivative, has not been obtained, where α and β are time and space derivative orders. It is essential to derive a general PFG signal attenuation expression including the FGPW effect for PFG anisotropic anomalous diffusion research. In this paper, two recently developed modified-Bloch equations, the fractal differential modified-Bloch equation and the fractional integral modified-Bloch equation, were extended to obtain general PFG signal attenuation expressions for anisotropic anomalous diffusion. Various cases of PFG anisotropic anomalous diffusion were investigated, including coupled and uncoupled anisotropic anomalous diffusion. The continuous-time random walk (CTRW) simulation was also carried out to support the theoretical results. The theory and the CTRW simulation agree with each other. The obtained signal attenuation expressions and the three-dimensional fractional modified-Bloch equations are important for analyzing PFG anisotropic anomalous diffusion in NMR and MRI.

Space radiation test model study. Report for 20 May 1985-20 February 1986

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

Nightingale, R.W.; Chiu, Y.T.; Davidson, G.T.

1986-03-14

Dynamic models of the energetic populations in the outer radiation belts are being developed to better understand the extreme variations of particle flux in response to magnetospheric and solar activity. The study utilizes the SCATHA SC3 high-energy electron data, covering energies from 47 keV to 5 MeV with fine pitch-angle measurements (3 deg field of view) over the L-shell range of 5.3 to 8.7. Butter-fly distributions in the dusk sector signify particle losses due to L shell splitting of the particle-drift orbits and the subsequent scattering of the particles from the orbits by the magnetopause. To model the temporal variationsmore » and diffusion procsses of the particle populations, the data were organized into phase-space distributions, binned according to altitude (L shell), energy, pitch angle, and time. These distributions can then be mapped to the equator and plotted for fixed first and second adiabatic invariants of the inherent particle motion. A new and efficient method for calculating the third adiabatic invariant using a line integral of the relevant magnetic potential at the particle mirror points has been developed and is undergoing testing. This method will provide a useful means of displaying the radial diffusion signatures of the outer radiation belts during the more-active periods when the L shell parameter is not a good concept due to severe drift-shell splitting. The first phase of fitting the energetic-electron phase-space distributions with a combined radial and pitch-angle diffusion formulation is well underway. Bessel functions are being fit to the data in an eigenmode expansion method to determine the diffusion coefficients.« less

A parallel algorithm for nonlinear convection-diffusion equations

NASA Technical Reports Server (NTRS)

Scroggs, Jeffrey S.

1990-01-01

A parallel algorithm for the efficient solution of nonlinear time-dependent convection-diffusion equations with small parameter on the diffusion term is presented. The method is based on a physically motivated domain decomposition that is dictated by singular perturbation analysis. The analysis is used to determine regions where certain reduced equations may be solved in place of the full equation. The method is suitable for the solution of problems arising in the simulation of fluid dynamics. Experimental results for a nonlinear equation in two-dimensions are presented.

Generalized fractional diffusion equations for accelerating subdiffusion and truncated Lévy flights

NASA Astrophysics Data System (ADS)

Chechkin, A. V.; Gonchar, V. Yu.; Gorenflo, R.; Korabel, N.; Sokolov, I. M.

2008-08-01

Fractional diffusion equations are widely used to describe anomalous diffusion processes where the characteristic displacement scales as a power of time. For processes lacking such scaling the corresponding description may be given by diffusion equations with fractional derivatives of distributed order. Such equations were introduced in A. V. Chechkin, R. Gorenflo, and I. Sokolov [Phys. Rev. E 66, 046129 (2002)] for the description of the processes getting more anomalous in the course of time (decelerating subdiffusion and accelerating superdiffusion). Here we discuss the properties of diffusion equations with fractional derivatives of the distributed order for the description of anomalous relaxation and diffusion phenomena getting less anomalous in the course of time, which we call, respectively, accelerating subdiffusion and decelerating superdiffusion. For the former process, by taking a relatively simple particular example with two fixed anomalous diffusion exponents we show that the proposed equation effectively describes the subdiffusion phenomenon with diffusion exponent varying in time. For the latter process we demonstrate by a particular example how the power-law truncated Lévy stable distribution evolves in time to the distribution with power-law asymptotics and Gaussian shape in the central part. The special case of two different orders is characteristic for the general situation in which the extreme orders dominate the asymptotics.

Symmetry classification of time-fractional diffusion equation

NASA Astrophysics Data System (ADS)

Naeem, I.; Khan, M. D.

2017-01-01

In this article, a new approach is proposed to construct the symmetry groups for a class of fractional differential equations which are expressed in the modified Riemann-Liouville fractional derivative. We perform a complete group classification of a nonlinear fractional diffusion equation which arises in fractals, acoustics, control theory, signal processing and many other applications. Introducing the suitable transformations, the fractional derivatives are converted to integer order derivatives and in consequence the nonlinear fractional diffusion equation transforms to a partial differential equation (PDE). Then the Lie symmetries are computed for resulting PDE and using inverse transformations, we derive the symmetries for fractional diffusion equation. All cases are discussed in detail and results for symmetry properties are compared for different values of α. This study provides a new way of computing symmetries for a class of fractional differential equations.

Heavy-tailed fractional Pearson diffusions.

PubMed

Leonenko, N N; Papić, I; Sikorskii, A; Šuvak, N

2017-11-01

We define heavy-tailed fractional reciprocal gamma and Fisher-Snedecor diffusions by a non-Markovian time change in the corresponding Pearson diffusions. Pearson diffusions are governed by the backward Kolmogorov equations with space-varying polynomial coefficients and are widely used in applications. The corresponding fractional reciprocal gamma and Fisher-Snedecor diffusions are governed by the fractional backward Kolmogorov equations and have heavy-tailed marginal distributions in the steady state. We derive the explicit expressions for the transition densities of the fractional reciprocal gamma and Fisher-Snedecor diffusions and strong solutions of the associated Cauchy problems for the fractional backward Kolmogorov equation.

Vertical distribution of ozone: a new method of determination using satellite measurements.

PubMed

Aruga, T; Igarashi, T

1976-01-01

A new method to determine the vertical distribution of atmospheric ozone over a wide range from the spectral measurement of backscattered solar uv radiation is proposed. Equations for the diffuse reflection in an inhomogeneous atmosphere are introduced, and some theoretical approximations are discussed. An inversion equation is formulated in such a way that the change of radiance at each wavelength, caused by the minute relative increment of ozone density at each altitude, is obtained exactly. The equation is solved by an iterative procedure using the weight function obtained in this work. The results of computer simulation indicate that the ozone distribution from the mesopause to the tropopause can be determined, and that although it is impossible to suggest exactly the complicated profile with fine structure, the smoothed ozone distribution and the total content can be determined with almost the same accuracy as the accuracies of measurement and theoretical calculation of the spectral intensity.

Double Diffusive Magnetohydrodynamic (MHD) Mixed Convective Slip Flow along a Radiating Moving Vertical Flat Plate with Convective Boundary Condition

PubMed Central

Rashidi, Mohammad M.; Kavyani, Neda; Abelman, Shirley; Uddin, Mohammed J.; Freidoonimehr, Navid

2014-01-01

In this study combined heat and mass transfer by mixed convective flow along a moving vertical flat plate with hydrodynamic slip and thermal convective boundary condition is investigated. Using similarity variables, the governing nonlinear partial differential equations are converted into a system of coupled nonlinear ordinary differential equations. The transformed equations are then solved using a semi-numerical/analytical method called the differential transform method and results are compared with numerical results. Close agreement is found between the present method and the numerical method. Effects of the controlling parameters, including convective heat transfer, magnetic field, buoyancy ratio, hydrodynamic slip, mixed convective, Prandtl number and Schmidt number are investigated on the dimensionless velocity, temperature and concentration profiles. In addition effects of different parameters on the skin friction factor, , local Nusselt number, , and local Sherwood number are shown and explained through tables. PMID:25343360

Understanding Coupling of Global and Diffuse Solar Radiation with Climatic Variability

NASA Astrophysics Data System (ADS)

Hamdan, Lubna

Global solar radiation data is very important for wide variety of applications and scientific studies. However, this data is not readily available because of the cost of measuring equipment and the tedious maintenance and calibration requirements. Wide variety of models have been introduced by researchers to estimate and/or predict the global solar radiations and its components (direct and diffuse radiation) using other readily obtainable atmospheric parameters. The goal of this research is to understand the coupling of global and diffuse solar radiation with climatic variability, by investigating the relationships between these radiations and atmospheric parameters. For this purpose, we applied multilinear regression analysis on the data of National Solar Radiation Database 1991--2010 Update. The analysis showed that the main atmospheric parameters that affect the amount of global radiation received on earth's surface are cloud cover and relative humidity. Global radiation correlates negatively with both variables. Linear models are excellent approximations for the relationship between atmospheric parameters and global radiation. A linear model with the predictors total cloud cover, relative humidity, and extraterrestrial radiation is able to explain around 98% of the variability in global radiation. For diffuse radiation, the analysis showed that the main atmospheric parameters that affect the amount received on earth's surface are cloud cover and aerosol optical depth. Diffuse radiation correlates positively with both variables. Linear models are very good approximations for the relationship between atmospheric parameters and diffuse radiation. A linear model with the predictors total cloud cover, aerosol optical depth, and extraterrestrial radiation is able to explain around 91% of the variability in diffuse radiation. Prediction analysis showed that the linear models we fitted were able to predict diffuse radiation with efficiency of test adjusted R2 values equal to 0.93, using the data of total cloud cover, aerosol optical depth, relative humidity and extraterrestrial radiation. However, for prediction purposes, using nonlinear terms or nonlinear models might enhance the prediction of diffuse radiation.

Feynman-Kac equations for reaction and diffusion processes

NASA Astrophysics Data System (ADS)

Hou, Ru; Deng, Weihua

2018-04-01

This paper provides a theoretical framework for deriving the forward and backward Feynman-Kac equations for the distribution of functionals of the path of a particle undergoing both diffusion and reaction processes. Once given the diffusion type and reaction rate, a specific forward or backward Feynman-Kac equation can be obtained. The results in this paper include those for normal/anomalous diffusions and reactions with linear/nonlinear rates. Using the derived equations, we apply our findings to compute some physical (experimentally measurable) statistics, including the occupation time in half-space, the first passage time, and the occupation time in half-interval with an absorbing or reflecting boundary, for the physical system with anomalous diffusion and spontaneous evanescence.

Effects of soot absorption coefficient-Planck function correlation on radiative heat transfer in oxygen-enriched propane turbulent diffusion flame

NASA Astrophysics Data System (ADS)

Consalvi, J. L.; Nmira, F.

2016-03-01

The main objective of this article is to quantify the influence of the soot absorption coefficient-Planck function correlation on radiative loss and flame structure in an oxygen-enhanced propane turbulent diffusion flame. Calculations were run with and without accounting for this correlation by using a standard k-ε model and the steady laminar flamelet model (SLF) coupled to a joint Probability Density Function (PDF) of mixture fraction, enthalpy defect, scalar dissipation rate, and soot quantities. The PDF transport equation is solved by using a Stochastic Eulerian Field (SEF) method. The modeling of soot production is carried out by using a flamelet-based semi-empirical acetylene/benzene soot model. Radiative heat transfer is modeled by using a wide band correlated-k model and turbulent radiation interactions (TRI) are accounted for by using the Optically-Thin Fluctuation Approximation (OTFA). Predicted soot volume fraction, radiant wall heat flux distribution and radiant fraction are in good agreement with the available experimental data. Model results show that soot absorption coefficient and Planck function are negatively correlated in the region of intense soot emission. Neglecting this correlation is found to increase significantly the radiative loss leading to a substantial impact on flame structure in terms of mean and rms values of temperature. In addition mean and rms values of soot volume fraction are found to be less sensitive to the correlation than temperature since soot formation occurs mainly in a region where its influence is low.

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

NASA Astrophysics Data System (ADS)

LöWe, H.; Helbig, N.

2012-10-01

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

A Green's Function Approach to Simulate DNA Damage by the Indirect Effect

NASA Technical Reports Server (NTRS)

Plante, Ianik; Cicinotta, Francis A.

2013-01-01

The DNA damage is of fundamental importance in the understanding of the effects of ionizing radiation. DNA is damaged by the direct effect of radiation (e.g. direct ionization) and by indirect effect (e.g. damage by.OH radicals created by the radiolysis of water). Despite years of research, many questions on the DNA damage by ionizing radiation remains. In the recent years, the Green's functions of the diffusion equation (GFDE) have been used extensively in biochemistry [1], notably to simulate biochemical networks in time and space [2]. In our future work on DNA damage, we wish to use an approach based on the GFDE to refine existing models on the indirect effect of ionizing radiation on DNA. To do so, we will use the code RITRACKS [3] developed at the NASA Johnson Space Center to simulate the radiation track structure and calculate the position of radiolytic species after irradiation. We have also recently developed an efficient Monte-Carlo sampling algorithm for the GFDE of reversible reactions with an intermediate state [4], which can be modified and adapted to simulate DNA damage by free radicals. To do so, we will use the known reaction rate constants between radicals (OH, eaq, H,...) and the DNA bases, sugars and phosphates and use the sampling algorithms to simulate the diffusion of free radicals and chemical reactions with DNA. These techniques should help the understanding of the contribution of the indirect effect in the formation of DNA damage and double-strand breaks.

Numerical approximations for fractional diffusion equations via a Chebyshev spectral-tau method

NASA Astrophysics Data System (ADS)

Doha, Eid H.; Bhrawy, Ali H.; Ezz-Eldien, Samer S.

2013-10-01

In this paper, a class of fractional diffusion equations with variable coefficients is considered. An accurate and efficient spectral tau technique for solving the fractional diffusion equations numerically is proposed. This method is based upon Chebyshev tau approximation together with Chebyshev operational matrix of Caputo fractional differentiation. Such approach has the advantage of reducing the problem to the solution of a system of algebraic equations, which may then be solved by any standard numerical technique. We apply this general method to solve four specific examples. In each of the examples considered, the numerical results show that the proposed method is of high accuracy and is efficient for solving the time-dependent fractional diffusion equations.

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

Wollaeger, Ryan T.; Van Rossum, Daniel R., E-mail: wollaeger@wisc.edu, E-mail: daan@flash.uchicago.edu

Implicit Monte Carlo (IMC) and Discrete Diffusion Monte Carlo (DDMC) are methods used to stochastically solve the radiative transport and diffusion equations, respectively. These methods combine into a hybrid transport-diffusion method we refer to as IMC-DDMC. We explore a multigroup IMC-DDMC scheme that in DDMC, combines frequency groups with sufficient optical thickness. We term this procedure ''opacity regrouping''. Opacity regrouping has previously been applied to IMC-DDMC calculations for problems in which the dependence of the opacity on frequency is monotonic. We generalize opacity regrouping to non-contiguous groups and implement this in SuperNu, a code designed to do radiation transport inmore » high-velocity outflows with non-monotonic opacities. We find that regrouping of non-contiguous opacity groups generally improves the speed of IMC-DDMC radiation transport. We present an asymptotic analysis that informs the nature of the Doppler shift in DDMC groups and summarize the derivation of the Gentile-Fleck factor for modified IMC-DDMC. We test SuperNu using numerical experiments including a quasi-manufactured analytic solution, a simple 10 group problem, and the W7 problem for Type Ia supernovae. We find that opacity regrouping is necessary to make our IMC-DDMC implementation feasible for the W7 problem and possibly Type Ia supernova simulations in general. We compare the bolometric light curves and spectra produced by the SuperNu and PHOENIX radiation transport codes for the W7 problem. The overall shape of the bolometric light curves are in good agreement, as are the spectra and their evolution with time. However, for the numerical specifications we considered, we find that the peak luminosity of the light curve calculated using SuperNu is ∼10% less than that calculated using PHOENIX.« less

The exit-time problem for a Markov jump process

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

Burch, N.; D'Elia, Marta; Lehoucq, Richard B.

2014-12-15

The purpose of our paper is to consider the exit-time problem for a finite-range Markov jump process, i.e, the distance the particle can jump is bounded independent of its location. Such jump diffusions are expedient models for anomalous transport exhibiting super-diffusion or nonstandard normal diffusion. We refer to the associated deterministic equation as a volume-constrained nonlocal diffusion equation. The volume constraint is the nonlocal analogue of a boundary condition necessary to demonstrate that the nonlocal diffusion equation is well-posed and is consistent with the jump process. A critical aspect of the analysis is a variational formulation and a recently developedmore » nonlocal vector calculus. Furthermore, this calculus allows us to pose nonlocal backward and forward Kolmogorov equations, the former equation granting the various moments of the exit-time distribution.« less

Nonlinear anomalous diffusion equation and fractal dimension: exact generalized Gaussian solution.

PubMed

Pedron, I T; Mendes, R S; Malacarne, L C; Lenzi, E K

2002-04-01

In this work we incorporate, in a unified way, two anomalous behaviors, the power law and stretched exponential ones, by considering the radial dependence of the N-dimensional nonlinear diffusion equation partial differential rho/ partial differential t=nabla.(Knablarho(nu))-nabla.(muFrho)-alpharho, where K=Dr(-theta), nu, theta, mu, and D are real parameters, F is the external force, and alpha is a time-dependent source. This equation unifies the O'Shaughnessy-Procaccia anomalous diffusion equation on fractals (nu=1) and the spherical anomalous diffusion for porous media (theta=0). An exact spherical symmetric solution of this nonlinear Fokker-Planck equation is obtained, leading to a large class of anomalous behaviors. Stationary solutions for this Fokker-Planck-like equation are also discussed by introducing an effective potential.

Diffusion approximations to the chemical master equation only have a consistent stochastic thermodynamics at chemical equilibrium

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.

Diffusion approximations to the chemical master equation only have a consistent stochastic thermodynamics at chemical equilibrium.

PubMed

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.

Green functions and Langevin equations for nonlinear diffusion equations: A comment on ‘Markov processes, Hurst exponents, and nonlinear diffusion equations’ by Bassler et al.

NASA Astrophysics Data System (ADS)

Frank, T. D.

2008-02-01

We discuss two central claims made in the study by Bassler et al. [K.E. Bassler, G.H. Gunaratne, J.L. McCauley, Physica A 369 (2006) 343]. Bassler et al. claimed that Green functions and Langevin equations cannot be defined for nonlinear diffusion equations. In addition, they claimed that nonlinear diffusion equations are linear partial differential equations disguised as nonlinear ones. We review bottom-up and top-down approaches that have been used in the literature to derive Green functions for nonlinear diffusion equations and, in doing so, show that the first claim needs to be revised. We show that the second claim as well needs to be revised. To this end, we point out similarities and differences between non-autonomous linear Fokker-Planck equations and autonomous nonlinear Fokker-Planck equations. In this context, we raise the question whether Bassler et al.’s approach to financial markets is physically plausible because it necessitates the introduction of external traders and causes. Such external entities can easily be eliminated when taking self-organization principles and concepts of nonextensive thermostatistics into account and modeling financial processes by means of nonlinear Fokker-Planck equations.

CFD simulation of simultaneous monotonic cooling and surface heat transfer coefficient

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

Mihálka, Peter, E-mail: usarmipe@savba.sk; Matiašovský, Peter, E-mail: usarmat@savba.sk

The monotonic heating regime method for determination of thermal diffusivity is based on the analysis of an unsteady-state (stabilised) thermal process characterised by an independence of the space-time temperature distribution on initial conditions. At the first kind of the monotonic regime a sample of simple geometry is heated / cooled at constant ambient temperature. The determination of thermal diffusivity requires the determination rate of a temperature change and simultaneous determination of the first eigenvalue. According to a characteristic equation the first eigenvalue is a function of the Biot number defined by a surface heat transfer coefficient and thermal conductivity ofmore » an analysed material. Knowing the surface heat transfer coefficient and the first eigenvalue the thermal conductivity can be determined. The surface heat transport coefficient during the monotonic regime can be determined by the continuous measurement of long-wave radiation heat flow and the photoelectric measurement of the air refractive index gradient in a boundary layer. CFD simulation of the cooling process was carried out to analyse local convective and radiative heat transfer coefficients more in detail. Influence of ambient air flow was analysed. The obtained eigenvalues and corresponding surface heat transfer coefficient values enable to determine thermal conductivity of the analysed specimen together with its thermal diffusivity during a monotonic heating regime.« less

Fully coupled simulation of cosmic reionization. I. numerical methods and tests

DOE PAGES

Norman, Michael L.; Reynolds, Daniel R.; So, Geoffrey C.; ...

2015-01-09

Here, we describe an extension of the Enzo code to enable fully coupled radiation hydrodynamical simulation of inhomogeneous reionization in large similar to(100 Mpc)(3) cosmological volumes with thousands to millions of point sources. We solve all dynamical, radiative transfer, thermal, and ionization processes self-consistently on the same mesh, as opposed to a postprocessing approach which coarse-grains the radiative transfer. But, we employ a simple subgrid model for star formation which we calibrate to observations. The numerical method presented is a modification of an earlier method presented in Reynolds et al. differing principally in the operator splitting algorithm we use tomore » advance the system of equations. Radiation transport is done in the gray flux-limited diffusion (FLD) approximation, which is solved by implicit time integration split off from the gas energy and ionization equations, which are solved separately. This results in a faster and more robust scheme for cosmological applications compared to the earlier method. The FLD equation is solved using the hypre optimally scalable geometric multigrid solver from LLNL. By treating the ionizing radiation as a grid field as opposed to rays, our method is scalable with respect to the number of ionizing sources, limited only by the parallel scaling properties of the radiation solver. We test the speed and accuracy of our approach on a number of standard verification and validation tests. We show by direct comparison with Enzo's adaptive ray tracing method Moray that the well-known inability of FLD to cast a shadow behind opaque clouds has a minor effect on the evolution of ionized volume and mass fractions in a reionization simulation validation test. Finally, we illustrate an application of our method to the problem of inhomogeneous reionization in a 80 Mpc comoving box resolved with 3200(3) Eulerian grid cells and dark matter particles.« less

Diffusion Influenced Adsorption Kinetics.

PubMed

Miura, Toshiaki; Seki, Kazuhiko

2015-08-27

When the kinetics of adsorption is influenced by the diffusive flow of solutes, the solute concentration at the surface is influenced by the surface coverage of solutes, which is given by the Langmuir-Hinshelwood adsorption equation. The diffusion equation with the boundary condition given by the Langmuir-Hinshelwood adsorption equation leads to the nonlinear integro-differential equation for the surface coverage. In this paper, we solved the nonlinear integro-differential equation using the Grünwald-Letnikov formula developed to solve fractional kinetics. Guided by the numerical results, analytical expressions for the upper and lower bounds of the exact numerical results were obtained. The upper and lower bounds were close to the exact numerical results in the diffusion- and reaction-controlled limits, respectively. We examined the validity of the two simple analytical expressions obtained in the diffusion-controlled limit. The results were generalized to include the effect of dispersive diffusion. We also investigated the effect of molecular rearrangement of anisotropic molecules on surface coverage.

Analytical solutions of the space-time fractional Telegraph and advection-diffusion equations

NASA Astrophysics Data System (ADS)

Tawfik, Ashraf M.; Fichtner, Horst; Schlickeiser, Reinhard; Elhanbaly, A.

2018-02-01

The aim of this paper is to develop a fractional derivative model of energetic particle transport for both uniform and non-uniform large-scale magnetic field by studying the fractional Telegraph equation and the fractional advection-diffusion equation. Analytical solutions of the space-time fractional Telegraph equation and space-time fractional advection-diffusion equation are obtained by use of the Caputo fractional derivative and the Laplace-Fourier technique. The solutions are given in terms of Fox's H function. As an illustration they are applied to the case of solar energetic particles.

Boundary value problems for multi-term fractional differential equations

NASA Astrophysics Data System (ADS)

Daftardar-Gejji, Varsha; Bhalekar, Sachin

2008-09-01

Multi-term fractional diffusion-wave equation along with the homogeneous/non-homogeneous boundary conditions has been solved using the method of separation of variables. It is observed that, unlike in the one term case, solution of multi-term fractional diffusion-wave equation is not necessarily non-negative, and hence does not represent anomalous diffusion of any kind.

A finite element formulation preserving symmetric and banded diffusion stiffness matrix characteristics for fractional differential equations

NASA Astrophysics Data System (ADS)

Lin, Zeng; Wang, Dongdong

2017-10-01

Due to the nonlocal property of the fractional derivative, the finite element analysis of fractional diffusion equation often leads to a dense and non-symmetric stiffness matrix, in contrast to the conventional finite element formulation with a particularly desirable symmetric and banded stiffness matrix structure for the typical diffusion equation. This work first proposes a finite element formulation that preserves the symmetry and banded stiffness matrix characteristics for the fractional diffusion equation. The key point of the proposed formulation is the symmetric weak form construction through introducing a fractional weight function. It turns out that the stiffness part of the present formulation is identical to its counterpart of the finite element method for the conventional diffusion equation and thus the stiffness matrix formulation becomes trivial. Meanwhile, the fractional derivative effect in the discrete formulation is completely transferred to the force vector, which is obviously much easier and efficient to compute than the dense fractional derivative stiffness matrix. Subsequently, it is further shown that for the general fractional advection-diffusion-reaction equation, the symmetric and banded structure can also be maintained for the diffusion stiffness matrix, although the total stiffness matrix is not symmetric in this case. More importantly, it is demonstrated that under certain conditions this symmetric diffusion stiffness matrix formulation is capable of producing very favorable numerical solutions in comparison with the conventional non-symmetric diffusion stiffness matrix finite element formulation. The effectiveness of the proposed methodology is illustrated through a series of numerical examples.

Albedo and flux extinction coefficient of impure snow for diffuse shortwave radiation

NASA Technical Reports Server (NTRS)

Choudhury, B. J.; Mo, T.; Wang, J. R.; Chang, A. T. C.

1981-01-01

Impurities enter a snowpack as a result of fallout of scavenging by falling snow crystals. Albedo and flux extinction coefficient of soot contaminated snowcovers were studied using a two stream approximation of the radiative transfer equation. The effect of soot was calculated by two methods: independent scattering by ice grains and impurities and average refractive index for ice grains. Both methods predict a qualitatively similar effect of soot; the albedo is decreased and the extinction coefficient is increased compared to that for pure snow in the visible region; the infrared properties are largely unaffected. Quantitatively, however, the effect of soot is more pronounced in the average refractive index method. Soot contamination provides a qualitative explanation for several snow observations.

Diffusion coefficients in organic-water solutions and comparison with Stokes-Einstein predictions

NASA Astrophysics Data System (ADS)

Evoy, E.; Kamal, S.; Bertram, A. K.

2017-12-01

Diffusion coefficients of organic species in particles containing secondary organic material (SOM) are necessary for predicting the growth and reactivity of these particles in the atmosphere. Previously, the Stokes-Einstein equation combined with viscosity measurements have been used to predict these diffusion coefficients. However, the accuracy of the Stokes-Einstein equation for predicting diffusion coefficients in SOM-water particles has not been quantified. To test the Stokes-Einstein equation, diffusion coefficients of fluorescent organic probe molecules were measured in citric acid-water and sorbitol-water solutions. These solutions were used as proxies for SOM-water particles found in the atmosphere. Measurements were performed as a function of water activity, ranging from 0.26-0.86, and as a function of viscosity ranging from 10-3 to 103 Pa s. Diffusion coefficients were measured using fluorescence recovery after photobleaching. The measured diffusion coefficients were compared with predictions made using the Stokes-Einstein equation combined with literature viscosity data. Within the uncertainties of the measurements, the measured diffusion coefficients agreed with the predicted diffusion coefficients, in all cases.

Diffusion approximations to the chemical master equation only have a consistent stochastic thermodynamics at chemical equilibrium

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

Horowitz, Jordan M., E-mail: jordan.horowitz@umb.edu

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 stochasticmore » thermodynamics, we show that one can still use the diffusive solutions to approximate the underlying thermodynamics of the chemical master equation.« less

S{sub 2}SA preconditioning for the S{sub n} equations with strictly non negative spatial discretization

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

Bruss, D. E.; Morel, J. E.; Ragusa, J. C.

2013-07-01

Preconditioners based upon sweeps and diffusion-synthetic acceleration have been constructed and applied to the zeroth and first spatial moments of the 1-D S{sub n} transport equation using a strictly non negative nonlinear spatial closure. Linear and nonlinear preconditioners have been analyzed. The effectiveness of various combinations of these preconditioners are compared. In one dimension, nonlinear sweep preconditioning is shown to be superior to linear sweep preconditioning, and DSA preconditioning using nonlinear sweeps in conjunction with a linear diffusion equation is found to be essentially equivalent to nonlinear sweeps in conjunction with a nonlinear diffusion equation. The ability to use amore » linear diffusion equation has important implications for preconditioning the S{sub n} equations with a strictly non negative spatial discretization in multiple dimensions. (authors)« less

Influence of Diffused Solar Radiation on the Solar Concentrating System of a Plant Shoot Configuration

NASA Astrophysics Data System (ADS)

Obara, Shin'ya

Investigation of a plant shoot configuration is used to obtain valuable information concerning the received light system. Additionally, analysis results concerning a plant shoot configuration interaction with direct solar radiation were taken from a past study. However, in order to consider a plant shoot as a received sunlight system, it is necessary to understand the received light characteristics of both direct solar radiation and diffused solar radiation. Under a clear sky, the ratio of direct solar radiation to diffused solar radiation is large. However, under a clouded sky, the amount of diffused solar radiation becomes larger. Therefore, in this paper, we investigate the received light characteristics of a plant shoot configuration under the influence of diffused solar radiation. As a result, we clarify the relationship between the amount of diffused solar radiation and the amount of received light as a function of the characteristics of the plant shoot configuration. In order to obtain diffused solar radiation, it is necessary to correspond to the radiation of the multi-directions. In the analysis, the characteristic of the difference in arrangement of the top leaf and the other leaf was obtained. Therefore, in analysis, leaves other than the top were distributed in the wide range.

Radiant exchange in partially specular architectural environments

NASA Astrophysics Data System (ADS)

Beamer, C. Walter; Muehleisen, Ralph T.

2003-10-01

The radiant exchange method, also known as radiosity, was originally developed for thermal radiative heat transfer applications. Later it was used to model architectural lighting systems, and more recently it has been extended to model acoustic systems. While there are subtle differences in these applications, the basic method is based on solving a system of energy balance equations, and it is best applied to spaces with mainly diffuse reflecting surfaces. The obvious drawback to this method is that it is based around the assumption that all surfaces in the system are diffuse reflectors. Because almost all architectural systems have at least some partially specular reflecting surfaces in the system it is important to extend the radiant exchange method to deal with this type of surface reflection. [Work supported by NSF.

Unsteady magnetohydrodynamics mixed convection flow in a rotating medium with double diffusion

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

Jiann, Lim Yeou; Ismail, Zulkhibri; Khan, Ilyas

2015-05-15

Exact solutions of an unsteady Magnetohydrodynamics (MHD) flow over an impulsively started vertical plate in a rotating medium are presented. The effects of thermal radiative and thermal diffusion on the fluid flow are also considered. The governing equations are modelled and solved for velocity, temperature and concentration using Laplace transforms technique. Expressions of velocity, temperature and concentration profiles are obtained and their numerical results are presented graphically. Skin friction, Sherwood number and Nusselt number are also computed and presented in tabular forms. The determined solutions can generate a large class of solutions as special cases corresponding to different motions withmore » technical relevance. The results obtained herein may be used to verify the validation of obtained numerical solutions for more complicated fluid flow problems.« less

Lévy/Anomalous Diffusion as a Mean-Field Theory for 3D Cloud Effects in Shortwave Radiative Transfer: Empirical Support, New Analytical Formulation, and Impact on Atmospheric Absorption

NASA Astrophysics Data System (ADS)

Buldyrev, S.; Davis, A.; Marshak, A.; Stanley, H. E.

2001-12-01

Two-stream radiation transport models, as used in all current GCM parameterization schemes, are mathematically equivalent to ``standard'' diffusion theory where the physical picture is a slow propagation of the diffuse radiation by Gaussian random walks. The space/time spread (technically, the Green function) of this diffusion process is described exactly by a Gaussian distribution; from the statistical physics viewpoint, this follows from the convergence of the sum of many (rescaled) steps between scattering events with a finite variance. This Gaussian picture follows directly from first principles (the radiative transfer equation) under the assumptions of horizontal uniformity and large optical depth, i.e., there is a homogeneous plane-parallel cloud somewhere in the column. The first-order effect of 3D variability of cloudiness, the main source of scattering, is to perturb the distribution of single steps between scatterings which, modulo the ``1-g'' rescaling, can be assumed effectively isotropic. The most natural generalization of the Gaussian distribution is the 1-parameter family of symmetric Lévy-stable distributions because the sum of many zero-mean random variables with infinite variance, but finite moments of order q < α (0 < α < 2), converge to them. It has been shown on heuristic grounds that for these Lévy-based random walks the typical number of scatterings is now (1-g)τ α for transmitted light. The appearance of a non-rational exponent is why this is referred to as ``anomalous'' diffusion. Note that standard/Gaussian diffusion is retrieved in the limit α = 2-. Lévy transport theory has been successfully used in the statistical physics literature to investigate a wide variety of systems with strongly nonlinear dynamics; these applications range from random advection in turbulent fluids to the erratic behavior of financial time-series and, most recently, self-regulating ecological systems. We will briefly survey the state-of-the-art observations that offer compelling empirical support for the Lévy/anomalous diffusion model in atmospheric radiation: (1) high-resolution spectroscopy of differential absorption in the O2 A-band from ground; (2) temporal transient records of lightning strokes transmitted through clouds to a sensitive detector in space; and (3) the Gamma-distributions of optical depths derived from Landsat cloud scenes at 30-m resolution. We will then introduce a rigorous analytical formulation of Lévy/anomalous transport through finite media based on fractional derivatives and Sonin calculus. A remarkable result from this new theoretical development is an extremal property of the α = 1+ case (divergent mean-free-path), as is observed in the cloudy atmosphere. Finally, we will discuss the implications of anomalous transport theory for bulk 3D effects on the current enhanced absorption problem as well as its role as the basis of a next-generation GCM radiation parameterization.

Detailed in situ laser calibration of the infrared imaging video bolometer for the JT-60U tokamak

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

Parchamy, H.; Peterson, B. J.; Konoshima, S.

2006-10-15

The infrared imaging video bolometer (IRVB) in JT-60U includes a single graphite-coated gold foil with an effective area of 9x7 cm{sup 2} and a thickness of 2.5 {mu}m. The thermal images of the foil resulting from the plasma radiation are provided by an IR camera. The calibration technique of the IRVB gives confidence in the absolute levels of the measured values of the plasma radiation. The in situ calibration is carried out in order to obtain local foil properties such as the thermal diffusivity {kappa} and the product of the thermal conductivity k and the thickness t{sub f} of themore » foil. These quantities are necessary for solving the two-dimensional heat diffusion equation of the foil which is used in the experiments. These parameters are determined by comparing the measured temperature profiles (for kt{sub f}) and their decays (for {kappa}) with the corresponding results of a finite element model using the measured HeNe laser power profile as a known radiation power source. The infrared camera (Indigo/Omega) is calibrated by fitting the temperature rise of a heated plate to the resulting camera data using the Stefan-Boltzmann law.« less

The Interface Influence in TiN/SiN x Multilayer Nanocomposite Under Irradiation

NASA Astrophysics Data System (ADS)

Uglov, V. V.; Safronov, I. V.; Kvasov, N. T.; Remnev, G. E.; Shimanski, V. I.

2018-01-01

The paper focuses on studying the kinetics of radiation-induced point defects formed in TiN/SiN x multilayer nanocomposites with account of their generation, diffusion recombination, and the influence of sinks functioning as interfaces. In order to describe the kinetics in nanocrystalline TiN and amorphous SiN x phases, a finite-difference method is used to solve the system of balance kinetic equations for absolute defect concentrations depending on the spatiotemporal variables. A model of the disclination-dislocation interface structure is used to study the absorption of radiation-induced point defects on the boundaries in created stress fields. It is shown that the interface effectively absorbs point defects in these phases of TiN/SiN x multilayer nanocomposite, thereby reducing their amount within the space between phases. This behavior of point defects partially explains a mechanism of the radiation resistance in this type of nanocomposites.

Group iterative methods for the solution of two-dimensional time-fractional diffusion equation

NASA Astrophysics Data System (ADS)

Balasim, Alla Tareq; Ali, Norhashidah Hj. Mohd.

2016-06-01

Variety of problems in science and engineering may be described by fractional partial differential equations (FPDE) in relation to space and/or time fractional derivatives. The difference between time fractional diffusion equations and standard diffusion equations lies primarily in the time derivative. Over the last few years, iterative schemes derived from the rotated finite difference approximation have been proven to work well in solving standard diffusion equations. However, its application on time fractional diffusion counterpart is still yet to be investigated. In this paper, we will present a preliminary study on the formulation and analysis of new explicit group iterative methods in solving a two-dimensional time fractional diffusion equation. These methods were derived from the standard and rotated Crank-Nicolson difference approximation formula. Several numerical experiments were conducted to show the efficiency of the developed schemes in terms of CPU time and iteration number. At the request of all authors of the paper an updated version of this article was published on 7 July 2016. The original version supplied to AIP Publishing contained an error in Table 1 and References 15 and 16 were incomplete. These errors have been corrected in the updated and republished article.

Equatorial superrotation in a thermally driven zonally symmetric circulation

NASA Technical Reports Server (NTRS)

Mayr, H. G.; Harris, I.

1981-01-01

Near the equator where the Coriolis force vanishes, the momentum balance for the axially symmetric circulation is established between horizontal and vertical diffusion, which, a priori, does not impose constraints on the direction or magnitude of the zonal winds. Solar radiation absorbed at low latitudes is a major force in driving large scale motions with air rising near the equator and falling at higher latitudes. In the upper leg of the meridional cell, angular momentum is redistributed so that the atmosphere tends to subrotate (or corotate) at low latitudes and superrotate at high latitudes. In the lower leg, however, the process is reversed and produces a tendency for the equatorial region to superrotate. The outcome depends on the energy budget which is closely coupled to the momentum budget through the thermal wind equation; a pressure (temperature) maximum is required to sustain equatorial superrotation. Such a condition arises in regions which are convectively unstable and the temperature lapse rate is superadiabatic. It should arise in the tropospheres of Jupiter and Saturn; planetary energy from the interior is carried to higher altitudes where radiation to space becomes important. Upward equatorial motions in the direct and indirect circulations (Ferrel-Thomson type) imposed by insolation can then trap dynamic energy for equatorial heating which can sustain the superrotation of the equatorial region.

Earth Reflected Solar Radiation Incident upon an Arbitrarily Oriented Spinning Flat Plate

NASA Technical Reports Server (NTRS)

Cunningham, Fred G.

1963-01-01

A general derivation is given for the earth reflected solar radiation input to a flat plate--a solar cell paddle, for example--which is spinning about an axis coincident with the axis of symmetry of the satellite to which it is affixed. The resulting equations are written for the general case so that arbitrary orientations of the spin axis with respect to the earth-satellite line and arbitrary orientations of the normal to the plate with respect to the spin axis can be treated. No attempt is made to perform the resulting integrations because of the complexity of the equations; nor is there any attempt to delineate the integration limits for the general case. However, the equations governing these limits are given. The appendixes contain: the results, in graphical form, of two representative examples; the general computer program for the calculation is given in Fortran notation; and the results of a calculation of the distribution of albedo energy on the proposed Echo II satellite. The value of the mean solar constant used is 1.395 times 10 (sup 4) ergs per centimeters-squared per second; the mean albedo of the earth is assumed to be 0.34; and the earth is assumed to be a diffuse reflector.

On the anisotropic advection-diffusion equation with time dependent coefficients

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

Hernandez-Coronado, Hector; Coronado, Manuel; Del-Castillo-Negrete, Diego B.

The advection-diffusion equation with time dependent velocity and anisotropic time dependent diffusion tensor is examined in regard to its non-classical transport features and to the use of a non-orthogonal coordinate system. Although this equation appears in diverse physical problems, particularly in particle transport in stochastic velocity fields and in underground porous media, a detailed analysis of its solutions is lacking. In order to study the effects of the time-dependent coefficients and the anisotropic diffusion on transport, we solve analytically the equation for an initial Dirac delta pulse. Here, we discuss the solutions to three cases: one based on power-law correlationmore » functions where the pulse diffuses faster than the classical rate ~t, a second case specically designed to display slower rate of diffusion than the classical one, and a third case to describe hydrodynamic dispersion in porous media« less

On the anisotropic advection-diffusion equation with time dependent coefficients

DOE PAGES

Hernandez-Coronado, Hector; Coronado, Manuel; Del-Castillo-Negrete, Diego B.

2017-02-01

The advection-diffusion equation with time dependent velocity and anisotropic time dependent diffusion tensor is examined in regard to its non-classical transport features and to the use of a non-orthogonal coordinate system. Although this equation appears in diverse physical problems, particularly in particle transport in stochastic velocity fields and in underground porous media, a detailed analysis of its solutions is lacking. In order to study the effects of the time-dependent coefficients and the anisotropic diffusion on transport, we solve analytically the equation for an initial Dirac delta pulse. Here, we discuss the solutions to three cases: one based on power-law correlationmore » functions where the pulse diffuses faster than the classical rate ~t, a second case specically designed to display slower rate of diffusion than the classical one, and a third case to describe hydrodynamic dispersion in porous media« less

General pulsed-field gradient signal attenuation expression based on a fractional integral modified-Bloch equation

NASA Astrophysics Data System (ADS)

Lin, Guoxing

2018-10-01

Anomalous diffusion has been investigated in many polymer and biological systems. The analysis of PFG anomalous diffusion relies on the ability to obtain the signal attenuation expression. However, the general analytical PFG signal attenuation expression based on the fractional derivative has not been previously reported. Additionally, the reported modified-Bloch equations for PFG anomalous diffusion in the literature yielded different results due to their different forms. Here, a new integral type modified-Bloch equation based on the fractional derivative for PFG anomalous diffusion is proposed, which is significantly different from the conventional differential type modified-Bloch equation. The merit of the integral type modified-Bloch equation is that the original properties of the contributions from linear or nonlinear processes remain unchanged at the instant of the combination. From the modified-Bloch equation, the general solutions are derived, which includes the finite gradient pulse width (FGPW) effect. The numerical evaluation of these PFG signal attenuation expressions can be obtained either by the Adomian decomposition, or a direct integration method that is fast and practicable. The theoretical results agree with the continuous-time random walk (CTRW) simulations performed in this paper. Additionally, the relaxation effect in PFG anomalous diffusion is found to be different from that in PFG normal diffusion. The new modified-Bloch equations and their solutions provide a fundamental tool to analyze PFG anomalous diffusion in nuclear magnetic resonance (NMR) and magnetic resonance imaging (MRI).

Analysis of activation energy in Couette-Poiseuille flow of nanofluid in the presence of chemical reaction and convective boundary conditions

NASA Astrophysics Data System (ADS)

Zeeshan, A.; Shehzad, N.; Ellahi, R.

2018-03-01

The motivation of the current article is to explore the energy activation in MHD radiative Couette-Poiseuille flow nanofluid in horizontal channel with convective boundary conditions. The mathematical model of Buongiorno [1] effectively describes the current flow analysis. Additionally, the impact of chemical reaction is also taken in account. The governing flow equations are simplified with the help of boundary layer approximations. Non-linear coupled equations for momentum, energy and mass transfer are tackled with analytical (HAM) technique. The influence of dimensionless convergence parameter like Brownian motion parameter, radiation parameter, buoyancy ratio parameter, dimensionless activation energy, thermophoresis parameter, temperature difference parameter, dimensionless reaction rate, Schmidt number, Brinkman number, Biot number and convection diffusion parameter on velocity, temperature and concentration profiles are discussed graphically and in tabular form. From the results, it is elaborate that the nanoparticle concentration is directly proportional to the chemical reaction with activation energy and the performance of Brownian motion on nanoparticle concentration gives reverse pattern to that of thermophoresis parameter.

Radiative transfer calculated from a Markov chain formalism

NASA Technical Reports Server (NTRS)

Esposito, L. W.; House, L. L.

1978-01-01

The theory of Markov chains is used to formulate the radiative transport problem in a general way by modeling the successive interactions of a photon as a stochastic process. Under the minimal requirement that the stochastic process is a Markov chain, the determination of the diffuse reflection or transmission from a scattering atmosphere is equivalent to the solution of a system of linear equations. This treatment is mathematically equivalent to, and thus has many of the advantages of, Monte Carlo methods, but can be considerably more rapid than Monte Carlo algorithms for numerical calculations in particular applications. We have verified the speed and accuracy of this formalism for the standard problem of finding the intensity of scattered light from a homogeneous plane-parallel atmosphere with an arbitrary phase function for scattering. Accurate results over a wide range of parameters were obtained with computation times comparable to those of a standard 'doubling' routine. The generality of this formalism thus allows fast, direct solutions to problems that were previously soluble only by Monte Carlo methods. Some comparisons are made with respect to integral equation methods.

Diffusion-enhanced Förster resonance energy transfer and the effects of external quenchers and the donor quantum yield.

PubMed

Jacob, Maik H; Dsouza, Roy N; Ghosh, Indrajit; Norouzy, Amir; Schwarzlose, Thomas; Nau, Werner M

2013-01-10

The structural and dynamic properties of a flexible peptidic chain codetermine its biological activity. These properties are imprinted in intrachain site-to-site distances as well as in diffusion coefficients of mutual site-to-site motion. Both distance distribution and diffusion determine the extent of Förster resonance energy transfer (FRET) between two chain sites labeled with a FRET donor and acceptor. Both could be obtained from time-resolved FRET measurements if their individual contributions to the FRET efficiency could be systematically varied. Because the FRET diffusion enhancement (FDE) depends on the donor-fluorescence lifetime, it has been proposed that the FDE can be reduced by shortening the donor lifetime through an external quencher. Benefiting from the high diffusion sensitivity of short-distance FRET, we tested this concept experimentally on a (Gly-Ser)(6) segment labeled with the donor/acceptor pair naphthylalanine/2,3-diazabicyclo[2.2.2]oct-2-ene (NAla/Dbo). Surprisingly, the very effective quencher potassium iodide (KI) had no effect at all on the average donor-acceptor distance, although the donor lifetime was shortened from ca. 36 ns in the absence of KI to ca. 3 ns in the presence of 30 mM KI. We show that the proposed approach had to fail because it is not the experimentally observed but the radiative donor lifetime that controls the FDE. Because of that, any FRET ensemble measurement can easily underestimate diffusion and might be misleading even if it employs the Haas-Steinberg diffusion equation (HSE). An extension of traditional FRET analysis allowed us to evaluate HSE simulations and to corroborate as well as generalize the experimental results. We demonstrate that diffusion-enhanced FRET depends on the radiative donor lifetime as it depends on the diffusion coefficient, a useful symmetry that can directly be applied to distinguish dynamic and structural effects of viscous cosolvents on the polymer chain. We demonstrate that the effective FRET rate and the recovered donor-acceptor distance depend on the quantum yield, most strongly in the absence of diffusion, which has to be accounted for in the interpretation of distance trends monitored by FRET.

ANALYTICAL SOLUTIONS OF THE ATMOSPHERIC DIFFUSION EQUATION WITH MULTIPLE SOURCES AND HEIGHT-DEPENDENT WIND SPEED AND EDDY DIFFUSIVITIES. (R825689C072)

EPA Science Inventory

Abstract

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...

ANALYTICAL SOLUTIONS OF THE ATMOSPHERIC DIFFUSION EQUATION WITH MULTIPLE SOURCES AND HEIGHT-DEPENDENT WIND SPEED AND EDDY DIFFUSIVITIES. (R825689C048)

EPA Science Inventory

Abstract

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...

A fast semi-discrete Kansa method to solve the two-dimensional spatiotemporal fractional diffusion equation

NASA Astrophysics Data System (ADS)

Sun, HongGuang; Liu, Xiaoting; Zhang, Yong; Pang, Guofei; Garrard, Rhiannon

2017-09-01

Fractional-order diffusion equations (FDEs) extend classical diffusion equations by quantifying anomalous diffusion frequently observed in heterogeneous media. Real-world diffusion can be multi-dimensional, requiring efficient numerical solvers that can handle long-term memory embedded in mass transport. To address this challenge, a semi-discrete Kansa method is developed to approximate the two-dimensional spatiotemporal FDE, where the Kansa approach first discretizes the FDE, then the Gauss-Jacobi quadrature rule solves the corresponding matrix, and finally the Mittag-Leffler function provides an analytical solution for the resultant time-fractional ordinary differential equation. Numerical experiments are then conducted to check how the accuracy and convergence rate of the numerical solution are affected by the distribution mode and number of spatial discretization nodes. Applications further show that the numerical method can efficiently solve two-dimensional spatiotemporal FDE models with either a continuous or discrete mixing measure. Hence this study provides an efficient and fast computational method for modeling super-diffusive, sub-diffusive, and mixed diffusive processes in large, two-dimensional domains with irregular shapes.

The First National Student Conference: NASA University Research Centers at Minority Institutions

NASA Technical Reports Server (NTRS)

Daso, Endwell O. (Editor); Mebane, Stacie (Editor)

1997-01-01

The conference includes contributions from 13 minority universities with NASA University Research Centers. Topics discussed include: leadership, survival strategies, life support systems, food systems, simulated hypergravity, chromium diffusion doping, radiation effects on dc-dc converters, metal oxide glasses, crystal growth of Bil3, science and communication on wheels, semiconductor thin films, numerical solution of random algebraic equations, fuzzy logic control, spatial resolution of satellite images, programming language development, nitric oxide in the thermosphere and mesosphere, high performance polyimides, crossover control in genetic algorithms, hyperthermal ion scattering, etc.

Verification of reflectance models in turbid waters

NASA Technical Reports Server (NTRS)

Tanis, F. J.; Lyzenga, D. R.

1981-01-01

Inherent optical parameters of very turbid waters were used to evaluate existing water reflectance models. Measured upwelling radiance spectra and Monte Carlo simulations of the radiative transfer equations were compared with results from models based upon two flow, quasi-single scattering, augmented isotropic scattering, and power series approximation. Each model was evaluated for three separate components of upwelling radiance: (1) direct sunlight; (2) diffuse skylight; and (3) internally reflected light. Limitations of existing water reflectance models as applied to turbid waters and possible applications to the extraction of water constituent information are discussed.

Non-equilibrium radiation from viscous chemically reacting two-phase exhaust plumes

NASA Technical Reports Server (NTRS)

Penny, M. M.; Smith, S. D.; Mikatarian, R. R.; Ring, L. R.; Anderson, P. G.

1976-01-01

A knowledge of the structure of the rocket exhaust plumes is necessary to solve problems involving plume signatures, base heating, plume/surface interactions, etc. An algorithm is presented which treats the viscous flow of multiphase chemically reacting fluids in a two-dimensional or axisymmetric supersonic flow field. The gas-particle flow solution is fully coupled with the chemical kinetics calculated using an implicit scheme to calculate chemical production rates. Viscous effects include chemical species diffusion with the viscosity coefficient calculated using a two-equation turbulent kinetic energy model.

A new paradigm for predicting zonal-mean climate and climate change

NASA Astrophysics Data System (ADS)

Armour, K.; Roe, G.; Donohoe, A.; Siler, N.; Markle, B. R.; Liu, X.; Feldl, N.; Battisti, D. S.; Frierson, D. M.

2016-12-01

How will the pole-to-equator temperature gradient, or large-scale patterns of precipitation, change under global warming? Answering such questions typically involves numerical simulations with comprehensive general circulation models (GCMs) that represent the complexities of climate forcing, radiative feedbacks, and atmosphere and ocean dynamics. Yet, our understanding of these predictions hinges on our ability to explain them through the lens of simple models and physical theories. Here we present evidence that zonal-mean climate, and its changes, can be understood in terms of a moist energy balance model that represents atmospheric heat transport as a simple diffusion of latent and sensible heat (as a down-gradient transport of moist static energy, with a diffusivity coefficient that is nearly constant with latitude). We show that the theoretical underpinnings of this model derive from the principle of maximum entropy production; that its predictions are empirically supported by atmospheric reanalyses; and that it successfully predicts the behavior of a hierarchy of climate models - from a gray radiation aquaplanet moist GCM, to comprehensive GCMs participating in CMIP5. As an example of the power of this paradigm, we show that, given only patterns of local radiative feedbacks and climate forcing, the moist energy balance model accurately predicts the evolution of zonal-mean temperature and atmospheric heat transport as simulated by the CMIP5 ensemble. These results suggest that, despite all of its dynamical complexity, the atmosphere essentially responds to energy imbalances by simply diffusing latent and sensible heat down-gradient; this principle appears to explain zonal-mean climate and its changes under global warming.

Existence and Stability of Traveling Waves for Degenerate Reaction-Diffusion Equation with Time Delay

NASA Astrophysics Data System (ADS)

Huang, Rui; Jin, Chunhua; Mei, Ming; Yin, Jingxue

2018-01-01

This paper deals with the existence and stability of traveling wave solutions for a degenerate reaction-diffusion equation with time delay. The degeneracy of spatial diffusion together with the effect of time delay causes us the essential difficulty for the existence of the traveling waves and their stabilities. In order to treat this case, we first show the existence of smooth- and sharp-type traveling wave solutions in the case of c≥c^* for the degenerate reaction-diffusion equation without delay, where c^*>0 is the critical wave speed of smooth traveling waves. Then, as a small perturbation, we obtain the existence of the smooth non-critical traveling waves for the degenerate diffusion equation with small time delay τ >0 . Furthermore, we prove the global existence and uniqueness of C^{α ,β } -solution to the time-delayed degenerate reaction-diffusion equation via compactness analysis. Finally, by the weighted energy method, we prove that the smooth non-critical traveling wave is globally stable in the weighted L^1 -space. The exponential convergence rate is also derived.

Existence and Stability of Traveling Waves for Degenerate Reaction-Diffusion Equation with Time Delay

NASA Astrophysics Data System (ADS)

Huang, Rui; Jin, Chunhua; Mei, Ming; Yin, Jingxue

2018-06-01

This paper deals with the existence and stability of traveling wave solutions for a degenerate reaction-diffusion equation with time delay. The degeneracy of spatial diffusion together with the effect of time delay causes us the essential difficulty for the existence of the traveling waves and their stabilities. In order to treat this case, we first show the existence of smooth- and sharp-type traveling wave solutions in the case of c≥c^* for the degenerate reaction-diffusion equation without delay, where c^*>0 is the critical wave speed of smooth traveling waves. Then, as a small perturbation, we obtain the existence of the smooth non-critical traveling waves for the degenerate diffusion equation with small time delay τ >0. Furthermore, we prove the global existence and uniqueness of C^{α ,β }-solution to the time-delayed degenerate reaction-diffusion equation via compactness analysis. Finally, by the weighted energy method, we prove that the smooth non-critical traveling wave is globally stable in the weighted L^1-space. The exponential convergence rate is also derived.

Diffusion phenomenon for linear dissipative wave equations in an exterior domain

NASA Astrophysics Data System (ADS)

Ikehata, Ryo

Under the general condition of the initial data, we will derive the crucial estimates which imply the diffusion phenomenon for the dissipative linear wave equations in an exterior domain. In order to derive the diffusion phenomenon for dissipative wave equations, the time integral method which was developed by Ikehata and Matsuyama (Sci. Math. Japon. 55 (2002) 33) plays an effective role.

Progress Towards a Rad-Hydro Code for Modern Computing Architectures LA-UR-10-02825

NASA Astrophysics Data System (ADS)

Wohlbier, J. G.; Lowrie, R. B.; Bergen, B.; Calef, M.

2010-11-01

We are entering an era of high performance computing where data movement is the overwhelming bottleneck to scalable performance, as opposed to the speed of floating-point operations per processor. All multi-core hardware paradigms, whether heterogeneous or homogeneous, be it the Cell processor, GPGPU, or multi-core x86, share this common trait. In multi-physics applications such as inertial confinement fusion or astrophysics, one may be solving multi-material hydrodynamics with tabular equation of state data lookups, radiation transport, nuclear reactions, and charged particle transport in a single time cycle. The algorithms are intensely data dependent, e.g., EOS, opacity, nuclear data, and multi-core hardware memory restrictions are forcing code developers to rethink code and algorithm design. For the past two years LANL has been funding a small effort referred to as Multi-Physics on Multi-Core to explore ideas for code design as pertaining to inertial confinement fusion and astrophysics applications. The near term goals of this project are to have a multi-material radiation hydrodynamics capability, with tabular equation of state lookups, on cartesian and curvilinear block structured meshes. In the longer term we plan to add fully implicit multi-group radiation diffusion and material heat conduction, and block structured AMR. We will report on our progress to date.

A Three-Wave Model of the Stratosphere with Coupled Dynamics, Radiation and Photochemistry. Appendix M

NASA Technical Reports Server (NTRS)

Shia, Run-Lie; Zhou, Shuntai; Ko, Malcolm K. W.; Sze, Nien-Dak; Salstein, David; Cady-Pereira, Karen

1997-01-01

A zonal mean chemistry transport model (2-D CTM) coupled with a semi-spectral dynamical model is used to simulate the distributions of trace gases in the present day atmosphere. The zonal-mean and eddy equations for the velocity and the geopotential height are solved in the semi-spectral dynamical model. The residual mean circulation is derived from these dynamical variables and used to advect the chemical species in the 2- D CTM. Based on a linearized wave transport equation, the eddy diffusion coefficients for chemical tracers are expressed in terms of the amplitude, frequency and growth rate of dynamical waves; local chemical loss rates; and a time constant parameterizing small scale mixing. The contributions to eddy flux are from the time varying wave amplitude (transient eddy), chemical reactions (chemical eddy) and small scale mixing. In spite of the high truncation in the dynamical module (only three longest waves are resolved), the model has simulated many observed characteristics of stratospheric dynamics and distribution of chemical species including ozone. Compared with the values commonly used in 2-D CTMs, the eddy diffusion coefficients for chemical species calculated in this model are smaller, especially in the subtropics. It is also found that the chemical eddy diffusion has only a small effects in determining the distribution of most slow species, including ozone in the stratosphere.

An enriched finite element method to fractional advection-diffusion equation

NASA Astrophysics Data System (ADS)

Luan, Shengzhi; Lian, Yanping; Ying, Yuping; Tang, Shaoqiang; Wagner, Gregory J.; Liu, Wing Kam

2017-08-01

In this paper, an enriched finite element method with fractional basis [ 1,x^{α }] for spatial fractional partial differential equations is proposed to obtain more stable and accurate numerical solutions. For pure fractional diffusion equation without advection, the enriched Galerkin finite element method formulation is demonstrated to simulate the exact solution successfully without any numerical oscillation, which is advantageous compared to the traditional Galerkin finite element method with integer basis [ 1,x] . For fractional advection-diffusion equation, the oscillatory behavior becomes complex due to the introduction of the advection term which can be characterized by a fractional element Peclet number. For the purpose of addressing the more complex numerical oscillation, an enriched Petrov-Galerkin finite element method is developed by using a dimensionless fractional stabilization parameter, which is formulated through a minimization of the residual of the nodal solution. The effectiveness and accuracy of the enriched finite element method are demonstrated by a series of numerical examples of fractional diffusion equation and fractional advection-diffusion equation, including both one-dimensional and two-dimensional, steady-state and time-dependent cases.

A Robust and Efficient Method for Steady State Patterns in Reaction-Diffusion Systems

PubMed Central

Lo, Wing-Cheong; Chen, Long; Wang, Ming; Nie, Qing

2012-01-01

An inhomogeneous steady state pattern of nonlinear reaction-diffusion equations with no-flux boundary conditions is usually computed by solving the corresponding time-dependent reaction-diffusion equations using temporal schemes. Nonlinear solvers (e.g., Newton’s method) take less CPU time in direct computation for the steady state; however, their convergence is sensitive to the initial guess, often leading to divergence or convergence to spatially homogeneous solution. Systematically numerical exploration of spatial patterns of reaction-diffusion equations under different parameter regimes requires that the numerical method be efficient and robust to initial condition or initial guess, with better likelihood of convergence to an inhomogeneous pattern. Here, a new approach that combines the advantages of temporal schemes in robustness and Newton’s method in fast convergence in solving steady states of reaction-diffusion equations is proposed. In particular, an adaptive implicit Euler with inexact solver (AIIE) method is found to be much more efficient than temporal schemes and more robust in convergence than typical nonlinear solvers (e.g., Newton’s method) in finding the inhomogeneous pattern. Application of this new approach to two reaction-diffusion equations in one, two, and three spatial dimensions, along with direct comparisons to several other existing methods, demonstrates that AIIE is a more desirable method for searching inhomogeneous spatial patterns of reaction-diffusion equations in a large parameter space. PMID:22773849

Parallel Solver for Diffuse Optical Tomography on Realistic Head Models With Scattering and Clear Regions.

PubMed

Placati, Silvio; Guermandi, Marco; Samore, Andrea; Scarselli, Eleonora Franchi; Guerrieri, Roberto

2016-09-01

Diffuse optical tomography is an imaging technique, based on evaluation of how light propagates within the human head to obtain the functional information about the brain. Precision in reconstructing such an optical properties map is highly affected by the accuracy of the light propagation model implemented, which needs to take into account the presence of clear and scattering tissues. We present a numerical solver based on the radiosity-diffusion model, integrating the anatomical information provided by a structural MRI. The solver is designed to run on parallel heterogeneous platforms based on multiple GPUs and CPUs. We demonstrate how the solver provides a 7 times speed-up over an isotropic-scattered parallel Monte Carlo engine based on a radiative transport equation for a domain composed of 2 million voxels, along with a significant improvement in accuracy. The speed-up greatly increases for larger domains, allowing us to compute the light distribution of a full human head ( ≈ 3 million voxels) in 116 s for the platform used.

Synchrotron radiation and diffusive shock acceleration - A short review and GRB perspective

NASA Astrophysics Data System (ADS)

Karlica, Mile

2015-12-01

In this talk we present the sponge" model and its possible implications on the GRB afterglow light curves. "Sponge" model describes source of GRB afterglow radiation as fragmented GRB ejecta where bubbles move through the rarefied medium. In the first part of the talk a short introduction to synchrotron radiation and Fermi acceleration was presented. In the assumption that X-ray luminosity of GRB afterglow phase comes from the kinetic energy losses of clouds in ejecta medium radiated as synchrotron radiation we solved currently very simple equation of motion to find which combination of cloud and medium regime describes the afterglow light curve the best. We proposed for the first step to watch simple combinations of expansion regimes for both bubbles and surrounding medium. The closest case to the numerical fit of GRB 150403A with time power law index k = 1.38 is the combination of constant bubbles and Sedov like expanding medium with time power law index k = 1.25. Of course the question of possible mixture of variuos regime combinations is still open within this model.

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

NASA Astrophysics Data System (ADS)

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

2017-10-01

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