Interfacial mixing in high-energy-density matter with a multiphysics kinetic model

NASA Astrophysics Data System (ADS)

Haack, Jeffrey R.; Hauck, Cory D.; Murillo, Michael S.

2017-12-01

We have extended a recently developed multispecies, multitemperature Bhatnagar-Gross-Krook model [Haack et al., J. Stat. Phys. 168, 822 (2017), 10.1007/s10955-017-1824-9], to include multiphysics capabilities that enable modeling of a wider range of physical conditions. In terms of geometry, we have extended from the spatially homogeneous setting to one spatial dimension. In terms of the physics, we have included an atomic ionization model, accurate collision physics across coupling regimes, self-consistent electric fields, and degeneracy in the electronic screening. We apply the model to a warm dense matter scenario in which the ablator-fuel interface of an inertial confinement fusion target is heated, but for larger length and time scales and for much higher temperatures than can be simulated using molecular dynamics. Relative to molecular dynamics, the kinetic model greatly extends the temperature regime and the spatiotemporal scales over which we are able to model. In our numerical results we observe hydrogen from the ablator material jetting into the fuel during the early stages of the implosion and compare the relative size of various diffusion components (Fickean diffusion, electrodiffusion, and barodiffusion) that drive this process. We also examine kinetic effects, such as anisotropic distributions and velocity separation, in order to determine when this problem can be described with a hydrodynamic model.

Parametrization of Drag and Turbulence for Urban Neighbourhoods with Trees

NASA Astrophysics Data System (ADS)

Krayenhoff, E. S.; Santiago, J.-L.; Martilli, A.; Christen, A.; Oke, T. R.

2015-08-01

Urban canopy parametrizations designed to be coupled with mesoscale models must predict the integrated effect of urban obstacles on the flow at each height in the canopy. To assess these neighbourhood-scale effects, results of microscale simulations may be horizontally-averaged. Obstacle-resolving computational fluid dynamics (CFD) simulations of neutrally-stratified flow through canopies of blocks (buildings) with varying distributions and densities of porous media (tree foliage) are conducted, and the spatially-averaged impacts on the flow of these building-tree combinations are assessed. The accuracy with which a one-dimensional (column) model with a one-equation (-) turbulence scheme represents spatially-averaged CFD results is evaluated. Individual physical mechanisms by which trees and buildings affect flow in the column model are evaluated in terms of relative importance. For the treed urban configurations considered, effects of buildings and trees may be considered independently. Building drag coefficients and length scale effects need not be altered due to the presence of tree foliage; therefore, parametrization of spatially-averaged flow through urban neighbourhoods with trees is greatly simplified. The new parametrization includes only source and sink terms significant for the prediction of spatially-averaged flow profiles: momentum drag due to buildings and trees (and the associated wake production of turbulent kinetic energy), modification of length scales by buildings, and enhanced dissipation of turbulent kinetic energy due to the small scale of tree foliage elements. Coefficients for the Santiago and Martilli (Boundary-Layer Meteorol 137: 417-439, 2010) parametrization of building drag coefficients and length scales are revised. Inclusion of foliage terms from the new parametrization in addition to the Santiago and Martilli building terms reduces root-mean-square difference (RMSD) of the column model streamwise velocity component and turbulent kinetic energy relative to the CFD model by 89 % in the canopy and 71 % above the canopy on average for the highest leaf area density scenarios tested: . RMSD values with the new parametrization are less than 20 % of mean layer magnitude for the streamwise velocity component within and above the canopy, and for above-canopy turbulent kinetic energy; RMSD values for within-canopy turbulent kinetic energy are negligible for most scenarios. The foliage-related portion of the new parametrization is required for scenarios with tree foliage of equal or greater height than the buildings, and for scenarios with foliage below roof height for building plan area densities less than approximately 0.25.

Frustration of resonant preheating by exotic kinetic terms

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

Rahmati, Shohreh; Seahra, Sanjeev S., E-mail: srahmati@unb.ca, E-mail: sseahra@unb.ca

2014-10-01

We study the effects of exotic kinetic terms on parametric resonance during the preheating epoch of the early universe. Specifically, we consider modifications to the action of ordinary matter fields motivated by generalized uncertainty principles, polymer quantization, as well as Dirac-Born-Infeld and k-essence models. To leading order in an ''exotic physics'' scale, the equations of motion derived from each of these models have the same algebraic form involving a nonlinear self-interaction in the matter sector. Neglecting spatial dependence, we show that the nonlinearity effectively shuts down the parametric resonance after a finite time period. We find numeric evidence that themore » frustration of parametric resonance persists to spatially inhomogenous matter in (1+1)-dimensions.« less

Processes of aggression described by kinetic method

NASA Astrophysics Data System (ADS)

Aristov, V. V.; Ilyin, O.

2014-12-01

In the last decades many investigations have been devoted to theoretical models in new areas concerning description of different biological, sociological and historical processes. In the present paper we suggest a model of the Nazi Germany invasion of Poland, France and USSR based on the kinetic theory. We model this process with the Cauchy boundary problem for the two-element kinetic equations with spatial initial conditions. The solution of the problem is given in the form of traveling wave. The propagation velocity of a frontline depends on the quotient between initial forces concentrations. Moreover it is obtained that the general solution of the model can be expressed in terms of quadratures and elementary functions. Finally it is shown that the frontline velocities are complied with the historical data.

Chemical Continuous Time Random Walks

NASA Astrophysics Data System (ADS)

Aquino, T.; Dentz, M.

2017-12-01

Traditional methods for modeling solute transport through heterogeneous media employ Eulerian schemes to solve for solute concentration. More recently, Lagrangian methods have removed the need for spatial discretization through the use of Monte Carlo implementations of Langevin equations for solute particle motions. While there have been recent advances in modeling chemically reactive transport with recourse to Lagrangian methods, these remain less developed than their Eulerian counterparts, and many open problems such as efficient convergence and reconstruction of the concentration field remain. We explore a different avenue and consider the question: In heterogeneous chemically reactive systems, is it possible to describe the evolution of macroscopic reactant concentrations without explicitly resolving the spatial transport? Traditional Kinetic Monte Carlo methods, such as the Gillespie algorithm, model chemical reactions as random walks in particle number space, without the introduction of spatial coordinates. The inter-reaction times are exponentially distributed under the assumption that the system is well mixed. In real systems, transport limitations lead to incomplete mixing and decreased reaction efficiency. We introduce an arbitrary inter-reaction time distribution, which may account for the impact of incomplete mixing. This process defines an inhomogeneous continuous time random walk in particle number space, from which we derive a generalized chemical Master equation and formulate a generalized Gillespie algorithm. We then determine the modified chemical rate laws for different inter-reaction time distributions. We trace Michaelis-Menten-type kinetics back to finite-mean delay times, and predict time-nonlocal macroscopic reaction kinetics as a consequence of broadly distributed delays. Non-Markovian kinetics exhibit weak ergodicity breaking and show key features of reactions under local non-equilibrium.

Self-consistent multidimensional electron kinetic model for inductively coupled plasma sources

NASA Astrophysics Data System (ADS)

Dai, Fa Foster

Inductively coupled plasma (ICP) sources have received increasing interest in microelectronics fabrication and lighting industry. In 2-D configuration space (r, z) and 2-D velocity domain (νθ,νz), a self- consistent electron kinetic analytic model is developed for various ICP sources. The electromagnetic (EM) model is established based on modal analysis, while the kinetic analysis gives the perturbed Maxwellian distribution of electrons by solving Boltzmann-Vlasov equation. The self- consistent algorithm combines the EM model and the kinetic analysis by updating their results consistently until the solution converges. The closed-form solutions in the analytical model provide rigorous and fast computing for the EM fields and the electron kinetic behavior. The kinetic analysis shows that the RF energy in an ICP source is extracted by a collisionless dissipation mechanism, if the electron thermovelocity is close to the RF phase velocities. A criterion for collisionless damping is thus given based on the analytic solutions. To achieve uniformly distributed plasma for plasma processing, we propose a novel discharge structure with both planar and vertical coil excitations. The theoretical results demonstrate improved uniformity for the excited azimuthal E-field in the chamber. Non-monotonic spatial decay in electric field and space current distributions was recently observed in weakly- collisional plasmas. The anomalous skin effect is found to be responsible for this phenomenon. The proposed model successfully models the non-monotonic spatial decay effect and achieves good agreements with the measurements for different applied RF powers. The proposed analytical model is compared with other theoretical models and different experimental measurements. The developed model is also applied to two kinds of ICP discharges used for electrodeless light sources. One structure uses a vertical internal coil antenna to excite plasmas and another has a metal shield to prevent the electromagnetic radiation. The theoretical results delivered by the proposed model agree quite well with the experimental measurements in many aspects. Therefore, the proposed self-consistent model provides an efficient and reliable means for designing ICP sources in various applications such as VLSI fabrication and electrodeless light sources.

WE-FG-206-06: Dual-Input Tracer Kinetic Modeling and Its Analog Implementation for Dynamic Contrast-Enhanced (DCE-) MRI of Malignant Mesothelioma (MPM)

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

Lee, S; Rimner, A; Hayes, S

Purpose: To use dual-input tracer kinetic modeling of the lung for mapping spatial heterogeneity of various kinetic parameters in malignant MPM Methods: Six MPM patients received DCE-MRI as part of their radiation therapy simulation scan. 5 patients had the epitheloid subtype of MPM, while one was biphasic. A 3D fast-field echo sequence with TR/TE/Flip angle of 3.62ms/1.69ms/15° was used for DCE-MRI acquisition. The scan was collected for 5 minutes with a temporal resolution of 5-9 seconds depending on the spatial extent of the tumor. A principal component analysis-based groupwise deformable registration was used to co-register all the DCE-MRI series formore » motion compensation. All the images were analyzed using five different dual-input tracer kinetic models implemented in analog continuous-time formalism: the Tofts-Kety (TK), extended TK (ETK), two compartment exchange (2CX), adiabatic approximation to the tissue homogeneity (AATH), and distributed parameter (DP) models. The following parameters were computed for each model: total blood flow (BF), pulmonary flow fraction (γ), pulmonary blood flow (BF-pa), systemic blood flow (BF-a), blood volume (BV), mean transit time (MTT), permeability-surface area product (PS), fractional interstitial volume (vi), extraction fraction (E), volume transfer constant (Ktrans) and efflux rate constant (kep). Results: Although the majority of patients had epitheloid histologies, kinetic parameter values varied across different models. One patient showed a higher total BF value in all models among the epitheloid histologies, although the γ value was varying among these different models. In one tumor with a large area of necrosis, the TK and ETK models showed higher E, Ktrans, and kep values and lower interstitial volume as compared to AATH and DP and 2CX models. Kinetic parameters such as BF-pa, BF-a, PS, Ktrans values were higher in surviving group compared to non-surviving group across most models. Conclusion: Dual-input tracer kinetic modeling is feasible in determining micro-vascular characteristics of MPM. This project was supported from Cycle for Survival and MSK Imaging and radiation science (IMRAS) grants.« less

Composition/Structure/Dynamics of comet and planetary satellite atmospheres

NASA Technical Reports Server (NTRS)

Combi, Michael R. (Principal Investigator)

1995-01-01

This research program addresses two cases of tenuous planetary atmospheres: comets and Io. The comet atmospheric research seeks to analyze a set of spatial profiles of CN in comet Halley taken in a 7.4-day period in April 1986; to apply a new dust coma model to various observations; and to analyze observations of the inner hydrogen coma, which can be optically thick to the resonance scattering of Lyman-alpha radiation, with the newly developed approach that combines a spherical radiative transfer model with our Monte Carlo H coma model. The Io research seeks to understand the atmospheric escape from Io with a hybrid-kinetic model for neutral gases and plasma given methods and algorithms developed for the study of neutral gas cometary atmospheres and the earth's polar wind and plasmasphere. Progress is reported on cometary Hydrogen Lyman-alpha studies; time-series analysis of cometary spatial profiles; model analysis of the dust comae of comets; and a global kinetic atmospheric model of Io.

MR-guided dynamic PET reconstruction with the kernel method and spectral temporal basis functions

NASA Astrophysics Data System (ADS)

Novosad, Philip; Reader, Andrew J.

2016-06-01

Recent advances in dynamic positron emission tomography (PET) reconstruction have demonstrated that it is possible to achieve markedly improved end-point kinetic parameter maps by incorporating a temporal model of the radiotracer directly into the reconstruction algorithm. In this work we have developed a highly constrained, fully dynamic PET reconstruction algorithm incorporating both spectral analysis temporal basis functions and spatial basis functions derived from the kernel method applied to a co-registered T1-weighted magnetic resonance (MR) image. The dynamic PET image is modelled as a linear combination of spatial and temporal basis functions, and a maximum likelihood estimate for the coefficients can be found using the expectation-maximization (EM) algorithm. Following reconstruction, kinetic fitting using any temporal model of interest can be applied. Based on a BrainWeb T1-weighted MR phantom, we performed a realistic dynamic [18F]FDG simulation study with two noise levels, and investigated the quantitative performance of the proposed reconstruction algorithm, comparing it with reconstructions incorporating either spectral analysis temporal basis functions alone or kernel spatial basis functions alone, as well as with conventional frame-independent reconstruction. Compared to the other reconstruction algorithms, the proposed algorithm achieved superior performance, offering a decrease in spatially averaged pixel-level root-mean-square-error on post-reconstruction kinetic parametric maps in the grey/white matter, as well as in the tumours when they were present on the co-registered MR image. When the tumours were not visible in the MR image, reconstruction with the proposed algorithm performed similarly to reconstruction with spectral temporal basis functions and was superior to both conventional frame-independent reconstruction and frame-independent reconstruction with kernel spatial basis functions. Furthermore, we demonstrate that a joint spectral/kernel model can also be used for effective post-reconstruction denoising, through the use of an EM-like image-space algorithm. Finally, we applied the proposed algorithm to reconstruction of real high-resolution dynamic [11C]SCH23390 data, showing promising results.

MR-guided dynamic PET reconstruction with the kernel method and spectral temporal basis functions.

PubMed

Novosad, Philip; Reader, Andrew J

2016-06-21

Recent advances in dynamic positron emission tomography (PET) reconstruction have demonstrated that it is possible to achieve markedly improved end-point kinetic parameter maps by incorporating a temporal model of the radiotracer directly into the reconstruction algorithm. In this work we have developed a highly constrained, fully dynamic PET reconstruction algorithm incorporating both spectral analysis temporal basis functions and spatial basis functions derived from the kernel method applied to a co-registered T1-weighted magnetic resonance (MR) image. The dynamic PET image is modelled as a linear combination of spatial and temporal basis functions, and a maximum likelihood estimate for the coefficients can be found using the expectation-maximization (EM) algorithm. Following reconstruction, kinetic fitting using any temporal model of interest can be applied. Based on a BrainWeb T1-weighted MR phantom, we performed a realistic dynamic [(18)F]FDG simulation study with two noise levels, and investigated the quantitative performance of the proposed reconstruction algorithm, comparing it with reconstructions incorporating either spectral analysis temporal basis functions alone or kernel spatial basis functions alone, as well as with conventional frame-independent reconstruction. Compared to the other reconstruction algorithms, the proposed algorithm achieved superior performance, offering a decrease in spatially averaged pixel-level root-mean-square-error on post-reconstruction kinetic parametric maps in the grey/white matter, as well as in the tumours when they were present on the co-registered MR image. When the tumours were not visible in the MR image, reconstruction with the proposed algorithm performed similarly to reconstruction with spectral temporal basis functions and was superior to both conventional frame-independent reconstruction and frame-independent reconstruction with kernel spatial basis functions. Furthermore, we demonstrate that a joint spectral/kernel model can also be used for effective post-reconstruction denoising, through the use of an EM-like image-space algorithm. Finally, we applied the proposed algorithm to reconstruction of real high-resolution dynamic [(11)C]SCH23390 data, showing promising results.

Stress enhanced calcium kinetics in a neuron.

PubMed

Kant, Aayush; Bhandakkar, Tanmay K; Medhekar, Nikhil V

2018-02-01

Accurate modeling of the mechanobiological response of a Traumatic Brain Injury is beneficial toward its effective clinical examination, treatment and prevention. Here, we present a stress history-dependent non-spatial kinetic model to predict the microscale phenomena of secondary insults due to accumulation of excess calcium ions (Ca[Formula: see text]) induced by the macroscale primary injuries. The model is able to capture the experimentally observed increase and subsequent partial recovery of intracellular Ca[Formula: see text] concentration in response to various types of mechanical impulses. We further establish the accuracy of the model by comparing our predictions with key experimental observations.

Local reaction kinetics by imaging☆

PubMed Central

Suchorski, Yuri; Rupprechter, Günther

2016-01-01

In the present contribution we present an overview of our recent studies using the “kinetics by imaging” approach for CO oxidation on heterogeneous model systems. The method is based on the correlation of the PEEM image intensity with catalytic activity: scaled down to the μm-sized surface regions, such correlation allows simultaneous local kinetic measurements on differently oriented individual domains of a polycrystalline metal-foil, including the construction of local kinetic phase diagrams. This allows spatially- and component-resolved kinetic studies and, e.g., a direct comparison of inherent catalytic properties of Pt(hkl)- and Pd(hkl)-domains or supported μm-sized Pd-powder agglomerates, studies of the local catalytic ignition and the role of defects and grain boundaries in the local reaction kinetics. PMID:26865736

Processes of aggression described by kinetic method

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

Aristov, V. V.; Ilyin, O.

In the last decades many investigations have been devoted to theoretical models in new areas concerning description of different biological, sociological and historical processes. In the present paper we suggest a model of the Nazi Germany invasion of Poland, France and USSR based on the kinetic theory. We model this process with the Cauchy boundary problem for the two-element kinetic equations with spatial initial conditions. The solution of the problem is given in the form of traveling wave. The propagation velocity of a frontline depends on the quotient between initial forces concentrations. Moreover it is obtained that the general solutionmore » of the model can be expressed in terms of quadratures and elementary functions. Finally it is shown that the frontline velocities are complied with the historical data.« less

Studies of ion kinetic effects in OMEGA shock-driven implosions using fusion burn imaging

NASA Astrophysics Data System (ADS)

Rosenberg, M. J.; Seguin, F. H.; Rinderknecht, H. G.; Sio, H.; Zylstra, A. B.; Gatu Johnson, M.; Frenje, J. A.; Li, C. K.; Petrasso, R. D.; Amendt, P. A.; Wilks, S. C.; Zimmerman, G.; Hoffman, N. M.; Kagan, G.; Molvig, K.; Glebov, V. Yu.; Stoeckl, C.; Marshall, F. J.; Seka, W.; Delettrez, J. A.; Sangster, T. C.; Betti, R.; Meyerhofer, D. D.; Atzeni, S.; Nikroo, A.

2014-10-01

Ion kinetic effects have been inferred in a series of shock-driven implosions at OMEGA from an increasing yield discrepancy between observations and hydrodynamic simulations as the ion-ion mean free path increases. To more precisely identify the nature and impact of ion kinetic effects, spatial burn profile measurements of DD and D3He reactions in these D3He-filled shock-driven implosions are presented and contrasted to both purely hydrodynamic models and models that include ion kinetic effects. It is shown that in implosions where the ion mean free path is equal to or greater than the size of the fuel region, purely hydrodynamic models fail to capture the observed burn profiles, while a model that includes ion diffusion is able to recover the observed burn profile shape. These results further elucidate the ion kinetic mechanisms that are present under long mean-free-path conditions after shock convergence in both shock-driven and ablatively-driven implosions. This work was supported in part by the U.S. DOE, NLUF, LLE, and LLNL.

Hydrolysis kinetics in anaerobic degradation of particulate organic material: an overview.

PubMed

Vavilin, V A; Fernandez, B; Palatsi, J; Flotats, X

2008-01-01

The applicability of different kinetics to the hydrolysis of particulate organic material in anaerobic digestion is discussed. Hydrolysis has traditionally been modelled according to the first-order kinetics. For complex substrate, the first-order kinetics should be modified in order to take into account hardly degradable material. It has been shown that models in which hydrolysis is coupled to the growth of hydrolytic bacteria work well at high or at fluctuant organic loading. In particular, the surface-related two-phase and the Contois models showed good fits to experimental data from a wide range of organic waste. Both models tend to the first-order kinetics at a high biomass-to-waste ratio and, for this reason, they can be considered as more general models. Examples on different inhibition processes that might affect the degradation of solid waste are reported. Acetogenesis or methanogenesis might be the rate-limiting stages in complex waste. In such cases, stimulation of hydrolysis (mechanically, chemically or biologically) may lead to a further inhibition of these stages, which ultimately affects hydrolysis as well. Since the hydrolysis process is characterized by surface and transport phenomena, new developments in spatially distributed models are considered fundamental to provide new insights in this complex process.

Kinetic Description of Ionospheric Outflows Based on the Exact Form of Fokker-Planck Collision Operator: Electrons

NASA Technical Reports Server (NTRS)

Khazanov, George V.; Khabibrakhmanov, Ildar K.; Glocer, Alex

2012-01-01

We present the results of a finite difference implementation of the kinetic Fokker-Planck model with an exact form of the nonlinear collisional operator, The model is time dependent and three-dimensional; one spatial dimension and two in velocity space. The spatial dimension is aligned with the local magnetic field, and the velocity space is defined by the magnitude of the velocity and the cosine of pitch angle. An important new feature of model, the concept of integration along the particle trajectories, is discussed in detail. Integration along the trajectories combined with the operator time splitting technique results in a solution scheme which accurately accounts for both the fast convection of the particles along the magnetic field lines and relatively slow collisional process. We present several tests of the model's performance and also discuss simulation results of the evolution of the plasma distribution for realistic conditions in Earth's plasmasphere under different scenarios.

A mathematical model for foreign body reactions in 2D.

PubMed

Su, Jianzhong; Gonzales, Humberto Perez; Todorov, Michail; Kojouharov, Hristo; Tang, Liping

2011-02-01

The foreign body reactions are commonly referred to the network of immune and inflammatory reactions of human or animals to foreign objects placed in tissues. They are basic biological processes, and are also highly relevant to bioengineering applications in implants, as fibrotic tissue formations surrounding medical implants have been found to substantially reduce the effectiveness of devices. Despite of intensive research on determining the mechanisms governing such complex responses, few mechanistic mathematical models have been developed to study such foreign body reactions. This study focuses on a kinetics-based predictive tool in order to analyze outcomes of multiple interactive complex reactions of various cells/proteins and biochemical processes and to understand transient behavior during the entire period (up to several months). A computational model in two spatial dimensions is constructed to investigate the time dynamics as well as spatial variation of foreign body reaction kinetics. The simulation results have been consistent with experimental data and the model can facilitate quantitative insights for study of foreign body reaction process in general.

Computational modeling of fully-ionized, magnetized plasmas using the fluid approximation

NASA Astrophysics Data System (ADS)

Schnack, Dalton

2005-10-01

Strongly magnetized plasmas are rich in spatial and temporal scales, making a computational approach useful for studying these systems. The most accurate model of a magnetized plasma is based on a kinetic equation that describes the evolution of the distribution function for each species in six-dimensional phase space. However, the high dimensionality renders this approach impractical for computations for long time scales in relevant geometry. Fluid models, derived by taking velocity moments of the kinetic equation [1] and truncating (closing) the hierarchy at some level, are an approximation to the kinetic model. The reduced dimensionality allows a wider range of spatial and/or temporal scales to be explored. Several approximations have been used [2-5]. Successful computational modeling requires understanding the ordering and closure approximations, the fundamental waves supported by the equations, and the numerical properties of the discretization scheme. We review and discuss several ordering schemes, their normal modes, and several algorithms that can be applied to obtain a numerical solution. The implementation of kinetic parallel closures is also discussed [6].[1] S. Chapman and T.G. Cowling, ``The Mathematical Theory of Non-Uniform Gases'', Cambridge University Press, Cambridge, UK (1939).[2] R.D. Hazeltine and J.D. Meiss, ``Plasma Confinement'', Addison-Wesley Publishing Company, Redwood City, CA (1992).[3] L.E. Sugiyama and W. Park, Physics of Plasmas 7, 4644 (2000).[4] J.J. Ramos, Physics of Plasmas, 10, 3601 (2003).[5] P.J. Catto and A.N. Simakov, Physics of Plasmas, 11, 90 (2004).[6] E.D. Held et al., Phys. Plasmas 11, 2419 (2004)

Power-law spatial dispersion from fractional Liouville equation

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

Tarasov, Vasily E.

2013-10-15

A microscopic model in the framework of fractional kinetics to describe spatial dispersion of power-law type is suggested. The Liouville equation with the Caputo fractional derivatives is used to obtain the power-law dependence of the absolute permittivity on the wave vector. The fractional differential equations for electrostatic potential in the media with power-law spatial dispersion are derived. The particular solutions of these equations for the electric potential of point charge in this media are considered.

Information hidden in the velocity distribution of ions and the exact kinetic Bohm criterion

NASA Astrophysics Data System (ADS)

Tsankov, Tsanko V.; Czarnetzki, Uwe

2017-05-01

Non-equilibrium distribution functions of electrons and ions play an important role in plasma physics. A prominent example is the kinetic Bohm criterion. Since its first introduction it has been controversial for theoretical reasons and due to the lack of experimental data, in particular on the ion distribution function. Here we resolve the theoretical as well as the experimental difficulties by an exact solution of the kinetic Boltzmann equation including charge exchange collisions and ionization. This also allows for the first time non-invasive measurement of spatially resolved ion velocity distributions, absolute values of the ion and electron densities, temperatures, and mean energies as well as the electric field and the plasma potential in the entire plasma. The non-invasive access to the spatially resolved distribution functions of electrons and ions is applied to the problem of the kinetic Bohm criterion. Theoretically a so far missing term in the criterion is derived and shown to be of key importance. With the new term the validity of the kinetic criterion at high collisionality and its agreement with the fluid picture are restored. All findings are supported by experimental data, theory and a numerical model with excellent agreement throughout.

Localized Models of Charged Particle Motion in Martian Crustal Magnetic Cusps

NASA Astrophysics Data System (ADS)

Brain, D. A.; Poppe, A. R.; Jarvinen, R.; Dong, Y.; Egan, H. L.; Fang, X.

2017-12-01

The induced magnetosphere of Mars is punctuated by localized but strong crustal magnetic fields that are observed to play host to a variety of phenomena typically associated with global magnetic fields, such as auroral processes and particle precipitation, field-aligned current systems, and ion outflow. Each of these phenomena occur on the night side, in small-scale magnetic `cusp' regions of vertically aligned field. Cusp regions are not yet capable of being spatially resolved in global scale models that include the ion kinetics necessary for simulating charged particle transport along cusps. Local models are therefore necessary if we are to understand how cusp processes operate at Mars. Here we present the first results of an effort to model the kinetic particle motion and electric fields in Martian cusps. We are adapting both a 1.5D Particle-in-Cell (PIC) model for lunar magnetic cusps regions to the Martian case and a hybrid model framework (used previously for the global Martian plasma interaction and for lunar magnetic anomaly regions) to cusps in 2D. By comparing the models we can asses the importance of electron kinetics in particle transport along cusp field lines. In this first stage of our study we model a moderately strong nightside cusp, with incident hot hydrogen plasma from above, and cold planetary (oxygen) plasma entering the simulation from below. We report on the spatial and temporal distribution of plasma along cusp field lines for this initial case.

Interfacial mixing in high energy-density matter with a multiphysics kinetic model

NASA Astrophysics Data System (ADS)

Haack, Jeff; Hauck, Cory; Murillo, Michael

2017-10-01

We have extended a recently-developed multispecies, multitemperature BGK model to include multiphysics capability that allows modeling of a wider range of plasma conditions. In particular, we have extended the model to describe one spatial dimension, and included a multispecies atomic ionization model, accurate collision physics across coupling regimes, self-consistent electric fields, and degeneracy in the electronic screening. We apply the new model to a warm dense matter scenario in which the ablator-fuel interface of an inertial confinement fusion target is heated, similar to a recent molecular dynamics study, but for larger length and time scales and for much higher temperatures. From our numerical results we are able to explore a variety of phenomena, including hydrogen jetting, kinetic effects (non-Maxwellian and anisotropic distributions), plasma physics (size, persistence and role of electric fields) and transport (relative sizes of Fickean diffision, electrodiffusion and barodiffusion). As compared with the recent molecular dynamics work the kinetic model greatly extends the accessible physical domains we are able to model.

Ion kinetic effects on the ignition and burn of inertial confinement fusion targets: A multi-scale approach

NASA Astrophysics Data System (ADS)

Peigney, B. E.; Larroche, O.; Tikhonchuk, V.

2014-12-01

In this article, we study the hydrodynamics and burn of the thermonuclear fuel in inertial confinement fusion pellets at the ion kinetic level. The analysis is based on a two-velocity-scale Vlasov-Fokker-Planck kinetic model that is specially tailored to treat fusion products (suprathermal α-particles) in a self-consistent manner with the thermal bulk. The model assumes spherical symmetry in configuration space and axial symmetry in velocity space around the mean flow velocity. A typical hot-spot ignition design is considered. Compared with fluid simulations where a multi-group diffusion scheme is applied to model α transport, the full ion-kinetic approach reveals significant non-local effects on the transport of energetic α-particles. This has a direct impact on hydrodynamic spatial profiles during combustion: the hot spot reactivity is reduced, while the inner dense fuel layers are pre-heated by the escaping α-suprathermal particles, which are transported farther out of the hot spot. We show how the kinetic transport enhancement of fusion products leads to a significant reduction of the fusion yield.

Ion kinetic effects on the ignition and burn of inertial confinement fusion targets: A multi-scale approach

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

Peigney, B. E.; Larroche, O.; Tikhonchuk, V.

2014-12-15

In this article, we study the hydrodynamics and burn of the thermonuclear fuel in inertial confinement fusion pellets at the ion kinetic level. The analysis is based on a two-velocity-scale Vlasov-Fokker-Planck kinetic model that is specially tailored to treat fusion products (suprathermal α-particles) in a self-consistent manner with the thermal bulk. The model assumes spherical symmetry in configuration space and axial symmetry in velocity space around the mean flow velocity. A typical hot-spot ignition design is considered. Compared with fluid simulations where a multi-group diffusion scheme is applied to model α transport, the full ion-kinetic approach reveals significant non-local effectsmore » on the transport of energetic α-particles. This has a direct impact on hydrodynamic spatial profiles during combustion: the hot spot reactivity is reduced, while the inner dense fuel layers are pre-heated by the escaping α-suprathermal particles, which are transported farther out of the hot spot. We show how the kinetic transport enhancement of fusion products leads to a significant reduction of the fusion yield.« less

A Kinetics Model for KrF Laser Amplifiers

NASA Astrophysics Data System (ADS)

Giuliani, J. L.; Kepple, P.; Lehmberg, R.; Obenschain, S. P.; Petrov, G.

1999-11-01

A computer kinetics code has been developed to model the temporal and spatial behavior of an e-beam pumped KrF laser amplifier. The deposition of the primary beam electrons is assumed to be spatially uniform and the energy distribution function of the nascent electron population is calculated to be near Maxwellian below 10 eV. For an initial Kr/Ar/F2 composition, the code calculates the densities of 24 species subject to over 100 reactions with 1-D spatial resolution (typically 16 zones) along the longitudinal lasing axis. Enthalpy accounting for each process is performed to partition the energy into internal, thermal, and radiative components. The electron as well as the heavy particle temperatures are followed for energy conservation and excitation rates. Transport of the lasing photons is performed along the axis on a dense subgrid using the method of characteristics. Amplified spontaneous emission is calculated using a discrete ordinates approach and includes contributions to the local intensity from the whole amplifier volume. Specular reflection off side walls and the rear mirror are included. Results of the model will be compared with data from the NRL NIKE laser and other published results.

Stochastic population dynamics in spatially extended predator-prey systems

NASA Astrophysics Data System (ADS)

Dobramysl, Ulrich; Mobilia, Mauro; Pleimling, Michel; Täuber, Uwe C.

2018-02-01

Spatially extended population dynamics models that incorporate demographic noise serve as case studies for the crucial role of fluctuations and correlations in biological systems. Numerical and analytic tools from non-equilibrium statistical physics capture the stochastic kinetics of these complex interacting many-particle systems beyond rate equation approximations. Including spatial structure and stochastic noise in models for predator-prey competition invalidates the neutral Lotka-Volterra population cycles. Stochastic models yield long-lived erratic oscillations stemming from a resonant amplification mechanism. Spatially extended predator-prey systems display noise-stabilized activity fronts that generate persistent correlations. Fluctuation-induced renormalizations of the oscillation parameters can be analyzed perturbatively via a Doi-Peliti field theory mapping of the master equation; related tools allow detailed characterization of extinction pathways. The critical steady-state and non-equilibrium relaxation dynamics at the predator extinction threshold are governed by the directed percolation universality class. Spatial predation rate variability results in more localized clusters, enhancing both competing species’ population densities. Affixing variable interaction rates to individual particles and allowing for trait inheritance subject to mutations induces fast evolutionary dynamics for the rate distributions. Stochastic spatial variants of three-species competition with ‘rock-paper-scissors’ interactions metaphorically describe cyclic dominance. These models illustrate intimate connections between population dynamics and evolutionary game theory, underscore the role of fluctuations to drive populations toward extinction, and demonstrate how space can support species diversity. Two-dimensional cyclic three-species May-Leonard models are characterized by the emergence of spiraling patterns whose properties are elucidated by a mapping onto a complex Ginzburg-Landau equation. Multiple-species extensions to general ‘food networks’ can be classified on the mean-field level, providing both fundamental understanding of ensuing cooperativity and profound insight into the rich spatio-temporal features and coarsening kinetics in the corresponding spatially extended systems. Novel space-time patterns emerge as a result of the formation of competing alliances; e.g. coarsening domains that each incorporate rock-paper-scissors competition games.

Temporal acceleration of spatially distributed kinetic Monte Carlo simulations

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

Chatterjee, Abhijit; Vlachos, Dionisios G.

The computational intensity of kinetic Monte Carlo (KMC) simulation is a major impediment in simulating large length and time scales. In recent work, an approximate method for KMC simulation of spatially uniform systems, termed the binomial {tau}-leap method, was introduced [A. Chatterjee, D.G. Vlachos, M.A. Katsoulakis, Binomial distribution based {tau}-leap accelerated stochastic simulation, J. Chem. Phys. 122 (2005) 024112], where molecular bundles instead of individual processes are executed over coarse-grained time increments. This temporal coarse-graining can lead to significant computational savings but its generalization to spatially lattice KMC simulation has not been realized yet. Here we extend the binomial {tau}-leapmore » method to lattice KMC simulations by combining it with spatially adaptive coarse-graining. Absolute stability and computational speed-up analyses for spatial systems along with simulations provide insights into the conditions where accuracy and substantial acceleration of the new spatio-temporal coarse-graining method are ensured. Model systems demonstrate that the r-time increment criterion of Chatterjee et al. obeys the absolute stability limit for values of r up to near 1.« less

Geometry Genetics and Evolution

NASA Astrophysics Data System (ADS)

Siggia, Eric

2011-03-01

Darwin argued that highly perfected organs such as the vertebrate eye could evolve by a series of small changes, each of which conferred a selective advantage. In the context of gene networks, this idea can be recast into a predictive algorithm, namely find networks that can be built by incremental adaptation (gradient search) to perform some task. It embodies a ``kinetic'' view of evolution where a solution that is quick to evolve is preferred over a global optimum. Examples of biochemical kinetic networks were evolved for temporal adaptation, temperature compensated entrainable clocks, explore-exploit trade off in signal discrimination, will be presented as well as networks that model the spatially periodic somites (vertebrae) and HOX gene expression in the vertebrate embryo. These models appear complex by the criterion of 19th century applied mathematics since there is no separation of time or spatial scales, yet they are all derivable by gradient optimization of simple functions (several in the Pareto evolution) often based on the Shannon entropy of the time or spatial response. Joint work with P. Francois, Physics Dept. McGill University. With P. Francois, Physics Dept. McGill University

Spatially organizing biochemistry: choosing a strategy to translate synthetic biology to the factory.

PubMed

Jakobson, Christopher M; Tullman-Ercek, Danielle; Mangan, Niall M

2018-05-29

Natural biochemical systems are ubiquitously organized both in space and time. Engineering the spatial organization of biochemistry has emerged as a key theme of synthetic biology, with numerous technologies promising improved biosynthetic pathway performance. One strategy, however, may produce disparate results for different biosynthetic pathways. We use a spatially resolved kinetic model to explore this fundamental design choice in systems and synthetic biology. We predict that two example biosynthetic pathways have distinct optimal organization strategies that vary based on pathway-dependent and cell-extrinsic factors. Moreover, we demonstrate that the optimal design varies as a function of kinetic and biophysical properties, as well as culture conditions. Our results suggest that organizing biosynthesis has the potential to substantially improve performance, but that choosing the appropriate strategy is key. The flexible design-space analysis we propose can be adapted to diverse biosynthetic pathways, and lays a foundation to rationally choose organization strategies for biosynthesis.

Phase-field Model for Interstitial Loop Growth Kinetics and Thermodynamic and Kinetic Models of Irradiated Fe-Cr Alloys

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

Li, Yulan; Hu, Shenyang Y.; Sun, Xin

2011-06-15

Microstructure evolution kinetics in irradiated materials has strongly spatial correlation. For example, void and second phases prefer to nucleate and grow at pre-existing defects such as dislocations, grain boundaries, and cracks. Inhomogeneous microstructure evolution results in inhomogeneity of microstructure and thermo-mechanical properties. Therefore, the simulation capability for predicting three dimensional (3-D) microstructure evolution kinetics and its subsequent impact on material properties and performance is crucial for scientific design of advanced nuclear materials and optimal operation conditions in order to reduce uncertainty in operational and safety margins. Very recently the meso-scale phase-field (PF) method has been used to predict gas bubblemore » evolution, void swelling, void lattice formation and void migration in irradiated materials,. Although most results of phase-field simulations are qualitative due to the lake of accurate thermodynamic and kinetic properties of defects, possible missing of important kinetic properties and processes, and the capability of current codes and computers for large time and length scale modeling, the simulations demonstrate that PF method is a promising simulation tool for predicting 3-D heterogeneous microstructure and property evolution, and providing microstructure evolution kinetics for higher scale level simulations of microstructure and property evolution such as mean field methods. This report consists of two parts. In part I, we will present a new phase-field model for predicting interstitial loop growth kinetics in irradiated materials. The effect of defect (vacancy/interstitial) generation, diffusion and recombination, sink strength, long-range elastic interaction, inhomogeneous and anisotropic mobility on microstructure evolution kinetics is taken into account in the model. The model is used to study the effect of elastic interaction on interstitial loop growth kinetics, the interstitial flux, and sink strength of interstitial loop for interstitials. In part II, we present a generic phase field model and discuss the thermodynamic and kinetic properties in phase-field models including the reaction kinetics of radiation defects and local free energy of irradiated materials. In particular, a two-sublattice thermodynamic model is suggested to describe the local free energy of alloys with irradiated defects. Fe-Cr alloy is taken as an example to explain the required thermodynamic and kinetic properties for quantitative phase-field modeling. Finally the great challenges in phase-field modeling will be discussed.« less

The influence of spatially and temporally high-resolution wind forcing on the power input to near-inertial waves in the ocean

NASA Astrophysics Data System (ADS)

Rimac, A.; Eden, C.; von Storch, J.

2012-12-01

Coexistence of stable stratification, the meridional overturning circulation and meso-scale eddies and their influence on the ocean's circulation still raise complex questions concerning the ocean energetics. Oceanic general circulation is mainly forced by the wind field and deep water tides. Its essential energetics are the conversion of kinetic energy of the winds and tides into oceanic potential and kinetic energy. Energy needed for the circulation is bound to internal wave fields. Direct internal wave generation by the wind at the sea surface is one of the sources of this energy. Previous studies using mixed-layer type of models and low frequency wind forcings (six-hourly and daily) left room for improvement. Using mixed-layer models it is not possible to assess the distribution of near-inertial energy into the deep ocean. Also, coarse temporal resolution of wind forcing strongly underestimates the near-inertial wave energy. To overcome this difficulty we use a high resolution ocean model with high frequency wind forcings. We establish the following model setup: We use the Max Planck Institute Ocean Model (MPIOM) on a tripolar grid with 45km horizontal resolution and 40 vertical levels. We run the model with wind forcings that vary in horizontal (250km versus 40km) and temporal resolution (six versus one-hourly). In our study we answer the following questions: How big is the wind kinetic energy input to the near-inertial waves? Is the kinetic energy of the near-inertial waves enhanced when high-frequency wind forcings are used? If so, by how much and why, due to higher level of temporal wind variability or due to better spatial representation of the near-inertial waves? How big is the total power of near-inertial waves generated by the wind at the surface of the ocean? We run the model for one year. Our model results show that the near-inertial waves are excited both using wind forcings of high and low horizontal and temporal resolution. Near-inertial energy is almost two times higher when we force the model with high frequency wind forcings. The influence on the energy mostly depends on the time difference between two forcing fields while the spatial difference has little influence.

Continuum-Kinetic Models and Numerical Methods for Multiphase Applications

NASA Astrophysics Data System (ADS)

Nault, Isaac Michael

This thesis presents a continuum-kinetic approach for modeling general problems in multiphase solid mechanics. In this context, a continuum model refers to any model, typically on the macro-scale, in which continuous state variables are used to capture the most important physics: conservation of mass, momentum, and energy. A kinetic model refers to any model, typically on the meso-scale, which captures the statistical motion and evolution of microscopic entitites. Multiphase phenomena usually involve non-negligible micro or meso-scopic effects at the interfaces between phases. The approach developed in the thesis attempts to combine the computational performance benefits of a continuum model with the physical accuracy of a kinetic model when applied to a multiphase problem. The approach is applied to modeling a single particle impact in Cold Spray, an engineering process that intimately involves the interaction of crystal grains with high-magnitude elastic waves. Such a situation could be classified a multiphase application due to the discrete nature of grains on the spatial scale of the problem. For this application, a hyper elasto-plastic model is solved by a finite volume method with approximate Riemann solver. The results of this model are compared for two types of plastic closure: a phenomenological macro-scale constitutive law, and a physics-based meso-scale Crystal Plasticity model.

A spatial focusing model for G protein signals. Regulator of G protein signaling (RGS) protien-mediated kinetic scaffolding.

PubMed

Zhong, Huailing; Wade, Susan M; Woolf, Peter J; Linderman, Jennifer J; Traynor, John R; Neubig, Richard R

2003-02-28

Regulators of G protein signaling (RGS) are GTPase-accelerating proteins (GAPs), which can inhibit heterotrimeric G protein pathways. In this study, we provide experimental and theoretical evidence that high concentrations of receptors (as at a synapse) can lead to saturation of GDP-GTP exchange making GTP hydrolysis rate-limiting. This results in local depletion of inactive heterotrimeric G-GDP, which is reversed by RGS GAP activity. Thus, RGS enhances receptor-mediated G protein activation even as it deactivates the G protein. Evidence supporting this model includes a GTP-dependent enhancement of guanosine 5'-3-O-(thio)triphosphate (GTPgammaS) binding to G(i) by RGS. The RGS domain of RGS4 is sufficient for this, not requiring the NH(2)- or COOH-terminal extensions. Furthermore, a kinetic model including only the GAP activity of RGS replicates the GTP-dependent enhancement of GTPgammaS binding observed experimentally. Finally in a Monte Carlo model, this mechanism results in a dramatic "spatial focusing" of active G protein. Near the receptor, G protein activity is maintained even with RGS due to the ability of RGS to reduce depletion of local Galpha-GDP levels permitting rapid recoupling to receptor and maintained G protein activation near the receptor. In contrast, distant signals are suppressed by the RGS, since Galpha-GDP is not depleted there. Thus, a novel RGS-mediated "kinetic scaffolding" mechanism is proposed which narrows the spatial range of active G protein around a cluster of receptors limiting the spill-over of G protein signals to more distant effector molecules, thus enhancing the specificity of G(i) protein signals.

Mound-Interface Kinetics in Dictyostelium Aggregation

NASA Astrophysics Data System (ADS)

Tutu, Hiroki

2002-09-01

The mound development of the cellular slime mold amoebae Dictyostelium discoideum is studied with an interface kinetic model for the height of cell layers. As a competitive role for the chemotaxis, we compare two types of curvature relaxations; the surface relaxation induced by cell-substrate affinity (model A), and that comes from a cell-cell adhesive effect (model B). It is found that both models are characterized by the growth law for the maximum mound height. Based on a self-similarity scaling hypothesis for the spatial structure of streaming pattern, we suggest a scaling law for the growth of mound-height hmax ˜ t1-1/α+β/α with α = 2 (4) for the model A (B) and a number 0 ≤ β < 1.

Scaling Effects of Cr(VI) Reduction Kinetics. The Role of Geochemical Heterogeneity

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

Wang, Li; Li, Li

2015-10-22

The natural subsurface is highly heterogeneous with minerals distributed in different spatial patterns. Fundamental understanding of how mineral spatial distribution patterns regulate sorption process is important for predicting the transport and fate of chemicals. Existing studies about the sorption was carried out in well-mixed batch reactors or uniformly packed columns, with few data available on the effects of spatial heterogeneities. As a result, there is a lack of data and understanding on how spatial heterogeneities control sorption processes. In this project, we aim to understand and develop modeling capabilities to predict the sorption of Cr(VI), an omnipresent contaminant in naturalmore » systems due to its natural occurrence and industrial utilization. We systematically examine the role of spatial patterns of illite, a common clay, in determining the extent of transport limitation and scaling effects associated with Cr(VI) sorption capacity and kinetics using column experiments and reactive transport modeling. Our results showed that the sorbed mass and rates can differ by an order of magnitude due to of the illite spatial heterogeneities and transport limitation. With constraints from data, we also developed the capabilities of modeling Cr(VI) in heterogeneous media. The developed model is then utilized to understand the general principles that govern the relationship between sorption and connectivity, a key measure of the spatial pattern characteristics. This correlation can be used to estimate Cr(VI) sorption characteristics in heterogeneous porous media. Insights gained here bridge gaps between laboratory and field application in hydrogeology and geochemical field, and advance predictive understanding of reactive transport processes in the natural heterogeneous subsurface. We believe that these findings will be of interest to a large number of environmental geochemists and engineers, hydrogeologists, and those interested in contaminant fate and transport, water quality and water composition, and natural attenuation processes in natural systems.« less

Applications of nonlinear science and kinetic equations to the spread of epidemics

NASA Astrophysics Data System (ADS)

Macinnis, David Robert

The study of the spread of epidemics is currently growing into a successful subfield of a combination of nonlinear science and statistical mechanics. Topics studied in this field include kinetic and mean field levels of epidemiological models. This thesis consists of the analysis of such topics and specifically directed at the Hantavirus, West Nile virus, and the Bubonic Plague. A successful reaction-diffusion equation approach developed recently by Abramson and Kenkre was able to describe spatiotemporal patterns of the Hantavirus model. From measurements of the parameters of their model it was found that the mice, the carriers of the infection, must be regarded as moving diffusively within attractive potentials representative of home ranges. Several attempts have been made to incorporate home ranges into their model. Two of these attempts are discussed within this thesis. A model to explain the transmission of the West Nile virus within bird and mosquito populations was recently developed by Kenkre, Parmenter, Peixoto, and Sadasiv who showed how spatially resolved issues could be discussed but restricted their analysis to mean field considerations. This thesis extends that study by investigating spatial resolution of the infected populations. Traveling waves of the bird and mosquito populations are found in the West Nile context. Infection control of various epidemics has become increasingly important to limit the potential force of infection into the human population. This thesis contains a quantitative attempt at a theory of such control (for the West Nile virus) via spraying of the mosquito population. Mean field and kinetic level models are proposed in this thesis to describe the transmission of the Bubonic Plague which involves flea and mammal populations. The various populations are found to undergo a variety of bifurcations as well as hysteresis in their steady state regime. Spatially resolved analysis of the populations is also presented.

Spatial Inhomogeneity of Kinetic and Magnetic Dissipations in Thermal Convection

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

Hotta, H.

We investigate the inhomogeneity of kinetic and magnetic dissipations in thermal convection using high-resolution calculations. In statistically steady turbulence, the injected and dissipated energies are balanced. This means that a large amount of energy is continuously converted into internal energy via dissipation. As in thermal convection, downflows are colder than upflows and the inhomogeneity of the dissipation potentially changes the convection structure. Our investigation of the inhomogeneity of the dissipation shows the following. (1) More dissipation is seen around the bottom of the calculation domain, and this tendency is promoted with the magnetic field. (2) The dissipation in the downflowmore » is much larger than that in the upflow. The dissipation in the downflow is more than 80% of the total at maximum. This tendency is also promoted with the magnetic field. (3) Although 2D probability density functions of the kinetic and magnetic dissipations versus the vertical velocity are similar, the kinetic and magnetic dissipations are not well correlated. Our result suggests that the spatial inhomogeneity of the dissipation is significant and should be considered when modeling a small-scale strong magnetic field generated with an efficient small-scale dynamo for low-resolution calculations.« less

TU-C-12A-11: Comparisons Between Cu-ATSM PET and DCE-CT Kinetic Parameters in Canine Sinonasal Tumors

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

La Fontaine, M; Bradshaw, T; Kubicek, L

2014-06-15

Purpose: Regions of poor perfusion within tumors may be associated with higher hypoxic levels. This study aimed to test this hypothesis by comparing measurements of hypoxia from Cu-ATSM PET to vasculature kinetic parameters from DCE-CT kinetic analysis. Methods: Ten canine patients with sinonasal tumors received one Cu-ATSM PET/CT scan and three DCE-CT scans prior to treatment. Cu-ATSM PET/CT and DCE-CT scans were registered and resampled to matching voxel dimensions. Kinetic analysis was performed on DCE-CT scans and for each patient, the resulting kinetic parameter values from the three DCE-CT scans were averaged together. Cu-ATSM SUVs were spatially correlated (r{sub spatial})more » on a voxel-to-voxel basis against the following DCE-CT kinetic parameters: transit time (t{sub 1}), blood flow (F), vasculature fraction (v{sub 1}), and permeability (PS). In addition, whole-tumor comparisons were performed by correlating (r{sub ROI}) the mean Cu-ATSM SUV (SUV{sub mean}) with median kinetic parameter values. Results: The spatial correlations (r{sub spatial}) were poor and ranged from -0.04 to 0.21 for all kinetic parameters. These low spatial correlations may be due to high variability in the DCE-CT kinetic parameter voxel values between scans. In our hypothesis, t{sub 1} was expected to have a positive correlation, while F was expected to have a negative correlation to hypoxia. However, in wholetumor analysis the opposite was found for both t{sub 1} (r{sub ROI} = -0.25) and F (r{sub ROI} = 0.56). PS and v{sub 1} may depict angiogenic responses to hypoxia and found positive correlations to Cu-ATSM SUV for PS (r{sub ROI} = 0.41), and v{sub 1} (r{sub ROI} = 0.57). Conclusion: Low spatial correlations were found between Cu-ATSM uptake and DCE-CT vasculature parameters, implying that poor perfusion is not associated with higher hypoxic regions. Across patients, the most hypoxic tumors tended to have higher blood flow values, which is contrary to our initial hypothesis. Funding: R01 CA136927.« less

Photocatalytic mineralization of commercial herbicides in a pilot-scale solar CPC reactor: photoreactor modeling and reaction kinetics constants independent of radiation field.

PubMed

Colina-Márquez, Jose; Machuca-Martínez, Fiderman; Li Puma, Gianluca

2009-12-01

The six-flux absorption-scattering model (SFM) of the radiation field in the photoreactor, combined with reaction kinetics and fluid-dynamic models, has proved to be suitable to describe the degradation of water pollutants in heterogeneous photocatalytic reactors, combining simplicity and accuracy. In this study, the above approach was extended to model the photocatalytic mineralization of a commercial herbicides mixture (2,4-D, diuron, and ametryne used in Colombian sugar cane crops) in a solar, pilot-scale, compound parabolic collector (CPC) photoreactor using a slurry suspension of TiO(2). The ray-tracing technique was used jointly with the SFM to determine the direction of both the direct and diffuse solar photon fluxes and the spatial profile of the local volumetric rate of photon absorption (LVRPA) in the CPC reactor. Herbicides mineralization kinetics with explicit photon absorption effects were utilized to remove the dependence of the observed rate constants from the reactor geometry and radiation field in the photoreactor. The results showed that the overall model fitted the experimental data of herbicides mineralization in the solar CPC reactor satisfactorily for both cloudy and sunny days. Using the above approach kinetic parameters independent of the radiation field in the reactor can be estimated directly from the results of experiments carried out in a solar CPC reactor. The SFM combined with reaction kinetics and fluid-dynamic models proved to be a simple, but reliable model, for solar photocatalytic applications.

Object Segmentation from Motion Discontinuities and Temporal Occlusions–A Biologically Inspired Model

PubMed Central

Beck, Cornelia; Ognibeni, Thilo; Neumann, Heiko

2008-01-01

Background Optic flow is an important cue for object detection. Humans are able to perceive objects in a scene using only kinetic boundaries, and can perform the task even when other shape cues are not provided. These kinetic boundaries are characterized by the presence of motion discontinuities in a local neighbourhood. In addition, temporal occlusions appear along the boundaries as the object in front covers the background and the objects that are spatially behind it. Methodology/Principal Findings From a technical point of view, the detection of motion boundaries for segmentation based on optic flow is a difficult task. This is due to the problem that flow detected along such boundaries is generally not reliable. We propose a model derived from mechanisms found in visual areas V1, MT, and MSTl of human and primate cortex that achieves robust detection along motion boundaries. It includes two separate mechanisms for both the detection of motion discontinuities and of occlusion regions based on how neurons respond to spatial and temporal contrast, respectively. The mechanisms are embedded in a biologically inspired architecture that integrates information of different model components of the visual processing due to feedback connections. In particular, mutual interactions between the detection of motion discontinuities and temporal occlusions allow a considerable improvement of the kinetic boundary detection. Conclusions/Significance A new model is proposed that uses optic flow cues to detect motion discontinuities and object occlusion. We suggest that by combining these results for motion discontinuities and object occlusion, object segmentation within the model can be improved. This idea could also be applied in other models for object segmentation. In addition, we discuss how this model is related to neurophysiological findings. The model was successfully tested both with artificial and real sequences including self and object motion. PMID:19043613

Modeling and simulation of CANDU reactor and its regulating system

NASA Astrophysics Data System (ADS)

Javidnia, Hooman

Analytical computer codes are indispensable tools in design, optimization, and control of nuclear power plants. Numerous codes have been developed to perform different types of analyses related to the nuclear power plants. A large number of these codes are designed to perform safety analyses. In the context of safety analyses, the control system is often neglected. Although there are good reasons for such a decision, that does not mean that the study of control systems in the nuclear power plants should be neglected altogether. In this thesis, a proof of concept code is developed as a tool that can be used in the design. optimization. and operation stages of the control system. The main objective in the design of this computer code is providing a tool that is easy to use by its target audience and is capable of producing high fidelity results that can be trusted to design the control system and optimize its performance. Since the overall plant control system covers a very wide range of processes, in this thesis the focus has been on one particular module of the the overall plant control system, namely, the reactor regulating system. The center of the reactor regulating system is the CANDU reactor. A nodal model for the reactor is used to represent the spatial neutronic kinetics of the core. The nodal model produces better results compared to the point kinetics model which is often used in the design and analysis of control system for nuclear reactors. The model can capture the spatial effects to some extent. although it is not as detailed as the finite difference methods. The criteria for choosing a nodal model of the core are: (1) the model should provide more detail than point kinetics and capture spatial effects, (2) it should not be too complex or overly detailed to slow down the simulation and provide details that are extraneous or unnecessary for a control engineer. Other than the reactor itself, there are auxiliary models that describe dynamics of different phenomena related to the transfer of the energy from the core. The main function of the reactor regulating system is to control the power of the reactor. This is achieved by using a set of detectors. reactivity devices. and digital control algorithms. Three main reactivity devices that are activated during short-term or intermediate-term transients are modeled in this thesis. The main elements of the digital control system are implemented in accordance to the program specifications for the actual control system in CANDU reactors. The simulation results are validated against requirements of the reactor regulating system. actual plant data. and pre-validated data from other computer codes. The validation process shows that the simulation results can be trusted in making engineering decisions regarding the reactor regulating system and prediction of the system performance in response to upset conditions or disturbances. KEYWORDS: CANDU reactors. reactor regulating system. nodal model. spatial kinetics. reactivity devices. simulation.

Spatial and temporal features of the growth of a bacterial species colonizing the zebrafish gut.

PubMed

Jemielita, Matthew; Taormina, Michael J; Burns, Adam R; Hampton, Jennifer S; Rolig, Annah S; Guillemin, Karen; Parthasarathy, Raghuveer

2014-12-16

The vertebrate intestine is home to microbial ecosystems that play key roles in host development and health. Little is known about the spatial and temporal dynamics of these microbial communities, limiting our understanding of fundamental properties, such as their mechanisms of growth, propagation, and persistence. To address this, we inoculated initially germ-free zebrafish larvae with fluorescently labeled strains of an Aeromonas species, representing an abundant genus in the zebrafish gut. Using light sheet fluorescence microscopy to obtain three-dimensional images spanning the gut, we quantified the entire bacterial load, as founding populations grew from tens to tens of thousands of cells over several hours. The data yield the first ever measurements of the growth kinetics of a microbial species inside a live vertebrate intestine and show dynamics that robustly fit a logistic growth model. Intriguingly, bacteria were nonuniformly distributed throughout the gut, and bacterial aggregates showed considerably higher growth rates than did discrete individuals. The form of aggregate growth indicates intrinsically higher division rates for clustered bacteria, rather than surface-mediated agglomeration onto clusters. Thus, the spatial organization of gut bacteria both relative to the host and to each other impacts overall growth kinetics, suggesting that spatial characterizations will be an important input to predictive models of host-associated microbial community assembly. Our intestines are home to vast numbers of microbes that influence many aspects of health and disease. Though we now know a great deal about the constituents of the gut microbiota, we understand very little about their spatial structure and temporal dynamics in humans or in any animal: how microbial populations establish themselves, grow, fluctuate, and persist. To address this, we made use of a model organism, the zebrafish, and a new optical imaging technique, light sheet fluorescence microscopy, to visualize for the first time the colonization of a live, vertebrate gut by specific bacteria with sufficient resolution to quantify the population over a range from a few individuals to tens of thousands of bacterial cells. Our results provide unprecedented measures of bacterial growth kinetics and also show the influence of spatial structure on bacterial populations, which can be revealed only by direct imaging. Copyright © 2014 Jemielita et al.

Transient Simulation of the Multi-SERTTA Experiment with MAMMOTH

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

Ortensi, Javier; Baker, Benjamin; Wang, Yaqi

This work details the MAMMOTH reactor physics simulations of the Static Environment Rodlet Transient Test Apparatus (SERTTA) conducted at Idaho National Laboratory in FY-2017. TREAT static-environment experiment vehicles are being developed to enable transient testing of Pressurized Water Reactor (PWR) type fuel specimens, including fuel concepts with enhanced accident tolerance (Accident Tolerant Fuels, ATF). The MAMMOTH simulations include point reactor kinetics as well as spatial dynamics for a temperature-limited transient. The strongly coupled multi-physics solutions of the neutron flux and temperature fields are second order accurate both in the spatial and temporal domains. MAMMOTH produces pellet stack powers that are within 1.5% of the Monte Carlo reference solutions. Some discrepancies between the MCNP model used in the design of the flux collars and the Serpent/MAMMOTH models lead to higher power and energy deposition values in Multi-SERTTA unit 1. The TREAT core results compare well with the safety case computed with point reactor kinetics in RELAP5-3D. The reactor period is 44 msec, which corresponds to a reactivity insertion of 2.685% delta k/kmore » $. The peak core power in the spatial dynamics simulation is 431 MW, which the point kinetics model over-predicts by 12%. The pulse width at half the maximum power is 0.177 sec. Subtle transient effects are apparent at the beginning insertion in the experimental samples due to the control rod removal. Additional difference due to transient effects are observed in the sample powers and enthalpy. The time dependence of the power coupling factor (PCF) is calculated for the various fuel stacks of the Multi-SERTTA vehicle. Sample temperatures in excess of 3100 K, the melting point UO$$_2$$, are computed with the adiabatic heat transfer model. The planned shaped-transient might introduce additional effects that cannot be predicted with PRK models. Future modeling will be focused on the shaped-transient by improving the control rod models in MAMMOTH and adding the BISON thermo-elastic models and thermal-fluids heat transfer.« less

Traveltime-based descriptions of transport and mixing in heterogeneous domains

NASA Astrophysics Data System (ADS)

Luo, Jian; Cirpka, Olaf A.

2008-09-01

Modeling mixing-controlled reactive transport using traditional spatial discretization of the domain requires identifying the spatial distributions of hydraulic and reactive parameters including mixing-related quantities such as dispersivities and kinetic mass transfer coefficients. In most applications, breakthrough curves (BTCs) of conservative and reactive compounds are measured at only a few locations and spatially explicit models are calibrated by matching these BTCs. A common difficulty in such applications is that the individual BTCs differ too strongly to justify the assumption of spatial homogeneity, whereas the number of observation points is too small to identify the spatial distribution of the decisive parameters. The key objective of the current study is to characterize physical transport by the analysis of conservative tracer BTCs and predict the macroscopic BTCs of compounds that react upon mixing from the interpretation of conservative tracer BTCs and reactive parameters determined in the laboratory. We do this in the framework of traveltime-based transport models which do not require spatially explicit, costly aquifer characterization. By considering BTCs of a conservative tracer measured on different scales, one can distinguish between mixing, which is a prerequisite for reactions, and spreading, which per se does not foster reactions. In the traveltime-based framework, the BTC of a solute crossing an observation plane, or ending in a well, is interpreted as the weighted average of concentrations in an ensemble of non-interacting streamtubes, each of which is characterized by a distinct traveltime value. Mixing is described by longitudinal dispersion and/or kinetic mass transfer along individual streamtubes, whereas spreading is characterized by the distribution of traveltimes, which also determines the weights associated with each stream tube. Key issues in using the traveltime-based framework include the description of mixing mechanisms and the estimation of the traveltime distribution. In this work, we account for both apparent longitudinal dispersion and kinetic mass transfer as mixing mechanisms, thus generalizing the stochastic-convective model with or without inter-phase mass transfer and the advective-dispersive streamtube model. We present a nonparametric approach of determining the traveltime distribution, given a BTC integrated over an observation plane and estimated mixing parameters. The latter approach is superior to fitting parametric models in cases wherein the true traveltime distribution exhibits multiple peaks or long tails. It is demonstrated that there is freedom for the combinations of mixing parameters and traveltime distributions to fit conservative BTCs and describe the tailing. A reactive transport case of a dual Michaelis-Menten problem demonstrates that the reactive mixing introduced by local dispersion and mass transfer may be described by apparent mean mass transfer with coefficients evaluated by local BTCs.

Eddy-driven low-frequency variability: physics and observability through altimetry

NASA Astrophysics Data System (ADS)

Penduff, Thierry; Sérazin, Guillaume; Arbic, Brian; Mueller, Malte; Richman, James G.; Shriver, Jay F.; Morten, Andrew J.; Scott, Robert B.

2015-04-01

Model studies have revealed the propensity of the eddying ocean circulation to generate strong low-frequency variability (LFV) intrinsically, i.e. without low-frequency atmospheric variability. In the present study, gridded satellite altimeter products, idealized quasi-geostrophic (QG) turbulent simulations, and realistic high-resolution global ocean simulations are used to study the spontaneous tendency of mesoscale (relatively high frequency and high wavenumber) kinetic energy to non-linearly cascade towards larger time and space scales. The QG model reveals that large-scale variability, arising from the well-known spatial inverse cascade, is associated with low frequencies. Low-frequency, low-wavenumber energy is maintained primarily by nonlinearities in the QG model, with forcing (by large-scale shear) and friction playing secondary roles. In realistic simulations, nonlinearities also generally drive kinetic energy to low frequencies and low wavenumbers. In some, but not all, regions of the gridded altimeter product, surface kinetic energy is also found to cascade toward low frequencies. Exercises conducted with the realistic model suggest that the spatial and temporal filtering inherent in the construction of gridded satellite altimeter maps may contribute to the discrepancies seen in some regions between the direction of frequency cascade in models versus gridded altimeter maps. Finally, the range of frequencies that are highly energized and engaged these cascades appears much greater than the range of highly energized and engaged wavenumbers. Global eddying simulations, performed in the context of the CHAOCEAN project in collaboration with the CAREER project, provide estimates of the range of timescales that these oceanic nonlinearities are likely to feed without external variability.

Class of self-limiting growth models in the presence of nonlinear diffusion

NASA Astrophysics Data System (ADS)

Kar, Sandip; Banik, Suman Kumar; Ray, Deb Shankar

2002-06-01

The source term in a reaction-diffusion system, in general, does not involve explicit time dependence. A class of self-limiting growth models dealing with animal and tumor growth and bacterial population in a culture, on the other hand, are described by kinetics with explicit functions of time. We analyze a reaction-diffusion system to study the propagation of spatial front for these models.

Fuzzy Performance between Surface Fitting and Energy Distribution in Turbulence Runner

PubMed Central

Liang, Zhongwei; Liu, Xiaochu; Ye, Bangyan; Brauwer, Richard Kars

2012-01-01

Because the application of surface fitting algorithms exerts a considerable fuzzy influence on the mathematical features of kinetic energy distribution, their relation mechanism in different external conditional parameters must be quantitatively analyzed. Through determining the kinetic energy value of each selected representative position coordinate point by calculating kinetic energy parameters, several typical algorithms of complicated surface fitting are applied for constructing microkinetic energy distribution surface models in the objective turbulence runner with those obtained kinetic energy values. On the base of calculating the newly proposed mathematical features, we construct fuzzy evaluation data sequence and present a new three-dimensional fuzzy quantitative evaluation method; then the value change tendencies of kinetic energy distribution surface features can be clearly quantified, and the fuzzy performance mechanism discipline between the performance results of surface fitting algorithms, the spatial features of turbulence kinetic energy distribution surface, and their respective environmental parameter conditions can be quantitatively analyzed in detail, which results in the acquirement of final conclusions concerning the inherent turbulence kinetic energy distribution performance mechanism and its mathematical relation. A further turbulence energy quantitative study can be ensured. PMID:23213287

Kinetic Models for Topological Nearest-Neighbor Interactions

NASA Astrophysics Data System (ADS)

Blanchet, Adrien; Degond, Pierre

2017-12-01

We consider systems of agents interacting through topological interactions. These have been shown to play an important part in animal and human behavior. Precisely, the system consists of a finite number of particles characterized by their positions and velocities. At random times a randomly chosen particle, the follower, adopts the velocity of its closest neighbor, the leader. We study the limit of a system size going to infinity and, under the assumption of propagation of chaos, show that the limit kinetic equation is a non-standard spatial diffusion equation for the particle distribution function. We also study the case wherein the particles interact with their K closest neighbors and show that the corresponding kinetic equation is the same. Finally, we prove that these models can be seen as a singular limit of the smooth rank-based model previously studied in Blanchet and Degond (J Stat Phys 163:41-60, 2016). The proofs are based on a combinatorial interpretation of the rank as well as some concentration of measure arguments.

Application and theoretical analysis of the flamelet model for supersonic turbulent combustion flows in the scramjet engine

NASA Astrophysics Data System (ADS)

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

2014-11-01

In the framework of Reynolds-averaged Navier-Stokes simulation, supersonic turbulent combustion flows at the German Aerospace Centre (DLR) combustor and Japan Aerospace Exploration Agency (JAXA) integrated scramjet engine are numerically simulated using the flamelet model. Based on the DLR combustor case, theoretical analysis and numerical experiments conclude that: the finite rate model only implicitly considers the large-scale turbulent effect and, due to the lack of the small-scale non-equilibrium effect, it would overshoot the peak temperature compared to the flamelet model in general. Furthermore, high-Mach-number compressibility affects the flamelet model mainly through two ways: the spatial pressure variation and the static enthalpy variation due to the kinetic energy. In the flamelet library, the mass fractions of the intermediate species, e.g. OH, are more sensible to the above two effects than the main species such as H2O. Additionally, in the combustion flowfield where the pressure is larger than the value adopted in the generation of the flamelet library or the conversion from the static enthalpy to the kinetic energy occurs, the temperature obtained by the flamelet model without taking compressibility effects into account would be undershot, and vice versa. The static enthalpy variation effect has only little influence on the temperature simulation of the flamelet model, while the effect of the spatial pressure variation may cause relatively large errors. From the JAXA case, it is found that the flamelet model cannot in general be used for an integrated scramjet engine. The existence of the inlet together with the transverse injection scheme could cause large spatial variations of pressure, so the pressure value adopted for the generation of a flamelet library should be fine-tuned according to a pre-simulation of pure mixing.

Pulsed IR Heating Studies of Single-Molecule DNA Duplex Dissociation Kinetics and Thermodynamics

PubMed Central

Holmstrom, Erik D.; Dupuis, Nicholas F.; Nesbitt, David J.

2014-01-01

Single-molecule fluorescence spectroscopy is a powerful technique that makes it possible to observe the conformational dynamics associated with biomolecular processes. The addition of precise temperature control to these experiments can yield valuable thermodynamic information about equilibrium and kinetic rate constants. To accomplish this, we have developed a microscopy technique based on infrared laser overtone/combination band absorption to heat small (≈10−11 liter) volumes of water. Detailed experimental characterization of this technique reveals three major advantages over conventional stage heating methods: 1), a larger range of steady-state temperatures (20–100°C); 2), substantially superior spatial (≤20 μm) control; and 3), substantially superior temporal (≈1 ms) control. The flexibility and breadth of this spatial and temporally resolved laser-heating approach is demonstrated in single-molecule fluorescence assays designed to probe the dissociation of a 21 bp DNA duplex. These studies are used to support a kinetic model based on nucleic acid end fraying that describes dissociation for both short (<10 bp) and long (>10 bp) DNA duplexes. These measurements have been extended to explore temperature-dependent kinetics for the 21 bp construct, which permit determination of single-molecule activation enthalpies and entropies for DNA duplex dissociation. PMID:24411254

Construction and comparison of parallel implicit kinetic solvers in three spatial dimensions

NASA Astrophysics Data System (ADS)

Titarev, Vladimir; Dumbser, Michael; Utyuzhnikov, Sergey

2014-01-01

The paper is devoted to the further development and systematic performance evaluation of a recent deterministic framework Nesvetay-3D for modelling three-dimensional rarefied gas flows. Firstly, a review of the existing discretization and parallelization strategies for solving numerically the Boltzmann kinetic equation with various model collision integrals is carried out. Secondly, a new parallelization strategy for the implicit time evolution method is implemented which improves scaling on large CPU clusters. Accuracy and scalability of the methods are demonstrated on a pressure-driven rarefied gas flow through a finite-length circular pipe as well as an external supersonic flow over a three-dimensional re-entry geometry of complicated aerodynamic shape.

Nonequilibrium transition and pattern formation in a linear reaction-diffusion system with self-regulated kinetics

NASA Astrophysics Data System (ADS)

Paul, Shibashis; Ghosh, Shyamolina; Ray, Deb Shankar

2018-02-01

We consider a reaction-diffusion system with linear, stochastic activator-inhibitor kinetics where the time evolution of concentration of a species at any spatial location depends on the relative average concentration of its neighbors. This self-regulating nature of kinetics brings in spatial correlation between the activator and the inhibitor. An interplay of this correlation in kinetics and disparity of diffusivities of the two species leads to symmetry breaking non-equilibrium transition resulting in stationary pattern formation. The role of initial noise strength and the linear reaction terms has been analyzed for pattern selection.

Predictive computation of genomic logic processing functions in embryonic development

PubMed Central

Peter, Isabelle S.; Faure, Emmanuel; Davidson, Eric H.

2012-01-01

Gene regulatory networks (GRNs) control the dynamic spatial patterns of regulatory gene expression in development. Thus, in principle, GRN models may provide system-level, causal explanations of developmental process. To test this assertion, we have transformed a relatively well-established GRN model into a predictive, dynamic Boolean computational model. This Boolean model computes spatial and temporal gene expression according to the regulatory logic and gene interactions specified in a GRN model for embryonic development in the sea urchin. Additional information input into the model included the progressive embryonic geometry and gene expression kinetics. The resulting model predicted gene expression patterns for a large number of individual regulatory genes each hour up to gastrulation (30 h) in four different spatial domains of the embryo. Direct comparison with experimental observations showed that the model predictively computed these patterns with remarkable spatial and temporal accuracy. In addition, we used this model to carry out in silico perturbations of regulatory functions and of embryonic spatial organization. The model computationally reproduced the altered developmental functions observed experimentally. Two major conclusions are that the starting GRN model contains sufficiently complete regulatory information to permit explanation of a complex developmental process of gene expression solely in terms of genomic regulatory code, and that the Boolean model provides a tool with which to test in silico regulatory circuitry and developmental perturbations. PMID:22927416

Nonequilibrium relaxations within the ground-state manifold in the antiferromagnetic Ising model on a triangular lattice.

PubMed

Kim, Eunhye; Lee, Sung Jong; Kim, Bongsoo

2007-02-01

We present an extensive Monte Carlo simulation study on the nonequilibrium kinetics of triangular antiferromagnetic Ising model within the ground state ensemble which consists of sectors, each of which is characterized by a unique value of the string density p through a dimer covering method. Building upon our recent work [Phys. Rev. E 68, 066127 (2003)] where we considered the nonequilibrium relaxation observed within the dominant sector with p=2/3, we here focus on the nonequilibrium kinetics within the minor sectors with p<2/3. The initial configurations are chosen as those in which the strings are straight and evenly distributed. In the minor sectors, we observe a characteristic spatial anisotropy in both equilibrium and nonequilibrium spatial correlations. We observe emergence of a critical relaxation region (in the spatial and temporal domain) which grows as p deviates from p=2/3. Spatial anisotropy appears in the equilibrium spatial correlation with the characteristic length scale xi(e,V)(p) diverging with vanishing string density as xi(e,V)(p) approximately p(-2) along the vertical direction, while along the horizontal direction the spatial length scale diverges as xi(e,H) approximately p(-1). Analytic forms for the anisotropic equilibrium correlation functions are given. We also find that the spin autocorrelation function A(t) shows a simple scaling behavior A(t)=A(t/tau(A)(p)), where the time scale tau(A)(p) shows a power-law divergence with vanishing p as tau(A)(p) approximately p(-phi) with phi approximately or equal to 4. These features can be understood in terms of random walk nature of the fluctuations of the strings within the typical separation between neighboring strings.

Assessment of mean-field microkinetic models for CO methanation on stepped metal surfaces using accelerated kinetic Monte Carlo

NASA Astrophysics Data System (ADS)

Andersen, Mie; Plaisance, Craig P.; Reuter, Karsten

2017-10-01

First-principles screening studies aimed at predicting the catalytic activity of transition metal (TM) catalysts have traditionally been based on mean-field (MF) microkinetic models, which neglect the effect of spatial correlations in the adsorbate layer. Here we critically assess the accuracy of such models for the specific case of CO methanation over stepped metals by comparing to spatially resolved kinetic Monte Carlo (kMC) simulations. We find that the typical low diffusion barriers offered by metal surfaces can be significantly increased at step sites, which results in persisting correlations in the adsorbate layer. As a consequence, MF models may overestimate the catalytic activity of TM catalysts by several orders of magnitude. The potential higher accuracy of kMC models comes at a higher computational cost, which can be especially challenging for surface reactions on metals due to a large disparity in the time scales of different processes. In order to overcome this issue, we implement and test a recently developed algorithm for achieving temporal acceleration of kMC simulations. While the algorithm overall performs quite well, we identify some challenging cases which may lead to a breakdown of acceleration algorithms and discuss possible directions for future algorithm development.

Limits of applicability of a two-temperature model under nonuniform heating of metal by an ultrashort laser pulse

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

Polyakov, D S; Yakovlev, E B

The heating of metals (silver and aluminium) by ultrashort laser pulses is analysed proceeding from a spatially nonuniform kinetic equation for the electron distribution function. The electron subsystem thermalisation is estimated in a wide range of absorbed pulse energy density. The limits of applicability are determined for the two-temperature model. (interaction of laser radiation with matter)

THE SPECTRAL AMPLITUDE OF STELLAR CONVECTION AND ITS SCALING IN THE HIGH-RAYLEIGH-NUMBER REGIME

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

Featherstone, Nicholas A.; Hindman, Bradley W., E-mail: feathern@colorado.edu

2016-02-10

Convection plays a central role in the dynamics of any stellar interior, and yet its operation remains largely hidden from direct observation. As a result, much of our understanding concerning stellar convection necessarily derives from theoretical and computational models. The Sun is, however, exceptional in that regard. The wealth of observational data afforded by its proximity provides a unique test bed for comparing convection models against observations. When such comparisons are carried out, surprising inconsistencies between those models and observations become apparent. Both photospheric and helioseismic measurements suggest that convection simulations may overestimate convective flow speeds on large spatial scales.more » Moreover, many solar convection simulations have difficulty reproducing the observed solar differential rotation owing to this apparent overestimation. We present a series of three-dimensional stellar convection simulations designed to examine how the amplitude and spectral distribution of convective flows are established within a star’s interior. While these simulations are nonmagnetic and nonrotating in nature, they demonstrate two robust phenomena. When run with sufficiently high Rayleigh number, the integrated kinetic energy of the convection becomes effectively independent of thermal diffusion, but the spectral distribution of that kinetic energy remains sensitive to both of these quantities. A simulation that has converged to a diffusion-independent value of kinetic energy will divide that energy between spatial scales such that low-wavenumber power is overestimated and high-wavenumber power is underestimated relative to a comparable system possessing higher Rayleigh number. We discuss the implications of these results in light of the current inconsistencies between models and observations.« less

Spatial Rule-Based Modeling: A Method and Its Application to the Human Mitotic Kinetochore

PubMed Central

Ibrahim, Bashar; Henze, Richard; Gruenert, Gerd; Egbert, Matthew; Huwald, Jan; Dittrich, Peter

2013-01-01

A common problem in the analysis of biological systems is the combinatorial explosion that emerges from the complexity of multi-protein assemblies. Conventional formalisms, like differential equations, Boolean networks and Bayesian networks, are unsuitable for dealing with the combinatorial explosion, because they are designed for a restricted state space with fixed dimensionality. To overcome this problem, the rule-based modeling language, BioNetGen, and the spatial extension, SRSim, have been developed. Here, we describe how to apply rule-based modeling to integrate experimental data from different sources into a single spatial simulation model and how to analyze the output of that model. The starting point for this approach can be a combination of molecular interaction data, reaction network data, proximities, binding and diffusion kinetics and molecular geometries at different levels of detail. We describe the technique and then use it to construct a model of the human mitotic inner and outer kinetochore, including the spindle assembly checkpoint signaling pathway. This allows us to demonstrate the utility of the procedure, show how a novel perspective for understanding such complex systems becomes accessible and elaborate on challenges that arise in the formulation, simulation and analysis of spatial rule-based models. PMID:24709796

Modeling Photo-Bleaching Kinetics to Create High Resolution Maps of Rod Rhodopsin in the Human Retina

PubMed Central

Ehler, Martin; Dobrosotskaya, Julia; Cunningham, Denise; Wong, Wai T.; Chew, Emily Y.; Czaja, Wojtek; Bonner, Robert F.

2015-01-01

We introduce and describe a novel non-invasive in-vivo method for mapping local rod rhodopsin distribution in the human retina over a 30-degree field. Our approach is based on analyzing the brightening of detected lipofuscin autofluorescence within small pixel clusters in registered imaging sequences taken with a commercial 488nm confocal scanning laser ophthalmoscope (cSLO) over a 1 minute period. We modeled the kinetics of rhodopsin bleaching by applying variational optimization techniques from applied mathematics. The physical model and the numerical analysis with its implementation are outlined in detail. This new technique enables the creation of spatial maps of the retinal rhodopsin and retinal pigment epithelium (RPE) bisretinoid distribution with an ≈ 50μm resolution. PMID:26196397

Assessment of ion kinetic effects in shock-driven inertial confinement fusion (ICF) implosions using fusion burn imaging

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

Rosenberg, M. J.; Séguin, F. H.; Amendt, P. A.

The significance and nature of ion kinetic effects in D³He-filled, shock-driven inertial confinement fusion implosions are assessed through measurements of fusion burn profiles. Over this series of experiments, the ratio of ion-ion mean free path to minimum shell radius (the Knudsen number, N K) was varied from 0.3 to 9 in order to probe hydrodynamic-like to strongly kinetic plasma conditions; as the Knudsen number increased, hydrodynamic models increasingly failed to match measured yields, while an empirically-tuned, first-step model of ion kinetic effects better captured the observed yield trends [Rosenberg et al., Phys. Rev. Lett. 112, 185001 (2014)]. Here, spatially resolvedmore » measurements of the fusion burn are used to examine kinetic ion transport effects in greater detail, adding an additional dimension of understanding that goes beyond zero-dimensional integrated quantities to one-dimensional profiles. In agreement with the previous findings, a comparison of measured and simulated burn profiles shows that models including ion transport effects are able to better match the experimental results. In implosions characterized by large Knudsen numbers (N K ~ 3), the fusion burn profiles predicted by hydrodynamics simulations that exclude ion mean free path effects are peaked far from the origin, in stark disagreement with the experimentally observed profiles, which are centrally peaked. In contrast, a hydrodynamics simulation that includes a model of ion diffusion is able to qualitatively match the measured profile shapes. Therefore, ion diffusion or diffusion-like processes are identified as a plausible explanation of the observed trends, though further refinement of the models is needed for a more complete and quantitative understanding of ion kinetic effects.« less

Assessment of ion kinetic effects in shock-driven inertial confinement fusion implosions using fusion burn imaging

NASA Astrophysics Data System (ADS)

Rosenberg, M. J.; Séguin, F. H.; Amendt, P. A.; Atzeni, S.; Rinderknecht, H. G.; Hoffman, N. M.; Zylstra, A. B.; Li, C. K.; Sio, H.; Gatu Johnson, M.; Frenje, J. A.; Petrasso, R. D.; Glebov, V. Yu.; Stoeckl, C.; Seka, W.; Marshall, F. J.; Delettrez, J. A.; Sangster, T. C.; Betti, R.; Wilks, S. C.; Pino, J.; Kagan, G.; Molvig, K.; Nikroo, A.

2015-06-01

The significance and nature of ion kinetic effects in D3He-filled, shock-driven inertial confinement fusion implosions are assessed through measurements of fusion burn profiles. Over this series of experiments, the ratio of ion-ion mean free path to minimum shell radius (the Knudsen number, NK) was varied from 0.3 to 9 in order to probe hydrodynamic-like to strongly kinetic plasma conditions; as the Knudsen number increased, hydrodynamic models increasingly failed to match measured yields, while an empirically-tuned, first-step model of ion kinetic effects better captured the observed yield trends [Rosenberg et al., Phys. Rev. Lett. 112, 185001 (2014)]. Here, spatially resolved measurements of the fusion burn are used to examine kinetic ion transport effects in greater detail, adding an additional dimension of understanding that goes beyond zero-dimensional integrated quantities to one-dimensional profiles. In agreement with the previous findings, a comparison of measured and simulated burn profiles shows that models including ion transport effects are able to better match the experimental results. In implosions characterized by large Knudsen numbers (NK ˜ 3), the fusion burn profiles predicted by hydrodynamics simulations that exclude ion mean free path effects are peaked far from the origin, in stark disagreement with the experimentally observed profiles, which are centrally peaked. In contrast, a hydrodynamics simulation that includes a model of ion diffusion is able to qualitatively match the measured profile shapes. Therefore, ion diffusion or diffusion-like processes are identified as a plausible explanation of the observed trends, though further refinement of the models is needed for a more complete and quantitative understanding of ion kinetic effects.

Assessment of ion kinetic effects in shock-driven inertial confinement fusion (ICF) implosions using fusion burn imaging

DOE PAGES

Rosenberg, M. J.; Séguin, F. H.; Amendt, P. A.; ...

2015-06-02

The significance and nature of ion kinetic effects in D³He-filled, shock-driven inertial confinement fusion implosions are assessed through measurements of fusion burn profiles. Over this series of experiments, the ratio of ion-ion mean free path to minimum shell radius (the Knudsen number, N K) was varied from 0.3 to 9 in order to probe hydrodynamic-like to strongly kinetic plasma conditions; as the Knudsen number increased, hydrodynamic models increasingly failed to match measured yields, while an empirically-tuned, first-step model of ion kinetic effects better captured the observed yield trends [Rosenberg et al., Phys. Rev. Lett. 112, 185001 (2014)]. Here, spatially resolvedmore » measurements of the fusion burn are used to examine kinetic ion transport effects in greater detail, adding an additional dimension of understanding that goes beyond zero-dimensional integrated quantities to one-dimensional profiles. In agreement with the previous findings, a comparison of measured and simulated burn profiles shows that models including ion transport effects are able to better match the experimental results. In implosions characterized by large Knudsen numbers (N K ~ 3), the fusion burn profiles predicted by hydrodynamics simulations that exclude ion mean free path effects are peaked far from the origin, in stark disagreement with the experimentally observed profiles, which are centrally peaked. In contrast, a hydrodynamics simulation that includes a model of ion diffusion is able to qualitatively match the measured profile shapes. Therefore, ion diffusion or diffusion-like processes are identified as a plausible explanation of the observed trends, though further refinement of the models is needed for a more complete and quantitative understanding of ion kinetic effects.« less

Crystallographic Lattice Boltzmann Method

PubMed Central

Namburi, Manjusha; Krithivasan, Siddharth; Ansumali, Santosh

2016-01-01

Current approaches to Direct Numerical Simulation (DNS) are computationally quite expensive for most realistic scientific and engineering applications of Fluid Dynamics such as automobiles or atmospheric flows. The Lattice Boltzmann Method (LBM), with its simplified kinetic descriptions, has emerged as an important tool for simulating hydrodynamics. In a heterogeneous computing environment, it is often preferred due to its flexibility and better parallel scaling. However, direct simulation of realistic applications, without the use of turbulence models, remains a distant dream even with highly efficient methods such as LBM. In LBM, a fictitious lattice with suitable isotropy in the velocity space is considered to recover Navier-Stokes hydrodynamics in macroscopic limit. The same lattice is mapped onto a cartesian grid for spatial discretization of the kinetic equation. In this paper, we present an inverted argument of the LBM, by making spatial discretization as the central theme. We argue that the optimal spatial discretization for LBM is a Body Centered Cubic (BCC) arrangement of grid points. We illustrate an order-of-magnitude gain in efficiency for LBM and thus a significant progress towards feasibility of DNS for realistic flows. PMID:27251098

Coulomb-driven energy boost of heavy ions for laser-plasma acceleration.

PubMed

Braenzel, J; Andreev, A A; Platonov, K; Klingsporn, M; Ehrentraut, L; Sandner, W; Schnürer, M

2015-03-27

An unprecedented increase of kinetic energy of laser accelerated heavy ions is demonstrated. Ultrathin gold foils have been irradiated by an ultrashort laser pulse at a peak intensity of 8×10^{19}  W/  cm^{2}. Highly charged gold ions with kinetic energies up to >200  MeV and a bandwidth limited energy distribution have been reached by using 1.3 J laser energy on target. 1D and 2D particle in cell simulations show how a spatial dependence on the ion's ionization leads to an enhancement of the accelerating electrical field. Our theoretical model considers a spatial distribution of the ionization inside the thin target, leading to a field enhancement for the heavy ions by Coulomb explosion. It is capable of explaining the energy boost of highly charged ions, enabling a higher efficiency for the laser-driven heavy ion acceleration.

A Spatially Continuous Model of Carbohydrate Digestion and Transport Processes in the Colon

PubMed Central

Moorthy, Arun S.; Brooks, Stephen P. J.; Kalmokoff, Martin; Eberl, Hermann J.

2015-01-01

A spatially continuous mathematical model of transport processes, anaerobic digestion and microbial complexity as would be expected in the human colon is presented. The model is a system of first-order partial differential equations with context determined number of dependent variables, and stiff, non-linear source terms. Numerical simulation of the model is used to elucidate information about the colon-microbiota complex. It is found that the composition of materials on outflow of the model does not well-describe the composition of material in other model locations, and inferences using outflow data varies according to model reactor representation. Additionally, increased microbial complexity allows the total microbial community to withstand major system perturbations in diet and community structure. However, distribution of strains and functional groups within the microbial community can be modified depending on perturbation length and microbial kinetic parameters. Preliminary model extensions and potential investigative opportunities using the computational model are discussed. PMID:26680208

Modeling of Ureolytic Calcite Precipitation for the Remediation of Sr-90 Using a Variable Velocity Streamtube Ensemble

NASA Astrophysics Data System (ADS)

Weathers, T. S.; Ginn, T. R.; Spycher, N.; Barkouki, T. H.; Fujita, Y.; Smith, R. W.

2009-12-01

Subsurface contamination is often mitigated with an injection/extraction well system. An understanding of heterogeneities within this radial flowfield is critical for modeling, prediction, and remediation of the subsurface. We address this using a Lagrangian approach: instead of depicting spatial extents of solutes in the subsurface we focus on their arrival distribution at the control well(s). A well-to-well treatment system that incorporates in situ microbially-mediated ureolysis to induce calcite precipitation for the immobilization of strontium-90 has been explored at the Vadose Zone Research Park (VZRP) near Idaho Falls, Idaho. PHREEQC2 is utilized to model the kinetically-controlled ureolysis and consequent calcite precipitation. PHREEQC2 provides a one-dimensional advective-dispersive transport option that can be and has been used in streamtube ensemble models. Traditionally, each streamtube maintains uniform velocity; however in radial flow in homogeneous media, the velocity within any given streamtube is variable in space, being highest at the input and output wells and approaching a minimum at the midpoint between the wells. This idealized velocity variability is of significance if kinetic reactions are present with multiple components, if kinetic reaction rates vary in space, if the reactions involve multiple phases (e.g. heterogeneous reactions), and/or if they impact physical characteristics (porosity/permeability), as does ureolytically driven calcite precipitation. Streamtube velocity patterns for any particular configuration of injection and withdrawal wells are available as explicit calculations from potential theory, and also from particle tracking programs. To approximate the actual spatial distribution of velocity along streamtubes, we assume idealized non-uniform velocity associated with homogeneous media. This is implemented in PHREEQC2 via a non-uniform spatial discretization within each streamtube that honors both the streamtube’s travel time and the idealized “fast-slow-fast” nonuniform velocity along the streamline. Breakthrough curves produced by each simulation are weighted by the path-respective flux fractions (obtained by deconvolution of tracer tests conducted at the VZRP) to obtain the flux-average of flow contributions to the observation well. Breakthrough data from urea injection experiments performed at the VZRP are compared to the model results from the PHREEQC2 variable velocity ensemble.

Exact solutions and low-frequency instability of the adiabatic auroral arc model

NASA Technical Reports Server (NTRS)

Cornwall, John M.

1988-01-01

The adiabatic auroral arc model couples a kinetic theory parallel current driven by mirror forces to horizontal ionospheric currents; the resulting equations are nonlinear. Some exact stationary solutions to these equations, some of them based on the Liouville equation, are developed, with both latitudinal and longitudinal spatial variations. These Liouville equation exact solutions are related to stability boundaries of low-frequency instabilities such as Kelvin-Helmholtz, as shown by a study of a simplified model.

Accounting for Space—Quantification of Cell-To-Cell Transmission Kinetics Using Virus Dynamics Models.

PubMed

Kumberger, Peter; Durso-Cain, Karina; Uprichard, Susan L; Dahari, Harel; Graw, Frederik

2018-04-17

Mathematical models based on ordinary differential equations (ODE) that describe the population dynamics of viruses and infected cells have been an essential tool to characterize and quantify viral infection dynamics. Although an important aspect of viral infection is the dynamics of viral spread, which includes transmission by cell-free virions and direct cell-to-cell transmission, models used so far ignored cell-to-cell transmission completely, or accounted for this process by simple mass-action kinetics between infected and uninfected cells. In this study, we show that the simple mass-action approach falls short when describing viral spread in a spatially-defined environment. Using simulated data, we present a model extension that allows correct quantification of cell-to-cell transmission dynamics within a monolayer of cells. By considering the decreasing proportion of cells that can contribute to cell-to-cell spread with progressing infection, our extension accounts for the transmission dynamics on a single cell level while still remaining applicable to standard population-based experimental measurements. While the ability to infer the proportion of cells infected by either of the transmission modes depends on the viral diffusion rate, the improved estimates obtained using our novel approach emphasize the need to correctly account for spatial aspects when analyzing viral spread.

A physiology-based parametric imaging method for FDG-PET data

NASA Astrophysics Data System (ADS)

Scussolini, Mara; Garbarino, Sara; Sambuceti, Gianmario; Caviglia, Giacomo; Piana, Michele

2017-12-01

Parametric imaging is a compartmental approach that processes nuclear imaging data to estimate the spatial distribution of the kinetic parameters governing tracer flow. The present paper proposes a novel and efficient computational method for parametric imaging which is potentially applicable to several compartmental models of diverse complexity and which is effective in the determination of the parametric maps of all kinetic coefficients. We consider applications to [18 F]-fluorodeoxyglucose positron emission tomography (FDG-PET) data and analyze the two-compartment catenary model describing the standard FDG metabolization by an homogeneous tissue and the three-compartment non-catenary model representing the renal physiology. We show uniqueness theorems for both models. The proposed imaging method starts from the reconstructed FDG-PET images of tracer concentration and preliminarily applies image processing algorithms for noise reduction and image segmentation. The optimization procedure solves pixel-wise the non-linear inverse problem of determining the kinetic parameters from dynamic concentration data through a regularized Gauss-Newton iterative algorithm. The reliability of the method is validated against synthetic data, for the two-compartment system, and experimental real data of murine models, for the renal three-compartment system.

Developments in FTIR spectroscopy of diamonds and better constraints on diamond thermal histories

NASA Astrophysics Data System (ADS)

Kohn, Simon; Speich, Laura; Smith, Christopher; Bulanova, Galina

2017-04-01

Fourier Transform Infrared (FTIR) spectroscopy is a commonly-used technique for investigating diamonds. It gives the most useful information if spatially-resolved measurements are used [1]. In this contribution we discuss the best way to acquire and present FTIR data from diamonds, using examples from Murowa (Zimbabwe), Argyle (Australia) and Machado River (Brazil). Examples of FTIR core-to-rim line scans, maps with high spatial resolution and maps with high spectral resolution that are fitted to extract the spatial variation of different nitrogen and hydrogen defects are presented. Model mantle residence temperatures are calculated from the concentration of A and B nitrogen-containing defects in the diamonds using known times of annealing in the mantle. A new, two-stage thermal annealing model is presented that better constrains the thermal history of the diamond and that of the mantle lithosphere in which the diamond resided. The effect of heterogeneity within the analysed FTIR volume is quantitatively assessed and errors in model temperatures that can be introduced by studying whole diamonds instead of thin plates are discussed. The kinetics of platelet growth and degradation will be discussed and the potential for two separate, kinetically-controlled defect reactions to be used to constrain a full thermal history of the diamond will be assessed. [1] Kohn, S.C., Speich, L., Smith, C.B. and Bulanova, G.P., 2016. FTIR thermochronometry of natural diamonds: A closer look. Lithos, 265, pp.148-158.

The scaling of contact rates with population density for the infectious disease models.

PubMed

Hu, Hao; Nigmatulina, Karima; Eckhoff, Philip

2013-08-01

Contact rates and patterns among individuals in a geographic area drive transmission of directly-transmitted pathogens, making it essential to understand and estimate contacts for simulation of disease dynamics. Under the uniform mixing assumption, one of two mechanisms is typically used to describe the relation between contact rate and population density: density-dependent or frequency-dependent. Based on existing evidence of population threshold and human mobility patterns, we formulated a spatial contact model to describe the appropriate form of transmission with initial growth at low density and saturation at higher density. We show that the two mechanisms are extreme cases that do not capture real population movement across all scales. Empirical data of human and wildlife diseases indicate that a nonlinear function may work better when looking at the full spectrum of densities. This estimation can be applied to large areas with population mixing in general activities. For crowds with unusually large densities (e.g., transportation terminals, stadiums, or mass gatherings), the lack of organized social contact structure deviates the physical contacts towards a special case of the spatial contact model - the dynamics of kinetic gas molecule collision. In this case, an ideal gas model with van der Waals correction fits well; existing movement observation data and the contact rate between individuals is estimated using kinetic theory. A complete picture of contact rate scaling with population density may help clarify the definition of transmission rates in heterogeneous, large-scale spatial systems. Copyright © 2013 The Authors. Published by Elsevier Inc. All rights reserved.

Spatial distribution of angular momentum inside the nucleon

NASA Astrophysics Data System (ADS)

Lorcé, Cédric; Mantovani, Luca; Pasquini, Barbara

2018-01-01

We discuss in detail the spatial distribution of angular momentum inside the nucleon. We show that the discrepancies between different definitions originate from terms that integrate to zero. Even though these terms can safely be dropped at the integrated level, they have to be taken into account when discussing distributions. Using the scalar diquark model, we illustrate our results and, for the first time, check explicitly that the equivalence between kinetic and canonical orbital angular momentum persists at the level of distributions, as expected in a system without gauge degrees of freedom.

Kinetics of radiation-induced precipitation at the alloy surface

NASA Astrophysics Data System (ADS)

Lam, N. Q.; Nguyen, T.; Leaf, G. K.; Yip, S.

1988-05-01

Radiation-induced precipitation of a new phase at the surface of an alloy during irradiation at elevated temperatures was studied with the aid of a kinetic model of segregation. The preferential coupling of solute atoms with the defect fluxes gives rise to a strong solute enrichment at the surface, which, if surpassing the solute solubility limit, leads to the formation of a precipitate layer. The moving precipitate/matrix interface was accommodated by means of a mathematical scheme that transforms spatial coordinates into a reference frame in which the boundaries are immobile. Sample calculations were performed for precipitation of the γ'-Ni 3Si layer on Ni-Si alloys undergoing electron irradiation. The dependences of the precipitation kinetics on the defect-production rate, irradiation temperature, internal defect sink concentration and alloy composition were investigated systematically.

Kinetics of the maintenance of the epidermis

NASA Astrophysics Data System (ADS)

Zhdanov, Vladimir P.; Cho, Nam-Joon

2013-08-01

The epidermis is the outermost layer of skin. It is comprised of keratin-containing cells called keratinocytes. Functionally, the epidermis serves as a physical barrier that can prevent infection and regulate body hydration. Maintenance and repair of the epidermis are important for human health. Mechanistically, these processes occur primarily via proliferation and differentiation of stem cells located in the basal monolayer. These processes are believed to depend on cell-cell communication and spatial constraints but existing kinetic models focus mainly on proliferation and differentiation. To address this issue, we present a mean-field kinetic model that takes these additional factors into account and describes the epidermis at a biosystem level. The corresponding equations operate with the populations of stem cells and differentiated cells in the basal layer. The keratinocytes located above the basal layer are treated at a more coarse-grained level by considering the thickness of the epidermis. The model clarifies the likely role of various negative feedbacks that may control the epidermis and, accordingly, provides insight into the cellular mechanisms underlying complex biological phenomena such as wound healing.

Heat and turbulent kinetic energy budgets for surface layer cooling induced by the passage of Hurricane Frances (2004)

NASA Astrophysics Data System (ADS)

Huang, Peisheng; Sanford, Thomas B.; Imberger, JöRg

2009-12-01

Heat and turbulent kinetic energy budgets of the ocean surface layer during the passage of Hurricane Frances were examined using a three-dimensional hydrodynamic model. In situ data obtained with the Electromagnetic-Autonomous Profiling Explorer (EM-APEX) floats were used to set up the initial conditions of the model simulation and to compare to the simulation results. The spatial heat budgets reveal that during the hurricane passage, not only the entrainment in the bottom of surface mixed layer but also the horizontal water advection were important factors determining the spatial pattern of sea surface temperature. At the free surface, the hurricane-brought precipitation contributed a negligible amount to the air-sea heat exchange, but the precipitation produced a negative buoyancy flux in the surface layer that overwhelmed the instability induced by the heat loss to the atmosphere. Integrated over the domain within 400 km of the hurricane eye on day 245.71 of 2004, the rate of heat anomaly in the surface water was estimated to be about 0.45 PW (1 PW = 1015 W), with about 20% (0.09 PW in total) of this was due to the heat exchange at the air-sea interface, and almost all the remainder (0.36 PW) was downward transported by oceanic vertical mixing. Shear production was the major source of turbulent kinetic energy amounting 88.5% of the source of turbulent kinetic energy, while the rest (11.5%) was attributed to the wind stirring at sea surface. The increase of ocean potential energy due to vertical mixing represented 7.3% of the energy deposited by wind stress.

KINETICS OF LOW SOURCE REACTOR STARTUPS. PART II

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

hurwitz, H. Jr.; MacMillan, D.B.; Smith, J.H.

1962-06-01

A computational technique is described for computation of the probability distribution of power level for a low source reactor startup. The technique uses a mathematical model, for the time-dependent probability distribution of neutron and precursor concentration, having finite neutron lifetime, one group of delayed neutron precursors, and no spatial dependence. Results obtained by the technique are given. (auth)

Decay and Spatial Diffusion of Turbulent Kinetic Energy In The Presence of a Linear Kinetic Energy Gradient

NASA Astrophysics Data System (ADS)

Meneveau, Charles

2015-11-01

A topic that elicited the interest of John Lumley is pressure transport in turbulence. In 1978 (JL, in Advances in Applied Mechanics, pages 123-176) he showed that pressure transport likely acts in the opposite direction to the spatial flux of kinetic energy due to triple velocity correlations. Here we examine a flow in which the interplay of turbulent decay and spatial transport is particularly relevant. Specifically, using a specially designed active grid and screens placed in the Corrsin wind tunnel, such a flow is realized. Data are acquired using X-wire thermal anemometry at different spanwise and downstream locations. In order to resolve the dissipation rate accurately, measurements are also acquired using the NSTAP probe developed and manufactured by Princeton researchers and kindly provided to us (M. Hultmark, Y. Fan, L. Smits). The results show power-law decay with downstream distance, with a decay exponent that becomes larger in the high kinetic energy side of the flow. Measurements of the dissipation enable us to obtain the spanwise gradient of the spatial flux. One possible explanation for the observations is upgrading transport of kinetic energy due to pressure-velocity correlations, although its magnitude required to close the budget appears very large. Absence of simultaneous pressure velocity measurement preclude us to fully elucidate the observed trends. In collaboration with Adrien Thormann, Johns Hopkins University. Financial support: National Science Foundation.

Analyzing and modeling the kinetics of amyloid beta pores associated with Alzheimer’s disease pathology

DOE PAGES

Ullah, Ghanim; Demuro, Angelo; Parker, Ian; ...

2015-09-08

Amyloid beta (Aβ) oligomers associated with Alzheimer’s disease (AD) form Ca 2+-permeable plasma membrane pores, leading to a disruption of the otherwise well-controlled intracellular calcium (Ca 2+) homeostasis. The resultant up-regulation of intracellular Ca 2+ concentration has detrimental implications for memory formation and cell survival. The gating kinetics and Ca 2+ permeability of Aβ pores are not well understood. We have used computational modeling in conjunction with the ability of optical patch-clamping for massively parallel imaging of Ca 2+ flux through thousands of pores in the cell membrane of Xenopus oocytes to elucidate the kinetic properties of Aβ pores. Themore » fluorescence time-series data from individual pores were idealized and used to develop data-driven Markov chain models for the kinetics of the Aβ pore at different stages of its evolution. Our study provides the first demonstration of developing Markov chain models for ion channel gating that are driven by optical-patch clamp data with the advantage of experiments being performed under close to physiological conditions. As a result, we demonstrate the up-regulation of gating of various Ca 2+ release channels due to Aβ pores and show that the extent and spatial range of such up-regulation increases as Aβ pores with low open probability and Ca 2+ permeability transition into those with high open probability and Ca 2+ permeability.« less

Bed turbulent kinetic energy boundary conditions for trapping efficiency and spatial distribution of sediments in basins.

PubMed

Isenmann, Gilles; Dufresne, Matthieu; Vazquez, José; Mose, Robert

2017-10-01

The purpose of this study is to develop and validate a numerical tool for evaluating the performance of a settling basin regarding the trapping of suspended matter. The Euler-Lagrange approach was chosen to model the flow and sediment transport. The numerical model developed relies on the open source library OpenFOAM ® , enhanced with new particle/wall interaction conditions to limit sediment deposition in zones with favourable hydrodynamic conditions (shear stress, turbulent kinetic energy). In particular, a new relation is proposed for calculating the turbulent kinetic energy threshold as a function of the properties of each particle (diameter and density). The numerical model is compared to three experimental datasets taken from the literature and collected for scale models of basins. The comparison of the numerical and experimental results permits concluding on the model's capacity to predict the trapping of particles in a settling basin with an absolute error in the region of 5% when the sediment depositions occur over the entire bed. In the case of sediment depositions localised in preferential zones, their distribution is reproduced well by the model and trapping efficiency is evaluated with an absolute error in the region of 10% (excluding cases of particles with very low density).

Unbiased Rare Event Sampling in Spatial Stochastic Systems Biology Models Using a Weighted Ensemble of Trajectories

PubMed Central

Donovan, Rory M.; Tapia, Jose-Juan; Sullivan, Devin P.; Faeder, James R.; Murphy, Robert F.; Dittrich, Markus; Zuckerman, Daniel M.

2016-01-01

The long-term goal of connecting scales in biological simulation can be facilitated by scale-agnostic methods. We demonstrate that the weighted ensemble (WE) strategy, initially developed for molecular simulations, applies effectively to spatially resolved cell-scale simulations. The WE approach runs an ensemble of parallel trajectories with assigned weights and uses a statistical resampling strategy of replicating and pruning trajectories to focus computational effort on difficult-to-sample regions. The method can also generate unbiased estimates of non-equilibrium and equilibrium observables, sometimes with significantly less aggregate computing time than would be possible using standard parallelization. Here, we use WE to orchestrate particle-based kinetic Monte Carlo simulations, which include spatial geometry (e.g., of organelles, plasma membrane) and biochemical interactions among mobile molecular species. We study a series of models exhibiting spatial, temporal and biochemical complexity and show that although WE has important limitations, it can achieve performance significantly exceeding standard parallel simulation—by orders of magnitude for some observables. PMID:26845334

A case study of the Weather Research and Forecasting model applied to the Joint Urban 2003 tracer field experiment. Part 2: Gas tracer dispersion

DOE PAGES

Nelson, Matthew A.; Brown, Michael J.; Halverson, Scot A.; ...

2016-07-28

Here, the Quick Urban & Industrial Complex (QUIC) atmospheric transport, and dispersion modelling, system was evaluated against the Joint Urban 2003 tracer-gas measurements. This was done using the wind and turbulence fields computed by the Weather Research and Forecasting (WRF) model. We compare the simulated and observed plume transport when using WRF-model-simulated wind fields, and local on-site wind measurements. Degradation of the WRF-model-based plume simulations was cased by errors in the simulated wind direction, and limitations in reproducing the small-scale wind-field variability. We explore two methods for importing turbulence from the WRF model simulations into the QUIC system. The firstmore » method uses parametrized turbulence profiles computed from WRF-model-computed boundary-layer similarity parameters; and the second method directly imports turbulent kinetic energy from the WRF model. Using the WRF model’s Mellor-Yamada-Janjic boundary-layer scheme, the parametrized turbulence profiles and the direct import of turbulent kinetic energy were found to overpredict and underpredict the observed turbulence quantities, respectively. Near-source building effects were found to propagate several km downwind. These building effects and the temporal/spatial variations in the observed wind field were often found to have a stronger influence over the lateral and vertical plume spread than the intensity of turbulence. Correcting the WRF model wind directions using a single observational location improved the performance of the WRF-model-based simulations, but using the spatially-varying flow fields generated from multiple observation profiles generally provided the best performance.« less

A case study of the Weather Research and Forecasting model applied to the Joint Urban 2003 tracer field experiment. Part 2: Gas tracer dispersion

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

Nelson, Matthew A.; Brown, Michael J.; Halverson, Scot A.

Here, the Quick Urban & Industrial Complex (QUIC) atmospheric transport, and dispersion modelling, system was evaluated against the Joint Urban 2003 tracer-gas measurements. This was done using the wind and turbulence fields computed by the Weather Research and Forecasting (WRF) model. We compare the simulated and observed plume transport when using WRF-model-simulated wind fields, and local on-site wind measurements. Degradation of the WRF-model-based plume simulations was cased by errors in the simulated wind direction, and limitations in reproducing the small-scale wind-field variability. We explore two methods for importing turbulence from the WRF model simulations into the QUIC system. The firstmore » method uses parametrized turbulence profiles computed from WRF-model-computed boundary-layer similarity parameters; and the second method directly imports turbulent kinetic energy from the WRF model. Using the WRF model’s Mellor-Yamada-Janjic boundary-layer scheme, the parametrized turbulence profiles and the direct import of turbulent kinetic energy were found to overpredict and underpredict the observed turbulence quantities, respectively. Near-source building effects were found to propagate several km downwind. These building effects and the temporal/spatial variations in the observed wind field were often found to have a stronger influence over the lateral and vertical plume spread than the intensity of turbulence. Correcting the WRF model wind directions using a single observational location improved the performance of the WRF-model-based simulations, but using the spatially-varying flow fields generated from multiple observation profiles generally provided the best performance.« less

Three-dimensional hybrid modeling of ion kinetic instabilities in space plasma

NASA Astrophysics Data System (ADS)

Ofman, L.

2017-12-01

Ion kinetic instabilities in space plasma are believed to play an imprortant role in energy transport, heating, dissipation of turbulence, as well as in generating of spectrum of magnetic fluctuations in the kinetic frequency range. The velocity distribution functions (VDFs) of unstable ion populations are generally non-Maxwellian and provide the free energy source that drives the waves. The VDFs were measured in-situ by satellites such as Helios, WIND, and would be obtained in the future Parkers' Solar Probe close to the Sun. In particular, temperature anisotropy provides a measure of VDF non-equilibroum structure, that together with parallel-beta determine the threshold of kinetic instabilities, such as mirror, ion-cyclotron, and firehose. Drifting population of alphas with respect to protons lead to the magnetosonic instability. So far, these isntabilities were studied primaraly using 1.5D or 2.5D particle-in-cell (PIC) or hybrid models (where electrons are modeled as a fluid), i.e., in 1 or 2 spatial dimensions with 3 components of velocity and magnetic field. I will present the results of recent full 3D hybrid models that studies these instabilities for heliospheric conditions and compare to previous modeling results. I will discuss the agreement and the differences between the 3D and more approximate models of the VDFs, the magnetic fluctuations spectra, and the temporal evolution of the anisotropy for typical instabilities relevant for space plasma. I will duscuss the use of the modeled VDFs for diagnostic of the physical processes that lead to space plasma energization from the observed VDFs in the heliospheric and magnetospheric plasma.

A third-order gas-kinetic CPR method for the Euler and Navier-Stokes equations on triangular meshes

NASA Astrophysics Data System (ADS)

Zhang, Chao; Li, Qibing; Fu, Song; Wang, Z. J.

2018-06-01

A third-order accurate gas-kinetic scheme based on the correction procedure via reconstruction (CPR) framework is developed for the Euler and Navier-Stokes equations on triangular meshes. The scheme combines the accuracy and efficiency of the CPR formulation with the multidimensional characteristics and robustness of the gas-kinetic flux solver. Comparing with high-order finite volume gas-kinetic methods, the current scheme is more compact and efficient by avoiding wide stencils on unstructured meshes. Unlike the traditional CPR method where the inviscid and viscous terms are treated differently, the inviscid and viscous fluxes in the current scheme are coupled and computed uniformly through the kinetic evolution model. In addition, the present scheme adopts a fully coupled spatial and temporal gas distribution function for the flux evaluation, achieving high-order accuracy in both space and time within a single step. Numerical tests with a wide range of flow problems, from nearly incompressible to supersonic flows with strong shocks, for both inviscid and viscous problems, demonstrate the high accuracy and efficiency of the present scheme.

Film growth kinetics and electric field patterning during electrospray deposition of block copolymer thin films

NASA Astrophysics Data System (ADS)

Toth, Kristof; Hu, Hanqiong; Choo, Youngwoo; Loewenberg, Michael; Osuji, Chinedum

The delivery of sub-micron droplets of dilute polymer solutions to a heated substrate by electrospray deposition (ESD) enables precisely controlled and continuous growth of block copolymer (BCP) thin films. Here we explore patterned deposition of BCP films by spatially varying the electric field at the substrate using an underlying charged grid, as well as film growth kinetics. Numerical analysis was performed to examine pattern fidelity by considering the trajectories of charged droplets during flight through imposed periodic field variations in the vicinity of the substrate. Our work uncovered an unexpected modality for improving the resolution of the patterning process via stronger field focusing through the use of a second oppositely charged grid beneath a primary focusing array, with an increase in highly localized droplet deposition on the intersecting nodes of the grid. Substrate coverage kinetics are considered for homopolymer deposition in the context of simple kinetic models incorporating temperature and molecular weight dependence of diffusivity. By contrast, film coverage kinetics for block copolymer depositions are additionally convoluted with preferential wetting and thickness-periodicity commensurability effects. NSF GRFP.

Particle-in-cell modeling for MJ scale dense plasma focus with varied anode shape

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

Link, A., E-mail: link6@llnl.gov; Halvorson, C., E-mail: link6@llnl.gov; Schmidt, A.

2014-12-15

Megajoule scale dense plasma focus (DPF) Z-pinches with deuterium gas fill are compact devices capable of producing 10{sup 12} neutrons per shot but past predictive models of large-scale DPF have not included kinetic effects such as ion beam formation or anomalous resistivity. We report on progress of developing a predictive DPF model by extending our 2D axisymmetric collisional kinetic particle-in-cell (PIC) simulations from the 4 kJ, 200 kA LLNL DPF to 1 MJ, 2 MA Gemini DPF using the PIC code LSP. These new simulations incorporate electrodes, an external pulsed-power driver circuit, and model the plasma from insulator lift-off throughmore » the pinch phase. To accommodate the vast range of relevant spatial and temporal scales involved in the Gemini DPF within the available computational resources, the simulations were performed using a new hybrid fluid-to-kinetic model. This new approach allows single simulations to begin in an electron/ion fluid mode from insulator lift-off through the 5-6 μs run-down of the 50+ cm anode, then transition to a fully kinetic PIC description during the run-in phase, when the current sheath is 2-3 mm from the central axis of the anode. Simulations are advanced through the final pinch phase using an adaptive variable time-step to capture the fs and sub-mm scales of the kinetic instabilities involved in the ion beam formation and neutron production. Validation assessments are being performed using a variety of different anode shapes, comparing against experimental measurements of neutron yield, neutron anisotropy and ion beam production.« less

The stability and slow dynamics of spot patterns in the 2D Brusselator model: The effect of open systems and heterogeneities

NASA Astrophysics Data System (ADS)

Tzou, J. C.; Ward, M. J.

2018-06-01

Spot patterns, whereby the activator field becomes spatially localized near certain dynamically-evolving discrete spatial locations in a bounded multi-dimensional domain, is a common occurrence for two-component reaction-diffusion (RD) systems in the singular limit of a large diffusivity ratio. In previous studies of 2-D localized spot patterns for various specific well-known RD systems, the domain boundary was assumed to be impermeable to both the activator and inhibitor, and the reaction-kinetics were assumed to be spatially uniform. As an extension of this previous theory, we use formal asymptotic methods to study the existence, stability, and slow dynamics of localized spot patterns for the singularly perturbed 2-D Brusselator RD model when the domain boundary is only partially impermeable, as modeled by an inhomogeneous Robin boundary condition, or when there is an influx of inhibitor across the domain boundary. In our analysis, we will also allow for the effect of a spatially variable bulk feed term in the reaction kinetics. By applying our extended theory to the special case of one-spot patterns and ring patterns of spots inside the unit disk, we provide a detailed analysis of the effect on spot patterns of these three different sources of heterogeneity. In particular, when there is an influx of inhibitor across the boundary of the unit disk, a ring pattern of spots can become pinned to a ring-radius closer to the domain boundary. Under a Robin condition, a quasi-equilibrium ring pattern of spots is shown to exhibit a novel saddle-node bifurcation behavior in terms of either the inhibitor diffusivity, the Robin constant, or the ambient background concentration. A spatially variable bulk feed term, with a concentrated source of "fuel" inside the domain, is shown to yield a saddle-node bifurcation structure of spot equilibria, which leads to qualitatively new spot-pinning behavior. Results from our asymptotic theory are validated from full numerical simulations of the Brusselator model.

Spatial averaging of a dissipative particle dynamics model for active suspensions

NASA Astrophysics Data System (ADS)

Panchenko, Alexander; Hinz, Denis F.; Fried, Eliot

2018-03-01

Starting from a fine-scale dissipative particle dynamics (DPD) model of self-motile point particles, we derive meso-scale continuum equations by applying a spatial averaging version of the Irving-Kirkwood-Noll procedure. Since the method does not rely on kinetic theory, the derivation is valid for highly concentrated particle systems. Spatial averaging yields stochastic continuum equations similar to those of Toner and Tu. However, our theory also involves a constitutive equation for the average fluctuation force. According to this equation, both the strength and the probability distribution vary with time and position through the effective mass density. The statistics of the fluctuation force also depend on the fine scale dissipative force equation, the physical temperature, and two additional parameters which characterize fluctuation strengths. Although the self-propulsion force entering our DPD model contains no explicit mechanism for aligning the velocities of neighboring particles, our averaged coarse-scale equations include the commonly encountered cubically nonlinear (internal) body force density.

Improved Enzyme Kinetic Model for Nitrification in Soils Amended with Ammonium. I. Literature Review,

DTIC Science & Technology

1980-01-01

ciple we can use these data in nitrification models, ON NITRIFICATION provided the concentration of oxygen in soil solution is known. This, however, is...Soil nitrifiers 30 0.8 (mixed culture) *As determined by Shah (1975). than within bulk soil solution , and McLaren and Packer these soils. The spatial...concentration concentration of 70 mg/l since this is the reported in soil solution depends on a soil-specific equilibrium influent concentration

Correlation functions in first-order phase transitions

NASA Astrophysics Data System (ADS)

Garrido, V.; Crespo, D.

1997-09-01

Most of the physical properties of systems underlying first-order phase transitions can be obtained from the spatial correlation functions. In this paper, we obtain expressions that allow us to calculate all the correlation functions from the droplet size distribution. Nucleation and growth kinetics is considered, and exact solutions are obtained for the case of isotropic growth by using self-similarity properties. The calculation is performed by using the particle size distribution obtained by a recently developed model (populational Kolmogorov-Johnson-Mehl-Avrami model). Since this model is less restrictive than that used in previously existing theories, the result is that the correlation functions can be obtained for any dependence of the kinetic parameters. The validity of the method is tested by comparison with the exact correlation functions, which had been obtained in the available cases by the time-cone method. Finally, the correlation functions corresponding to the microstructure developed in partitioning transformations are obtained.

Influence of the partial volume correction method on 18F-fluorodeoxyglucose brain kinetic modelling from dynamic PET images reconstructed with resolution model based OSEM

PubMed Central

Bowen, Spencer L.; Byars, Larry G.; Michel, Christian J.; Chonde, Daniel B.; Catana, Ciprian

2014-01-01

Kinetic parameters estimated from dynamic 18F-fluorodeoxyglucose PET acquisitions have been used frequently to assess brain function in humans. Neglecting partial volume correction (PVC) for a dynamic series has been shown to produce significant bias in model estimates. Accurate PVC requires a space-variant model describing the reconstructed image spatial point spread function (PSF) that accounts for resolution limitations, including non-uniformities across the field of view due to the parallax effect. For OSEM, image resolution convergence is local and influenced significantly by the number of iterations, the count density, and background-to-target ratio. As both count density and background-to-target values for a brain structure can change during a dynamic scan, the local image resolution may also concurrently vary. When PVC is applied post-reconstruction the kinetic parameter estimates may be biased when neglecting the frame-dependent resolution. We explored the influence of the PVC method and implementation on kinetic parameters estimated by fitting 18F-fluorodeoxyglucose dynamic data acquired on a dedicated brain PET scanner and reconstructed with and without PSF modelling in the OSEM algorithm. The performance of several PVC algorithms was quantified with a phantom experiment, an anthropomorphic Monte Carlo simulation, and a patient scan. Using the last frame reconstructed image only for regional spread function (RSF) generation, as opposed to computing RSFs for each frame independently, and applying perturbation GTM PVC with PSF based OSEM produced the lowest magnitude bias kinetic parameter estimates in most instances, although at the cost of increased noise compared to the PVC methods utilizing conventional OSEM. Use of the last frame RSFs for PVC with no PSF modelling in the OSEM algorithm produced the lowest bias in CMRGlc estimates, although by less than 5% in most cases compared to the other PVC methods. The results indicate that the PVC implementation and choice of PSF modelling in the reconstruction can significantly impact model parameters. PMID:24052021

Upscaling anomalous reactive kinetics (A+B-->C) from pore scale Lagrangian velocity analysis

NASA Astrophysics Data System (ADS)

De Anna, P.; Tartakovsky, A. M.; Le Borgne, T.; Dentz, M.

2011-12-01

Natural flow fields in porous media display a complex spatio-temporal organization due to heterogeneous geological structures at different scales. This multiscale disorder implies anomalous dispersion, mixing and reaction kinetics (Berkowitz et al. RG 2006, Tartakovsky PRE 2010). Here, we focus on the upscaling of anomalous kinetics arising from pore scale, non Gaussian and correlated, velocity distributions. We consider reactive front simulations, where a component A displaces a component B that saturates initially the porous domain. The reactive component C is produced at the dispersive front located at interface between the A and B domains. The simulations are performed with the SPH method. As the mixing zone grows, the total mass of C produced increases with time. The scaling of this evolution with time is different from that which would be obtained from the homogeneous advection dispersion reaction equation. This anomalous kinetics property is related to spatial structure of the reactive mixture, and its evolution with time under the combined action of advective and diffusive processes. We discuss the different scaling regimes arising depending on the dominant process that governs mixing. In order to upscale these processes, we analyze the Lagrangian velocity properties, which are characterized by the non Gaussian distributions and long range temporal correlation. The main origin of these properties is the existence of very low velocity regions where solute particles can remain trapped for a long time. Another source of strong correlation is the channeling of flow in localized high velocity regions, which created finger-like structures in the concentration field. We show the spatial Markovian, and temporal non Markovian, nature of the Lagrangian velocity field. Therefore, an upscaled model can be defined as a correlated Continuous Time Random Walk (Le Borgne et al. PRL 2008). A key feature of this model is the definition of a transition probability density for Lagrangian velocities across a characteristic correlation distance. We quantify this transition probability density from pore scale simulations and use it in the effective stochastic model. In this framework, we investigate the ability of this effective model to represent correctly dispersion and mixing.

Reductive dechlorination of trichloroethene DNAPL source zones: source zone architecture versus electron donor availability

NASA Astrophysics Data System (ADS)

Krol, M.; Kokkinaki, A.; Sleep, B.

2014-12-01

The persistence of dense-non-aqueous-phase liquids (DNAPLs) in the subsurface has led practitioners and regulatory agencies to turn towards low-maintenance, low-cost remediation methods. Biological degradation has been suggested as a possible solution, based on the well-proven ability of certain microbial species to break down dissolved chlorinated ethenes under favorable conditions. However, the biodegradation of pure phase chlorinated ethenes is subject to additional constraints: the continuous release of electron acceptor at a rate governed by mass transfer kinetics, and the temporal and spatial heterogeneity of DNAPL source zones which leads to spatially and temporally variable availability of the reactants for reductive dechlorination. In this work, we investigate the relationship between various DNAPL source zone characteristics and reaction kinetics using COMPSIM, a multiphase groundwater model that considers non-equilibrium mass transfer and Monod-type kinetics for reductive dechlorination. Numerical simulations are performed for simple, homogeneous trichloroethene DNAPL source zones to demonstrate the effect of single source zone characteristics, as well as for larger, more realistic heterogeneous source zones. It is shown that source zone size, and mass transfer kinetics may have a decisive effect on the predicted bio-enhancement. Finally, we evaluate the performance of DNAPL bioremediation for realistic, thermodynamically constrained, concentrations of electron donor. Our results indicate that the latter may be the most important limitation for the success of DNAPL bioremediation, leading to reduced bio-enhancement and, in many cases, comparable performance with water flooding.

On the Limiting Markov Process of Energy Exchanges in a Rarely Interacting Ball-Piston Gas

NASA Astrophysics Data System (ADS)

Bálint, Péter; Gilbert, Thomas; Nándori, Péter; Szász, Domokos; Tóth, Imre Péter

2017-02-01

We analyse the process of energy exchanges generated by the elastic collisions between a point-particle, confined to a two-dimensional cell with convex boundaries, and a `piston', i.e. a line-segment, which moves back and forth along a one-dimensional interval partially intersecting the cell. This model can be considered as the elementary building block of a spatially extended high-dimensional billiard modeling heat transport in a class of hybrid materials exhibiting the kinetics of gases and spatial structure of solids. Using heuristic arguments and numerical analysis, we argue that, in a regime of rare interactions, the billiard process converges to a Markov jump process for the energy exchanges and obtain the expression of its generator.

Phase-space analysis of the Schwinger effect in inhomogeneous electromagnetic fields

NASA Astrophysics Data System (ADS)

Kohlfürst, Christian

2018-05-01

Schwinger pair production in spatially and temporally inhomogeneous electric and magnetic fields is studied. The focus is on the particle phase-space distribution within a high-intensity few-cycle pulse. Accurate numerical solutions of a quantum kinetic theory (DHW formalism) are presented in momentum space and, with the aid of coarse-graining techniques, in a mixed spatial-momentum representation. Additionally, signatures of the carrier-envelope phase as well as spin-field interactions are discussed on the basis of a trajectory-based model taking into account instantaneous pair production and relativistic single-particle dynamics. Although our simple semi-classical single-particle model cannot describe every aspect of the particle production process (quantum interferences), essential features such as spin-field interactions are captured.

Orestes Kinetics Model for the Electra KrF Laser

NASA Astrophysics Data System (ADS)

Giuliani, J. L.; Kepple, P.; Lehmberg, R. H.; Myers, M. C.; Sethian, J. D.; Petrov, G.; Wolford, M.; Hegeler, F.

2003-10-01

Orestes is a first principles simulation code for the electron deposition, plasma chemistry, laser transport, and amplified spontaneous emission (ASE) in an e-beam pumped KrF laser. Orestes has been benchmarked against results from Nike at NRL and the Keio laser facility. The modeling tasks are to support ongoing oscillator experiments on the Electra laser ( 500 J), to predict performance of Electra as an amplifier, and to develop scaling relations for larger systems such as envisioned for an inertial fusion energy power plant. In Orestes the energy deposition of the primary beam electrons is assumed to be spatially uniform, but the excitation and ionization of the Ar/Kr/F2 target gas by the secondary electrons is determined from the energy distribution function as calculated by a Boltzmann code. The subsequent plasma kinetics of 23 species subject to over 100 reactions is followed with 1-D spatial resolution along the lasing axis. In addition, the vibrational relaxation among excited electronic states of the KrF molecule are included in the kinetics since lasing at 248 nm can occur from several vibrational lines of the B state. Transport of the lasing photons is solved by the method of characteristics. The time dependent ASE is calculated in 3-D using a ``local look-back'' scheme with discrete ordinates and includes specular reflection off the side walls and rear mirror. Gain narrowing is treated by multi-frequency transport of the ASE. Calculations for the gain, saturation intensity, extraction efficiency, and laser output from the Orestes model will be presented and compared with available data from Electra operated as an oscillator. Potential implications for the difference in optimal F2 concentration will be discussed along with the effects of window transmissivity at 248 nm.

Whole-body PET parametric imaging employing direct 4D nested reconstruction and a generalized non-linear Patlak model

NASA Astrophysics Data System (ADS)

Karakatsanis, Nicolas A.; Rahmim, Arman

2014-03-01

Graphical analysis is employed in the research setting to provide quantitative estimation of PET tracer kinetics from dynamic images at a single bed. Recently, we proposed a multi-bed dynamic acquisition framework enabling clinically feasible whole-body parametric PET imaging by employing post-reconstruction parameter estimation. In addition, by incorporating linear Patlak modeling within the system matrix, we enabled direct 4D reconstruction in order to effectively circumvent noise amplification in dynamic whole-body imaging. However, direct 4D Patlak reconstruction exhibits a relatively slow convergence due to the presence of non-sparse spatial correlations in temporal kinetic analysis. In addition, the standard Patlak model does not account for reversible uptake, thus underestimating the influx rate Ki. We have developed a novel whole-body PET parametric reconstruction framework in the STIR platform, a widely employed open-source reconstruction toolkit, a) enabling accelerated convergence of direct 4D multi-bed reconstruction, by employing a nested algorithm to decouple the temporal parameter estimation from the spatial image update process, and b) enhancing the quantitative performance particularly in regions with reversible uptake, by pursuing a non-linear generalized Patlak 4D nested reconstruction algorithm. A set of published kinetic parameters and the XCAT phantom were employed for the simulation of dynamic multi-bed acquisitions. Quantitative analysis on the Ki images demonstrated considerable acceleration in the convergence of the nested 4D whole-body Patlak algorithm. In addition, our simulated and patient whole-body data in the postreconstruction domain indicated the quantitative benefits of our extended generalized Patlak 4D nested reconstruction for tumor diagnosis and treatment response monitoring.

Front propagation in one-dimensional spatially periodic bistable media

NASA Astrophysics Data System (ADS)

Löber, Jakob; Bär, Markus; Engel, Harald

2012-12-01

Front propagation in heterogeneous bistable media is studied using the Schlögl model as a representative example. Spatially periodic modulations in the parameters of the bistable kinetics are taken into account perturbatively. Depending on the ratio L/l (L is the spatial period of the heterogeneity, l is the front width), appropriate singular perturbation techniques are applied to derive an ordinary differential equation for the position of the front in the presence of the heterogeneities. From this equation, the dependence of the average propagation speed on L/l as well as on the modulation amplitude is calculated. The analytical results obtained predict velocity overshoot, different cases of propagation failure, and the propagation speed for very large spatial periods in quantitative agreement with the results of direct numerical simulations of the underlying reaction-diffusion equation.

Modeling Np and Pu transport with a surface complexation model and spatially variant sorption capacities: Implications for reactive transport modeling and performance assessments of nuclear waste disposal sites

USGS Publications Warehouse

Glynn, P.D.

2003-01-01

One-dimensional (1D) geochemical transport modeling is used to demonstrate the effects of speciation and sorption reactions on the ground-water transport of Np and Pu, two redox-sensitive elements. Earlier 1D simulations (Reardon, 1981) considered the kinetically limited dissolution of calcite and its effect on ion-exchange reactions (involving 90Sr, Ca, Na, Mg and K), and documented the spatial variation of a 90Sr partition coefficient under both transient and steady-state chemical conditions. In contrast, the simulations presented here assume local equilibrium for all reactions, and consider sorption on constant potential, rather than constant charge, surfaces. Reardon's (1981) seminal findings on the spatial and temporal variability of partitioning (of 90Sr) are reexamined and found partially caused by his assumption of a kinetically limited reaction. In the present work, sorption is assumed the predominant retardation process controlling Pu and Np transport, and is simulated using a diffuse-double-layer-surface-complexation (DDLSC) model. Transport simulations consider the infiltration of Np- and Pu-contaminated waters into an initially uncontaminated environment, followed by the cleanup of the resultant contamination with uncontaminated water. Simulations are conducted using different spatial distributions of sorption capacities (with the same total potential sorption capacity, but with different variances and spatial correlation structures). Results obtained differ markedly from those that would be obtained in transport simulations using constant Kd, Langmuir or Freundlich sorption models. When possible, simulation results (breakthrough curves) are fitted to a constant K d advection-dispersion transport model and compared. Functional differences often are great enough that they prevent a meaningful fit of the simulation results with a constant K d (or even a Langmuir or Freundlich) model, even in the case of Np, a weakly sorbed radionuclide under the simulation conditions. Functional behaviors that cannot be fit include concentration trend reversals and radionuclide desorption spikes. Other simulation results are fit successfully but the fitted parameters (Kd and dispersivity) vary significantly depending on simulation conditions (e.g. "infiltration" vs. "cleanup" conditions). Notably, an increase in the variance of the specified sorption capacities results in a marked increase in the dispersion of the radionuclides. The results presented have implications for the simulation of radionuclide migration in performance assessments of nuclear waste-disposal sites, for the future monitoring of those sites, and more generally for modeling contaminant transport in ground-water environments. ?? 2003 Published by Elsevier Science Ltd.

Reliability of spatiotemporal and kinetic gait parameters determined by a new instrumented treadmill system.

PubMed

Reed, Lloyd F; Urry, Stephen R; Wearing, Scott C

2013-08-21

Despite the emerging use of treadmills integrated with pressure platforms as outcome tools in both clinical and research settings, published evidence regarding the measurement properties of these new systems is limited. This study evaluated the within- and between-day repeatability of spatial, temporal and vertical ground reaction force parameters measured by a treadmill system instrumented with a capacitance-based pressure platform. Thirty three healthy adults (mean age, 21.5 ± 2.8 years; height, 168.4 ± 9.9 cm; and mass, 67.8 ± 18.6 kg), walked barefoot on a treadmill system (FDM-THM-S, Zebris Medical GmbH) on three separate occasions. For each testing session, participants set their preferred pace but were blinded to treadmill speed. Spatial (foot rotation, step width, stride and step length), temporal (stride and step times, duration of stance, swing and single and double support) and peak vertical ground reaction force variables were collected over a 30-second capture period, equating to an average of 52 ± 5 steps of steady-state walking. Testing was repeated one week following the initial trial and again, for a third time, 20 minutes later. Repeated measures ANOVAs within a generalized linear modelling framework were used to assess between-session differences in gait parameters. Agreement between gait parameters measured within the same day (session 2 and 3) and between days (session 1 and 2; 1 and 3) were evaluated using the 95% repeatability coefficient. There were statistically significant differences in the majority (14/16) of temporal, spatial and kinetic gait parameters over the three test sessions (P < .01). The minimum change that could be detected with 95% confidence ranged between 3% and 17% for temporal parameters, 14% and 33% for spatial parameters, and 4% and 20% for kinetic parameters between days. Within-day repeatability was similar to that observed between days. Temporal and kinetic gait parameters were typically more consistent than spatial parameters. The 95% repeatability coefficient for vertical force peaks ranged between ± 53 and ± 63 N. The limits of agreement in spatial parameters and ground reaction forces for the treadmill system encompass previously reported changes with neuromuscular pathology and footwear interventions. These findings provide clinicians and researchers with an indication of the repeatability and sensitivity of the Zebris treadmill system to detect changes in common spatiotemporal gait parameters and vertical ground reaction forces.

A Point Kinetics Model for Estimating Neutron Multiplication of Bare Uranium Metal in Tagged Neutron Measurements

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

Tweardy, Matthew C.; McConchie, Seth; Hayward, Jason P.

An extension of the point kinetics model is developed in this paper to describe the neutron multiplicity response of a bare uranium object under interrogation by an associated particle imaging deuterium-tritium (D-T) measurement system. This extended model is used to estimate the total neutron multiplication of the uranium. Both MCNPX-PoliMi simulations and data from active interrogation measurements of highly enriched and depleted uranium geometries are used to evaluate the potential of this method and to identify the sources of systematic error. The detection efficiency correction for measured coincidence response is identified as a large source of systematic error. If themore » detection process is not considered, results suggest that the method can estimate total multiplication to within 13% of the simulated value. Values for multiplicity constants in the point kinetics equations are sensitive to enrichment due to (n, xn) interactions by D-T neutrons and can introduce another significant source of systematic bias. This can theoretically be corrected if isotopic composition is known a priori. Finally, the spatial dependence of multiplication is also suspected of introducing further systematic bias for high multiplication uranium objects.« less

A Point Kinetics Model for Estimating Neutron Multiplication of Bare Uranium Metal in Tagged Neutron Measurements

DOE PAGES

Tweardy, Matthew C.; McConchie, Seth; Hayward, Jason P.

2017-06-13

An extension of the point kinetics model is developed in this paper to describe the neutron multiplicity response of a bare uranium object under interrogation by an associated particle imaging deuterium-tritium (D-T) measurement system. This extended model is used to estimate the total neutron multiplication of the uranium. Both MCNPX-PoliMi simulations and data from active interrogation measurements of highly enriched and depleted uranium geometries are used to evaluate the potential of this method and to identify the sources of systematic error. The detection efficiency correction for measured coincidence response is identified as a large source of systematic error. If themore » detection process is not considered, results suggest that the method can estimate total multiplication to within 13% of the simulated value. Values for multiplicity constants in the point kinetics equations are sensitive to enrichment due to (n, xn) interactions by D-T neutrons and can introduce another significant source of systematic bias. This can theoretically be corrected if isotopic composition is known a priori. Finally, the spatial dependence of multiplication is also suspected of introducing further systematic bias for high multiplication uranium objects.« less

Evanescent field microscopy techniques for studying dynamics at the surface of living cells

NASA Astrophysics Data System (ADS)

Sund, Susan E.

This thesis presents two distinct optical microscopy techniques for applications in cell biophysics: (a)the extension to living cells of an established technique, total internal reflection/fluorescence recovery after photobleaching (TIR/FRAP) for the first time in imaging mode; and (b)the novel development of polarized total internal reflection fluorescence (p- TIRF) to study membrane orientation in living cells. Although reversible chemistry is crucial to dynamical processes in living cells, relatively little is known about the relevant chemical kinetic rates in vivo. TIR/FRAP, an established technique which can measure reversible biomolecular kinetic rates at surfaces, is extended here to measure kinetic parameters of microinjected rhodamine actin at the cytofacial surface of the plasma membrane of living cultured smooth muscle cells. For the first time, spatial imaging (with a CCD camera) is used in conjunction with TIR/FRAP. TIR/FRAP imaging allows production of spatially resolved images of kinetic data, and calculation of correlation distances, cell-wide gradients, and kinetic parameter dependence on initial fluorescence intensity. In living cells, membrane curvature occurs both in easily imaged large scale morphological features, and also in less visualizable submicroscopic regions of activity such as endocytosis, exocytosis, and cell surface ruffling. A fluorescence microscopic method, p-TIRF, is introduced here to visualize such regions. The method is based on fluorescence of the oriented membrane probe diI- C18-(3) (diI) excited by evanescent field light polarized either perpendicular or parallel to the plane of the substrate coverslip. The excitation efficiency from each polarization depends on the membrane orientation, and thus the ratio of the observed fluorescence excited by these two polarizations vividly shows regions of microscopic and submicroscopic curvature of the membrane. A theoretical background of the technique and experimental verifications are presented in samples of protein solutions, model lipid bilayers, and living cells. Sequential digital images of the polarized TIR fluorescence ratios show spatially-resolved time- course maps of membrane orientations on diI labeled macrophages from which low visibility membrane structures can be identified and quantified. The TIR images are sharpened and contrast-enhanced by deconvoluting them with an experimentally-measured point spread function.

Reply to 'Comments on upscaling geochemical reaction rates usingpore-scale network modeling' by Peter C. Lichtner and Qinjun Kang

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

Li, Li; Peters, Catherine A.; Celia, Michael A.

2006-05-03

Our paper "Upscaling geochemical reaction rates usingpore-scale network modeling" presents a novel application of pore-scalenetwork modeling to upscale mineral dissolution and precipitationreaction rates from the pore scale to the continuum scale, anddemonstrates the methodology by analyzing the scaling behavior ofanorthite and kaolinite reaction kinetics under conditions related to CO2sequestration. We conclude that under highly acidic conditions relevantto CO2 sequestration, the traditional continuum-based methodology may notcapture the spatial variation in concentrations from pore to pore, andscaling tools may be important in correctly modeling reactive transportprocesses in such systems. This work addresses the important butdifficult question of scaling mineral dissolution and precipitationreactionmore » kinetics, which is often ignored in fields such as geochemistry,water resources, and contaminant hydrology. Although scaling of physicalprocesses has been studied for almost three decades, very few studieshave examined the scaling issues related to chemical processes, despitetheir importance in governing the transport and fate of contaminants insubsurface systems.« less

Binomial tau-leap spatial stochastic simulation algorithm for applications in chemical kinetics.

PubMed

Marquez-Lago, Tatiana T; Burrage, Kevin

2007-09-14

In cell biology, cell signaling pathway problems are often tackled with deterministic temporal models, well mixed stochastic simulators, and/or hybrid methods. But, in fact, three dimensional stochastic spatial modeling of reactions happening inside the cell is needed in order to fully understand these cell signaling pathways. This is because noise effects, low molecular concentrations, and spatial heterogeneity can all affect the cellular dynamics. However, there are ways in which important effects can be accounted without going to the extent of using highly resolved spatial simulators (such as single-particle software), hence reducing the overall computation time significantly. We present a new coarse grained modified version of the next subvolume method that allows the user to consider both diffusion and reaction events in relatively long simulation time spans as compared with the original method and other commonly used fully stochastic computational methods. Benchmarking of the simulation algorithm was performed through comparison with the next subvolume method and well mixed models (MATLAB), as well as stochastic particle reaction and transport simulations (CHEMCELL, Sandia National Laboratories). Additionally, we construct a model based on a set of chemical reactions in the epidermal growth factor receptor pathway. For this particular application and a bistable chemical system example, we analyze and outline the advantages of our presented binomial tau-leap spatial stochastic simulation algorithm, in terms of efficiency and accuracy, in scenarios of both molecular homogeneity and heterogeneity.

Mathematical Modeling of Protein Misfolding Mechanisms in Neurological Diseases: A Historical Overview.

PubMed

Carbonell, Felix; Iturria-Medina, Yasser; Evans, Alan C

2018-01-01

Protein misfolding refers to a process where proteins become structurally abnormal and lose their specific 3-dimensional spatial configuration. The histopathological presence of misfolded protein (MP) aggregates has been associated as the primary evidence of multiple neurological diseases, including Prion diseases, Alzheimer's disease, Parkinson's disease, and Creutzfeldt-Jacob disease. However, the exact mechanisms of MP aggregation and propagation, as well as their impact in the long-term patient's clinical condition are still not well understood. With this aim, a variety of mathematical models has been proposed for a better insight into the kinetic rate laws that govern the microscopic processes of protein aggregation. Complementary, another class of large-scale models rely on modern molecular imaging techniques for describing the phenomenological effects of MP propagation over the whole brain. Unfortunately, those neuroimaging-based studies do not take full advantage of the tremendous capabilities offered by the chemical kinetics modeling approach. Actually, it has been barely acknowledged that the vast majority of large-scale models have foundations on previous mathematical approaches that describe the chemical kinetics of protein replication and propagation. The purpose of the current manuscript is to present a historical review about the development of mathematical models for describing both microscopic processes that occur during the MP aggregation and large-scale events that characterize the progression of neurodegenerative MP-mediated diseases.

Impulsive penetration : a viable mechanism for plasma entry across the magnetopause ?

NASA Astrophysics Data System (ADS)

De Keyser, Johan; Echim, Marius; Darrouzet, Fabien; Gunell, Herbert

Density inhomogeneities in the solar wind may cross the bow shock, and retain an excess earthward momentum in the magnetosheath upon approaching the magnetopause. Also, the bow shock dynamics as well as the behaviour of the magnetopause itself may introduce spatial inhomogeneities in the magnetosheath density and/or flow. Plasma entities with excess momentum may penetrate across the magnetopause, by the impulsive penetration mechanism. This plasma entry mechanism requires the existence of a polarization electric field in the moving blob, that is sustained by charge separation layers in the interfaces at the flanks of the blob. Both direct observation and simulation of plasma entry across the magnetopause following the impulsive penetration mechanism are hard. It is difficult to prove that observed plasma entry is really due to the impulsive penetration mechanism since the required charge separation layers or the resulting polarization electric field are hard to measure directly. Simply assessing the geometry is not easy, although multi-spacecraft missions like Cluster have resolved many of the ambiguities inherent in single-spacecraft measurements. Impulsive penetration is difficult to simulate as it operates on the fluid, the ion, and the electron scales simultaneously. It requires not only a high spatial resolution, but also a high precision to properly represent the charge imbalance in the flank interfaces. We have modelled impulsive penetration with a kinetic model, by simplifying the problem. The fully kinetic model is 3-dimensional in velocity space, but we consider spatial structure only along a single spatial dimension, namely the coordinate transverse to the blob’s direction of motion. We thereby assume that the blob is elongated both along the magnetic field and in the direction of motion. The model is semi-analytic and is able to represent the charge imbalance in the blob edges very well. In a second modelling step, we consider a slow, quasi-static change of this structure as the blob penetrates deeper into the magnetosphere, resulting in a description of the evolution of the penetrating plasma blob as a consequence of both adiabatic and non-adiabatic deceleration. Although the simulation considers this a simplified geometry, it sheds some light on some fundamental aspects of this plasma entry mechanism.

Stimulated luminescence emission from localized recombination in randomly distributed defects.

PubMed

Jain, Mayank; Guralnik, Benny; Andersen, Martin Thalbitzer

2012-09-26

We present a new kinetic model describing localized electronic recombination through the excited state of the donor (d) to an acceptor (a) centre in luminescent materials. In contrast to the existing models based on the localized transition model (LTM) of Halperin and Braner (1960 Phys. Rev. 117 408-15) which assumes a fixed d → a tunnelling probability for the entire crystal, our model is based on nearest-neighbour recombination within randomly distributed centres. Such a random distribution can occur through the entire volume or within the defect complexes of the dosimeter, and implies that the tunnelling probability varies with the donor-acceptor (d-a) separation distance. We first develop an 'exact kinetic model' that incorporates this variation in tunnelling probabilities, and evolves both in spatial as well as temporal domains. We then develop a simplified one-dimensional, semi-analytical model that evolves only in the temporal domain. An excellent agreement is observed between thermally and optically stimulated luminescence (TL and OSL) results produced from the two models. In comparison to the first-order kinetic behaviour of the LTM of Halperin and Braner (1960 Phys. Rev. 117 408-15), our model results in a highly asymmetric TL peak; this peak can be understood to derive from a continuum of several first-order TL peaks. Our model also shows an extended power law behaviour for OSL (or prompt luminescence), which is expected from localized recombination mechanisms in materials with random distribution of centres.

PHAST Version 2-A Program for Simulating Groundwater Flow, Solute Transport, and Multicomponent Geochemical Reactions

USGS Publications Warehouse

Parkhurst, David L.; Kipp, Kenneth L.; Charlton, Scott R.

2010-01-01

The computer program PHAST (PHREEQC And HST3D) simulates multicomponent, reactive solute transport in three-dimensional saturated groundwater flow systems. PHAST is a versatile groundwater flow and solute-transport simulator with capabilities to model a wide range of equilibrium and kinetic geochemical reactions. The flow and transport calculations are based on a modified version of HST3D that is restricted to constant fluid density and constant temperature. The geochemical reactions are simulated with the geochemical model PHREEQC, which is embedded in PHAST. Major enhancements in PHAST Version 2 allow spatial data to be defined in a combination of map and grid coordinate systems, independent of a specific model grid (without node-by-node input). At run time, aquifer properties are interpolated from the spatial data to the model grid; regridding requires only redefinition of the grid without modification of the spatial data. PHAST is applicable to the study of natural and contaminated groundwater systems at a variety of scales ranging from laboratory experiments to local and regional field scales. PHAST can be used in studies of migration of nutrients, inorganic and organic contaminants, and radionuclides; in projects such as aquifer storage and recovery or engineered remediation; and in investigations of the natural rock/water interactions in aquifers. PHAST is not appropriate for unsaturated-zone flow, multiphase flow, or density-dependent flow. A variety of boundary conditions are available in PHAST to simulate flow and transport, including specified-head, flux (specified-flux), and leaky (head-dependent) conditions, as well as the special cases of rivers, drains, and wells. Chemical reactions in PHAST include (1) homogeneous equilibria using an ion-association or Pitzer specific interaction thermodynamic model; (2) heterogeneous equilibria between the aqueous solution and minerals, ion exchange sites, surface complexation sites, solid solutions, and gases; and (3) kinetic reactions with rates that are a function of solution composition. The aqueous model (elements, chemical reactions, and equilibrium constants), minerals, exchangers, surfaces, gases, kinetic reactants, and rate expressions may be defined or modified by the user. A number of options are available to save results of simulations to output files. The data may be saved in three formats: a format suitable for viewing with a text editor; a format suitable for exporting to spreadsheets and postprocessing programs; and in Hierarchical Data Format (HDF), which is a compressed binary format. Data in the HDF file can be visualized on Windows computers with the program Model Viewer and extracted with the utility program PHASTHDF; both programs are distributed with PHAST.

Imaging Fourier transform spectroscopy of the boundary layer plume from laser irradiated polymers and carbon materials

NASA Astrophysics Data System (ADS)

Acosta, Roberto I.

The high-energy laser (HEL) lethality community needs for enhanced laser weapons systems requires a better understanding of a wide variety of emerging threats. In order to reduce the dimensionality of laser-materials interaction it is necessary to develop novel predictive capabilities of these events. The objective is to better understand the fundamentals of laser lethality testing by developing empirical models from hyperspectral imagery, enabling a robust library of experiments for vulnerability assessments. Emissive plumes from laser irradiated fiberglass reinforced polymers (FRP), poly(methyl methacrylate) (PMMA) and porous graphite targets were investigated primarily using a mid-wave infrared (MWIR) imaging Fourier transform spectrometer (FTS). Polymer and graphite targets were irradiated with a continuous wave (cw) fiber lasers. Data was acquired with a spectral resolution of 2 cm-1 and spatial resolution as high as 0.52 mm2 per pixel. Strong emission from H2O, CO, CO2 and hydrocarbons were observed in the MWIR between 1900-4000 cm-1. A single-layer radiative transfer model was developed to estimate spatial maps of temperature and column densities of CO and CO2 from the hyperspectral imagery of the boundary layer plume. The spectral model was used to compute the absorption cross sections of CO and CO2, using spectral line parameters from the high temperature extension of the HITRAN. Also, spatial maps of gas-phase temperature and methyl methacrylate (MMA) concentration were developed from laser irradiated carbon black-pigmented PMMA at irradiances of 4-22 W/cm2. Global kinetics interplay between heterogeneous and homogeneous combustion kinetics are shown from experimental observations at high spatial resolutions. Overall the boundary layer profile at steady-state is consistent with CO being mainly produced at the surface by heterogeneous reactions followed by a rapid homogeneous combustion in the boundary layer towards buoyancy.

Anomalous electron transport in Hall-effect thrusters: Comparison between quasi-linear kinetic theory and particle-in-cell simulations

NASA Astrophysics Data System (ADS)

Lafleur, T.; Martorelli, R.; Chabert, P.; Bourdon, A.

2018-06-01

Kinetic drift instabilities have been implicated as a possible mechanism leading to anomalous electron cross-field transport in E × B discharges, such as Hall-effect thrusters. Such instabilities, which are driven by the large disparity in electron and ion drift velocities, present a significant challenge to modelling efforts without resorting to time-consuming particle-in-cell (PIC) simulations. Here, we test aspects of quasi-linear kinetic theory with 2D PIC simulations with the aim of developing a self-consistent treatment of these instabilities. The specific quantities of interest are the instability growth rate (which determines the spatial and temporal evolution of the instability amplitude), and the instability-enhanced electron-ion friction force (which leads to "anomalous" electron transport). By using the self-consistently obtained electron distribution functions from the PIC simulations (which are in general non-Maxwellian), we find that the predictions of the quasi-linear kinetic theory are in good agreement with the simulation results. By contrast, the use of Maxwellian distributions leads to a growth rate and electron-ion friction force that is around 2-4 times higher, and consequently significantly overestimates the electron transport. A possible method for self-consistently modelling the distribution functions without requiring PIC simulations is discussed.

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

Perry, William L; Gunderson, Jake A; Dickson, Peter M

There has been a long history of interest in the decomposition kinetics of HMX and HMX-based formulations due to the widespread use of this explosive in high performance systems. The kinetics allow us to predict, or attempt to predict, the behavior of the explosive when subjected to thermal hazard scenarios that lead to ignition via impact, spark, friction or external heat. The latter, commonly referred to as 'cook off', has been widely studied and contemporary kinetic and transport models accurately predict time and location of ignition for simple geometries. However, there has been relatively little attention given to the problemmore » of localized ignition that results from the first three ignition sources of impact, spark and friction. The use of a zero-order single-rate expression describing the exothermic decomposition of explosives dates to the early work of Frank-Kamanetskii in the late 1930s and continued through the 60's and 70's. This expression provides very general qualitative insight, but cannot provide accurate spatial or timing details of slow cook off ignition. In the 70s, Catalano, et al., noted that single step kinetics would not accurately predict time to ignition in the one-dimensional time to explosion apparatus (ODTX). In the early 80s, Tarver and McGuire published their well-known three step kinetic expression that included an endothermic decomposition step. This scheme significantly improved the accuracy of ignition time prediction for the ODTX. However, the Tarver/McGuire model could not produce the internal temperature profiles observed in the small-scale radial experiments nor could it accurately predict the location of ignition. Those factors are suspected to significantly affect the post-ignition behavior and better models were needed. Brill, et al. noted that the enthalpy change due to the beta-delta crystal phase transition was similar to the assumed endothermic decomposition step in the Tarver/McGuire model. Henson, et al., deduced the kinetics and thermodynamics of the phase transition, providing Dickson, et al. with the information necessary to develop a four-step model that included a two-step nucleation and growth mechanism for the {beta}-{delta} phase transition. Initially, an irreversible scheme was proposed. That model accurately predicted the spatial and temporal cook off behavior of the small-scale radial experiment under slow heating conditions, but did not accurately capture the endothermic phase transition at a faster heating rate. The current version of the four-step model includes reversibility and accurately describes the small-scale radial experiment over a wide range of heating rates. We have observed impact-induced friction ignition of PBX 9501 with grit embedded between the explosive and the lower anvil surface. Observation was done using an infrared camera looking through the sapphire bottom anvil. Time to ignition and temperature-time behavior were recorded. The time to ignition was approximately 500 microseconds and the temperature was approximately 1000 K. The four step reversible kinetic scheme was previously validated for slow cook off scenarios. Our intention was to test the validity for significantly faster hot-spot processes, such as the impact-induced grit friction process studied here. We found the model predicted the ignition time within experimental error. There are caveats to consider when evaluating the agreement. The primary input to the model was friction work over an area computed by a stress analysis. The work rate itself, and the relative velocity of the grit and substrate both have a strong dependence on the initial position of the grit. Any errors in the analysis or the initial grit position would affect the model results. At this time, we do not know the sensitivity to these issues. However, the good agreement does suggest the four step kinetic scheme may have universal applicability for HMX systems.« less

Emergence of chaos in a spatially confined reactive system

NASA Astrophysics Data System (ADS)

Voorsluijs, Valérie; De Decker, Yannick

2016-11-01

In spatially restricted media, interactions between particles and local fluctuations of density can lead to important deviations of the dynamics from the unconfined, deterministic picture. In this context, we investigated how molecular crowding can affect the emergence of chaos in small reactive systems. We developed to this end an amended version of the Willamowski-Rössler model, where we account for the impenetrability of the reactive species. We analyzed the deterministic kinetics of this model and studied it with spatially-extended stochastic simulations in which the mobility of particles is included explicitly. We show that homogeneous fluctuations can lead to a destruction of chaos through a fluctuation-induced collision between chaotic trajectories and absorbing states. However, an interplay between the size of the system and the mobility of particles can counterbalance this effect so that chaos can indeed be found when particles diffuse slowly. This unexpected effect can be traced back to the emergence of spatial correlations which strongly affect the dynamics. The mobility of particles effectively acts as a new bifurcation parameter, enabling the system to switch from stationary states to absorbing states, oscillations or chaos.

In silico study on the effects of matrix structure in controlled drug release

NASA Astrophysics Data System (ADS)

Villalobos, Rafael; Cordero, Salomón; Maria Vidales, Ana; Domínguez, Armando

2006-07-01

Purpose: To study the effects of drug concentration and spatial distribution of the medicament, in porous solid dosage forms, on the kinetics and total yield of drug release. Methods: Cubic networks are used as models of drug release systems. They were constructed by means of the dual site-bond model framework, which allows a substrate to have adequate geometrical and topological distribution of its pore elements. Drug particles can move inside the networks by following a random walk model with excluded volume interactions between the particles. The drug release time evolution for different drug concentration and different initial drug spatial distribution has been monitored. Results: The numerical results show that in all the studied cases, drug release presents an anomalous behavior, and the consequences of the matrix structural properties, i.e., drug spatial distribution and drug concentration, on the drug release profile have been quantified. Conclusions: The Weibull function provides a simple connection between the model parameters and the microstructure of the drug release device. A critical modeling of drug release from matrix-type delivery systems is important in order to understand the transport mechanisms that are implicated, and to predict the effect of the device design parameters on the release rate.

Kinetics of the Tau PET Tracer 18F-AV-1451 (T807) in Subjects with Normal Cognitive Function, Mild Cognitive Impairment, and Alzheimer Disease.

PubMed

Shcherbinin, Sergey; Schwarz, Adam J; Joshi, Abhinay; Navitsky, Michael; Flitter, Matthew; Shankle, William R; Devous, Michael D; Mintun, Mark A

2016-10-01

We report kinetic modeling results of dynamic acquisition data from 0 to 100 min after injection with the tau PET tracer 18 F-AV-1451 in 19 subjects. Subjects were clinically diagnosed as 4 young cognitively normal, 5 old cognitively normal, 5 mild cognitive impairment, and 5 Alzheimer disease (AD). Kinetic modeling was performed using Logan graphical analysis with the cerebellum crus as a reference region. Voxelwise binding potential ([Formula: see text]) and SUV ratio ([Formula: see text]) images were compared. In AD subjects, slower and spatially nonuniform clearance from cortical regions was observed as compared with the controls, which led to focal uptake and elevated retention in the imaging data from 80 to 100 min after injection. BP from the dynamic data from 0 to 100 min correlated strongly (R 2 > 0.86) with corresponding regional [Formula: see text] values. In the putamen, the observed kinetics (positive [Formula: see text] at the tracer delivery stage and plateauing time-SUVR curves for all diagnostic categories) may suggest either additional off-target binding or a second binding site with different kinetics. The kinetics of the 18 F-AV-1451 tracer in cortical areas, as examined in this small group of subjects, differed by diagnostic stage. A delayed 80- to 100-min scan provided a reasonable substitute for a dynamic 0- to 100-min acquisition for cortical regions although other windows (e.g., 75-105 min) may be useful to evaluate. © 2016 by the Society of Nuclear Medicine and Molecular Imaging, Inc.

Spatial Frequency Responses of Anisotropic Refractive Index Gratings Formed in Holographic Polymer Dispersed Liquid Crystals

PubMed Central

Fukuda, Yoshiaki; Tomita, Yasuo

2016-01-01

We report on an experimental investigation of spatial frequency responses of anisotropic transmission refractive index gratings formed in holographic polymer dispersed liquid crystals (HPDLCs). We studied two different types of HPDLC materials employing two different monomer systems: one with acrylate monomer capable of radical mediated chain-growth polymerizations and the other with thiol-ene monomer capable of step-growth polymerizations. It was found that the photopolymerization kinetics of the two HPDLC materials could be well explained by the autocatalytic model. We also measured grating-spacing dependences of anisotropic refractive index gratings at a recording wavelength of 532 nm. It was found that the HPDLC material with the thiol-ene monomer gave higher spatial frequency responses than that with the acrylate monomer. Statistical thermodynamic simulation suggested that such a spatial frequency dependence was attributed primarily to a difference in the size of formed liquid crystal droplets due to different photopolymerization mechanisms. PMID:28773314

Spatial Frequency Responses of Anisotropic Refractive Index Gratings Formed in Holographic Polymer Dispersed Liquid Crystals.

PubMed

Fukuda, Yoshiaki; Tomita, Yasuo

2016-03-10

We report on an experimental investigation of spatial frequency responses of anisotropic transmission refractive index gratings formed in holographic polymer dispersed liquid crystals (HPDLCs). We studied two different types of HPDLC materials employing two different monomer systems: one with acrylate monomer capable of radical mediated chain-growth polymerizations and the other with thiol-ene monomer capable of step-growth polymerizations. It was found that the photopolymerization kinetics of the two HPDLC materials could be well explained by the autocatalytic model. We also measured grating-spacing dependences of anisotropic refractive index gratings at a recording wavelength of 532 nm. It was found that the HPDLC material with the thiol-ene monomer gave higher spatial frequency responses than that with the acrylate monomer. Statistical thermodynamic simulation suggested that such a spatial frequency dependence was attributed primarily to a difference in the size of formed liquid crystal droplets due to different photopolymerization mechanisms.

Birefringent breakup of Dirac fermions on a square optical lattice

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

Kennett, Malcolm P.; Komeilizadeh, Nazanin; Kaveh, Kamran

2011-05-15

We introduce a lattice model for fermions in a spatially periodic magnetic field that also has spatially periodic hopping amplitudes. We discuss how this model might be realized with cold atoms in an artificial magnetic field on a square optical lattice. When there is an average flux of half a flux quantum per plaquette, the spectrum of low-energy excitations can be described by massless Dirac fermions in which the usually doubly degenerate Dirac cones split into cones with different ''speeds of light.'' These gapless birefringent Dirac fermions arise because of broken chiral symmetry in the kinetic energy term of themore » effective low-energy Hamiltonian. We characterize the effects of various perturbations to the low-energy spectrum, including staggered potentials, interactions, and domain-wall topological defects.« less

Connection between two statistical approaches for the modelling of particle velocity and concentration distributions in turbulent flow: The mesoscopic Eulerian formalism and the two-point probability density function method

NASA Astrophysics Data System (ADS)

Simonin, Olivier; Zaichik, Leonid I.; Alipchenkov, Vladimir M.; Février, Pierre

2006-12-01

The objective of the paper is to elucidate a connection between two approaches that have been separately proposed for modelling the statistical spatial properties of inertial particles in turbulent fluid flows. One of the approaches proposed recently by Février, Simonin, and Squires [J. Fluid Mech. 533, 1 (2005)] is based on the partitioning of particle turbulent velocity field into spatially correlated (mesoscopic Eulerian) and random-uncorrelated (quasi-Brownian) components. The other approach stems from a kinetic equation for the two-point probability density function of the velocity distributions of two particles [Zaichik and Alipchenkov, Phys. Fluids 15, 1776 (2003)]. Comparisons between these approaches are performed for isotropic homogeneous turbulence and demonstrate encouraging agreement.

The effective compliance of spatially evolving planar wing-cracks

NASA Astrophysics Data System (ADS)

Ayyagari, R. S.; Daphalapurkar, N. P.; Ramesh, K. T.

2018-02-01

We present an analytic closed form solution for anisotropic change in compliance due to the spatial evolution of planar wing-cracks in a material subjected to largely compressive loading. A fully three-dimensional anisotropic compliance tensor is defined and evaluated considering the wing-crack mechanism, using a mixed-approach based on kinematic and energetic arguments to derive the coefficients in incremental compliance. Material, kinematic and kinetic parametric influences on the increments in compliance are studied in order to understand their physical implications on material failure. Model verification is carried out through comparisons to experimental uniaxial compression results to showcase the predictive capabilities of the current study.

Computational Study of Thrombus Formation and Clotting Factor Effects under Venous Flow Conditions

PubMed Central

Govindarajan, Vijay; Rakesh, Vineet; Reifman, Jaques; Mitrophanov, Alexander Y.

2016-01-01

A comprehensive understanding of thrombus formation as a physicochemical process that has evolved to protect the integrity of the human vasculature is critical to our ability to predict and control pathological states caused by a malfunctioning blood coagulation system. Despite numerous investigations, the spatial and temporal details of thrombus growth as a multicomponent process are not fully understood. Here, we used computational modeling to investigate the temporal changes in the spatial distributions of the key enzymatic (i.e., thrombin) and structural (i.e., platelets and fibrin) components within a growing thrombus. Moreover, we investigated the interplay between clot structure and its mechanical properties, such as hydraulic resistance to flow. Our model relied on the coupling of computational fluid dynamics and biochemical kinetics, and was validated using flow-chamber data from a previous experimental study. The model allowed us to identify the distinct patterns characterizing the spatial distributions of thrombin, platelets, and fibrin accumulating within a thrombus. Our modeling results suggested that under the simulated conditions, thrombin kinetics was determined predominantly by prothrombinase. Furthermore, our simulations showed that thrombus resistance imparted by fibrin was ∼30-fold higher than that imparted by platelets. Yet, thrombus-mediated bloodflow occlusion was driven primarily by the platelet deposition process, because the height of the platelet accumulation domain was approximately twice that of the fibrin accumulation domain. Fibrinogen supplementation in normal blood resulted in a nonlinear increase in thrombus resistance, and for a supplemented fibrinogen level of 48%, the thrombus resistance increased by ∼2.7-fold. Finally, our model predicted that restoring the normal levels of clotting factors II, IX, and X while simultaneously restoring fibrinogen (to 88% of its normal level) in diluted blood can restore fibrin generation to ∼78% of its normal level and hence improve clot formation under dilution. PMID:27119646

A passive physical model for DnaK chaperoning

NASA Astrophysics Data System (ADS)

Uhl, Lionel; Dumont, Audrey; Dukan, Sam

2018-03-01

Almost all living organisms use protein chaperones with a view to preventing proteins from misfolding or aggregation either spontaneously or during cellular stress. This work uses a reaction-diffusion stochastic model to describe the dynamic localization of the Hsp70 chaperone DnaK in Escherichia coli cells during transient proteotoxic collapse characterized by the accumulation of insoluble proteins. In the model, misfolded (‘abnormal’) proteins are produced during alcoholic stress and have the propensity to aggregate with a polymerization-like kinetics. When aggregates diffuse more slowly they grow larger. According to Michaelis-Menten-type kinetics, DnaK has the propensity to bind with misfolded proteins or aggregates in order to catalyse refolding. To match experimental fluorescence microscopy data showing clusters of DnaK-GFP localized in multiple foci, the model includes spatial zones with local reduced diffusion rates to generate spontaneous assemblies of DnaK called ‘foci’. Numerical simulations of our model succeed in reproducing the kinetics of DnaK localization experimentally observed. DnaK starts from foci, moves to large aggregates during acute stress, resolves those aggregates during recovery and finally returns to its initial punctate localization pattern. Finally, we compare real biological events with hypothetical repartitions of the protein aggregates or DnaK. We then notice that DnaK action is more efficient on protein aggregates than on protein homogeneously distributed.

Experimental and Computational Study fo CH, CH*, and OH* in an Axisymmetric Laminar Diffusion Flame

NASA Technical Reports Server (NTRS)

Walsh, K. T.

1998-01-01

In this study, we extend the results of previous combined numerical and experimental investigations of an axisymmetric laminar diffusion flame in which difference Raman spectroscopy, laser-induced fluorescence (LIF), and a multidimensional flame model were used to generate profiles of the temperature and major and minor species. A procedure is outlined by which the number densities of ground-state CH (X(sup 2)II) excited-state CH (A(sup 2)Delta, denoted CH*), and excited-state OH (A(sup 2)Sigma, denoted OH*) are measured and modeled. CH* and OH* number densities are deconvoluted from line-of-sight flame-emission measurements. Ground-state CH is measured using linear LIF. The computations are done with GRI Mech 2.11 as well as an alternate hydrocarbon mechanism. In both cases, additional reactions for the production and consumption of CH* and OH* are added from recent kinetic studies. Collisional quenching and spontaneous emission are responsible for the de-excitation of the excited-state radicals. As with our previous investigations, GRI Mech 2.11 continues to produce very good agreement with the overall flame length observed in the experiments, while significantly under predicting the flame lift-off height. The alternate kinetic scheme is much more accurate in predicting lift-off height but overpredicts the over-all flame length. Ground-state CH profiles predicted with GRI Mech 2.11 are in excellent agreement with the corresponding measurements, regarding both spatial distribution and absolute concentration (measured at 4 ppm) of the CH radical. Calculations of the excited-state species show reasonable agreement with the measurements as far as spatial distribution and overall characteristics are concerned. For OH*, the measured peak mole fraction, 1.3 x 10(exp -8), compared well with computed peaks, while the measured peak level for CH*, 2 x 10(exp -9), was severely underpredicted by both kinetic schemes, indicating that the formation and destruction kinetics associated with excited-state species in flames require further research.

Implementation of Solute Transport in the Vadose Zone into the `HYDRUS Package for MODFLOW'

NASA Astrophysics Data System (ADS)

Simunek, J.; Beegum, S.; Szymkiewicz, A.; Sudheer, K. P.

2017-12-01

The 'HYDRUS package for MODFLOW' was developed by Seo et al. (2007) and Twarakavi et al. (2008) to simultaneously evaluate transient water flow in both unsaturated and saturated zones. The package, which is based on the HYDRUS-1D model (Šimůnek et al., 2016) simulating unsaturated water flow in the vadose zone, was incorporated into MODFLOW (Harbaugh et al., 2000) simulating saturated groundwater flow. The HYDRUS package in the coupled model can be used to represent the effects of various unsaturated zone processes, including infiltration, evaporation, root water uptake, capillary rise, and recharge in homogeneous or layered soil profiles. The coupled model is effective in addressing spatially-variable saturated-unsaturated hydrological processes at the regional scale, allowing for complex layering in the unsaturated zone, spatially and temporarily variable water fluxes at the soil surface and in the root zone, and with alternating recharge and discharge fluxes (Twarakavi et al., 2008). One of the major limitations of the coupled model was that it could not be used to simulate at the same time solute transport. However, solute transport is highly dependent on water table fluctuations due to temporal and spatial variations in groundwater recharge. This is an important concern when the coupled model is used for analyzing groundwater contamination due to transport through the unsaturated zone. The objective of this study is to integrate the solute transport model (the solute transport part of HYDRUS-1D for the unsaturated zone and MT3DMS (Zheng and Wang, 1999; Zheng, 2009) for the saturated zone) into an existing coupled water flow model. The unsaturated zone component of the coupled model can consider solute transport involving many biogeochemical processes and reactions, including first-order degradation, volatilization, linear or nonlinear sorption, one-site kinetic sorption, two-site sorption, and two-kinetic sites sorption (Šimůnek and van Genuchten, 2008). Due to complex interactions at the groundwater table, certain modifications of the pressure head (compared to the original coupling) and solute concentration profiles were incorporated into the HYDRUS package. The developed integrated model is verified using HYDRUS-2D and analyzed for its computational time requirements.

Kinetics of particle deposition at heterogeneous surfaces

NASA Astrophysics Data System (ADS)

Stojiljković, D. Lj.; Vrhovac, S. B.

2017-12-01

The random sequential adsorption (RSA) approach is used to analyze adsorption of spherical particles of fixed diameter d0 on nonuniform surfaces covered by square cells arranged in a square lattice pattern. To characterize such pattern two dimensionless parameters are used: the cell size α and the cell-cell separation β, measured in terms of the particle diameter d0. Adsorption is assumed to occur if the particle (projected) center lies within a cell area. We focus on the kinetics of deposition process in the case when no more than a single disk can be placed onto any square cell (α < 1 /√{ 2 } ≈ 0 . 707). We find that the asymptotic approach of the coverage fraction θ(t) to the jamming limit θJ is algebraic if the parameters α and β satisfy the simple condition, β + α / 2 < 1. If this condition is not satisfied, the late time kinetics of deposition process is not consistent with the power law behavior. However, if the geometry of the pattern approaches towards ;noninteracting conditions; (β > 1), when adsorption on each cell can be decoupled, approach of the coverage fraction θ(t) to θJ becomes closer to the exponential law. Consequently, changing the pattern parameters in the present model allows to interpolate the deposition kinetics between the continuum limit and the lattice-like behavior. Structural properties of the jammed-state coverings are studied in terms of the radial distribution function g(r) and spatial distribution of particles inside the cell. Various, non-trivial spatial distributions are observed depending on the geometry of the pattern.

Longitudinal patterns and response lengths of algae in riverine ecosystems: A model analysis emphasising benthic-pelagic interactions.

PubMed

Jäger, Christoph G; Borchardt, Dietrich

2018-04-07

In riverine ecosystems primary production is principally possible in two habitats: in the benthic layer by sessile algae and in the surface water by planktonic algae being transported downstream. The relevance of these two habitats generally changes along the rivers' continuum. However, analyses of the interaction of algae in these two habitats and their controlling factors in riverine ecosystems are, so far, very rare. We use a simplified advection-diffusion model system combined with ecological process kinetics to analyse the interaction of benthic and planktonic algae and nutrients along idealised streams and rivers at regional to large scales. Because many of the underlying processes affecting algal dynamics are influenced by depth, we focus particularly on the impact of river depth on this interaction. At constant environmental conditions all state variables approach stable spatial equilibria along the river, independent of the boundary conditions at the upstream end. Because our model is very robust against changes of turbulent diffusion and stream velocity, these spatial equilibria can be analysed by a simplified ordinary differential equation (ode) version of our model. This model variant reveals that at shallower river depths, phytoplankton can exist only when it is subsidised by detaching benthic algae, and in turn, at deeper river depths, benthic algae can exist only in low biomasses which are subsidised by sinking planktonic algae. We generalise the spatial dynamics of the model system using different conditions at the upstream end of the model, which mimic various natural or anthropogenic factors (pristine source, dam, inflow of a waste water treatment plant, and dilution from e.g. a tributary) and analyse how these scenarios influence different aspects of the longitudinal spatial dynamics of the full spatial model: the relation of spatial equilibrium to spatial maximum, the distance to the spatial maximum, and the response length. Generally, our results imply that shallow systems recover within significantly shorter distances from spatially distinct disturbances when compared to deep systems, independent of the type of disturbance. Copyright © 2018 The Authors. Published by Elsevier Ltd.. All rights reserved.

ATLASGAL-selected massive clumps in the inner Galaxy. VI. Kinetic temperature and spatial density measured with formaldehyde

NASA Astrophysics Data System (ADS)

Tang, X. D.; Henkel, C.; Wyrowski, F.; Giannetti, A.; Menten, K. M.; Csengeri, T.; Leurini, S.; Urquhart, J. S.; König, C.; Güsten, R.; Lin, Y. X.; Zheng, X. W.; Esimbek, J.; Zhou, J. J.

2018-03-01

Context. Formaldehyde (H2CO) is a reliable tracer to accurately measure the physical parameters of dense gas in star-forming regions. Aim. We aim to determine directly the kinetic temperature and spatial density with formaldehyde for the 100 brightest ATLASGAL-selected clumps (the TOP100 sample) at 870 μm representing various evolutionary stages of high-mass star formation. Methods: Ten transitions (J = 3-2 and 4-3) of ortho- and para-H2CO near 211, 218, 225, and 291 GHz were observed with the Atacama Pathfinder EXperiment (APEX) 12 m telescope. Results: Using non-LTE models with RADEX, we derived the gas kinetic temperature and spatial density with the measured para-H2CO 321-220/303-202, 422-321/404-303, and 404-303/303-202 ratios. The gas kinetic temperatures derived from the para-H2CO 321-220/303-202 and 422-321/404-303 line ratios are high, ranging from 43 to >300 K with an unweighted average of 91 ± 4 K. Deduced Tkin values from the J = 3-2 and 4-3 transitions are similar. Spatial densities of the gas derived from the para-H2CO 404-303/303-202 line ratios yield 0.6-8.3 × 106 cm-3 with an unweighted average of 1.5 (±0.1) × 106 cm-3. A comparison of kinetic temperatures derived from para-H2CO, NH3, and dust emission indicates that para-H2CO traces a distinctly higher temperature than the NH3 (2, 2)/(1, 1) transitions and the dust, tracing heated gas more directly associated with the star formation process. The H2CO line widths are found to be correlated with bolometric luminosity and increase with the evolutionary stage of the clumps, which suggests that higher luminosities tend to be associated with a more turbulent molecular medium. It seems that the spatial densities measured with H2CO do not vary significantly with the evolutionary stage of the clumps. However, averaged gas kinetic temperatures derived from H2CO increase with time through the evolution of the clumps. The high temperature of the gas traced by H2CO may be mainly caused by radiation from embedded young massive stars and the interaction of outflows with the ambient medium. For Lbol/Mclump ≳ 10 L⊙/M⊙, we find a rough correlation between gas kinetic temperature and this ratio, which is indicative of the evolutionary stage of the individual clumps. The strong relationship between H2CO line luminosities and clump masses is apparently linear during the late evolutionary stages of the clumps, indicating that LH_2CO does reliably trace the mass of warm dense molecular gas. In our massive clumps H2CO line luminosities are approximately linearly correlated with bolometric luminosities over about four orders of magnitude in Lbol, which suggests that the mass of dense molecular gas traced by the H2CO line luminosity is well correlated with star formation. Source and H2CO parameters (Tables A.1-A.7) are only available at the CDS via anonymous ftp to http://cdsarc.u-strasbg.fr (http://130.79.128.5) or via http://cdsarc.u-strasbg.fr/viz-bin/qcat?J/A+A/611/A6

Resolving runaway electron distributions in space, time, and energy

NASA Astrophysics Data System (ADS)

Paz-Soldan, C.; Cooper, C. M.; Aleynikov, P.; Eidietis, N. W.; Lvovskiy, A.; Pace, D. C.; Brennan, D. P.; Hollmann, E. M.; Liu, C.; Moyer, R. A.; Shiraki, D.

2018-05-01

Areas of agreement and disagreement with present-day models of runaway electron (RE) evolution are revealed by measuring MeV-level bremsstrahlung radiation from runaway electrons (REs) with a pinhole camera. Spatially resolved measurements localize the RE beam, reveal energy-dependent RE transport, and can be used to perform full two-dimensional (energy and pitch-angle) inversions of the RE phase-space distribution. Energy-resolved measurements find qualitative agreement with modeling on the role of collisional and synchrotron damping in modifying the RE distribution shape. Measurements are consistent with predictions of phase-space attractors that accumulate REs, with non-monotonic features observed in the distribution. Temporally resolved measurements find qualitative agreement with modeling on the impact of collisional and synchrotron damping in varying the RE growth and decay rate. Anomalous RE loss is observed and found to be largest at low energy. Possible roles for kinetic instability or spatial transport to resolve these anomalies are discussed.

The influence of spatially and temporally high-resolution wind forcing on the power input to near-inertial waves in the ocean

NASA Astrophysics Data System (ADS)

Rimac, Antonija; von Storch, Jin-Song; Eden, Carsten

2013-04-01

The estimated power required to sustain global general circulation in the ocean is about 2 TW. This power is supplied with wind stress and tides. Energy spectrum shows pronounced maxima at near-inertial frequency. Near-inertial waves excited by high-frequency winds represent an important source for deep ocean mixing since they can propagate into the deep ocean and dissipate far away from the generation sites. The energy input by winds to near-inertial waves has been studied mostly using slab ocean models and wind stress forcing with coarse temporal resolution (e.g. 6-hourly). Slab ocean models lack the ability to reproduce fundamental aspects of kinetic energy balance and systematically overestimate the wind work. Also, slab ocean models do not account the energy used for the mixed layer deepening or the energy radiating downward into the deep ocean. Coarse temporal resolution of the wind forcing strongly underestimates the near-inertial energy. To overcome this difficulty we use an eddy permitting ocean model with high-frequency wind forcing. We establish the following model setup: We use the Max Planck Institute Ocean Model (MPIOM) on a tripolar grid with 45 km horizontal resolution and 40 vertical levels. We run the model with wind forcings that vary in horizontal and temporal resolution. We use high-resolution (1-hourly with 35 km horizontal resolution) and low-resolution winds (6-hourly with 250 km horizontal resolution). We address the following questions: Is the kinetic energy of near-inertial waves enhanced when high-resolution wind forcings are used? If so, is this due to higher level of overall wind variability or higher spatial or temporal resolution of wind forcing? How large is the power of near-inertial waves generated by winds? Our results show that near-inertial waves are enhanced and the near-inertial kinetic energy is two times higher (in the storm track regions 3.5 times higher) when high-resolution winds are used. Filtering high-resolution winds in space and time, the near-inertial kinetic energy reduces. The reduction is faster when a temporal filter is used suggesting that the high-frequency wind forcing is more efficient in generating near-inertial wave energy than the small-scale wind forcing. Using low-resolution wind forcing the wind generated power to near-inertial waves is 0.55 TW. When we use high-resolution wind forcing the result is 1.6 TW meaning that the result increases by 300%.

Solute segregation kinetics and dislocation depinning in a binary alloy

NASA Astrophysics Data System (ADS)

Dontsova, E.; Rottler, J.; Sinclair, C. W.

2015-06-01

Static strain aging, a phenomenon caused by diffusion of solute atoms to dislocations, is an important contributor to the strength of substitutional alloys. Accurate modeling of this complex process requires both atomic spatial resolution and diffusional time scales, which is very challenging to achieve with commonly used atomistic computational methods. In this paper, we use the recently developed "diffusive molecular dynamics" (DMD) method that is capable of describing the kinetics of the solute segregation process at the atomic level while operating on diffusive time scales in a computationally efficient way. We study static strain aging in the Al-Mg system and calculate the depinning shear stress between edge and screw dislocations and their solute atmospheres formed for various waiting times with different solute content and for a range of temperatures. A simple phenomenological model is also proposed that describes the observed behavior of the critical shear stress as a function of segregation level.

Consistent kinetic simulation of plasma and sputtering in low temperature plasmas

NASA Astrophysics Data System (ADS)

Schmidt, Frederik; Trieschmann, Jan; Mussenbrock, Thomas

2016-09-01

Plasmas are commonly used in sputtering applications for the deposition of thin films. Although magnetron sources are a prominent choice, capacitively coupled plasmas have certain advantages (e.g., sputtering of non-conducting and/or ferromagnetic materials, aside of excellent control of the ion energy distribution). In order to understand the collective plasma and sputtering dynamics, a kinetic simulation model is helpful. Particle-in-Cell has been proven to be successful in simulating the plasma dynamics, while the Test-Multi-Particle-Method can be used to describe the sputtered neutral species. In this talk a consistent combination of these methods is presented by consistently coupling the simulated ion flux as input to a neutral particle transport model. The combined model is used to simulate and discuss the spatially dependent densities, fluxes and velocity distributions of all particles. This work is supported by the German Research Foundation (DFG) in the frame of Transregional Collaborative Research Center (SFB) TR-87.

Mixing-controlled reactive transport on travel times in heterogeneous media

NASA Astrophysics Data System (ADS)

Luo, J.; Cirpka, O.

2008-05-01

Modeling mixing-controlled reactive transport using traditional spatial discretization of the domain requires identifying the spatial distributions of hydraulic and reactive parameters including mixing-related quantities such as dispersivities and kinetic mass-transfer coefficients. In most applications, breakthrough curves of conservative and reactive compounds are measured at only a few locations and models are calibrated by matching these breakthrough curves, which is an ill posed inverse problem. By contrast, travel-time based transport models avoid costly aquifer characterization. By considering breakthrough curves measured on different scales, one can distinguish between mixing, which is a prerequisite for reactions, and spreading, which per se does not foster reactions. In the travel-time based framework, the breakthrough curve of a solute crossing an observation plane, or ending in a well, is interpreted as the weighted average of concentrations in an ensemble of non-interacting streamtubes, each of which is characterized by a distinct travel-time value. Mixing is described by longitudinal dispersion and/or kinetic mass transfer along individual streamtubes, whereas spreading is characterized by the distribution of travel times which also determines the weights associated to each stream tube. Key issues in using the travel-time based framework include the description of mixing mechanisms and the estimation of the travel-time distribution. In this work, we account for both apparent longitudinal dispersion and kinetic mass transfer as mixing mechanisms, thus generalizing the stochastic-convective model with or without inter-phase mass transfer and the advective-dispersive streamtube model. We present a nonparametric approach of determining the travel-time distribution, given a breakthrough curve integrated over an observation plane and estimated mixing parameters. The latter approach is superior to fitting parametric models in cases where the true travel-time distribution exhibits multiple peaks or long tails. It is demonstrated that there is freedom for the combinations of mixing parameters and travel-time distributions to fit conservative breakthrough curves and describe the tailing. Reactive transport cases with a bimolecular instantaneous irreversible reaction and a dual Michaelis-Menten problem demonstrate that the mixing introduced by local dispersion and mass transfer may be described by apparent mean mass transfer with coefficients evaluated by local breakthrough curves.

The mosaic structure of plasma bulk flows in the Earth's magnetotail

NASA Technical Reports Server (NTRS)

Ashour-Abdalla, M.; Richard, R. L.; Zelenyi, L. M.; Peroomian, V.; Bosqued, J. M.

1995-01-01

Moments of plasma distributions observed in the magnetotail vary with different time scales. In this paper we attempt to explain the observed variability on intermediate timescales of approximately 10-20 min that result from the simultaneous energization and spatial structuring of solar wind plasma in the distant magnetotail. These processes stimulate the formation of a system of spatially disjointed. highly accelerated filaments (beamlets) in the tail. We use the results from large-scale kinetic modeling of magnetotail formation from a plasma mantle source to calculate moments of ion distribution functions throughout the tail. Statistical restrictions related to the limited number of particles in our system naturally reduce the spatial resolution of our results, but we show that our model is valid on intermediate spatial scales Delta(x) x Delta(z) equal to approximately 1 R(sub E) x 1000 km. For these spatial scales the resulting pattern, which resembles a mosaic, appears to be quite variable. The complexity of the pattern is related to the spatial interference between beamlets accelerated at various locations within the distant tail which mirror in the strong near-Earth magnetic field. Global motion of the magnetotail results in the displacement of spacecraft with respect to this mosaic pattern and can produce variations in all of the moments (especially the x-component of the bulk velocity) on intermediate timescales. The results obtained enable us to view the magnetotail plasma as consisting of two different populations: a tailward-Earthward system of highly accelerated beamlets interfering with each other, and an energized quasithermal population which gradually builds as the Earth is approached. In the near-Earth tail, these populations merge into a hot quasi-isotropic ion population typical of the near-Earth plasma sheet. The transformation of plasma sheet boundary layer (PSBL) beam energy into central plasma sheet (CPS) quasi-thermal energy occurs in the absence of collisions or noise. This paper also clarifies the relationship between the global scale where an MHD description might be appropriate and the lower intermediate scales where MHD fails and large-scale kinetic theory should be used.

Goal-oriented sensitivity analysis for lattice kinetic Monte Carlo simulations

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

Arampatzis, Georgios, E-mail: garab@math.uoc.gr; Department of Mathematics and Statistics, University of Massachusetts, Amherst, Massachusetts 01003; Katsoulakis, Markos A., E-mail: markos@math.umass.edu

2014-03-28

In this paper we propose a new class of coupling methods for the sensitivity analysis of high dimensional stochastic systems and in particular for lattice Kinetic Monte Carlo (KMC). Sensitivity analysis for stochastic systems is typically based on approximating continuous derivatives with respect to model parameters by the mean value of samples from a finite difference scheme. Instead of using independent samples the proposed algorithm reduces the variance of the estimator by developing a strongly correlated-“coupled”- stochastic process for both the perturbed and unperturbed stochastic processes, defined in a common state space. The novelty of our construction is that themore » new coupled process depends on the targeted observables, e.g., coverage, Hamiltonian, spatial correlations, surface roughness, etc., hence we refer to the proposed method as goal-oriented sensitivity analysis. In particular, the rates of the coupled Continuous Time Markov Chain are obtained as solutions to a goal-oriented optimization problem, depending on the observable of interest, by considering the minimization functional of the corresponding variance. We show that this functional can be used as a diagnostic tool for the design and evaluation of different classes of couplings. Furthermore, the resulting KMC sensitivity algorithm has an easy implementation that is based on the Bortz–Kalos–Lebowitz algorithm's philosophy, where events are divided in classes depending on level sets of the observable of interest. Finally, we demonstrate in several examples including adsorption, desorption, and diffusion Kinetic Monte Carlo that for the same confidence interval and observable, the proposed goal-oriented algorithm can be two orders of magnitude faster than existing coupling algorithms for spatial KMC such as the Common Random Number approach. We also provide a complete implementation of the proposed sensitivity analysis algorithms, including various spatial KMC examples, in a supplementary MATLAB source code.« less

Numerical analysis of flows of rarefied gases in long channels with octagonal cross section shapes

NASA Astrophysics Data System (ADS)

Szalmas, L.

2014-12-01

Isothermal, pressure driven rarefied gas flows through long channels with octagonal cross section shapes are analyzed computationally. The capillary is between inlet and outlet reservoirs. The cross section is constant along the axial direction. The boundary condition at the solid-gas interface is assumed to be diffuse reflection. Since the channel is long, the gaseous velocity is small compared to the average molecular speed. Consequently, a linearized description can be used. The flow is described by the linearized Bhatnagar-Gross-Krook kinetic model. The solution of the problem is divided into two stages. First, the local flow field is determined by assuming the local pressure gradient. Secondly, the global flow behavior is deduced by the consideration of the conservation of the mass along the axis of the capillary. The kinetic equation is solved by the discrete velocity method on the cross section. Both spatial and velocity spaces are discretized. A body fitted rectangular grid is used for the spatial space. Near the boundary, first-order, while in the interior part of the flow domain, second-order finite-differences are applied to approximate the spatial derivatives. This combination results into an efficient and straightforward numerical treatment. The velocity space is represented by a Gauss-Legendre quadrature. The kinetic equation is solved in an iterative manner. The local dimensionless flow rate is calculated and tabulated for a wide range of the gaseous rarefaction for octagonal cross sections with various geometrical parameters. It exhibits the Knudsen minimum phenomenon. The flow rates in the octagonal channel are compared to those through capillaries with circular and square cross sections. Typical velocity profiles are also shown. The mass flow rate and the distribution of the pressure are determined and presented for global pressure driven flows.

MRI of perfluorocarbon emulsion kinetics in rodent mammary tumours

NASA Astrophysics Data System (ADS)

Fan, Xiaobing; River, Jonathan N.; Muresan, Adrian S.; Popescu, Carmen; Zamora, Marta; Culp, Rita M.; Karczmar, Gregory S.

2006-01-01

Perfluorocarbon (PFC) emulsions can be imaged directly by fluorine-19 MRI. We developed an optimized protocol for preparing PFC droplets of uniform size, evaluated use of the resulting droplets as blood pool contrast agents, studied their uptake by tumours and determined the spatial resolution with which they can be imaged at 4.7 T. Perfluorocarbon droplets of three different average sizes (324, 293 and 225 nm) were prepared using a microemulsifier. Images of PFC droplets with good signal-to-noise ratio were acquired with 625 µm in-plane resolution, 3 mm slice thickness and acquisition time of ~4.5 min per image. Kinetics of washout were determined using a simple mathematical model. The maximum uptake of the PFC droplets was three times greater at the tumour rim than in muscle, but the washout rate was two to three times slower in the tumour. The results are consistent with leakage of the droplets into the tumour extravascular space due to the hyper-permeability of tumour capillaries. PFC droplets may allow practical and quantitative measurements of blood volume and capillary permeability in tumours with reasonable spatial resolution.

Effects of Time-Dependent Inflow Perturbations on Turbulent Flow in a Street Canyon

NASA Astrophysics Data System (ADS)

Duan, G.; Ngan, K.

2017-12-01

Urban flow and turbulence are driven by atmospheric flows with larger horizontal scales. Since building-resolving computational fluid dynamics models typically employ steady Dirichlet boundary conditions or forcing, the accuracy of numerical simulations may be limited by the neglect of perturbations. We investigate the sensitivity of flow within a unit-aspect-ratio street canyon to time-dependent perturbations near the inflow boundary. Using large-eddy simulation, time-periodic perturbations to the streamwise velocity component are incorporated via the nudging technique. Spatial averages of pointwise differences between unperturbed and perturbed velocity fields (i.e., the error kinetic energy) show a clear dependence on the perturbation period, though spatial structures are largely insensitive to the time-dependent forcing. The response of the error kinetic energy is maximized for perturbation periods comparable to the time scale of the mean canyon circulation. Frequency spectra indicate that this behaviour arises from a resonance between the inflow forcing and the mean motion around closed streamlines. The robustness of the results is confirmed using perturbations derived from measurements of roof-level wind speed.

The Formation and Spatiotemporal Progress of the pH Wave Induced by the Temperature Gradient in the Thin-Layer H2O2-Na2S2O3-H2SO4-CuSO4 Dynamical System.

PubMed

Jędrusiak, Mikołaj; Orlik, Marek

2016-03-31

The H2O2-S2O3(2-)-H(+)-Cu(2+) dynamical system exhibits sustained oscillations under flow conditions but reveals only a single initial peak of the indicator electrode potential and pH variation under batch isothermal conditions. Thus, in the latter case, there is no possibility of the coupling of the oscillations and diffusion which could lead to formation of sustained spatiotemporal patterns in this process. However, in the inhomogeneous temperature field, due to dependence of the local reaction kinetics on temperature, spatial inhomogeneities of pH distribution can develop which, in the presence of an appropriate indicator, thymol blue, manifest themselves as the color front traveling along the quasi-one-dimensional reactor. In this work, we describe the experimental conditions under which the above-mentioned phenomena can be observed and present their numerical model based on thermokinetic coupling and spatial coordinate introduced to earlier isothermal homogeneous kinetic mechanism.

VOXEL-LEVEL MAPPING OF TRACER KINETICS IN PET STUDIES: A STATISTICAL APPROACH EMPHASIZING TISSUE LIFE TABLES.

PubMed

O'Sullivan, Finbarr; Muzi, Mark; Mankoff, David A; Eary, Janet F; Spence, Alexander M; Krohn, Kenneth A

2014-06-01

Most radiotracers used in dynamic positron emission tomography (PET) scanning act in a linear time-invariant fashion so that the measured time-course data are a convolution between the time course of the tracer in the arterial supply and the local tissue impulse response, known as the tissue residue function. In statistical terms the residue is a life table for the transit time of injected radiotracer atoms. The residue provides a description of the tracer kinetic information measurable by a dynamic PET scan. Decomposition of the residue function allows separation of rapid vascular kinetics from slower blood-tissue exchanges and tissue retention. For voxel-level analysis, we propose that residues be modeled by mixtures of nonparametrically derived basis residues obtained by segmentation of the full data volume. Spatial and temporal aspects of diagnostics associated with voxel-level model fitting are emphasized. Illustrative examples, some involving cancer imaging studies, are presented. Data from cerebral PET scanning with 18 F fluoro-deoxyglucose (FDG) and 15 O water (H2O) in normal subjects is used to evaluate the approach. Cross-validation is used to make regional comparisons between residues estimated using adaptive mixture models with more conventional compartmental modeling techniques. Simulations studies are used to theoretically examine mean square error performance and to explore the benefit of voxel-level analysis when the primary interest is a statistical summary of regional kinetics. The work highlights the contribution that multivariate analysis tools and life-table concepts can make in the recovery of local metabolic information from dynamic PET studies, particularly ones in which the assumptions of compartmental-like models, with residues that are sums of exponentials, might not be certain.

Kinetics of propagation of the lattice excitation in a swift heavy ion track

NASA Astrophysics Data System (ADS)

Lipp, V. P.; Volkov, A. E.; Sorokin, M. V.; Rethfeld, B.

2011-05-01

In this research we verify the applicability of the temperature and heat diffusion conceptions for the description of subpicosecond lattice excitations in nanometric tracks of swift heavy ions (SHI) decelerated in solids in the electronic stopping regime. The method is based on the molecular dynamics (MD) analysis of temporal evolutions of the local kinetic and configurational temperatures of a lattice. We used solid argon as the model system. MD simulations demonstrated that in a SHI track (a) thermalization of lattice excitations takes time of several picoseconds, and (b) application of the parabolic heat diffusion equations for the description of spatial and temporal propagation of lattice excitations is questionable at least up to 10 ps after the ion passage.

Influence of the partial volume correction method on 18F-fluorodeoxyglucose brain kinetic modelling from dynamic PET images reconstructed with resolution model based OSEM

NASA Astrophysics Data System (ADS)

Bowen, Spencer L.; Byars, Larry G.; Michel, Christian J.; Chonde, Daniel B.; Catana, Ciprian

2013-10-01

Kinetic parameters estimated from dynamic 18F-fluorodeoxyglucose (18F-FDG) PET acquisitions have been used frequently to assess brain function in humans. Neglecting partial volume correction (PVC) for a dynamic series has been shown to produce significant bias in model estimates. Accurate PVC requires a space-variant model describing the reconstructed image spatial point spread function (PSF) that accounts for resolution limitations, including non-uniformities across the field of view due to the parallax effect. For ordered subsets expectation maximization (OSEM), image resolution convergence is local and influenced significantly by the number of iterations, the count density, and background-to-target ratio. As both count density and background-to-target values for a brain structure can change during a dynamic scan, the local image resolution may also concurrently vary. When PVC is applied post-reconstruction the kinetic parameter estimates may be biased when neglecting the frame-dependent resolution. We explored the influence of the PVC method and implementation on kinetic parameters estimated by fitting 18F-FDG dynamic data acquired on a dedicated brain PET scanner and reconstructed with and without PSF modelling in the OSEM algorithm. The performance of several PVC algorithms was quantified with a phantom experiment, an anthropomorphic Monte Carlo simulation, and a patient scan. Using the last frame reconstructed image only for regional spread function (RSF) generation, as opposed to computing RSFs for each frame independently, and applying perturbation geometric transfer matrix PVC with PSF based OSEM produced the lowest magnitude bias kinetic parameter estimates in most instances, although at the cost of increased noise compared to the PVC methods utilizing conventional OSEM. Use of the last frame RSFs for PVC with no PSF modelling in the OSEM algorithm produced the lowest bias in cerebral metabolic rate of glucose estimates, although by less than 5% in most cases compared to the other PVC methods. The results indicate that the PVC implementation and choice of PSF modelling in the reconstruction can significantly impact model parameters.

Tracer kinetic modelling for DCE-MRI quantification of subtle blood-brain barrier permeability.

PubMed

Heye, Anna K; Thrippleton, Michael J; Armitage, Paul A; Valdés Hernández, Maria Del C; Makin, Stephen D; Glatz, Andreas; Sakka, Eleni; Wardlaw, Joanna M

2016-01-15

There is evidence that subtle breakdown of the blood-brain barrier (BBB) is a pathophysiological component of several diseases, including cerebral small vessel disease and some dementias. Dynamic contrast-enhanced MRI (DCE-MRI) combined with tracer kinetic modelling is widely used for assessing permeability and perfusion in brain tumours and body tissues where contrast agents readily accumulate in the extracellular space. However, in diseases where leakage is subtle, the optimal approach for measuring BBB integrity is likely to differ since the magnitude and rate of enhancement caused by leakage are extremely low; several methods have been reported in the literature, yielding a wide range of parameters even in healthy subjects. We hypothesised that the Patlak model is a suitable approach for measuring low-level BBB permeability with low temporal resolution and high spatial resolution and brain coverage, and that normal levels of scanner instability would influence permeability measurements. DCE-MRI was performed in a cohort of mild stroke patients (n=201) with a range of cerebral small vessel disease severity. We fitted these data to a set of nested tracer kinetic models, ranking their performance according to the Akaike information criterion. To assess the influence of scanner drift, we scanned 15 healthy volunteers that underwent a "sham" DCE-MRI procedure without administration of contrast agent. Numerical simulations were performed to investigate model validity and the effect of scanner drift. The Patlak model was found to be most appropriate for fitting low-permeability data, and the simulations showed vp and K(Trans) estimates to be reasonably robust to the model assumptions. However, signal drift (measured at approximately 0.1% per minute and comparable to literature reports in other settings) led to systematic errors in calculated tracer kinetic parameters, particularly at low permeabilities. Our findings justify the growing use of the Patlak model in low-permeability states, which has the potential to provide valuable information regarding BBB integrity in a range of diseases. However, absolute values of the resulting tracer kinetic parameters should be interpreted with extreme caution, and the size and influence of signal drift should be measured where possible. Copyright © 2015 The Authors. Published by Elsevier Inc. All rights reserved.

Influence of the partial volume correction method on (18)F-fluorodeoxyglucose brain kinetic modelling from dynamic PET images reconstructed with resolution model based OSEM.

PubMed

Bowen, Spencer L; Byars, Larry G; Michel, Christian J; Chonde, Daniel B; Catana, Ciprian

2013-10-21

Kinetic parameters estimated from dynamic (18)F-fluorodeoxyglucose ((18)F-FDG) PET acquisitions have been used frequently to assess brain function in humans. Neglecting partial volume correction (PVC) for a dynamic series has been shown to produce significant bias in model estimates. Accurate PVC requires a space-variant model describing the reconstructed image spatial point spread function (PSF) that accounts for resolution limitations, including non-uniformities across the field of view due to the parallax effect. For ordered subsets expectation maximization (OSEM), image resolution convergence is local and influenced significantly by the number of iterations, the count density, and background-to-target ratio. As both count density and background-to-target values for a brain structure can change during a dynamic scan, the local image resolution may also concurrently vary. When PVC is applied post-reconstruction the kinetic parameter estimates may be biased when neglecting the frame-dependent resolution. We explored the influence of the PVC method and implementation on kinetic parameters estimated by fitting (18)F-FDG dynamic data acquired on a dedicated brain PET scanner and reconstructed with and without PSF modelling in the OSEM algorithm. The performance of several PVC algorithms was quantified with a phantom experiment, an anthropomorphic Monte Carlo simulation, and a patient scan. Using the last frame reconstructed image only for regional spread function (RSF) generation, as opposed to computing RSFs for each frame independently, and applying perturbation geometric transfer matrix PVC with PSF based OSEM produced the lowest magnitude bias kinetic parameter estimates in most instances, although at the cost of increased noise compared to the PVC methods utilizing conventional OSEM. Use of the last frame RSFs for PVC with no PSF modelling in the OSEM algorithm produced the lowest bias in cerebral metabolic rate of glucose estimates, although by less than 5% in most cases compared to the other PVC methods. The results indicate that the PVC implementation and choice of PSF modelling in the reconstruction can significantly impact model parameters.

Tracer kinetic modelling for DCE-MRI quantification of subtle blood–brain barrier permeability

PubMed Central

Heye, Anna K.; Thrippleton, Michael J.; Armitage, Paul A.; Valdés Hernández, Maria del C.; Makin, Stephen D.; Glatz, Andreas; Sakka, Eleni; Wardlaw, Joanna M.

2016-01-01

There is evidence that subtle breakdown of the blood–brain barrier (BBB) is a pathophysiological component of several diseases, including cerebral small vessel disease and some dementias. Dynamic contrast-enhanced MRI (DCE-MRI) combined with tracer kinetic modelling is widely used for assessing permeability and perfusion in brain tumours and body tissues where contrast agents readily accumulate in the extracellular space. However, in diseases where leakage is subtle, the optimal approach for measuring BBB integrity is likely to differ since the magnitude and rate of enhancement caused by leakage are extremely low; several methods have been reported in the literature, yielding a wide range of parameters even in healthy subjects. We hypothesised that the Patlak model is a suitable approach for measuring low-level BBB permeability with low temporal resolution and high spatial resolution and brain coverage, and that normal levels of scanner instability would influence permeability measurements. DCE-MRI was performed in a cohort of mild stroke patients (n = 201) with a range of cerebral small vessel disease severity. We fitted these data to a set of nested tracer kinetic models, ranking their performance according to the Akaike information criterion. To assess the influence of scanner drift, we scanned 15 healthy volunteers that underwent a “sham” DCE-MRI procedure without administration of contrast agent. Numerical simulations were performed to investigate model validity and the effect of scanner drift. The Patlak model was found to be most appropriate for fitting low-permeability data, and the simulations showed vp and KTrans estimates to be reasonably robust to the model assumptions. However, signal drift (measured at approximately 0.1% per minute and comparable to literature reports in other settings) led to systematic errors in calculated tracer kinetic parameters, particularly at low permeabilities. Our findings justify the growing use of the Patlak model in low-permeability states, which has the potential to provide valuable information regarding BBB integrity in a range of diseases. However, absolute values of the resulting tracer kinetic parameters should be interpreted with extreme caution, and the size and influence of signal drift should be measured where possible. PMID:26477653

Investigation of plasma-sheath resonances in low pressure discharges

NASA Astrophysics Data System (ADS)

Naggary, Schabnam; Kemaneci, Efe; Brinkmann, Ralf Peter; Megahed, Mustafa

2016-09-01

Plasma sheath resonances (PSR) arise from a periodic exchange between the kinetic electron energy in the plasma bulk and the electric field energy in the sheath and can easily be excited by the sheath-generated harmonics of the applied RF. In this contribution, we employ a series of models to obtain a well-defined description of these phenomena. In the first part, we use a global model to study the influence of the nonlinear charge-voltage characteristics on the electron dynamics. However, the global model is restricted to the assumption of spatially constant potential at each driven and grounded electrode and thus delivers only the fundamental mode of the current. In order to remedy the deficiency, we introduce a spatially resolved model for arbitrary reactor geometries with no assumptions on the homogeneity of the plasma. An exact evaluation of the analytical solution is realized on the assumption of a cylinderical plasma reactor geometry with uniform conductance. Furthermore, the spatially resolved model is capable of being utilized for a more realistic CCP reactor geometry and non homogeneous plasma provided the conductance distribution is known. For this purpose, we use the CFD-ACE+ tool. The results show that the proposed multi-mode model provides a significant improvement. The authors gratefully acknowledge the financial support by the ESI Group and the SFB- TR 87.

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

Ullah, Ghanim; Demuro, Angelo; Parker, Ian

Amyloid beta (Aβ) oligomers associated with Alzheimer’s disease (AD) form Ca 2+-permeable plasma membrane pores, leading to a disruption of the otherwise well-controlled intracellular calcium (Ca 2+) homeostasis. The resultant up-regulation of intracellular Ca 2+ concentration has detrimental implications for memory formation and cell survival. The gating kinetics and Ca 2+ permeability of Aβ pores are not well understood. We have used computational modeling in conjunction with the ability of optical patch-clamping for massively parallel imaging of Ca 2+ flux through thousands of pores in the cell membrane of Xenopus oocytes to elucidate the kinetic properties of Aβ pores. Themore » fluorescence time-series data from individual pores were idealized and used to develop data-driven Markov chain models for the kinetics of the Aβ pore at different stages of its evolution. Our study provides the first demonstration of developing Markov chain models for ion channel gating that are driven by optical-patch clamp data with the advantage of experiments being performed under close to physiological conditions. As a result, we demonstrate the up-regulation of gating of various Ca 2+ release channels due to Aβ pores and show that the extent and spatial range of such up-regulation increases as Aβ pores with low open probability and Ca 2+ permeability transition into those with high open probability and Ca 2+ permeability.« less

Quantifying fat, oil, and grease deposit formation kinetics.

PubMed

Iasmin, Mahbuba; Dean, Lisa O; Ducoste, Joel J

2016-01-01

Fat, oil, and grease (FOG) deposits formed in sanitary sewers are calcium-based saponified solids that are responsible for a significant number of nationwide sanitary sewer overflows (SSOs) across United States. In the current study, the kinetics of lab-based saponified solids were determined to understand the kinetics of FOG deposit formation in sewers for two types of fat (Canola and Beef Tallow) and two types of calcium sources (calcium chloride and calcium sulfate) under three pH (7 ± 0.5, 10 ± 0.5, and ≈14) and two temperature conditions (22 ± 0.5 and 45 ± 0.5 °C). The results of this study displayed quick reactions of a fraction of fats with calcium ions to form calcium based saponified solids. Results further showed that increased palmitic fatty acid content in source fats, the magnitude of the pH, and temperature significantly affect the FOG deposit formation and saponification rates. The experimental data of the kinetics were compared with two empirical models: a) Cotte saponification model and b) Foubert crystallization model and a mass-action based mechanistic model that included alkali driven hydrolysis of triglycerides. Results showed that the mass action based mechanistic model was able to predict changes in the rate of formation of saponified solids under the different experimental conditions compared to both empirical models. The mass-action based saponification model also revealed that the hydrolysis of Beef Tallow was slower compared to liquid Canola fat resulting in smaller quantities of saponified solids. This mechanistic saponification model, with its ability to track the saponified solids chemical precursors, may provide an initial framework to predict the spatial formation of FOG deposits in municipal sewers using system wide sewer collection modeling software. Copyright © 2015 Elsevier Ltd. All rights reserved.

Comparison of kinetic and extended magnetohydrodynamics computational models for the linear ion temperature gradient instability in slab geometry

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

Schnack, D. D.; Department of Physics, University of Wisconsin-Madison, Madison, Wisconsin 53706; Cheng, J.

We perform linear stability studies of the ion temperature gradient (ITG) instability in unsheared slab geometry using kinetic and extended magnetohydrodynamics (MHD) models, in the regime k{sub ∥}/k{sub ⊥}≪1. The ITG is a parallel (to B) sound wave that may be destabilized by finite ion Larmor radius (FLR) effects in the presence of a gradient in the equilibrium ion temperature. The ITG is stable in both ideal and resistive MHD; for a given temperature scale length L{sub Ti0}, instability requires that either k{sub ⊥}ρ{sub i} or ρ{sub i}/L{sub Ti0} be sufficiently large. Kinetic models capture FLR effects to all ordersmore » in either parameter. In the extended MHD model, these effects are captured only to lowest order by means of the Braginskii ion gyro-viscous stress tensor and the ion diamagnetic heat flux. We present the linear electrostatic dispersion relations for the ITG for both kinetic Vlasov and extended MHD (two-fluid) models in the local approximation. In the low frequency fluid regime, these reduce to the same cubic equation for the complex eigenvalue ω=ω{sub r}+iγ. An explicit solution is derived for the growth rate and real frequency in this regime. These are found to depend on a single non-dimensional parameter. We also compute the eigenvalues and the eigenfunctions with the extended MHD code NIMROD, and a hybrid kinetic δf code that assumes six-dimensional Vlasov ions and isothermal fluid electrons, as functions of k{sub ⊥}ρ{sub i} and ρ{sub i}/L{sub Ti0} using a spatially dependent equilibrium. These solutions are compared with each other, and with the predictions of the local kinetic and fluid dispersion relations. Kinetic and fluid calculations agree well at and near the marginal stability point, but diverge as k{sub ⊥}ρ{sub i} or ρ{sub i}/L{sub Ti0} increases. There is good qualitative agreement between the models for the shape of the unstable global eigenfunction for L{sub Ti0}/ρ{sub i}=30 and 20. The results quantify how far fluid calculations can be extended accurately into the kinetic regime. We conclude that for the linear ITG problem in slab geometry with unsheared magnetic field when k{sub ∥}/k{sub ⊥}≪1, the extended MHD model may be a reliable physical model for this problem when ρ{sub i}/L{sub Ti0}<10{sup −2} and k{sub ⊥}ρ{sub i}<0.2.« less

Generalized fluid theory including non-Maxwellian kinetic effects

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

Izacard, Olivier

The results obtained by the plasma physics community for the validation and the prediction of turbulence and transport in magnetized plasmas come mainly from the use of very central processing unit (CPU)-consuming particle-in-cell or (gyro)kinetic codes which naturally include non-Maxwellian kinetic effects. To date, fluid codes are not considered to be relevant for the description of these kinetic effects. Here, after revisiting the limitations of the current fluid theory developed in the 19th century, we generalize the fluid theory including kinetic effects such as non-Maxwellian super-thermal tails with as few fluid equations as possible. The collisionless and collisional fluid closuresmore » from the nonlinear Landau Fokker–Planck collision operator are shown for an arbitrary collisionality. Indeed, the first fluid models associated with two examples of collisionless fluid closures are obtained by assuming an analytic non-Maxwellian distribution function. One of the main differences with the literature is our analytic representation of the distribution function in the velocity phase space with as few hidden variables as possible thanks to the use of non-orthogonal basis sets. These new non-Maxwellian fluid equations could initiate the next generation of fluid codes including kinetic effects and can be expanded to other scientific disciplines such as astrophysics, condensed matter or hydrodynamics. As a validation test, we perform a numerical simulation based on a minimal reduced INMDF fluid model. The result of this test is the discovery of the origin of particle and heat diffusion. The diffusion is due to the competition between a growing INMDF on short time scales due to spatial gradients and the thermalization on longer time scales. Here, the results shown here could provide the insights to break some of the unsolved puzzles of turbulence.« less

Generalized fluid theory including non-Maxwellian kinetic effects

DOE PAGES

Izacard, Olivier

2017-03-29

The results obtained by the plasma physics community for the validation and the prediction of turbulence and transport in magnetized plasmas come mainly from the use of very central processing unit (CPU)-consuming particle-in-cell or (gyro)kinetic codes which naturally include non-Maxwellian kinetic effects. To date, fluid codes are not considered to be relevant for the description of these kinetic effects. Here, after revisiting the limitations of the current fluid theory developed in the 19th century, we generalize the fluid theory including kinetic effects such as non-Maxwellian super-thermal tails with as few fluid equations as possible. The collisionless and collisional fluid closuresmore » from the nonlinear Landau Fokker–Planck collision operator are shown for an arbitrary collisionality. Indeed, the first fluid models associated with two examples of collisionless fluid closures are obtained by assuming an analytic non-Maxwellian distribution function. One of the main differences with the literature is our analytic representation of the distribution function in the velocity phase space with as few hidden variables as possible thanks to the use of non-orthogonal basis sets. These new non-Maxwellian fluid equations could initiate the next generation of fluid codes including kinetic effects and can be expanded to other scientific disciplines such as astrophysics, condensed matter or hydrodynamics. As a validation test, we perform a numerical simulation based on a minimal reduced INMDF fluid model. The result of this test is the discovery of the origin of particle and heat diffusion. The diffusion is due to the competition between a growing INMDF on short time scales due to spatial gradients and the thermalization on longer time scales. Here, the results shown here could provide the insights to break some of the unsolved puzzles of turbulence.« less

Autocatalytic polymerization generates persistent random walk of crawling cells.

PubMed

Sambeth, R; Baumgaertner, A

2001-05-28

The autocatalytic polymerization kinetics of the cytoskeletal actin network provides the basic mechanism for a persistent random walk of a crawling cell. It is shown that network remodeling by branching processes near the cell membrane is essential for the bimodal spatial stability of the network which induces a spontaneous breaking of isotropic cell motion. Details of the phenomena are analyzed using a simple polymerization model studied by analytical and simulation methods.

Unchanged mitochondrial organization and compartmentation of high-energy phosphates in creatine-deficient GAMT−/− mouse hearts

PubMed Central

Branovets, Jelena; Sepp, Mervi; Kotlyarova, Svetlana; Jepihhina, Natalja; Sokolova, Niina; Aksentijevic, Dunja; Lygate, Craig A.; Neubauer, Stefan; Birkedal, Rikke

2013-01-01

Disruption of the creatine kinase (CK) system in hearts of CK-deficient mice leads to changes in the ultrastructure and regulation of mitochondrial respiration. We expected to see similar changes in creatine-deficient mice, which lack the enzyme guanidinoacetate methyltransferase (GAMT) to produce creatine. The aim of this study was to characterize the changes in cardiomyocyte mitochondrial organization, regulation of respiration, and intracellular compartmentation associated with GAMT deficiency. Three-dimensional mitochondrial organization was assessed by confocal microscopy. On populations of permeabilized cardiomyocytes, we recorded ADP and ATP kinetics of respiration, competition between mitochondria and pyruvate kinase for ADP produced by ATPases, ADP kinetics of endogenous pyruvate kinase, and ATP kinetics of ATPases. These data were analyzed by mathematical models to estimate intracellular compartmentation. Quantitative analysis of morphological and kinetic data as well as derived model fits showed no difference between GAMT-deficient and wild-type mice. We conclude that inactivation of the CK system by GAMT deficiency does not alter mitochondrial organization and intracellular compartmentation in relaxed cardiomyocytes. Thus, our results suggest that the healthy heart is able to preserve cardiac function at a basal level in the absence of CK-facilitated energy transfer without compromising intracellular organization and the regulation of mitochondrial energy homeostasis. This raises questions on the importance of the CK system as a spatial energy buffer in unstressed cardiomyocytes. PMID:23792673

Quantitative modeling of the reaction/diffusion kinetics of two-chemistry photopolymers

NASA Astrophysics Data System (ADS)

Kowalski, Benjamin Andrew

Optically driven diffusion in photopolymers is an appealing material platform for a broad range of applications, in which the recorded refractive index patterns serve either as images (e.g. data storage, display holography) or as optical elements (e.g. custom GRIN components, integrated optical devices). A quantitative understanding of the reaction/diffusion kinetics is difficult to obtain directly, but is nevertheless necessary in order to fully exploit the wide array of design freedoms in these materials. A general strategy for characterizing these kinetics is proposed, in which key processes are decoupled and independently measured. This strategy enables prediction of a material's potential refractive index change, solely on the basis of its chemical components. The degree to which a material does not reach this potential reveals the fraction of monomer that has participated in unwanted reactions, reducing spatial resolution and dynamic range. This approach is demonstrated for a model material similar to commercial media, achieving quantitative predictions of index response over three orders of exposure dose (~1 to ~103 mJ cm-2) and three orders of feature size (0.35 to 500 microns). The resulting insights enable guided, rational design of new material formulations with demonstrated performance improvement.

Spatio-temporal diffusion of dynamic PET images

NASA Astrophysics Data System (ADS)

Tauber, C.; Stute, S.; Chau, M.; Spiteri, P.; Chalon, S.; Guilloteau, D.; Buvat, I.

2011-10-01

Positron emission tomography (PET) images are corrupted by noise. This is especially true in dynamic PET imaging where short frames are required to capture the peak of activity concentration after the radiotracer injection. High noise results in a possible bias in quantification, as the compartmental models used to estimate the kinetic parameters are sensitive to noise. This paper describes a new post-reconstruction filter to increase the signal-to-noise ratio in dynamic PET imaging. It consists in a spatio-temporal robust diffusion of the 4D image based on the time activity curve (TAC) in each voxel. It reduces the noise in homogeneous areas while preserving the distinct kinetics in regions of interest corresponding to different underlying physiological processes. Neither anatomical priors nor the kinetic model are required. We propose an automatic selection of the scale parameter involved in the diffusion process based on a robust statistical analysis of the distances between TACs. The method is evaluated using Monte Carlo simulations of brain activity distributions. We demonstrate the usefulness of the method and its superior performance over two other post-reconstruction spatial and temporal filters. Our simulations suggest that the proposed method can be used to significantly increase the signal-to-noise ratio in dynamic PET imaging.

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

Nguyen, Thao; Luscher, D. J.; Wilkerson, J. W.

We developed a framework for dislocation-based viscoplasticity and dynamic ductile failure to model high strain rate deformation and damage in single crystals. The rate-dependence of the crystal plasticity formulation is based on the physics of relativistic dislocation kinetics suited for extremely high strain rates. The damage evolution is based on the dynamics of void growth, which are governed by both micro-inertia as well as dislocation kinetics and dislocation substructure evolution. Furthermore, an averaging scheme is proposed in order to approximate the evolution of the dislocation substructure in both the macroscale as well as its spatial distribution at the microscale. Inmore » addition, a concept of a single equivalent dislocation density that effectively captures the collective influence of dislocation density on all active slip systems is proposed here. Together, these concepts and approximations enable the use of semi-analytic solutions for void growth dynamics developed in [J. Wilkerson and K. Ramesh. A dynamic void growth model governed by dislocation kinetics. J. Mech. Phys. Solids, 70:262–280, 2014.], which greatly reduce the computational overhead that would otherwise be required. The resulting homogenized framework has been implemented into a commercially available finite element package, and a validation study against a suite of direct numerical simulations was carried out.« less

Three-dimensional spectral-spatial EPR imaging of free radicals in the heart: a technique for imaging tissue metabolism and oxygenation.

PubMed Central

Kuppusamy, P; Chzhan, M; Vij, K; Shteynbuk, M; Lefer, D J; Giannella, E; Zweier, J L

1994-01-01

It has been hypothesized that free radical metabolism and oxygenation in living organs and tissues such as the heart may vary over the spatially defined tissue structure. In an effort to study these spatially defined differences, we have developed electron paramagnetic resonance imaging instrumentation enabling the performance of three-dimensional spectral-spatial images of free radicals infused into the heart and large vessels. Using this instrumentation, high-quality three-dimensional spectral-spatial images of isolated perfused rat hearts and rabbit aortas are obtained. In the isolated aorta, it is shown that spatially and spectrally accurate images of the vessel lumen and wall could be obtained in this living vascular tissue. In the isolated rat heart, imaging experiments were performed to determine the kinetics of radical clearance at different spatial locations within the heart during myocardial ischemia. The kinetic data show the existence of regional and transmural differences in myocardial free radical clearance. It is further demonstrated that EPR imaging can be used to noninvasively measure spatially localized oxygen concentrations in the heart. Thus, the technique of spectral-spatial EPR imaging is shown to be a powerful tool in providing spatial information regarding the free radical distribution, metabolism, and tissue oxygenation in living biological organs and tissues. Images PMID:8159757

RockFall analyst: A GIS extension for three-dimensional and spatially distributed rockfall hazard modeling

NASA Astrophysics Data System (ADS)

Lan, Hengxing; Derek Martin, C.; Lim, C. H.

2007-02-01

Geographic information system (GIS) modeling is used in combination with three-dimensional (3D) rockfall process modeling to assess rockfall hazards. A GIS extension, RockFall Analyst (RA), which is capable of effectively handling large amounts of geospatial information relative to rockfall behaviors, has been developed in ArcGIS using ArcObjects and C#. The 3D rockfall model considers dynamic processes on a cell plane basis. It uses inputs of distributed parameters in terms of raster and polygon features created in GIS. Two major components are included in RA: particle-based rockfall process modeling and geostatistics-based rockfall raster modeling. Rockfall process simulation results, 3D rockfall trajectories and their velocity features either for point seeders or polyline seeders are stored in 3D shape files. Distributed raster modeling, based on 3D rockfall trajectories and a spatial geostatistical technique, represents the distribution of spatial frequency, the flying and/or bouncing height, and the kinetic energy of falling rocks. A distribution of rockfall hazard can be created by taking these rockfall characteristics into account. A barrier analysis tool is also provided in RA to aid barrier design. An application of these modeling techniques to a case study is provided. The RA has been tested in ArcGIS 8.2, 8.3, 9.0 and 9.1.

Dynamics of Reactive Microbial Hotspots in Concentration Gradient.

NASA Astrophysics Data System (ADS)

Hubert, A.; Farasin, J.; Tabuteau, H.; Dufresne, A.; Meheust, Y.; Le Borgne, T.

2017-12-01

In subsurface environments, bacteria play a major role in controlling the kinetics of a broad range of biogeochemical reactions. In such environments, nutrients fluxes and solute concentrations needed for bacteria metabolism may be highly variable in space and intermittent in time. This can lead to the formation of reactive hotspots where and when conditions are favorable to particular microorganisms, hence inducing biogeochemical reaction kinetics that differ significantly from those measured in homogeneous model environments. To investigate the impact of chemical gradients on the spatial structure and temporal dynamics of subsurface microorganism populations, we develop microfluidic cells allowing for a precise control of flow and chemical gradient conditions, as well as quantitative monitoring of the bacteria's spatial distribution and biofilm development. Using the non-motile Escherichia coli JW1908-1 strain and Gallionella capsiferriformans ES-2 as model organisms, we investigate the behavior and development of bacteria over a range of single and double concentration gradients in the concentrations of nutrients, electron donors and electron acceptors. We measure bacterial activity and population growth locally in precisely known hydrodynamic and chemical environments. This approach allows time-resolved monitoring of the location and intensity of reactive hotspots in micromodels as a function of the flow and chemical gradient conditions. We compare reactive microbial hotspot dynamics in our micromodels to classic growth laws and well-known growth parameters for the laboratory model bacteria Escherichia coli.We also discuss consequences for the formation and temporal dynamics of biofilms in the subsurface.

Foundations of modelling of nonequilibrium low-temperature plasmas

NASA Astrophysics Data System (ADS)

Alves, L. L.; Bogaerts, A.; Guerra, V.; Turner, M. M.

2018-02-01

This work explains the need for plasma models, introduces arguments for choosing the type of model that better fits the purpose of each study, and presents the basics of the most common nonequilibrium low-temperature plasma models and the information available from each one, along with an extensive list of references for complementary in-depth reading. The paper presents the following models, organised according to the level of multi-dimensional description of the plasma: kinetic models, based on either a statistical particle-in-cell/Monte-Carlo approach or the solution to the Boltzmann equation (in the latter case, special focus is given to the description of the electron kinetics); multi-fluid models, based on the solution to the hydrodynamic equations; global (spatially-average) models, based on the solution to the particle and energy rate-balance equations for the main plasma species, usually including a very complete reaction chemistry; mesoscopic models for plasma-surface interaction, adopting either a deterministic approach or a stochastic dynamical Monte-Carlo approach. For each plasma model, the paper puts forward the physics context, introduces the fundamental equations, presents advantages and limitations, also from a numerical perspective, and illustrates its application with some examples. Whenever pertinent, the interconnection between models is also discussed, in view of multi-scale hybrid approaches.

Modeling photo-bleaching kinetics to map local variations in rod rhodopsin density

NASA Astrophysics Data System (ADS)

Ehler, M.; Dobrosotskaya, J.; King, E. J.; Czaja, W.; Bonner, R. F.

2011-03-01

Localized rod photoreceptor and rhodopsin losses have been observed in post mortem histology both in normal aging and in age-related maculopathy. We propose to noninvasively map local rod rhodopsin density through analysis of the brightening of the underlying lipofuscin autofluorescence (LAF) in confocal scanning laser ophthalmoscopy (cSLO) imaging sequences starting in the dark adapted eye. The detected LAF increases as rhodopsin is bleached (time constant ~ 25sec) by the average retinal irradiance of the cSLO 488nm laser beam. We fit parameters of analytical expressions for the kinetics of rhodopsin bleaching that Lamb validated using electroretinogram recordings in human. By performing localized (~ 100μm) kinetic analysis, we create high resolution maps of the rhodopsin density. This new noninvasive imaging and analysis approach appears well-suited for measuring localized changes in the rod photoreceptors and correlating them at high spatial resolution with localized pathological changes of the retinal pigment epithelium (RPE) seen in steady-state LAF images.

A model of early formation of uranium molecular oxides in laser-ablated plasmas

NASA Astrophysics Data System (ADS)

Finko, Mikhail; Curreli, Davide; Azer, Magdi; Weisz, David; Crowhurst, Jonathan; Rose, Timothy; Koroglu, Batikan; Radousky, Harry; Zaug, Joseph; Armstrong, Mike

2017-10-01

An important problem within the field of nuclear forensics is fractionation: the formation of post-detonation nuclear debris whose composition does not reflect that of the source weapon. We are investigating uranium fractionation in rapidly cooling plasma using a combined experimental and modeling approach. In particular, we use laser ablation of uranium metal samples to produce a low-temperature plasma with physical conditions similar to a condensing nuclear fireball. Here we present a first plasma-chemistry model of uranium molecular species formation during the early stage of laser ablated plasma evolution in atmospheric oxygen. The system is simulated using a global kinetic model with rate coefficients calculated according to literature data and the application of reaction rate theory. The model allows for a detailed analysis of the evolution of key uranium molecular species and represents the first step in producing a uranium fireball model that is kinetically validated against spatially and temporally resolved spectroscopy measurements. This project was sponsored by the DoD, Defense Threat Reduction Agency, Grant HDTRA1-16- 1-0020. This work was performed in part under the auspices of the U.S. DoE by Lawrence Livermore National Laboratory under Contract DE-AC52- 07NA27344.

Implementation of logic functions and computations by chemical kinetics

NASA Astrophysics Data System (ADS)

Hjelmfelt, A.; Ross, J.

We review our work on the computational functions of the kinetics of chemical networks. We examine spatially homogeneous networks which are based on prototypical reactions occurring in living cells and show the construction of logic gates and sequential and parallel networks. This work motivates the study of an important biochemical pathway, glycolysis, and we demonstrate that the switch that controls the flux in the direction of glycolysis or gluconeogenesis may be described as a fuzzy AND operator. We also study a spatially inhomogeneous network which shares features of theoretical and biological neural networks.

High fidelity kinetic modeling of magnetic reconnection in laboratory plasma

NASA Astrophysics Data System (ADS)

Stanier, A.; Daughton, W. S.

2017-12-01

Over the past decade, a great deal of progress has been made towards understanding the physics of magnetic reconnection in weakly collisional regimes of relevance to both fusion devices, and to space and astrophysical plasmas. However, there remain some outstanding unsolved problems in reconnection physics, such as the generation and influence of plasmoids (flux ropes) within reconnection layers, the development of magnetic turbulence, the role of current driven and streaming instabilities, and the influence of electron pressure anisotropy on the layer structure. Due to the importance of these questions, new laboratory reconnection experiments are being built to allow controlled and reproducible study of such questions with the simultaneous acquisition of high time resolution measurements at a large number of spatial points. These experiments include the FLARE facility at Princeton University and the T-REX experiment at the University of Wisconsin. To guide and interpret these new experiments, and to extrapolate the results to space applications, new investments in kinetic modeling tools are required. We have recently developed a cylindrical version of the VPIC Particle-In-Cell code with the capability to perform first-principles kinetic simulations that approach experimental device size with more realistic geometry and drive coils. This cylindrical version inherits much of the optimization work that has been done recently for the next generation many-cores architectures with wider vector registers, and achieves comparable conservation properties as the Cartesian code. Namely it features exact discrete charge conservation, and a so-called "energy-conserving" scheme where the energy is conserved in the limit of continuous time, i.e. without contribution from spatial discretization (Lewis, 1970). We will present initial results of modeling magnetic reconnection in the experiments mentioned above. Since the VPIC code is open source (https://github.com/losalamos/vpic), this new cylindrical version will also be freely available to the community.

Modeling Dynamic Contrast-Enhanced MRI Data with a Constrained Local AIF.

PubMed

Duan, Chong; Kallehauge, Jesper F; Pérez-Torres, Carlos J; Bretthorst, G Larry; Beeman, Scott C; Tanderup, Kari; Ackerman, Joseph J H; Garbow, Joel R

2018-02-01

This study aims to develop a constrained local arterial input function (cL-AIF) to improve quantitative analysis of dynamic contrast-enhanced (DCE)-magnetic resonance imaging (MRI) data by accounting for the contrast-agent bolus amplitude error in the voxel-specific AIF. Bayesian probability theory-based parameter estimation and model selection were used to compare tracer kinetic modeling employing either the measured remote-AIF (R-AIF, i.e., the traditional approach) or an inferred cL-AIF against both in silico DCE-MRI data and clinical, cervical cancer DCE-MRI data. When the data model included the cL-AIF, tracer kinetic parameters were correctly estimated from in silico data under contrast-to-noise conditions typical of clinical DCE-MRI experiments. Considering the clinical cervical cancer data, Bayesian model selection was performed for all tumor voxels of the 16 patients (35,602 voxels in total). Among those voxels, a tracer kinetic model that employed the voxel-specific cL-AIF was preferred (i.e., had a higher posterior probability) in 80 % of the voxels compared to the direct use of a single R-AIF. Maps of spatial variation in voxel-specific AIF bolus amplitude and arrival time for heterogeneous tissues, such as cervical cancer, are accessible with the cL-AIF approach. The cL-AIF method, which estimates unique local-AIF amplitude and arrival time for each voxel within the tissue of interest, provides better modeling of DCE-MRI data than the use of a single, measured R-AIF. The Bayesian-based data analysis described herein affords estimates of uncertainties for each model parameter, via posterior probability density functions, and voxel-wise comparison across methods/models, via model selection in data modeling.

A generalized electron energy probability function for inductively coupled plasmas under conditions of nonlocal electron kinetics

NASA Astrophysics Data System (ADS)

Mouchtouris, S.; Kokkoris, G.

2018-01-01

A generalized equation for the electron energy probability function (EEPF) of inductively coupled Ar plasmas is proposed under conditions of nonlocal electron kinetics and diffusive cooling. The proposed equation describes the local EEPF in a discharge and the independent variable is the kinetic energy of electrons. The EEPF consists of a bulk and a depleted tail part and incorporates the effect of the plasma potential, Vp, and pressure. Due to diffusive cooling, the break point of the EEPF is eVp. The pressure alters the shape of the bulk and the slope of the tail part. The parameters of the proposed EEPF are extracted by fitting to measure EEPFs (at one point in the reactor) at different pressures. By coupling the proposed EEPF with a hybrid plasma model, measurements in the gaseous electronics conference reference reactor concerning (a) the electron density and temperature and the plasma potential, either spatially resolved or at different pressure (10-50 mTorr) and power, and (b) the ion current density of the electrode, are well reproduced. The effect of the choice of the EEPF on the results is investigated by a comparison to an EEPF coming from the Boltzmann equation (local electron kinetics approach) and to a Maxwellian EEPF. The accuracy of the results and the fact that the proposed EEPF is predefined renders its use a reliable alternative with a low computational cost compared to stochastic electron kinetic models at low pressure conditions, which can be extended to other gases and/or different electron heating mechanisms.

GoPhast: a graphical user interface for PHAST

USGS Publications Warehouse

Winston, Richard B.

2006-01-01

GoPhast is a graphical user interface (GUI) for the USGS model PHAST. PHAST simulates multicomponent, reactive solute transport in three-dimensional, saturated, ground-water flow systems. PHAST can model both equilibrium and kinetic geochemical reactions. PHAST is derived from HST3D (flow and transport) and PHREEQC (geochemical calculations). The flow and transport calculations are restricted to constant fluid density and constant temperature. The complexity of the input required by PHAST makes manual construction of its input files tedious and error-prone. GoPhast streamlines the creation of the input file and helps reduce errors. GoPhast allows the user to define the spatial input for the PHAST flow and transport data file by drawing points, lines, or polygons on top, front, and side views of the model domain. These objects can have up to two associated formulas that define their extent perpendicular to the view plane, allowing the objects to be three-dimensional. Formulas are also used to specify the values of spatial data (data sets) both globally and for individual objects. Objects can be used to specify the values of data sets independent of the spatial and temporal discretization of the model. Thus, the grid and simulation periods for the model can be changed without respecifying spatial data pertaining to the hydrogeologic framework and boundary conditions. This report describes the operation of GoPhast and demonstrates its use with examples. GoPhast runs on Windows 2000, Windows XP, and Linux operating systems.

Actin dynamics at the living cell submembrane imaged by total internal reflection fluorescence photobleaching.

PubMed Central

Sund, S E; Axelrod, D

2000-01-01

Although reversible chemistry is crucial to dynamical processes in living cells, relatively little is known about relevant chemical kinetic rates in vivo. Total internal reflection/fluorescence recovery after photobleaching (TIR/FRAP), an established technique previously demonstrated to measure reversible biomolecular kinetic rates at surfaces in vitro, is extended here to measure reversible biomolecular kinetic rates of actin at the cytofacial (subplasma membrane) surface of living cells. For the first time, spatial imaging (with a charge-coupled device camera) is used in conjunction with TIR/FRAP. TIR/FRAP imaging produces both spatial maps of kinetic parameters (off-rates and mobile fractions) and estimates of kinetic correlation distances, cell-wide kinetic gradients, and dependences of kinetic parameters on initial fluorescence intensity. For microinjected rhodamine actin in living cultured smooth muscle (BC3H1) cells, the unbinding rate at or near the cytofacial surface of the plasma membrane (averaged over the entire cell) is measured at 0.032 +/- 0.007 s(-1). The corresponding rate for actin marked by microinjected rhodamine phalloidin is very similar, 0.033 +/- 0.013 s(-1), suggesting that TIR/FRAP is reporting the dynamics of entire filaments or protofilaments. For submembrane fluorescence-marked actin, the intensity, off-rate, and mobile fraction show a positive correlation over a characteristic distance of 1-3 microm and a negative correlation over larger distances greater than approximately 7-14 microm. Furthermore, the kinetic parameters display a statistically significant cell-wide gradient, with the cell having a "fast" and "slow" end with respect to actin kinetics. PMID:10969025

Nanoscale dynamics of Joule heating and bubble nucleation in a solid-state nanopore.

PubMed

Levine, Edlyn V; Burns, Michael M; Golovchenko, Jene A

2016-01-01

We present a mathematical model for Joule heating of an electrolytic solution in a nanopore. The model couples the electrical and thermal dynamics responsible for rapid and extreme superheating of the electrolyte within the nanopore. The model is implemented numerically with a finite element calculation, yielding a time and spatially resolved temperature distribution in the nanopore region. Temperatures near the thermodynamic limit of superheat are predicted to be attained just before the explosive nucleation of a vapor bubble is observed experimentally. Knowledge of this temperature distribution enables the evaluation of related phenomena including bubble nucleation kinetics, relaxation oscillation, and bubble dynamics.

Optimization and universality of Brownian search in a basic model of quenched heterogeneous media

NASA Astrophysics Data System (ADS)

Godec, Aljaž; Metzler, Ralf

2015-05-01

The kinetics of a variety of transport-controlled processes can be reduced to the problem of determining the mean time needed to arrive at a given location for the first time, the so-called mean first-passage time (MFPT) problem. The occurrence of occasional large jumps or intermittent patterns combining various types of motion are known to outperform the standard random walk with respect to the MFPT, by reducing oversampling of space. Here we show that a regular but spatially heterogeneous random walk can significantly and universally enhance the search in any spatial dimension. In a generic minimal model we consider a spherically symmetric system comprising two concentric regions with piecewise constant diffusivity. The MFPT is analyzed under the constraint of conserved average dynamics, that is, the spatially averaged diffusivity is kept constant. Our analytical calculations and extensive numerical simulations demonstrate the existence of an optimal heterogeneity minimizing the MFPT to the target. We prove that the MFPT for a random walk is completely dominated by what we term direct trajectories towards the target and reveal a remarkable universality of the spatially heterogeneous search with respect to target size and system dimensionality. In contrast to intermittent strategies, which are most profitable in low spatial dimensions, the spatially inhomogeneous search performs best in higher dimensions. Discussing our results alongside recent experiments on single-particle tracking in living cells, we argue that the observed spatial heterogeneity may be beneficial for cellular signaling processes.

Hydrodynamic effects of kinetic power extraction by in-stream tidal turbines

NASA Astrophysics Data System (ADS)

Polagye, Brian L.

The hydrodynamic effects of extracting kinetic power from tidal streams presents unique challenges to the development of in-stream tidal power. In-stream tidal turbines superficially resemble wind turbines and extract kinetic power from the ebb and flood of strong tidal currents. Extraction increases the resistance to flow, leading to changes in tidal range, transport, mixing, and the kinetic resource itself. These far-field changes have environmental, social, and economic implications that must be understood to develop the in-stream resource. This dissertation describes the development of a one-dimensional numerical channel model and its application to the study of these effects. The model is applied to determine the roles played by site geometry, network topology, tidal regime, and device dynamics. A comparison is also made between theoretical and modeled predictions for the maximum amount of power which could be extracted from a tidal energy site. The model is extended to a simulation of kinetic power extraction from Puget Sound, Washington. In general, extracting tidal energy will have a number of far-field effects, in proportion to the level of power extraction. At the theoretical limit, these effects can be very significant (e.g., 50% reduction in transport), but are predicted to be immeasurably small for pilot-scale projects. Depending on the specifics of the site, far-field effects may either augment or reduce the existing tidal regime. Changes to the tide, in particular, have significant spatial variability. Since tidal streams are generally subcritical, effects are felt throughout the estuary, not just at the site of extraction. The one dimensional numerical modeling is supported by a robust theory for predicting the performance characteristics of in-stream devices. The far-field effects of tidal power depend on the total power dissipated by turbines, rather than the power extracted. When the low-speed wake downstream of a turbine mixes with the free-stream, power is lost, such that the total power dissipated by the turbine is significantly greater than the power extracted. This dissertation concludes with a framework for three-dimensional numerical modeling of near-field extraction effects.

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

Ginn, Timothy R.; Weathers, Tess

Biogeochemical modeling using PHREEQC2 and a streamtube ensemble approach is utilized to understand a well-to-well subsurface treatment system at the Vadose Zone Research Park (VZRP) near Idaho Falls, Idaho. Treatment involves in situ microbially-mediated ureolysis to induce calcite precipitation for the immobilization of strontium-90. PHREEQC2 is utilized to model the kinetically-controlled ureolysis and consequent calcite precipitation. Reaction kinetics, equilibrium phases, and cation exchange are used within PHREEQC2 to track pH and levels of calcium, ammonium, urea, and calcite precipitation over time, within a series of one-dimensional advective-dispersive transport paths creating a streamtube ensemble representation of the well-to-well transport. An understandingmore » of the impact of physical heterogeneities within this radial flowfield is critical for remediation design; we address this via the streamtube approach: instead of depicting spatial extents of solutes in the subsurface we focus on their arrival distribution at the control well(s). Traditionally, each streamtube maintains uniform velocity; however in radial flow in homogeneous media, the velocity within any given streamtube is spatially-variable in a common way, being highest at the input and output wells and approaching a minimum at the midpoint between the wells. This idealized velocity variability is of significance in the case of ureolytically driven calcite precipitation. Streamtube velocity patterns for any particular configuration of injection and withdrawal wells are available as explicit calculations from potential theory, and also from particle tracking programs. To approximate the actual spatial distribution of velocity along streamtubes, we assume idealized radial non-uniform velocity associated with homogeneous media. This is implemented in PHREEQC2 via a non-uniform spatial discretization within each streamtube that honors both the streamtube’s travel time and the idealized “fast-slow-fast” pattern of non-uniform velocity along the streamline. Breakthrough curves produced by each simulation are weighted by the path-respective flux fractions (obtained by deconvolution of tracer tests conducted at the VZRP) to obtain the flux-average of flow contributions to the observation well.« less

Joint reconstruction of dynamic PET activity and kinetic parametric images using total variation constrained dictionary sparse coding

NASA Astrophysics Data System (ADS)

Yu, Haiqing; Chen, Shuhang; Chen, Yunmei; Liu, Huafeng

2017-05-01

Dynamic positron emission tomography (PET) is capable of providing both spatial and temporal information of radio tracers in vivo. In this paper, we present a novel joint estimation framework to reconstruct temporal sequences of dynamic PET images and the coefficients characterizing the system impulse response function, from which the associated parametric images of the system macro parameters for tracer kinetics can be estimated. The proposed algorithm, which combines statistical data measurement and tracer kinetic models, integrates a dictionary sparse coding (DSC) into a total variational minimization based algorithm for simultaneous reconstruction of the activity distribution and parametric map from measured emission sinograms. DSC, based on the compartmental theory, provides biologically meaningful regularization, and total variation regularization is incorporated to provide edge-preserving guidance. We rely on techniques from minimization algorithms (the alternating direction method of multipliers) to first generate the estimated activity distributions with sub-optimal kinetic parameter estimates, and then recover the parametric maps given these activity estimates. These coupled iterative steps are repeated as necessary until convergence. Experiments with synthetic, Monte Carlo generated data, and real patient data have been conducted, and the results are very promising.

Turbulent transport across an interface between dry and humid air in a stratified environment

NASA Astrophysics Data System (ADS)

Gallana, Luca; de Santi, Francesca; di Savino, Silvio; Iovieno, Michele; Ricchiardone, Renzo; Tordella, Daniela

2014-11-01

The transport of energy and water vapor across a thin layer which separates two decaying isotropic turbulent flows with different kinetic energy and humidity is considered. The interface is placed in a shearless stratified environment in temporal decay. This system reproduces a few aspects of small scale turbulent transport across a dry air/moist air interface in an atmospheric like context. In our incompressible DNS at Reλ = 250 , Boussinesq's approximation is used for momentum and energy transport while the vapor is modeled as a passive scalar (Kumar, Schumacher & Shaw 2014). We investigated different stratification levels with an initial Fr between 0.8 and 8 in presence of a kinetic energy ratio equal to 7. As the buoyancy term becomes of the same order of the inertial ones, a spatial redistribution of kinetic energy, dissipation and vapor concentration is observed. This eventually leads to the onset of a well of kinetic energy in the low energy side of the mixing layer which blocks the entrainment of dry air. Results are discussed and compared with laboratory and numerical experiments. A posteriori estimates of the eventual compression/expansion of fluid particles inside the interfacial mixing layer are given (Nance & Durran 1994).

Mesoscale Thermodynamically motivated Statistical Mechanics based Kinetic Model for Sintering monoliths

NASA Astrophysics Data System (ADS)

Mohan, Nisha

Modeling the evolution of microstructure during sintering is a persistent challenge in ceramics science, although needed as the microstructure impacts properties of an engineered material. Bridging the gap between microscopic and continuum models, kinetic Monte Carlo (kMC) methods provide a stochastic approach towards sintering and microstructure evolution. These kMC models work at the mesoscale, with length and time-scales between those of atomistic and continuum approaches. We develop a sintering/compacting model for the two-phase sintering of boron nitride ceramics and allotropes alike. Our formulation includes mechanisms for phase transformation between h-BN and c-BN and takes into account thermodynamics of pressure and temperature on interaction energies and mechanism rates. In addition to replicating the micro-structure evolution observed in experiments, it also captures the phase diagram of Boron Nitride materials. Results have been analyzed in terms of phase diagrams and crystal growth. It also serves with insights to guide the choice of additives and conditions for the sintering process.While detailed time and spatial resolutions are lost in any MC, the progression of stochastic events still captures plausible local energy minima and long-time temporal developments. DARPA.

Imaging burst kinetics and spatial coordination during serial killing by single natural killer cells

PubMed Central

Choi, Paul J.; Mitchison, Timothy J.

2013-01-01

Cytotoxic lymphocytes eliminate virus-infected and cancerous cells by immune recognition and killing through the perforin-granzyme pathway. Traditional killing assays measure average target cell lysis at fixed times and high effector:target ratios. Such assays obscure kinetic details that might reveal novel physiology. We engineered target cells to report on granzyme activity, used very low effector:target ratios to observe potential serial killing, and performed low magnification time-lapse imaging to reveal time-dependent statistics of natural killer (NK) killing at the single-cell level. Most kills occurred during serial killing, and a single NK cell killed up to 10 targets over a 6-h assay. The first kill was slower than subsequent kills, especially on poor targets, or when NK signaling pathways were partially inhibited. Spatial analysis showed that sequential kills were usually adjacent. We propose that NK cells integrate signals from the previous and current target, possibly by simultaneous contact. The resulting burst kinetics and spatial coordination may control the activity of NK cells in tissues. PMID:23576740

Adaptive Nodal Transport Methods for Reactor Transient Analysis

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

Thomas Downar; E. Lewis

2005-08-31

Develop methods for adaptively treating the angular, spatial, and time dependence of the neutron flux in reactor transient analysis. These methods were demonstrated in the DOE transport nodal code VARIANT and the US NRC spatial kinetics code, PARCS.

Interwell Connectivity Evaluation Using Injection and Production Fluctuation Data

NASA Astrophysics Data System (ADS)

Shang, Barry Zhongqi

The development of multiscale methods for computational simulation of biophysical systems represents a significant challenge. Effective computational models that bridge physical insights obtained from atomistic simulations and experimental findings are lacking. An accurate passing of information between these scales would enable: (1) an improved physical understanding of structure-function relationships, and (2) enhanced rational strategies for molecular engineering and materials design. Two approaches are described in this dissertation to facilitate these multiscale goals. In Part I, we develop a lattice kinetic Monte Carlo model to simulate cellulose decomposition by cellulase enzymes and to understand the effects of spatial confinement on enzyme kinetics. An enhanced mechanistic understanding of this reaction system could enhance the design of cellulose bioconversion technologies for renewable and sustainable energy. Using our model, we simulate the reaction up to experimental conversion times of days, while simultaneously capturing the microscopic kinetic behaviors. Therefore, the influence of molecular-scale kinetics on the macroscopic conversion rate is made transparent. The inclusion of spatial constraints in the kinetic model represents a significant advance over classical mass-action models commonly used to describe this reaction system. We find that restrictions due to enzyme jamming and substrate heterogeneity at the molecular level play a dominate role in limiting cellulose conversion. We identify that the key rate limitations are the slow rates of enzyme complexation with glucan chains and the competition between enzyme processivity and jamming. We show that the kinetics of complexation, which involves extraction of a glucan chain end from the cellulose surface and threading through the enzyme active site, occurs slowly on the order of hours, while intrinsic hydrolytic bond cleavage occurs on the order of seconds. We also elucidate the subtle trade-off between processivity and jamming. Highly processive enzymes cleave a large fraction of a glucan chain during each processive run but are prone to jamming at obstacles. Less processive enzymes avoid jamming but cleave only a small fraction of a chain. Optimizing this trade-off maximizes the cellulose conversion rate. We also elucidate the molecular-scale kinetic origins for synergy among cellulases in enzyme mixtures. In contrast to the currently accepted theory, we show that the ability of an endoglucanase to increase the concentration of chain ends for exoglucanases is insufficient for synergy to occur. Rather, endoglucanases must enhance the rate of complexation between exoglucanases and the newly created chain ends. This enhancement occurs when the endoglucanase is able to partially decrystallize the cellulose surface. We show generally that the driving forces for complexation and jamming, which govern the kinetics of pure exoglucanases, also control the degree of synergy in endo-exo mixtures. In Part II, we focus our attention on a different multiscale problem. This challenge is the development of coarse-grained models from atomistic models to access larger length- and time-scales in a simulation. This problem is difficult because it requires a delicate balance between maintaining (1) physical simplicity in the coarse-grained model and (2) physical consistency with the atomistic model. To achieve these goals, we develop a scheme to coarse-grain an atomistic fluid model into a fluctuating hydrodynamics (FHD) model. The FHD model describes the solvent as a field of fluctuating mass, momentum, and energy densities. The dynamics of the fluid are governed by continuum balance equations and fluctuation-dissipation relations based on the constitutive transport laws. The incorporation of both macroscopic transport and microscopic fluctuation phenomena could provide richer physical insight into the behaviors of biophysical systems driven by hydrodynamic fluctuations, such as hydrophobic assembly and crystal nucleation. We further extend our coarse-graining method by developing an interfacial FHD model using information obtained from simulations of an atomistic liquid-vapor interface. We illustrate that a phenomenological Ginzburg-Landau free energy employed in the FHD model can effectively represent the attractive molecular interactions of the atomistic model, which give rise to phase separation. For argon and water, we show that the interfacial FHD model can reproduce the compressibility, surface tension, and capillary wave spectrum of the atomistic model. Via this approach, simulations that explore the coupling between hydrodynamic fluctuations and phase equilibria with molecular-scale consistency are now possible. In both Parts I and II, the emerging theme is that the combination of bottom-up coarse graining and top-down phenomenology is essential for enabling a multiscale approach to remain physically consistent with molecular-scale interactions while simultaneously capturing the collective macroscopic behaviors. This hybrid strategy enables the resulting computational models to be both physically insightful and practically meaningful. (Abstract shortened by UMI.).

Discrete-State Stochastic Models of Calcium-Regulated Calcium Influx and Subspace Dynamics Are Not Well-Approximated by ODEs That Neglect Concentration Fluctuations

PubMed Central

Weinberg, Seth H.; Smith, Gregory D.

2012-01-01

Cardiac myocyte calcium signaling is often modeled using deterministic ordinary differential equations (ODEs) and mass-action kinetics. However, spatially restricted “domains” associated with calcium influx are small enough (e.g., 10−17 liters) that local signaling may involve 1–100 calcium ions. Is it appropriate to model the dynamics of subspace calcium using deterministic ODEs or, alternatively, do we require stochastic descriptions that account for the fundamentally discrete nature of these local calcium signals? To address this question, we constructed a minimal Markov model of a calcium-regulated calcium channel and associated subspace. We compared the expected value of fluctuating subspace calcium concentration (a result that accounts for the small subspace volume) with the corresponding deterministic model (an approximation that assumes large system size). When subspace calcium did not regulate calcium influx, the deterministic and stochastic descriptions agreed. However, when calcium binding altered channel activity in the model, the continuous deterministic description often deviated significantly from the discrete stochastic model, unless the subspace volume is unrealistically large and/or the kinetics of the calcium binding are sufficiently fast. This principle was also demonstrated using a physiologically realistic model of calmodulin regulation of L-type calcium channels introduced by Yue and coworkers. PMID:23509597

The effect of the electric wind on the spatial distribution of chemical species in the positive corona discharge

NASA Astrophysics Data System (ADS)

Yanallah, K.; Pontiga, F.; Bouazza, M. R.; Chen, J. H.

2017-08-01

The electrohydrodynamic air flow generated by a positive corona discharge, and its effect on the spatial distribution of chemical species within a wire-plate corona reactor, have been numerically simulated. The computational model is based on the solutions of the Navier-Stokes equation and the continuity equation of each chemical species generated by the electrical discharge. A simplified analytical expression of the electric force density, which only requires the current density as the input parameter, has been used in the Navier-Stokes equation to obtain the velocity field. For the solution of the continuity equations, a plasma chemistry model that includes the most important reactions between electrons, atoms and molecules in air has been used. Similar to the electric force, the electron density distribution has been approximated by using a semi-analytical expression appropriate for the electrode geometry. The results of the study show that the spatial distribution of chemical species can be very different, and depends on the interplay between the electrohydrodynamic flow, the chemical kinetics of the species and its characteristic lifetime.

Mars-solar wind interaction: LatHyS, an improved parallel 3-D multispecies hybrid model

NASA Astrophysics Data System (ADS)

Modolo, Ronan; Hess, Sebastien; Mancini, Marco; Leblanc, Francois; Chaufray, Jean-Yves; Brain, David; Leclercq, Ludivine; Esteban-Hernández, Rosa; Chanteur, Gerard; Weill, Philippe; González-Galindo, Francisco; Forget, Francois; Yagi, Manabu; Mazelle, Christian

2016-07-01

In order to better represent Mars-solar wind interaction, we present an unprecedented model achieving spatial resolution down to 50 km, a so far unexplored resolution for global kinetic models of the Martian ionized environment. Such resolution approaches the ionospheric plasma scale height. In practice, the model is derived from a first version described in Modolo et al. (2005). An important effort of parallelization has been conducted and is presented here. A better description of the ionosphere was also implemented including ionospheric chemistry, electrical conductivities, and a drag force modeling the ion-neutral collisions in the ionosphere. This new version of the code, named LatHyS (Latmos Hybrid Simulation), is here used to characterize the impact of various spatial resolutions on simulation results. In addition, and following a global model challenge effort, we present the results of simulation run for three cases which allow addressing the effect of the suprathermal corona and of the solar EUV activity on the magnetospheric plasma boundaries and on the global escape. Simulation results showed that global patterns are relatively similar for the different spatial resolution runs, but finest grid runs provide a better representation of the ionosphere and display more details of the planetary plasma dynamic. Simulation results suggest that a significant fraction of escaping O+ ions is originated from below 1200 km altitude.

Formation of Polar Stratospheric Clouds in the Atmosphere

NASA Astrophysics Data System (ADS)

Aloyan, Artash; Yermakov, Alex; Arutyunyan, Vardan; Larin, Igor

2014-05-01

A new mathematical model of the global transport of gaseous species and aerosols in the atmosphere and the formation of polar stratospheric clouds (PSCs) in both hemispheres was constructed. PSCs play a significant role in ozone chemistry since heterogeneous reactions proceed on their particle surfaces and in the bulk, affecting the gas composition of the atmosphere, specifically, the content of chlorine and nitrogen compounds, which are actively involved in the destruction of ozone. Stratospheric clouds are generated by co-condensation of water vapor and nitric acid on sulfate particles and in some cases during the freezing of supercooled water as well as when nitric acid vapors are dissolved in sulfate aerosol particles [1]. These clouds differ in their chemical composition and microphysics [2]. In this study, we propose new kinetic equations describing the variability of species in the gas and condensed phases to simulate the formation of PSCs. Most models for the formation of PSCs use constant background values of sulfate aerosols in the lower stratosphere. This approach is too simplistic since sulfate aerosols in the stratosphere are characterized by considerably nonuniform spatial and temporal variations. Two PSC types are considered: Type 1 refers to the formation of nitric acid trihydrate (NAT) and Type 2 refers to the formation of particles composed of different proportions of H2SO4/HNO3/H2O. Their formation is coupled with the spatial problem of sulfate aerosol generation in the upper troposphere and lower stratosphere incorporating the chemical and kinetic transformation processes (photochemistry, nucleation, condensation/evaporation, and coagulation) and using a non-equilibrium particle-size distribution [3]. In this formulation, the system of equations is closed and allows an adequate description of the PSC dynamics in the stratosphere. Using the model developed, numerical experiments were performed to reproduce the spatial and temporal variability of polar clouds in both hemispheres for the winter time period. The numerical experiments were performed in the following sequence. In the first stage, we address the transport of multicomponent gaseous species, the formation of sulfate aerosols in the troposphere and lower stratosphere (spherical atmosphere), the chemical and kinetic transformations, and the biogenic and anthropogenic emissions of related chemical components [3]. This model makes it possible to reproduce the distribution of sulfate particles in the size range from 3 nm to 1 mcm. Next, the base model was improved by using a new module describing the dynamics of phase transition of substances in gaseous and condensed phases that are typical for different types of PSCs. Here, we used the methods of thermodynamics. Conclusions •The model developed allow us to reproduce the size distribution of sulfate particles generated from precursor gases in the troposphere and stratosphere; •The numerical experiments show that the model adequately reproduces the spatial characteristics of the PSC formation in the atmosphere. References 1.Carslaw K.S., Peter T., Clegg S.L. Modeling the composition of liquid stratospheric clouds. Rev. Geophys. 35, 125, 1997 2.Drdla, K., Shoeberl, M.R., and Browell, E.V., Microphysical modeling of the 1999-2000 Arctic winter. J. Geophys. Res., 2003, vol. 108, No. D5, p. 8312. 3.Aloyan, A.E., Yermakov, A.N., Arutyunyan, V.O., Sulfate aerosol formation in the troposphere and lower stratosphere, in Possibilities of Climate Stabilization by Using Novel Technologies, Moscow: Rosgidromet, 2012, pp. 75-98.

The spatial-temporal characteristics of type I collagen-based extracellular matrix.

PubMed

Jones, Christopher Allen Rucksack; Liang, Long; Lin, Daniel; Jiao, Yang; Sun, Bo

2014-11-28

Type I collagen abounds in mammalian extracellular matrix (ECM) and is crucial to many biophysical processes. While previous studies have mostly focused on bulk averaged properties, here we provide a comprehensive and quantitative spatial-temporal characterization of the microstructure of type I collagen-based ECM as the gelation temperature varies. The structural characteristics including the density and nematic correlation functions are obtained by analyzing confocal images of collagen gels prepared at a wide range of gelation temperatures (from 16 °C to 36 °C). As temperature increases, the gel microstructure varies from a "bundled" network with strong orientational correlation between the fibers to an isotropic homogeneous network with no significant orientational correlation, as manifested by the decaying of length scales in the correlation functions. We develop a kinetic Monte-Carlo collagen growth model to better understand how ECM microstructure depends on various environmental or kinetic factors. We show that the nucleation rate, growth rate, and an effective hydrodynamic alignment of collagen fibers fully determines the spatiotemporal fluctuations of the density and orientational order of collagen gel microstructure. Also the temperature dependence of the growth rate and nucleation rate follow the prediction of classical nucleation theory.

Magnetohydrodynamics with Embedded Particle-in-Cell Simulation of Mercury's Magnetosphere

NASA Astrophysics Data System (ADS)

Chen, Y.; Toth, G.; Jia, X.; Gombosi, T. I.; Markidis, S.

2015-12-01

Mercury's magnetosphere is much more dynamic than other planetary magnetospheres because of Mercury's weak intrinsic magnetic field and its proximity to the Sun. Magnetic reconnection and Kelvin-Helmholtz phenomena occur in Mercury's magnetopause and magnetotail at higher frequencies than in other planetary magnetosphere. For instance, chains of flux transfer events (FTEs) on the magnetopause, have been frequentlyobserved by the the MErcury Surface, Space ENvironment, GEochemistry and Ranging (MESSENGER) spacecraft (Slavin et al., 2012). Because ion Larmor radius is comparable to typical spatial scales in Mercury's magnetosphere, finite Larmor radius effects need to be accounted for. In addition, it is important to take in account non-ideal dissipation mechanisms to accurately describe magnetic reconnection. A kinetic approach allows us to model these phenomena accurately. However, kinetic global simulations, even for small-size magnetospheres like Mercury's, are currently unfeasible because of the high computational cost. In this work, we carry out global simulations of Mercury's magnetosphere with the recently developed MHD-EPIC model, which is a two-way coupling of the extended magnetohydrodynamic (XMHD) code BATS-R-US with the implicit Particle-in-Cell (PIC) model iPIC3D. The PIC model can cover the regions where kinetic effects are most important, such as reconnection sites. The BATS-R-US code, on the other hand, can efficiently handle the rest of the computational domain where the MHD or Hall MHD description is sufficient. We will present our preliminary results and comparison with MESSENGER observations.

Epidermal Homeostasis and Radiation Responses in a Multiscale Tissue Modeling Framework

NASA Technical Reports Server (NTRS)

Hu, Shaowen; Cucinotta, Francis A.

2013-01-01

The surface of skin is lined with several thin layers of epithelial cells that are maintained throughout life time by a small population of stem cells. High dose radiation exposures could injure and deplete the underlying proliferative cells and induce cutaneous radiation syndrome. In this work we propose a multiscale computational model for skin epidermal dynamics that links phenomena occurring at the subcellular, cellular, and tissue levels of organization, to simulate the experimental data of the radiation response of swine epidermis, which is closely similar to human epidermis. Incorporating experimentally measured histological and cell kinetic parameters, we obtain results of population kinetics and proliferation indexes comparable to observations in unirradiated and acutely irradiated swine experiments. At the sub-cellular level, several recently published Wnt signaling controlled cell-cycle models are applied and the roles of key components and parameters are analyzed. Based on our simulation results, we demonstrate that a moderate increase of proliferation rate for the survival proliferative cells is sufficient to fully repopulate the area denuded by high dose radiation, as long as the integrity of underlying basement membrane is maintained. Our work highlights the importance of considering proliferation kinetics as well as the spatial organization of tissues when conducting in vivo investigations of radiation responses. This integrated model allow us to test the validity of several basic biological rules at the cellular level and sub-cellular mechanisms by qualitatively comparing simulation results with published research, and enhance our understanding of the pathophysiological effects of ionizing radiation on skin.

Trivial constraints on orbital-free kinetic energy density functionals

NASA Astrophysics Data System (ADS)

Luo, Kai; Trickey, S. B.

2018-03-01

Approximate kinetic energy density functionals (KEDFs) are central to orbital-free density functional theory. Limitations on the spatial derivative dependencies of KEDFs have been claimed from differential virial theorems. We identify a central defect in the argument: the relationships are not true for an arbitrary density but hold only for the minimizing density and corresponding chemical potential. Contrary to the claims therefore, the relationships are not constraints and provide no independent information about the spatial derivative dependencies of approximate KEDFs. A simple argument also shows that validity for arbitrary v-representable densities is not restored by appeal to the density-potential bijection.

Energy and water vapor transport across a simplified cloud-clear air interface

NASA Astrophysics Data System (ADS)

Gallana, L.; Di Savino, S.; De Santi, F.; Iovieno, M.; Tordella, D.

2014-11-01

We consider a simplified physics of the could interface where condensation, evaporation and radiation are neglected and momentum, thermal energy and water vapor transport is represented in terms of the Boussinesq model coupled to a passive scalar transport equation for the vapor. The interface is modeled as a layer separating two isotropic turbulent regions with different kinetic energy and vapor concentration. In particular, we focus on the small scale part of the inertial range of the atmospheric boundary layer as well as on the dissipative range of scales which are important to the micro-physics of warm clouds. We have numerically investigated stably stratified interfaces by locally perturbing at an initial instant the standard temperature lapse rate at the cloud interface and then observing the temporal evolution of the system. When the buoyancy term becomes of the same order of the inertial one, we observe a spatial redistribution of the kinetic energy which produce a concomitant pit of kinetic energy within the mixing layer. In this situation, the mixing layer contains two interfacial regions with opposite kinetic energy gradient, which in turn produces two intermittent sublayers in the velocity fluctuations field. This changes the structure of the field with respect to the corresponding non-stratified shearless mixing: the communication between the two turbulent region is weak, and the growth of the mixing layer stops. These results are discussed with respect to Large Eddy Simulations data for the Planetary Boundary Layers.

Global Particle-in-Cell Simulations of Mercury's Magnetosphere

NASA Astrophysics Data System (ADS)

Schriver, D.; Travnicek, P. M.; Lapenta, G.; Amaya, J.; Gonzalez, D.; Richard, R. L.; Berchem, J.; Hellinger, P.

2017-12-01

Spacecraft observations of Mercury's magnetosphere have shown that kinetic ion and electron particle effects play a major role in the transport, acceleration, and loss of plasma within the magnetospheric system. Kinetic processes include reconnection, the breakdown of particle adiabaticity and wave-particle interactions. Because of the vast range in spatial scales involved in magnetospheric dynamics, from local electron Debye length scales ( meters) to solar wind/planetary magnetic scale lengths (tens to hundreds of planetary radii), fully self-consistent kinetic simulations of a global planetary magnetosphere remain challenging. Most global simulations of Earth's and other planet's magnetosphere are carried out using MHD, enhanced MHD (e.g., Hall MHD), hybrid, or a combination of MHD and particle in cell (PIC) simulations. Here, 3D kinetic self-consistent hybrid (ion particle, electron fluid) and full PIC (ion and electron particle) simulations of the solar wind interaction with Mercury's magnetosphere are carried out. Using the implicit PIC and hybrid simulations, Mercury's relatively small, but highly kinetic magnetosphere will be examined to determine how the self-consistent inclusion of electrons affects magnetic reconnection, particle transport and acceleration of plasma at Mercury. Also the spatial and energy profiles of precipitating magnetospheric ions and electrons onto Mercury's surface, which can strongly affect the regolith in terms of space weathering and particle outflow, will be examined with the PIC and hybrid codes. MESSENGER spacecraft observations are used both to initiate and validate the global kinetic simulations to achieve a deeper understanding of the role kinetic physics play in magnetospheric dynamics.

Measurement of turbulent spatial structure and kinetic energy spectrum by exact temporal-to-spatial mapping

NASA Astrophysics Data System (ADS)

Buchhave, Preben; Velte, Clara M.

2017-08-01

We present a method for converting a time record of turbulent velocity measured at a point in a flow to a spatial velocity record consisting of consecutive convection elements. The spatial record allows computation of dynamic statistical moments such as turbulent kinetic wavenumber spectra and spatial structure functions in a way that completely bypasses the need for Taylor's hypothesis. The spatial statistics agree with the classical counterparts, such as the total kinetic energy spectrum, at least for spatial extents up to the Taylor microscale. The requirements for applying the method are access to the instantaneous velocity magnitude, in addition to the desired flow quantity, and a high temporal resolution in comparison to the relevant time scales of the flow. We map, without distortion and bias, notoriously difficult developing turbulent high intensity flows using three main aspects that distinguish these measurements from previous work in the field: (1) The measurements are conducted using laser Doppler anemometry and are therefore not contaminated by directional ambiguity (in contrast to, e.g., frequently employed hot-wire anemometers); (2) the measurement data are extracted using a correctly and transparently functioning processor and are analysed using methods derived from first principles to provide unbiased estimates of the velocity statistics; (3) the exact mapping proposed herein has been applied to the high turbulence intensity flows investigated to avoid the significant distortions caused by Taylor's hypothesis. The method is first confirmed to produce the correct statistics using computer simulations and later applied to measurements in some of the most difficult regions of a round turbulent jet—the non-equilibrium developing region and the outermost parts of the developed jet. The proposed mapping is successfully validated using corresponding directly measured spatial statistics in the fully developed jet, even in the difficult outer regions of the jet where the average convection velocity is negligible and turbulence intensities increase dramatically. The measurements in the developing region reveal interesting features of an incomplete Richardson-Kolmogorov cascade under development.

Predicting the response of the deep-ocean microbiome to geochemical perturbations by hydrothermal vents.

PubMed

Reed, Daniel C; Breier, John A; Jiang, Houshuo; Anantharaman, Karthik; Klausmeier, Christopher A; Toner, Brandy M; Hancock, Cathrine; Speer, Kevin; Thurnherr, Andreas M; Dick, Gregory J

2015-08-01

Submarine hydrothermal vents perturb the deep-ocean microbiome by injecting reduced chemical species into the water column that act as an energy source for chemosynthetic organisms. These systems thus provide excellent natural laboratories for studying the response of microbial communities to shifts in marine geochemistry. The present study explores the processes that regulate coupled microbial-geochemical dynamics in hydrothermal plumes by means of a novel mathematical model, which combines thermodynamics, growth and reaction kinetics, and transport processes derived from a fluid dynamics model. Simulations of a plume located in the ABE vent field of the Lau basin were able to reproduce metagenomic observations well and demonstrated that the magnitude of primary production and rate of autotrophic growth are largely regulated by the energetics of metabolisms and the availability of electron donors, as opposed to kinetic parameters. Ambient seawater was the dominant source of microbes to the plume and sulphur oxidisers constituted almost 90% of the modelled community in the neutrally-buoyant plume. Data from drifters deployed in the region allowed the different time scales of metabolisms to be cast in a spatial context, which demonstrated spatial succession in the microbial community. While growth was shown to occur over distances of tens of kilometers, microbes persisted over hundreds of kilometers. Given that high-temperature hydrothermal systems are found less than 100 km apart on average, plumes may act as important vectors between different vent fields and other environments that are hospitable to similar organisms, such as oil spills and oxygen minimum zones.

Predicting the response of the deep-ocean microbiome to geochemical perturbations by hydrothermal vents

PubMed Central

Reed, Daniel C; Breier, John A; Jiang, Houshuo; Anantharaman, Karthik; Klausmeier, Christopher A; Toner, Brandy M; Hancock, Cathrine; Speer, Kevin; Thurnherr, Andreas M; Dick, Gregory J

2015-01-01

Submarine hydrothermal vents perturb the deep-ocean microbiome by injecting reduced chemical species into the water column that act as an energy source for chemosynthetic organisms. These systems thus provide excellent natural laboratories for studying the response of microbial communities to shifts in marine geochemistry. The present study explores the processes that regulate coupled microbial-geochemical dynamics in hydrothermal plumes by means of a novel mathematical model, which combines thermodynamics, growth and reaction kinetics, and transport processes derived from a fluid dynamics model. Simulations of a plume located in the ABE vent field of the Lau basin were able to reproduce metagenomic observations well and demonstrated that the magnitude of primary production and rate of autotrophic growth are largely regulated by the energetics of metabolisms and the availability of electron donors, as opposed to kinetic parameters. Ambient seawater was the dominant source of microbes to the plume and sulphur oxidisers constituted almost 90% of the modelled community in the neutrally-buoyant plume. Data from drifters deployed in the region allowed the different time scales of metabolisms to be cast in a spatial context, which demonstrated spatial succession in the microbial community. While growth was shown to occur over distances of tens of kilometers, microbes persisted over hundreds of kilometers. Given that high-temperature hydrothermal systems are found less than 100 km apart on average, plumes may act as important vectors between different vent fields and other environments that are hospitable to similar organisms, such as oil spills and oxygen minimum zones. PMID:25658053

Turbulent kinetic energy of the ocean winds over the Kuroshio Extension from QuikSCAT winds (1999-2009)

NASA Astrophysics Data System (ADS)

Yu, Kai; Dong, Changming; King, Gregory P.

2017-06-01

We investigate mesoscale turbulence (10-1000 km) in the ocean winds over the Kuroshio Extension (28°N-40°N, 140°E-180°E) using the QuikSCAT data set (November 1999 to October 2009). We calculate the second (Djj) and third-order structure functions (Djjj) and the spatial variance (Vj) as a function of scale r (j=L,T denotes, respectively, the longitudinal (divergent) and transverse (vortical) component). The most interesting results of the analysis follow. Although both Vj>(r>) and Djj>(r>) measure the turbulent kinetic energy (TKE), we find that Vj>(r>) is the more robust measure. The spatial variance density (dVj/dr) has a broad peak near 450 km (close to the midlatitude Rossby radius of deformation). On interannual time scales, TKE correlates well with the El Niño 3.4 index. According to turbulence theory, the kinetic energy cascades downscale (upscale) if DLLL>(r>) (also skewness SL=DLLL/DLL3/2) is negative (positive). Our results for the Kuroshio Extension are consistent with a downscale cascade (indicating convergence dominates). Furthermore, classical turbulence theory predicts that SL=-0.3 and independent of r; however, we find SL varies strongly with r, from -4 at small scales to -0.3 at large scales. This nonclassical behavior implies strong-scale interaction, which we attribute to the rapid, and sometimes explosive, growth of storms in the region through baroclinic instability. Finally, we find that ST (a measure of cyclonic/anticyclonic asymmetry) is positive (cyclonic) and also varies strongly with r, from 4 at small scales to 0.5 at large scales. New turbulence models are needed to explain these results, and that will benefit Weather Prediction and climate modeling.