Science.gov

Sample records for kinetic equations

  1. Kinetic Equations for Economic Sciences

    NASA Astrophysics Data System (ADS)

    Bisi, M.; Brugna, C.

    2010-04-01

    We discuss, both from the analytical and the numerical point of view, a kinetic model for wealth distribution in a simple market economy which models, besides binary trade interactions, also taxation and redistribution of collected wealth.

  2. Boltzmann kinetic equation for filtered fluid turbulence.

    PubMed

    Girimaji, Sharath S

    2007-07-20

    We develop a kinetic Boltzmann equation for describing filtered fluid turbulence applicable for continuum and noncontinuum effects. The effect of unresolved turbulent motion on the resolved distribution function is elucidated and closure modeling issues of kinetic Boltzmann and Navier-Stokes descriptions are reconciled. This could pave the way for unifying turbulence modeling at kinetic and continuum levels and the development of numerical methods that are valid over a wide range of flow physics. PMID:17678288

  3. Turbulence kinetic energy equation for dilute suspensions

    NASA Technical Reports Server (NTRS)

    Abou-Arab, T. W.; Roco, M. C.

    1989-01-01

    A multiphase turbulence closure model is presented which employs one transport equation, namely the turbulence kinetic energy equation. The proposed form of this equation is different from the earlier formulations in some aspects. The power spectrum of the carrier fluid is divided into two regions, which interact in different ways and at different rates with the suspended particles as a function of the particle-eddy size ratio and density ratio. The length scale is described algebraically. A mass/time averaging procedure for the momentum and kinetic energy equations is adopted. The resulting turbulence correlations are modeled under less retrictive assumptions comparative to previous work. The closures for the momentum and kinetic energy equations are given. Comparisons of the predictions with experimental results on liquid-solid jet and gas-solid pipe flow show satisfactory agreement.

  4. Spectrum Analysis of Some Kinetic Equations

    NASA Astrophysics Data System (ADS)

    Yang, Tong; Yu, Hongjun

    2016-05-01

    We analyze the spectrum structure of some kinetic equations qualitatively by using semigroup theory and linear operator perturbation theory. The models include the classical Boltzmann equation for hard potentials with or without angular cutoff and the Landau equation with {γ≥q-2} . As an application, we show that the solutions to these two fundamental equations are asymptotically equivalent (mod time decay rate {t^{-5/4}} ) as {tto∞} to that of the compressible Navier-Stokes equations for initial data around an equilibrium state.

  5. Stochastic thermodynamics for linear kinetic equations

    NASA Astrophysics Data System (ADS)

    Van den Broeck, C.; Toral, R.

    2015-07-01

    Stochastic thermodynamics is formulated for variables that are odd under time reversal. The invariance under spatial rotation of the collision rates due to the isotropy of the heat bath is shown to be a crucial ingredient. An alternative detailed fluctuation theorem is derived, expressed solely in terms of forward statistics. It is illustrated for a linear kinetic equation with kangaroo rates.

  6. On a Kinetic Equation for Coalescing Particles

    NASA Astrophysics Data System (ADS)

    Escobedo, Miguel; Laurençot, Philippe; Mischler, Stéphane

    Existence of global weak solutions to a spatially inhomogeneous kinetic model for coalescing particles is proved, each particle being identified by its mass, momentum and position. The large time convergence to zero is also shown. The cornestone of our analysis is that, for any nonnegative and convex function, the associated Orlicz norm is a Liapunov functional. Existence and asymptotic behaviour then rely on weak and strong compactness methods in L1 in the spirit of the DiPerna-Lions theory for the Boltzmann equation.

  7. Benchmarks for the point kinetics equations

    SciTech Connect

    Ganapol, B.; Picca, P.; Previti, A.; Mostacci, D.

    2013-07-01

    A new numerical algorithm is presented for the solution to the point kinetics equations (PKEs), whose accurate solution has been sought for over 60 years. The method couples the simplest of finite difference methods, a backward Euler, with Richardsons extrapolation, also called an acceleration. From this coupling, a series of benchmarks have emerged. These include cases from the literature as well as several new ones. The novelty of this presentation lies in the breadth of reactivity insertions considered, covering both prescribed and feedback reactivities, and the extreme 8- to 9- digit accuracy achievable. The benchmarks presented are to provide guidance to those who wish to develop further numerical improvements. (authors)

  8. Hypocoercivity of linear degenerately dissipative kinetic equations

    NASA Astrophysics Data System (ADS)

    Duan, Renjun

    2011-08-01

    In this paper we develop a general approach of studying the hypocoercivity for a class of linear kinetic equations with both transport and degenerately dissipative terms. As concrete examples, the relaxation operator, Fokker-Planck operator and linearized Boltzmann operator are considered when the spatial domain takes the whole space or torus and when there is a confining force or not. The key part of the developed approach is to construct some equivalent temporal energy functionals for obtaining time rates of the solution trending towards equilibrium in some Hilbert spaces. The result in the case of the linear Boltzmann equation with confining forces is new. The proof mainly makes use of the macro-micro decomposition combined with Kawashima's argument on dissipation of the hyperbolic-parabolic system. At the end, a Korn-type inequality with probability measure is provided to deal with dissipation of momentum components.

  9. Test of Relativistic Kinetic Energy Equation

    NASA Astrophysics Data System (ADS)

    Chaudhary, Bharat

    2014-03-01

    Kinetic energy of a body equals the work done on it by a force, constant or variable. Force is the time rate of change of momentum. Momentum is mass times velocity. According to special relativity mass and velocity both are variables. Therefore, the differentiation of their product (momentum) has two terms, both are variables. One term is the product of mass and acceleration. The other is of velocity and the rate of change of mass. They together equal the applied force. Since the force equals the sum of two variable terms, it also becomes a variable even if it was a constant earlier. Therefore it is a flaw. There are two more flaws in the force equation. They are found by putting the force equal to zero. When this is done, the acceleration doesn't become zero. This is physically incompatible and is therefore a flaw. The other flaw in the equation is found by integrating the right side terms and evaluating the constant of integration from the initial conditions. Then we get a term containing logarithm of zero that is undefined, therefore the expression so obtained is meaningless. Since it comes from the relativistic definition of force, therefore we conclude that this definition is wrong. Thus we find that there are three flaws in the relativistic definition of force. They all make the relativistic equation of force wrong.

  10. Neutrino quantum kinetic equations: The collision term

    DOE PAGESBeta

    Blaschke, Daniel N.; Cirigliano, Vincenzo

    2016-08-25

    We derive the collision term relevant for neutrino quantum kinetic equations in the early universe and compact astrophysical objects, displaying its full matrix structure in both flavor and spin degrees of freedom. We include in our analysis neutrino-neutrino processes, scattering and annihilation with electrons and positrons, and neutrino scattering off nucleons (the latter in the low-density limit). After presenting the general structure of the collision terms, we take two instructive limiting cases. The one-flavor limit highlights the structure in helicity space and allows for a straightforward interpretation of the off-diagonal entries in terms of the product of scattering amplitudes ofmore » the two helicity states. As a result, the isotropic limit is relevant for studies of the early universe: in this case the terms involving spin coherence vanish and the collision term can be expressed in terms of two-dimensional integrals, suitable for computational implementation.« less

  11. Kinetic equation for spin-polarized plasmas

    SciTech Connect

    Cowley, S.C.; Kulsrud, R.M.; Valeo, E.

    1984-07-01

    The usual kinetic description of a plasma is extended to include variables to describe the spin. The distribution function, over phase-space and the new spin variables, provides a sufficient description of a spin-polarized plasma. The evolution equation for the distribution function is given. The equations derived are used to calculate depolarization due to four processes, inhomogeneous fields, collisions, collisions in inhomogeneous fields, and waves. It is found that depolarization by field inhomogeneity on scales large compared with the gyroradius is totally negligible. The same is true for collisional depolarization. Collisions in inhomogeneous fields yield a depolarization rate of order 10/sup -4/S/sup -1/ for deuterons and a negligible rate for tritons in a typical fusion reactor design. This is still sufficiently small on reactor time scales. However, small amplitude magnetic fluctuations (of order one gauss) resonant with the spin precession frequency can lead to significant depolarization (depolarises triton in ten seconds and deuteron in a hundred seconds.)

  12. Moment equations for chromatography based on Langmuir type reaction kinetics.

    PubMed

    Miyabe, Kanji

    2014-08-22

    Moment equations were derived for chromatography, in which the reaction kinetics between solute molecules and functional ligands on the stationary phase was represented by the Langmuir type rate equation. A set of basic equations of the general rate model of chromatography representing the mass balance, mass transfer rate, and reaction kinetics in the column were analytically solved in the Laplace domain. The moment equations for the first absolute moment and the second central moment in the real time domain were derived from the analytical solution in the Laplace domain. The moment equations were used for predicting the chromatographic behavior under hypothetical HPLC conditions. The influence of the parameters relating to the adsorption equilibrium and to the reaction kinetics on the chromatographic behavior was quantitatively evaluated. It is expected that the moment equations are effective for a detailed analysis of the influence of the mass transfer rates and of the Langmuir type reaction kinetics on the column efficiency. PMID:24999066

  13. Drift-free kinetic equations for turbulent dispersion.

    PubMed

    Bragg, A; Swailes, D C; Skartlien, R

    2012-11-01

    The dispersion of passive scalars and inertial particles in a turbulent flow can be described in terms of probability density functions (PDFs) defining the statistical distribution of relevant scalar or particle variables. The construction of transport equations governing the evolution of such PDFs has been the subject of numerous studies, and various authors have presented formulations for this type of equation, usually referred to as a kinetic equation. In the literature it is often stated, and widely assumed, that these PDF kinetic equation formulations are equivalent. In this paper it is shown that this is not the case, and the significance of differences among the various forms is considered. In particular, consideration is given to which form of equation is most appropriate for modeling dispersion in inhomogeneous turbulence and most consistent with the underlying particle equation of motion. In this regard the PDF equations for inertial particles are considered in the limit of zero particle Stokes number and assessed against the fully mixed (zero-drift) condition for fluid points. A long-standing question regarding the validity of kinetic equations in the fluid-point limit is answered; it is demonstrated formally that one version of the kinetic equation (derived using the Furutsu-Novikov method) provides a model that satisfies this zero-drift condition exactly in both homogeneous and inhomogeneous systems. In contrast, other forms of the kinetic equation do not satisfy this limit or apply only in a limited regime. PMID:23214875

  14. Calculation of coupling coefficients for equations of multipoint kinetics

    NASA Astrophysics Data System (ADS)

    Ioannisian, M. V.

    2013-12-01

    The multipoint kinetics equations for fission reaction rate are developed. The algorithm for computation of coupling coefficients is implemented within the MCU-5 code. Results from approbation of the method using the model problem and experimental data are presented.

  15. Fractional kinetic equation for Hamiltonian chaos

    NASA Astrophysics Data System (ADS)

    Zaslavsky, G. M.

    1994-09-01

    Hamiltonian chaotic dynamics of particles (or passive particles in fluids) can be described by a fractional generalization of the Fokker-Planck-Kolmogorov equation (FFPK) which is defined by two fractional critical exponents (α, β) responsible for the space and time derivatives of the distribution function correspondingly. A renormalization method has been proposed to determine (α, β) from the first principles (ie. from the Hamiltonian). The anomalous transport exponent μ is derived as μ = β/α or μ = β/2α for the first order mean displacement in self-similar transport.

  16. Turbulent kinetic energy equation and free mixing

    NASA Technical Reports Server (NTRS)

    Morel, T.; Torda, T. P.; Bradshaw, P.

    1973-01-01

    Calculation of free shear flows was carried out to investigate the usefulness of several concepts which were previously successfully applied to wall flows. The method belongs to the class of differential approaches. The turbulence is taken into account by the introduction of one additional partial differential equation, the transport equation for the turbulent shear stress. The structure of turbulence is modeled after Bradshaw et al. This model was used successfully in boundary layers and its applicability to other flows is demonstrated. The work reported differs substantially from that of an earlier attempt to use this approach for calculation of free flows. The most important difference is that the region around the center line is treated by invoking the interaction hypothesis (concerning the structure of turbulence in the regions separated by the velocity extrema). The compressibility effects on shear layer spreading at low and moderate Mach numbers were investigated. In the absence of detailed experiments in free flows, the evidence from boundary layers that at low Mach numbers the structure of turbulence is unaffected by the compressibility was relied on. The present model was tested over a range of self-preserving and developing flows including pressure gradients using identical empirical input. The dependence of the structure of turbulence on the spreading rate of the shear layer was established.

  17. The Linearized Kinetic Equation -- A Functional Analytic Approach

    NASA Astrophysics Data System (ADS)

    Brinkmann, Ralf Peter

    2009-10-01

    Kinetic models of plasma phenomena are difficult to address for two reasons. They i) are given as systems of nonlinear coupled integro-differential equations, and ii) involve generally six-dimensional distribution functions f(r,v,t). In situations which can be addressed in a linear regime, the first difficulty disappears, but the second one still poses considerable practical problems. This contribution presents an abstract approach to linearized kinetic theory which employs the methods of functional analysis. A kinetic electron equation with elastic electron-neutral interaction is studied in the electrostatic approximation. Under certain boundary conditions, a nonlinear functional, the kinetic free energy, exists which has the properties of a Lyapunov functional. In the linear regime, the functional becomes a quadratic form which motivates the definition of a bilinear scalar product, turning the space of all distribution functions into a Hilbert space. The linearized kinetic equation can then be described in terms of dynamical operators with well-defined properties. Abstract solutions can be constructed which have mathematically plausible properties. As an example, the formalism is applied to the example of the multipole resonance probe (MRP). Under the assumption of a Maxwellian background distribution, the kinetic model of that diagnostics device is compared to a previously investigated fluid model.

  18. Drift kinetic equation exact through second order in gyroradius expansion

    SciTech Connect

    Simakov, Andrei N.; Catto, Peter J.

    2005-01-01

    The drift kinetic equation of Hazeltine [R. D. Hazeltine, Plasma Phys. 15, 77 (1973)] for a magnetized plasma of arbitrary collisionality is widely believed to be exact through the second order in the gyroradius expansion. It is demonstrated that this equation is only exact through the first order. The reason is that when evaluating the second-order gyrophase dependent distribution function, Hazeltine neglected contributions from the first-order gyrophase dependent distribution function, and then used this incomplete expression to derive the equation for the gyrophase independent distribution function. Consequently, the second-order distribution function and the stress tensor derived by this approach are incomplete. By relaxing slightly Hazeltine's orderings one is able to obtain a drift kinetic equation accurate through the second order in the gyroradius expansion. In addition, the gyroviscous stress tensor for plasmas of arbitrary collisionality is obtained.

  19. Consistent description of kinetic equation with triangle anomaly

    SciTech Connect

    Pu Shi; Gao Jianhua; Wang Qun

    2011-05-01

    We provide a consistent description of the kinetic equation with a triangle anomaly which is compatible with the entropy principle of the second law of thermodynamics and the charge/energy-momentum conservation equations. In general an anomalous source term is necessary to ensure that the equations for the charge and energy-momentum conservation are satisfied and that the correction terms of distribution functions are compatible to these equations. The constraining equations from the entropy principle are derived for the anomaly-induced leading order corrections to the particle distribution functions. The correction terms can be determined for the minimum number of unknown coefficients in one charge and two charge cases by solving the constraining equations.

  20. Enzyme Kinetics and the Michaelis-Menten Equation

    ERIC Educational Resources Information Center

    Biaglow, Andrew; Erickson, Keith; McMurran, Shawnee

    2010-01-01

    The concepts presented in this article represent the cornerstone of classical mathematical biology. The central problem of the article relates to enzyme kinetics, which is a biochemical system. However, the theoretical underpinnings that lead to the formation of systems of time-dependent ordinary differential equations have been applied widely to…

  1. General Entropic Approximations for Canonical Systems Described by Kinetic Equations

    NASA Astrophysics Data System (ADS)

    Pavan, V.

    2011-02-01

    In this paper we extend the general construction of entropic approximation for kinetic operators modelling canonical systems. More precisely, this paper aims at pursuing to thermalized systems the works of Levermore, Schneider and Junk on moments problems relying on entropy minimization in order to construct BGK approximations and moments based equations.

  2. Kinetic schemes for the ultra-relativistic Euler equations

    NASA Astrophysics Data System (ADS)

    Kunik, Matthias; Qamar, Shamsul; Warnecke, Gerald

    2003-05-01

    We present a kinetic numerical scheme for the relativistic Euler equations, which describe the flow of a perfect fluid in terms of the particle density n, the spatial part of the four-velocity u and the pressure p. The kinetic approach is very simple in the ultra-relativistic limit, but may also be applied to more general cases. The basic ingredients of the kinetic scheme are the phase-density in equilibrium and the free flight. The phase-density generalizes the non-relativistic Maxwellian for a gas in local equilibrium. The free flight is given by solutions of a collision free kinetic transport equation. The scheme presented here is an explicit method and unconditionally stable. We establish that the conservation laws of mass, momentum and energy as well as the entropy inequality are everywhere exactly satisfied by the solution of the kinetic scheme. For that reason we obtain weak admissible Euler solutions including arbitrarily complicated shock interactions. In the numerical case studies the results obtained from the kinetic scheme are compared with the first order upwind and centered schemes.

  3. Bounce-averaged Kinetic Equations and Neoclassical Polarization Density

    SciTech Connect

    First Author = B.H. Fong; T.S. Hahm

    1998-07-01

    The rigorous formulation of the bounce-averaged equations is presented based upon the Poincare-Cartan one-form andLie perturbation methods. The resulting bounce-averaged Vlasov equation is Hamiltonian, thus suitable for theself-consistent simulation of low-frequency electrostatic turbulence in the trapped ion mode regime. In the bounce-kineticPoisson equation, the "neoclassical polarization density" arises from the difference between bounce-averaged banana centerand real trapped particle densities across a field line. This representation of the neoclassical polarization drift as ashielding term provides a systematic way to study the long-term behavior of the turbulence-driven E x B flow.

  4. Spectral function and kinetic equation for a normal Fermi liquid

    SciTech Connect

    Arshad, M.; Siddique, I.; Kondratyev, A. S.

    2007-08-01

    On the basis of the Kadanoff-Baym (KB) version of the time-dependent Green's function method, an Ansatz for the approximation of a spectral function is offered. The Ansatz possesses all the advantages of quasiparticle and extended quasiparticle approximations and satisfies the KB equation for a spectral function in the case of slightly nonequilibrium system when disturbances in space and time are taken into consideration in the gradient approximation. This feature opens opportunities for the microscopic derivation of the Landau kinetic equation for the quasiparticle distribution function of the normal Fermi liquid and provides the widening of these equations' temperature range of validity.

  5. Kinetic equations for diffusion in the presence of entropic barriers.

    PubMed

    Reguera, D; Rubí, J M

    2001-12-01

    We use the mesoscopic nonequilibrium thermodynamics theory to derive the general kinetic equation of a system in the presence of potential barriers. The result is applied to a description of the evolution of systems whose dynamics is influenced by entropic barriers. We analyze in detail the case of diffusion in a domain of irregular geometry in which the presence of the boundaries induces an entropy barrier when approaching the exact dynamics by a coarsening of the description. The corresponding kinetic equation, named the Fick-Jacobs equation, is obtained, and its validity is generalized through the formulation of a scaling law for the diffusion coefficient which depends on the shape of the boundaries. The method we propose can be useful to analyze the dynamics of systems at the nanoscale where the presence of entropy barriers is a common feature. PMID:11736170

  6. A kinetic equation with kinetic entropy functions for scalar conservation laws

    NASA Technical Reports Server (NTRS)

    Perthame, Benoit; Tadmor, Eitan

    1990-01-01

    A nonlinear kinetic equation is constructed and proved to be well-adapted to describe general multidimensional scalar conservation laws. In particular, it is proved to be well-posed uniformly in epsilon - the microscopic scale. It is also shown that the proposed kinetic equation is equipped with a family of kinetic entropy functions - analogous to Boltzmann's microscopic H-function, such that they recover Krushkov-type entropy inequality on the macroscopic scale. Finally, it is proved by both - BV compactness arguments in the one-dimensional case, that the local density of kinetic particles admits a continuum limit, as it converges strongly with epsilon below 0 to the unique entropy solution of the corresponding conservation law.

  7. Derivation of the Drude conductivity from quantum kinetic equations

    NASA Astrophysics Data System (ADS)

    Kitamura, Hikaru

    2015-11-01

    The Drude formula of ac (frequency-dependent) electric conductivity has been established as a simple and practically useful model to understand the electromagnetic response of simple free-electron-like metals. In most textbooks of solid-state physics, the Drude formula is derived from either a classical equation of motion or the semiclassical Boltzmann transport equation. On the other hand, quantum-mechanical derivation of the Drude conductivity, which requires an appropriate treatment of phonon-assisted intraband transitions with small momentum transfer, has not been well documented except for the zero- or high-frequency case. Here, a lucid derivation of the Drude conductivity that covers the entire frequency range is presented by means of quantum kinetic equations in the density-matrix formalism. The derivation is straightforward so that advanced undergraduate students or early-year graduate students will be able to gain insight into the link between the microscopic Schrödinger equation and macroscopic transport.

  8. Analytical Derivation of Moment Equations in Stochastic Chemical Kinetics

    PubMed Central

    Sotiropoulos, Vassilios; Kaznessis, Yiannis N.

    2011-01-01

    The master probability equation captures the dynamic behavior of a variety of stochastic phenomena that can be modeled as Markov processes. Analytical solutions to the master equation are hard to come by though because they require the enumeration of all possible states and the determination of the transition probabilities between any two states. These two tasks quickly become intractable for all but the simplest of systems. Instead of determining how the probability distribution changes in time, we can express the master probability distribution as a function of its moments, and, we can then write transient equations for the probability distribution moments. In 1949, Moyal defined the derivative, or jump, moments of the master probability distribution. These are measures of the rate of change in the probability distribution moment values, i.e. what the impact is of any given transition between states on the moment values. In this paper we present a general scheme for deriving analytical moment equations for any N-dimensional Markov process as a function of the jump moments. Importantly, we propose a scheme to derive analytical expressions for the jump moments for any N-dimensional Markov process. To better illustrate the concepts, we focus on stochastic chemical kinetics models for which we derive analytical relations for jump moments of arbitrary order. Chemical kinetics models are widely used to capture the dynamic behavior of biological systems. The elements in the jump moment expressions are a function of the stoichiometric matrix and the reaction propensities, i.e the probabilistic reaction rates. We use two toy examples, a linear and a non-linear set of reactions, to demonstrate the applicability and limitations of the scheme. Finally, we provide an estimate on the minimum number of moments necessary to obtain statistical significant data that would uniquely determine the dynamics of the underlying stochastic chemical kinetic system. The first two moments

  9. Immersed boundary method for Boltzmann model kinetic equations

    NASA Astrophysics Data System (ADS)

    Pekardan, Cem; Chigullapalli, Sruti; Sun, Lin; Alexeenko, Alina

    2012-11-01

    Three different immersed boundary method formulations are presented for Boltzmann model kinetic equations such as Bhatnagar-Gross-Krook (BGK) and Ellipsoidal statistical Bhatnagar-Gross-Krook (ESBGK) model equations. 1D unsteady IBM solution for a moving piston is compared with the DSMC results and 2D quasi-steady microscale gas damping solutions are verified by a conformal finite volume method solver. Transient analysis for a sinusoidally moving beam is also carried out for the different pressure conditions (1 atm, 0.1 atm and 0.01 atm) corresponding to Kn=0.05,0.5 and 5. Interrelaxation method (Method 2) is shown to provide a faster convergence as compared to the traditional interpolation scheme used in continuum IBM formulations. Unsteady damping in rarefied regime is characterized by a significant phase-lag which is not captured by quasi-steady approximations.

  10. Coarse-grained kinetic equations for quantum systems

    NASA Astrophysics Data System (ADS)

    Petrov, E. G.

    2013-01-01

    The nonequilibrium density matrix method is employed to derive a master equation for the averaged state populations of an open quantum system subjected to an external high frequency stochastic field. It is shown that if the characteristic time τstoch of the stochastic process is much lower than the characteristic time τsteady of the establishment of the system steady state populations, then on the time scale Δ t ˜ τsteady, the evolution of the system populations can be described by the coarse-grained kinetic equations with the averaged transition rates. As an example, the exact averaging is carried out for the dichotomous Markov process of the kangaroo type.

  11. Hamiltonian fluid reductions of drift-kinetic equations and the link with water-bags

    NASA Astrophysics Data System (ADS)

    Perin, M.; Chandre, C.; Tassi, E.

    2016-07-01

    Hamiltonian models for the first three moments of the drift-kinetic distribution function, namely the density, the fluid velocity and the parallel pressure, are derived from the Hamiltonian structure of the drift-kinetic equations. The link with the water-bag closure is established, showing that, unlike the one-dimensional Vlasov equations, these solutions are the only Hamiltonian fluid reductions for the drift-kinetic equation. These models are discussed through their equations of motion and their Casimir invariants.

  12. Pauli equation for semiconductor quantum dot photoluminescence kinetics investigation

    NASA Astrophysics Data System (ADS)

    Turkov, Vadim K.; Leonov, Mikhail Y.; Rukhlenko, Ivan D.; Fedorov, Anatoly V.

    2012-11-01

    We develop a theory of secondary emission from a single quantum dot, when the lowest-energy states of its electron-hole pairs are involved in the photoluminescence process. For the sake of definiteness, our model allows for two states contributing to the luminescence. We analyze the dependency of secondary emission intensity on the energy gap between the states, while considering that the gap is determined by the quantum dot's size. An analytical expression for the time-dependent signal of thermalized luminescence is obtained using an analytical solution to the kinetic Pauli equation. This expression yields the signal of stationary luminescence as the spectral width of the excitation pulse tends to zero.

  13. Kinetic equations for baryogenesis via sterile neutrino oscillation

    NASA Astrophysics Data System (ADS)

    Asaka, Takehiko; Eijima, Sintaro; Ishida, Hiroyuki

    2012-02-01

    We investigate baryogenesis in the νMSM (neutrino Minimal Standard Model), which is the SM extended by three right-handed neutrinos with masses below the electroweak scale. The baryon asymmetry of the universe can be generated by the mechanism via flavor oscillation of right-handed (sterile) neutrinos which are responsible to masses of active neutrinos confirmed by various experiments. We present the kinetic equations for the matrix of densities of leptons which describe the generation of asymmetries. Especially, the momentum dependence of the matrix of densities is taken into account. By solving these equations numerically, it is found that the momentum distribution is significantly distorted from the equilibrium one, since the production for the modes with lower momenta k << T (T is the temperature of the universe) is enhanced, while suppressed for higher modes. As a result, the most important mode for the yields of sterile neutrinos as well as the baryon asymmetry is k simeq 2T, which is smaller than langlekrangle inferred from the thermal average. The comparison with the previous works is also discussed.

  14. Verification of continuum drift kinetic equation solvers in NIMROD

    NASA Astrophysics Data System (ADS)

    Held, E. D.; Kruger, S. E.; Ji, J.-Y.; Belli, E. A.; Lyons, B. C.

    2015-03-01

    Verification of continuum solutions to the electron and ion drift kinetic equations (DKEs) in NIMROD [C. R. Sovinec et al., J. Comp. Phys. 195, 355 (2004)] is demonstrated through comparison with several neoclassical transport codes, most notably NEO [E. A. Belli and J. Candy, Plasma Phys. Controlled Fusion 54, 015015 (2012)]. The DKE solutions use NIMROD's spatial representation, 2D finite-elements in the poloidal plane and a 1D Fourier expansion in toroidal angle. For 2D velocity space, a novel 1D expansion in finite elements is applied for the pitch angle dependence and a collocation grid is used for the normalized speed coordinate. The full, linearized Coulomb collision operator is kept and shown to be important for obtaining quantitative results. Bootstrap currents, parallel ion flows, and radial particle and heat fluxes show quantitative agreement between NIMROD and NEO for a variety of tokamak equilibria. In addition, velocity space distribution function contours for ions and electrons show nearly identical detailed structure and agree quantitatively. A Θ-centered, implicit time discretization and a block-preconditioned, iterative linear algebra solver provide efficient electron and ion DKE solutions that ultimately will be used to obtain closures for NIMROD's evolving fluid model.

  15. Verification of continuum drift kinetic equation solvers in NIMROD

    SciTech Connect

    Held, E. D.; Ji, J.-Y.; Kruger, S. E.; Belli, E. A.; Lyons, B. C.

    2015-03-15

    Verification of continuum solutions to the electron and ion drift kinetic equations (DKEs) in NIMROD [C. R. Sovinec et al., J. Comp. Phys. 195, 355 (2004)] is demonstrated through comparison with several neoclassical transport codes, most notably NEO [E. A. Belli and J. Candy, Plasma Phys. Controlled Fusion 54, 015015 (2012)]. The DKE solutions use NIMROD's spatial representation, 2D finite-elements in the poloidal plane and a 1D Fourier expansion in toroidal angle. For 2D velocity space, a novel 1D expansion in finite elements is applied for the pitch angle dependence and a collocation grid is used for the normalized speed coordinate. The full, linearized Coulomb collision operator is kept and shown to be important for obtaining quantitative results. Bootstrap currents, parallel ion flows, and radial particle and heat fluxes show quantitative agreement between NIMROD and NEO for a variety of tokamak equilibria. In addition, velocity space distribution function contours for ions and electrons show nearly identical detailed structure and agree quantitatively. A Θ-centered, implicit time discretization and a block-preconditioned, iterative linear algebra solver provide efficient electron and ion DKE solutions that ultimately will be used to obtain closures for NIMROD's evolving fluid model.

  16. A quasi-linear kinetic equation for cosmic rays in the interplanetary medium

    NASA Technical Reports Server (NTRS)

    Luhmann, J. G.

    1976-01-01

    A kinetic equation for interplanetary cosmic rays is set up with the aid of weak-plasma-turbulence theory for an idealized radially symmetric model of the interplanetary magnetic field. As a starting point, this treatment invokes the Vlasov equation instead of the traditional Fokker-Planck equation. Quasi-linear theory is applied to obtain a momentum diffusion equation for the heliocentric frame of reference which describes the interaction of cosmic rays with convecting magnetic irregularities in the solar-wind plasma. Under restricted conditions, the well-known equation of solar modulation can be obtained from this kinetic equation.

  17. A comparison of the efficiency of numerical methods for integrating chemical kinetic rate equations

    NASA Technical Reports Server (NTRS)

    Radhakrishnan, K.

    1984-01-01

    The efficiency of several algorithms used for numerical integration of stiff ordinary differential equations was compared. The methods examined included two general purpose codes EPISODE and LSODE and three codes (CHEMEQ, CREK1D and GCKP84) developed specifically to integrate chemical kinetic rate equations. The codes were applied to two test problems drawn from combustion kinetics. The comparisons show that LSODE is the fastest code available for the integration of combustion kinetic rate equations. It is shown that an iterative solution of the algebraic energy conservation equation to compute the temperature can be more efficient then evaluating the temperature by integrating its time-derivative.

  18. Effective computation of stochastic protein kinetic equation by reducing stiffness via variable transformation

    NASA Astrophysics Data System (ADS)

    Wang, Lijin

    2016-06-01

    The stochastic protein kinetic equations can be stiff for certain parameters, which makes their numerical simulation rely on very small time step sizes, resulting in large computational cost and accumulated round-off errors. For such situation, we provide a method of reducing stiffness of the stochastic protein kinetic equation by means of a kind of variable transformation. Theoretical and numerical analysis show effectiveness of this method. Its generalization to a more general class of stochastic differential equation models is also discussed.

  19. Kinetic equation for strongly interacting dense Fermi systems

    NASA Astrophysics Data System (ADS)

    Lipavský, P.; Morawetz, K.; Špička, V.

    We review the non-relativistic Green's-function approach to the kinetic equations for Fermi liquids far from equilibrium. The emphasis is on the consistent treatment of the off-shell motion between collisions and on the non-instant and non-local picture of binary collisions.The resulting kinetic equation is of the Boltzmann type, and it represents an interpolation between the theory of transport in metals and the theory of moderately dense gases. The free motion of particles is renormalised by various mean field and mass corrections in the spirit of Landau's quasiparticles in metals. The collisions are non-local in the spirit of Enskog's theory of non-ideal gases. The collisions are moreover non-instant, a feature which is absent in the theory of gases, but which is shown to be important for dense Fermi systems.In spite of its formal complexity, the presented theory has a simple implementation within the Monte-Carlo simulation schemes. Applications in nuclear physics are given for heavy-ion reactions and the results are compared with the former theory and recent experimental data.The effect of the off-shell motion and the non-local and non-instant collisions on the dynamics of the system can be characterised in terms of thermodynamic functions such as the energy density or the pressure tensor. Non-equilibrium counterparts of these functions and the corresponding balance equations are derived and discussed from two points of view. Firstly, they are used to prove the conservation laws. Secondly, the role of individual microscopic mechanisms in fluxes of particles and momenta and in transformations of the energy is clarified. Nous examinons la technique des fonctions de Green non relativistes appliquée aux équations cinétiques pour les liquides de Fermi hors équilibre. L'accent est mis sur le traitement cohérent des effets hors couche entre les collisions ainsi que sur l'aspect non-local et non-instantané des collisions binaires.L'équation cinétique r

  20. Parameter Estimates in Differential Equation Models for Chemical Kinetics

    ERIC Educational Resources Information Center

    Winkel, Brian

    2011-01-01

    We discuss the need for devoting time in differential equations courses to modelling and the completion of the modelling process with efforts to estimate the parameters in the models using data. We estimate the parameters present in several differential equation models of chemical reactions of order n, where n = 0, 1, 2, and apply more general…

  1. Formulation and closure of compressible turbulence equations in the light of kinetic theory

    NASA Technical Reports Server (NTRS)

    Tsuge, S.; Sagara, K.

    1976-01-01

    Fluid-dynamic moment equations, based on a kinetic hierarchy system, are derived governing the interaction between turbulent and thermal fluctuations. The kinetic theory is shown to reduce the inherent complexity of the conventional formalism of compressible turbulence theory and to minimize arbitrariness in formulating the closure condition.

  2. Entropy Analysis of Kinetic Flux Vector Splitting Schemes for the Compressible Euler Equations

    NASA Technical Reports Server (NTRS)

    Shiuhong, Lui; Xu, Jun

    1999-01-01

    Flux Vector Splitting (FVS) scheme is one group of approximate Riemann solvers for the compressible Euler equations. In this paper, the discretized entropy condition of the Kinetic Flux Vector Splitting (KFVS) scheme based on the gas-kinetic theory is proved. The proof of the entropy condition involves the entropy definition difference between the distinguishable and indistinguishable particles.

  3. Kinetic Equation for Two-Particle Distribution Function in Boltzmann Gas Mixtures and Equation of Motion for Quasiparticle Pairs

    NASA Astrophysics Data System (ADS)

    Saveliev, V. L.

    2011-05-01

    Pair collisions is the main interaction process in the Boltzmann gas dynamics. By making use of exactly the same physical assumptions as was used by Ludwig Boltzmann we write the kinetic equation for two-particle distribution function of molecules in the gas mixtures. Instead of the collision integral, there are the linear scattering operator and the chaos projector in the right part of this equation. Because the scattering operator is more simple then Boltzmann collision integral this equation opens new opportunities for mathematical description of the Boltzmann gas dynamics.

  4. Integrating chemical kinetic rate equations by selective use of stiff and nonstiff methods

    NASA Technical Reports Server (NTRS)

    Radhakrishnan, K.

    1985-01-01

    The effect of switching between nonstiff and stiff methods on the efficiency of algorithms for integrating chemical kinetic rate equations is presented. Different integration methods are tested by application of the packaged code LSODE to four practical combustion kinetics problems. The problems describe adiabatic, homogeneous gas-phase combustion reactions. It is shown that selective use of nonstiff and stiff methods in different regimes of a typical batch combustion problem is faster than the use of either method for the entire problem. The implications of this result to the development of fast integration techniques for combustion kinetic rate equations are discussed.

  5. Integrating chemical kinetic rate equations by selective use of stiff and nonstiff methods

    NASA Technical Reports Server (NTRS)

    Radhakrishnan, K.

    1985-01-01

    The effect of switching between nonstiff and stiff methods on the efficiency of algorithms for integrating chemical kinetic rate equations was examined. Different integration methods were tested by application of the packaged code LSODE to four practical combustion kinetics problems. The problems describe adiabatic, and homogeneous gas phase combustion reactions. It is shown that selective use of nonstiff and stiff methods in different regimes of a typical batch combustion problem is faster than the use of either method for the entire problem. The implications which result in the development of fast integration techniques for combustion kinetic rate equations are discussed.

  6. Nonlinear Drift-Kinetic Equation in the Presence of a Circularly Polarized Wave

    NASA Technical Reports Server (NTRS)

    Khazanov, G. V.; Krivorutsky, E. N.; Whitaker, Ann F. (Technical Monitor)

    2001-01-01

    Equations of the single particle motion and nonlinear kinetic equation for plasma in the presence of a circularly polarized wave of arbitrary frequency in the drift approximation are presented. The nonstationarity and inhomogeneity of the plasma-wave system are taken into account.

  7. Nonlinear Drift-Kinetic Equation in the Presence of a Circularly Polarized Wave

    NASA Technical Reports Server (NTRS)

    Khazanov, G. V.; Krivorutsky, E. N.; Six, N. Frank (Technical Monitor)

    2002-01-01

    Equations of the single particle motion and nonlinear kinetic equation for plasma in the presence of a circularly polarized wave of arbitrary frequency in the drift approximation are presented. The nonstationarity and inhomogeneity of the plasma-wave system are taken into account. The time dependent part of the ponderomotive force is discussed.

  8. A comparison of the efficiency of numerical methods for integrating chemical kinetic rate equations

    NASA Technical Reports Server (NTRS)

    Radhakrishnan, K.

    1984-01-01

    A comparison of the efficiency of several algorithms recently developed for the efficient numerical integration of stiff ordinary differential equations is presented. The methods examined include two general-purpose codes EPISODE and LSODE and three codes (CHEMEQ, CREK1D, and GCKP84) developed specifically to integrate chemical kinetic rate equations. The codes are applied to two test problems drawn from combustion kinetics. The comparisons show that LSODE is the fastest code currently available for the integration of combustion kinetic rate equations. An important finding is that an iterative solution of the algebraic energy conservation equation to compute the temperature can be more efficient than evaluating the temperature by integrating its time-derivative.

  9. Mapping of the classical kinetic balance equations onto the Pauli equation

    NASA Astrophysics Data System (ADS)

    Pesci, Adriana I.; Goldstein, Raymond E.; Uys, Hermann

    2005-01-01

    Here we find a mapping onto the Pauli equation of the first two balance equations derived from the classical Boltzmann equation. The essence of this mapping, which we previously used to obtain the particular case of the Sturm-Liouville operator known as Schrödinger's equation, consists of applying a Fourier transform to the momentum coordinate of the distribution function. This procedure introduces a natural parameter η with units of angular momentum. The main difference between the two cases is the conditions imposed on the probability distribution function, a difference most clearly understood at the level of the hydrodynamic equations generated in the first steps of the mapping. The case leading to the Sturm-Liouville operator corresponds to an irrotational flow, while here the ansatz leading to the Pauli equation corresponds to a fluid with non-zero vorticity. In the context of fluid dynamics, the magnitude of the angular momentum associated with the vorticity is η/2. To perform the mapping we follow the standard technique common in hydrodynamic problems, namely writing the Lagrangian for the Euler equations with the corresponding constraints expressed in terms of the Clebsch variables.

  10. New integration techniques for chemical kinetic rate equations. II - Accuracy comparison

    NASA Technical Reports Server (NTRS)

    Radhakrishnan, K.

    1986-01-01

    A comparison of the accuracy of several techniques recently developed for solving stiff differential equations is presented. The techniques examined include two general purpose codes EEPISODE and LSODE developed for an arbitrary system of ordinary differential equations, and three specialized codes CHEMEQ, CREKID, and GCKP84 developed specifically to solve chemical kinetic rate equations. The accuracy comparisons are made by applying these solution procedures to two practical combustion kinetics problems. Both problems describe adiabatic, homogeneous, gas phase chemical reactions at constant pressure, and include all three combustion regimes: induction heat release, and equilibration. The comparisons show that LSODE is the most efficient code - in the sense that it requires the least computational work to attain a specified accuracy level. An important finding is that an iterative solution of the algebraic enthalpy conservation equation for the temperature can be more accurate and efficient than computing the temperature by integrating its time derivative.

  11. New integration techniques for chemical kinetic rate equations. 2: Accuracy comparison

    NASA Technical Reports Server (NTRS)

    Radhakrishnan, K.

    1985-01-01

    A comparison of the accuracy of several techniques recently developed for solving stiff differential equations is presented. The techniques examined include two general purpose codes EEPISODE and LSODE developed for an arbitrary system of ordinary differential equations, and three specialized codes CHEMEQ, CREKID, and GCKP84 developed specifically to solve chemical kinetic rate equations. The accuracy comparisons are made by applying these solution procedures to two practical combustion kinetics problems. Both problems describe adiabatic, homogeneous, gas phase chemical reactions at constant pressure, and include all three combustion regimes: induction, heat release, and equilibration. The comparisons show that LSODE is the most efficient code - in the sense that it requires the least computational work to attain a specified accuracy level. An important finding is that an iterative solution of the algebraic enthalpy conservation equation for the temperature can be more accurate and efficient than computing the temperature by integrating its time derivative.

  12. Analytical solution of the kinetic equation for a uniform plasma in a magnetic field

    SciTech Connect

    Ji, Jeong-Young; Held, Eric D.

    2010-07-15

    The kinetic equation for a single-component uniform plasma in a magnetic field is analytically solved by the moment method. The linear system of ordinary differential equations for the moments is decomposed into subsystems of lower dimensions by a geometric method. The eigensystem of each subsystem shows that parallel moments decay monotonically, but perpendicular lth harmonic moments decay while oscillating with the l,l-2,...,-th harmonics of gyrofrequency. A generalization to a multicomponent plasma is discussed.

  13. BRIEF COMMUNICATION: On the drift kinetic equation driven by plasma flows

    NASA Astrophysics Data System (ADS)

    Shaing, K. C.

    2010-07-01

    A drift kinetic equation that is driven by plasma flows has previously been derived by Shaing and Spong 1990 (Phys. Fluids B 2 1190). The terms that are driven by particle speed that is parallel to the magnetic field B have been neglected. Here, such terms are discussed to examine their importance to the equation and to show that these terms do not contribute to the calculations of plasma viscosity in large aspect ratio toroidal plasmas, e.g. tokamaks and stellarators.

  14. Kinetic equations for a density matrix describing nonlinear effects in spectral line wings

    SciTech Connect

    Parkhomenko, A. I. Shalagin, A. M.

    2011-11-15

    Kinetic quantum equations are derived for a density matrix with collision integrals describing nonlinear effects in spectra line wings. These equations take into account the earlier established inequality of the spectral densities of Einstein coefficients for absorption and stimulated radiation emission by a two-level quantum system in the far wing of a spectral line in the case of frequent collisions. The relationship of the absorption and stimulated emission probabilities with the characteristics of radiation and an elementary scattering event is found.

  15. Second-order accurate kinetic schemes for the ultra-relativistic Euler equations

    NASA Astrophysics Data System (ADS)

    Kunik, Matthias; Qamar, Shamsul; Warnecke, Gerald

    2003-12-01

    A second-order accurate kinetic scheme for the numerical solution of the relativistic Euler equations is presented. These equations describe the flow of a perfect fluid in terms of the particle density n, the spatial part of the four-velocity u and the pressure p. The kinetic scheme, is based on the well-known fact that the relativistic Euler equations are the moments of the relativistic Boltzmann equation of the kinetic theory of gases when the distribution function is a relativistic Maxwellian. The kinetic scheme consists of two phases, the convection phase (free-flight) and collision phase. The velocity distribution function at the end of the free-flight is the solution of the collisionless transport equation. The collision phase instantaneously relaxes the distribution to the local Maxwellian distribution. The fluid dynamic variables of density, velocity, and internal energy are obtained as moments of the velocity distribution function at the end of the free-flight phase. The scheme presented here is an explicit method and unconditionally stable. The conservation laws of mass, momentum and energy as well as the entropy inequality are everywhere exactly satisfied by the solution of the kinetic scheme. The scheme also satisfies positivity and L1-stability. The scheme can be easily made into a total variation diminishing method for the distribution function through a suitable choice of the interpolation strategy. In the numerical case studies the results obtained from the first- and second-order kinetic schemes are compared with the first- and second-order upwind and central schemes. We also calculate the experimental order of convergence and numerical L1-stability of the scheme for smooth initial data.

  16. Fundamental equations of a mixture of gas and small spherical solid particles from simple kinetic theory.

    NASA Technical Reports Server (NTRS)

    Pai, S. I.

    1973-01-01

    The fundamental equations of a mixture of a gas and pseudofluid of small spherical solid particles are derived from the Boltzmann equation of two-fluid theory. The distribution function of the gas molecules is defined in the same manner as in the ordinary kinetic theory of gases, but the distribution function for the solid particles is different from that of the gas molecules, because it is necessary to take into account the different size and physical properties of solid particles. In the proposed simple kinetic theory, two additional parameters are introduced: one is the radius of the spheres and the other is the instantaneous temperature of the solid particles in the distribution of the solid particles. The Boltzmann equation for each species of the mixture is formally written, and the transfer equations of these Boltzmann equations are derived and compared to the well-known fundamental equations of the mixture of a gas and small solid particles from continuum theory. The equations obtained reveal some insight into various terms in the fundamental equations. For instance, the partial pressure of the pseudofluid of solid particles is not negligible if the volume fraction of solid particles is not negligible as in the case of lunar ash flow.

  17. Numerical methods for solving moment equations in kinetic theory of neuronal network dynamics

    NASA Astrophysics Data System (ADS)

    Rangan, Aaditya V.; Cai, David; Tao, Louis

    2007-02-01

    Recently developed kinetic theory and related closures for neuronal network dynamics have been demonstrated to be a powerful theoretical framework for investigating coarse-grained dynamical properties of neuronal networks. The moment equations arising from the kinetic theory are a system of (1 + 1)-dimensional nonlinear partial differential equations (PDE) on a bounded domain with nonlinear boundary conditions. The PDEs themselves are self-consistently specified by parameters which are functions of the boundary values of the solution. The moment equations can be stiff in space and time. Numerical methods are presented here for efficiently and accurately solving these moment equations. The essential ingredients in our numerical methods include: (i) the system is discretized in time with an implicit Euler method within a spectral deferred correction framework, therefore, the PDEs of the kinetic theory are reduced to a sequence, in time, of boundary value problems (BVPs) with nonlinear boundary conditions; (ii) a set of auxiliary parameters is introduced to recast the original BVP with nonlinear boundary conditions as BVPs with linear boundary conditions - with additional algebraic constraints on the auxiliary parameters; (iii) a careful combination of two Newton's iterates for the nonlinear BVP with linear boundary condition, interlaced with a Newton's iterate for solving the associated algebraic constraints is constructed to achieve quadratic convergence for obtaining the solutions with self-consistent parameters. It is shown that a simple fixed-point iteration can only achieve a linear convergence for the self-consistent parameters. The practicability and efficiency of our numerical methods for solving the moment equations of the kinetic theory are illustrated with numerical examples. It is further demonstrated that the moment equations derived from the kinetic theory of neuronal network dynamics can very well capture the coarse-grained dynamical properties of

  18. New integration techniques for chemical kinetic rate equations. I - Efficiency comparison

    NASA Technical Reports Server (NTRS)

    Radhakrishnan, K.

    1986-01-01

    A comparison of the efficiency of several recently developed numerical techniques for solving chemical kinetic rate equations is presented. The solution procedures examined include two general-purpose codes, EPISODE and LSODE, developed as multipurpose differential equation solvers, and three specialzed codes, CHEMEQ, CREK1D, and GCKP84, developed specifically for chemical kinetics. The efficiency comparison is made by applying these codes to two practical combustion kinetics problems. Both problems describe adiabatic, constant-pressure, gas-phase chemical reactions and include all three combustion regimes: induction, heat release, and equilibration. The comparison shows that LSODE is the fastest routine currently available for solving chemical kinetic rate equations. An important finding is that an iterative solution of the algebraic enthalpy conservation equation for temperature can be significantly faster than evaluation of the temperature by integration of its time derivative. Significant increases in computational speed are realized by updating the reaction rate constants only when the temperature change exceeds an amount Delta-T that is problem dependent. An approximate expression for the automatic evaluation of Delta-T is presented and is shown to result in increased computational speed.

  19. The Master kinetic equation for the statistical treatment of the Boltzmann-Sinai classical dynamical system

    NASA Astrophysics Data System (ADS)

    Tessarotto, Massimo; Cremaschini, Claudio

    2014-07-01

    In this investigation, exact particular realizations are sought for the microscopic statistical description which is associated with the classical dynamical system (CDS) formed by N identical smooth hard spheres subject to elastic collisions ( S N -CDS). The problem is posed in the framework of the ab initio statistical description of S N -CDS recently developed. It is shown that the Liouville equation associated with SN-CDS admits an exact particular solution for the N-body probability density function (PDF). This is factorized in terms of the i-th particle 1-body PDF (for all i = 1, N) via suitable weighting factors, which are denoted here as particle occupation coefficients. The latter are found to depend functionally only on the 1-body PDFs which are associated with each of the remaining particles belonging to S N -CDS. Furthermore, the 1-body PDF is proved to obey a well-defined statistical equation, referred to here as Master kinetic equation. This is an exact kinetic equation which takes into account the occurrence of configuration-space correlations due to the finite size of the extended particles, while depending functionally on the same 1-body PDF only. The asymptotic approximation of the Master equation, which holds in validity of the Boltzmann-Grad limit, is shown to recover in a suitable asymptotic sense the customary Boltzmann equation. Finally, a critical analysis is presented of the original and modified versions of the Enskog kinetic equation, as well as of some of the non-linear kinetic approaches formulated in the past for dense granular gases. Their conditions of validity and main differences with respect to the present theory are pointed out.

  20. Geodesic acoustic mode in anisotropic plasmas using double adiabatic model and gyro-kinetic equation

    SciTech Connect

    Ren, Haijun; Cao, Jintao

    2014-12-15

    Geodesic acoustic mode in anisotropic tokamak plasmas is theoretically analyzed by using double adiabatic model and gyro-kinetic equation. The bi-Maxwellian distribution function for guiding-center ions is assumed to obtain a self-consistent form, yielding pressures satisfying the magnetohydrodynamic (MHD) anisotropic equilibrium condition. The double adiabatic model gives the dispersion relation of geodesic acoustic mode (GAM), which agrees well with the one derived from gyro-kinetic equation. The GAM frequency increases with the ratio of pressures, p{sub ⊥}/p{sub ∥}, and the Landau damping rate is dramatically decreased by p{sub ⊥}/p{sub ∥}. MHD result shows a low-frequency zonal flow existing for all p{sub ⊥}/p{sub ∥}, while according to the kinetic dispersion relation, no low-frequency branch exists for p{sub ⊥}/p{sub ∥}≳ 2.

  1. General theory of multistage geminate reactions of isolated pairs of reactants. I. Kinetic equations

    SciTech Connect

    Doktorov, Alexander B.; Kipriyanov, Alexey A.

    2014-05-14

    General matrix approach to the consideration of multistage geminate reactions of isolated pairs of reactants depending on reactant mobility is formulated on the basis of the concept of “effective” particles. Various elementary reactions (stages of multistage reaction including physicochemical processes of internal quantum state changes) proceeding with the participation of isolated pairs of reactants (or isolated reactants) are taken into account. Investigation has been made in terms of kinetic approach implying the derivation of general (matrix) kinetic equations for local and mean probabilities of finding any of the reaction species in the sample under study (or for local and mean concentrations). The recipes for the calculation of kinetic coefficients of the equations for mean quantities in terms of relative coordinates of reactants have been formulated in the general case of inhomogeneous reacting systems. Important specific case of homogeneous reacting systems is considered.

  2. Mapping of the relativistic kinetic balance equations onto the Klein Gordon and second-order Dirac equations

    NASA Astrophysics Data System (ADS)

    Pesci, Adriana I.; Goldstein, Raymond E.; Uys, Hermann

    2005-05-01

    In previous work we have shown that the quantum potential can be derived from the classical kinetic equations both for particles with and without spin. Here, we extend these mappings to the relativistic case. The essence of the analysis consists of Fourier transforming the momentum coordinate of the distribution function. This procedure introduces a natural parameter η with units of angular momentum. In the non-relativistic case the ansatz of either separability, or separability and additivity, imposed on the probability distribution function produces mappings onto the Schrödinger equation and the Pauli equation, respectively. The former corresponds to an irrotational flow, the latter to a fluid with non-zero vorticity. In this work we show that the relativistic mappings lead to the Klein-Gordon equation in the irrotational case and to the second-order Dirac equation in the rotational case. These mappings are irreversible; an approximate inverse is the Wigner function. Taken together, these results provide a statistical interpretation of quantum mechanics.

  3. Functional integral derivation of the kinetic equation of two-dimensional point vortices

    NASA Astrophysics Data System (ADS)

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

    2016-08-01

    We present a brief derivation of the kinetic equation describing the secular evolution of point vortices in two-dimensional hydrodynamics, by relying on a functional integral formalism. We start from Liouville's equation which describes the exact dynamics of a two-dimensional system of point vortices. At the order 1 / N, the evolution of the system is characterised by the first two equations of the BBGKY hierarchy involving the system's 1-body distribution function and its 2-body correlation function. Thanks to the introduction of auxiliary fields, these two evolution constraints may be rewritten as a functional integral. When functionally integrated over the 2-body correlation function, this rewriting leads to a new constraint coupling the 1-body distribution function and the two auxiliary fields. Once inverted, this constraint provides, through a new route, the closed non-linear kinetic equation satisfied by the 1-body distribution function. Such a method sheds new lights on the origin of these kinetic equations complementing the traditional derivation methods.

  4. Kinetic effects on Alfven wave nonlinearity. II - The modified nonlinear wave equation

    NASA Technical Reports Server (NTRS)

    Spangler, Steven R.

    1990-01-01

    A previously developed Vlasov theory is used here to study the role of resonant particle and other kinetic effects on Alfven wave nonlinearity. A hybrid fluid-Vlasov equation approach is used to obtain a modified version of the derivative nonlinear Schroedinger equation. The differences between a scalar model for the plasma pressure and a tensor model are discussed. The susceptibilty of the modified nonlinear wave equation to modulational instability is studied. The modulational instability normally associated with the derivative nonlinear Schroedinger equation will, under most circumstances, be restricted to left circularly polarized waves. The nonlocal term in the modified nonlinear wave equation engenders a new modulational instability that is independent of beta and the sense of circular polarization. This new instability may explain the occurrence of wave packet steepening for all values of the plasma beta in the vicinity of the earth's bow shock.

  5. Modified Gompertz equation for electrotherapy murine tumor growth kinetics: predictions and new hypotheses

    PubMed Central

    2010-01-01

    Background Electrotherapy effectiveness at different doses has been demonstrated in preclinical and clinical studies; however, several aspects that occur in the tumor growth kinetics before and after treatment have not yet been revealed. Mathematical modeling is a useful instrument that can reveal some of these aspects. The aim of this paper is to describe the complete growth kinetics of unperturbed and perturbed tumors through use of the modified Gompertz equation in order to generate useful insight into the mechanisms that underpin this devastating disease. Methods The complete tumor growth kinetics for control and treated groups are obtained by interpolation and extrapolation methods with different time steps, using experimental data of fibrosarcoma Sa-37. In the modified Gompertz equation, a delay time is introduced to describe the tumor's natural history before treatment. Different graphical strategies are used in order to reveal new information in the complete kinetics of this tumor type. Results The first stage of complete tumor growth kinetics is highly non linear. The model, at this stage, shows different aspects that agree with those reported theoretically and experimentally. Tumor reversibility and the proportionality between regions before and after electrotherapy are demonstrated. In tumors that reach partial remission, two antagonistic post-treatment processes are induced, whereas in complete remission, two unknown antitumor mechanisms are induced. Conclusion The modified Gompertz equation is likely to lead to insights within cancer research. Such insights hold promise for increasing our understanding of tumors as self-organizing systems and, the possible existence of phase transitions in tumor growth kinetics, which, in turn, may have significant impacts both on cancer research and on clinical practice. PMID:21029411

  6. Kinetic analysis of thermally relativistic flow with dissipation. II. Relativistic Boltzmann equation versus its kinetic models

    NASA Astrophysics Data System (ADS)

    Yano, Ryosuke; Matsumoto, Jun; Suzuki, Kojiro

    2011-06-01

    Thermally relativistic flow with dissipation was analyzed by solving the rarefied supersonic flow of thermally relativistic matter around a triangle prism by Yano and Suzuki [Phys. Rev. DPRVDAQ1550-7998 83, 023517 (2011)10.1103/PhysRevD.83.023517], where the Anderson-Witting (AW) model was used as a solver. In this paper, we solve the same problem, which was analyzed by Yano and Suzuki, using the relativistic Boltzmann equation (RBE). To solve the RBE, the conventional direct simulation Monte Carlo method for the nonrelativistic Boltzmann equation is extended to a new direct simulation Monte Carlo method for the RBE. Additionally, we solve the modified Marle (MM) model proposed by Yano-Suzuki-Kuroda for comparisons. The solution of the thermally relativistic shock layer around the triangle prism obtained using the relativistic Boltzmann equation is considered by focusing on profiles of macroscopic quantities, such as the density, velocity, temperature, heat flux and dynamic pressure along the stagnation streamline (SSL). Differences among profiles of the number density, velocity and temperature along the SSL obtained using the RBE, the AW and MM. models are described in the framework of the relativistic Navier-Stokes-Fourier law. Finally, distribution functions on the SSL obtained using the RBE are compared with those obtained using the AW and MM models. The distribution function inside the shock wave obtained using the RBE does not indicate a bimodal form, which is obtained using the AW and MM models, but a smooth deceleration of thermally relativistic matter inside a shock wave.

  7. Towards an ultra efficient kinetic scheme. Part I: Basics on the BGK equation

    NASA Astrophysics Data System (ADS)

    Dimarco, Giacomo; Loubere, Raphaël

    2013-12-01

    In this paper we present a new ultra efficient numerical method for solving kinetic equations. In this preliminary work, we present the scheme in the case of the BGK relaxation operator. The scheme, being based on a splitting technique between transport and collision, can be easily extended to other collisional operators as the Boltzmann collision integral or to other kinetic equations such as the Vlasov equation. The key idea, on which the method relies, is to solve the collision part on a grid and then to solve exactly the transport linear part by following the characteristics backward in time. The main difference between the method proposed and semi-Lagrangian methods is that here we do not need to reconstruct the distribution function at each time step. This allows to tremendously reduce the computational cost of the method and it permits for the first time, to the author's knowledge, to compute solutions of full six dimensional kinetic equations on a single processor laptop machine. Numerical examples, up to the full three dimensional case, are presented which validate the method and assess its efficiency in 1D, 2D and 3D.

  8. Equations for the kinetic modeling of supersonically flowing electrically excited lasers

    NASA Technical Reports Server (NTRS)

    Lind, R. C.

    1973-01-01

    The equations for the kinetic modeling of a supersonically flowing electrically excited laser system are presented. The work focuses on the use of diatomic gases, in particular carbon monoxide mixtures. The equations presented include the vibrational rate equation which describes the vibrational population distribution, the electron, ion and electronic level rate equations, the gasdynamic equations for an ionized gas in the presence of an applied electric field, and the free electron Boltzmann equation including flow and gradient coupling terms. The model developed accounts for vibration-vibration collisions, vibration-translation collisions, electron-molecule inelastic excitation and superelastic de-excitation collisions, charge particle collisions, ionization and three body recombination collisions, elastic collisions, and radiative decay, all of which take place in such a system. A simplified form of the free electron Boltzmann equation is developed and discussed with emphasis placed on its coupling with the supersonic flow. A brief description of a possible solution procedure for the set of coupled equations is then discussed.

  9. A gas-kinetic BGK scheme for the compressible Navier-Stokes equations

    NASA Technical Reports Server (NTRS)

    Xu, Kun

    2000-01-01

    This paper presents an improved gas-kinetic scheme based on the Bhatnagar-Gross-Krook (BGK) model for the compressible Navier-Stokes equations. The current method extends the previous gas-kinetic Navier-Stokes solver developed by Xu and Prendergast by implementing a general nonequilibrium state to represent the gas distribution function at the beginning of each time step. As a result, the requirement in the previous scheme, such as the particle collision time being less than the time step for the validity of the BGK Navier-Stokes solution, is removed. Therefore, the applicable regime of the current method is much enlarged and the Navier-Stokes solution can be obtained accurately regardless of the ratio between the collision time and the time step. The gas-kinetic Navier-Stokes solver developed by Chou and Baganoff is the limiting case of the current method, and it is valid only under such a limiting condition. Also, in this paper, the appropriate implementation of boundary condition for the kinetic scheme, different kinetic limiting cases, and the Prandtl number fix are presented. The connection among artificial dissipative central schemes, Godunov-type schemes, and the gas-kinetic BGK method is discussed. Many numerical tests are included to validate the current method.

  10. New forms of two-particle and one-particle kinetic equations

    NASA Astrophysics Data System (ADS)

    Saveliev, V. L.; Yonemura, S.

    2012-11-01

    Pair collisions are the main interaction process in the Boltzmann gas dynamics. By making use of exactly the same physical assumptions as was done by Ludwig Boltzmann we wrote the kinetic equation for two-particle distribution function of molecules in gas mixtures. Instead of the collision integral, there are the linear scattering operator and the chaos projector in the right part of this equation. We developed a new technique for factorization of the scattering operator on the bases of right inverses to the Casimir operator of the group of rotations. We exactly transformed the Boltzmann collision integral to the Landau-Fokker-Planck like form.

  11. Solving the master equation without kinetic Monte Carlo: Tensor train approximations for a CO oxidation model

    NASA Astrophysics Data System (ADS)

    Gelß, Patrick; Matera, Sebastian; Schütte, Christof

    2016-06-01

    In multiscale modeling of heterogeneous catalytic processes, one crucial point is the solution of a Markovian master equation describing the stochastic reaction kinetics. Usually, this is too high-dimensional to be solved with standard numerical techniques and one has to rely on sampling approaches based on the kinetic Monte Carlo method. In this study we break the curse of dimensionality for the direct solution of the Markovian master equation by exploiting the Tensor Train Format for this purpose. The performance of the approach is demonstrated on a first principles based, reduced model for the CO oxidation on the RuO2(110) surface. We investigate the complexity for increasing system size and for various reaction conditions. The advantage over the stochastic simulation approach is illustrated by a problem with increased stiffness.

  12. An asymptotic-preserving scheme for linear kinetic equation with fractional diffusion limit

    NASA Astrophysics Data System (ADS)

    Wang, Li; Yan, Bokai

    2016-05-01

    We present a new asymptotic-preserving scheme for the linear Boltzmann equation which, under appropriate scaling, leads to a fractional diffusion limit. Our scheme rests on novel micro-macro decomposition to the distribution function, which splits the original kinetic equation following a reshuffled Hilbert expansion. As opposed to classical diffusion limit, a major difficulty comes from the fat tail in the equilibrium which makes the truncation in velocity space depending on the small parameter. Our idea is, while solving the macro-micro part in a truncated velocity domain (truncation only depends on numerical accuracy), to incorporate an integrated tail over the velocity space that is beyond the truncation, and its major component can be precomputed once with any accuracy. Such an addition is essential to drive the solution to the correct asymptotic limit. Numerical experiments validate its efficiency in both kinetic and fractional diffusive regimes.

  13. A Steady-State Approximation to the Two-Dimensional Master Equation for Chemical Kinetics Calculations.

    PubMed

    Nguyen, Thanh Lam; Stanton, John F

    2015-07-16

    In the field of chemical kinetics, the solution of a two-dimensional master equation that depends explicitly on both total internal energy (E) and total angular momentum (J) is a challenging problem. In this work, a weak-E/fixed-J collisional model (i.e., weak-collisional internal energy relaxation/free-collisional angular momentum relaxation) is used along with the steady-state approach to solve the resulting (simplified) two-dimensional (E,J)-grained master equation. The corresponding solutions give thermal rate constants and product branching ratios as functions of both temperature and pressure. We also have developed a program that can be used to predict and analyze experimental chemical kinetics results. This expedient technique, when combined with highly accurate potential energy surfaces, is cable of providing results that may be meaningfully compared to experiments. The reaction of singlet oxygen with methane proceeding through vibrationally excited methanol is used as an illustrative example. PMID:25815602

  14. Analysis of atmospheric flow over a surface protrusion using the turbulence kinetic energy equation

    NASA Technical Reports Server (NTRS)

    Frost, W.; Harper, W. L.; Fichtl, G. H.

    1975-01-01

    Atmospheric flow fields resulting from a semi-elliptical surface obstruction in an otherwise horizontally homogeneous statistically stationary flow are modelled with the boundary-layer/Boussinesq-approximation of the governing equation of fluid mechanics. The turbulence kinetic energy equation is used to determine the dissipative effects of turbulent shear on the mean flow. Mean-flow results are compared with those given in a previous paper where the same problem was attacked using a Prandtl mixing-length hypothesis. Iso-lines of turbulence kinetic energy and turbulence intensity are plotted in the plane of the flow. They highlight regions of high turbulence intensity in the stagnation zone and sharp gradients in intensity along the transition from adverse to favourable pressure gradient.

  15. A critical analysis of the accuracy of several numerical techniques for combustion kinetic rate equations

    NASA Technical Reports Server (NTRS)

    Radhadrishnan, Krishnan

    1993-01-01

    A detailed analysis of the accuracy of several techniques recently developed for integrating stiff ordinary differential equations is presented. The techniques include two general-purpose codes EPISODE and LSODE developed for an arbitrary system of ordinary differential equations, and three specialized codes CHEMEQ, CREK1D, and GCKP4 developed specifically to solve chemical kinetic rate equations. The accuracy study is made by application of these codes to two practical combustion kinetics problems. Both problems describe adiabatic, homogeneous, gas-phase chemical reactions at constant pressure, and include all three combustion regimes: induction, heat release, and equilibration. To illustrate the error variation in the different combustion regimes the species are divided into three types (reactants, intermediates, and products), and error versus time plots are presented for each species type and the temperature. These plots show that CHEMEQ is the most accurate code during induction and early heat release. During late heat release and equilibration, however, the other codes are more accurate. A single global quantity, a mean integrated root-mean-square error, that measures the average error incurred in solving the complete problem is used to compare the accuracy of the codes. Among the codes examined, LSODE is the most accurate for solving chemical kinetics problems. It is also the most efficient code, in the sense that it requires the least computational work to attain a specified accuracy level. An important finding is that use of the algebraic enthalpy conservation equation to compute the temperature can be more accurate and efficient than integrating the temperature differential equation.

  16. A coarse-grained kinetic equation for neutral particles in turbulent fusion plasmas

    SciTech Connect

    Mekkaoui, A.; Marandet, Y.; Genesio, P.; Rosato, J.; Stamm, R.; Capes, H.; Koubiti, M.; Godbert-Mouret, L.; Catoire, F.

    2012-06-15

    A coarse-grained kinetic equation for neutral particles (atoms, molecules) in magnetized fusion plasmas, valid on time scales large compared to the turbulence correlation time, is presented. This equation includes the effects of plasma density fluctuations, described by gamma statistics, on the transport of neutral particles. These effects have so far been neglected in plasma edge modeling, in spite of the fact that the amplitude of fluctuations can be of order unity. Density fluctuations are shown to have a marked effect on the screening of neutrals and on the spatial localization of the ionization source, in particular at high density. The coarse-grained equations obtained in this work are readily implemented in edge code suites currently used for fusion plasma analysis and future divertor design (ITER, DEMO).

  17. Kinetic equations for hopping transport and spin relaxation in a random magnetic field

    NASA Astrophysics Data System (ADS)

    Shumilin, A. V.; Kabanov, V. V.

    2015-07-01

    We derive the kinetic equations for a hopping transport that take into account an electron spin and the possibility of double occupation. In the Ohmic regime, the equations are reduced to the generalized Miller-Abrahams resistor network. We apply these equations to the problem of the magnetic moment relaxation due to the interaction with the random hyperfine fields. It is shown that in a wide range of parameters the relaxation rate is governed by the hops with the similar rates as spin precession frequency. It is demonstrated that at the large time scale spin relaxation is nonexponential. We argue that the nonexponential relaxation of the magnetic moment is related to the spin of electrons in the slow-relaxing traps. Interestingly, the traps can significantly influence the spin relaxation in the infinite conducting cluster at large times.

  18. Solving kinetic equations with adaptive mesh in phase space for rarefied gas dynamics and plasma physics (Invited)

    SciTech Connect

    Kolobov, Vladimir; Arslanbekov, Robert; Frolova, Anna

    2014-12-09

    The paper describes an Adaptive Mesh in Phase Space (AMPS) technique for solving kinetic equations with deterministic mesh-based methods. The AMPS technique allows automatic generation of adaptive Cartesian mesh in both physical and velocity spaces using a Tree-of-Trees data structure. We illustrate advantages of AMPS for simulations of rarefied gas dynamics and electron kinetics on low temperature plasmas. In particular, we consider formation of the velocity distribution functions in hypersonic flows, particle kinetics near oscillating boundaries, and electron kinetics in a radio-frequency sheath. AMPS provide substantial savings in computational cost and increased efficiency of the mesh-based kinetic solvers.

  19. Solving kinetic equations with adaptive mesh in phase space for rarefied gas dynamics and plasma physics (Invited)

    NASA Astrophysics Data System (ADS)

    Kolobov, Vladimir; Arslanbekov, Robert; Frolova, Anna

    2014-12-01

    The paper describes an Adaptive Mesh in Phase Space (AMPS) technique for solving kinetic equations with deterministic mesh-based methods. The AMPS technique allows automatic generation of adaptive Cartesian mesh in both physical and velocity spaces using a Tree-of-Trees data structure. We illustrate advantages of AMPS for simulations of rarefied gas dynamics and electron kinetics on low temperature plasmas. In particular, we consider formation of the velocity distribution functions in hypersonic flows, particle kinetics near oscillating boundaries, and electron kinetics in a radio-frequency sheath. AMPS provide substantial savings in computational cost and increased efficiency of the mesh-based kinetic solvers.

  20. Kinetic Thomas-Fermi solutions of the Gross-Pitaevskii equation

    NASA Astrophysics Data System (ADS)

    Ölschläger, M.; Wirth, G.; Smith, C. Morais; Hemmerich, A.

    2009-04-01

    Approximate solutions of the Gross-Pitaevskii (GP) equation, obtained upon neglection of the kinetic energy, are well known as Thomas-Fermi solutions. They are characterized by the compensation of the local potential by the collisional energy. In this article we consider exact solutions of the GP-equation with this property and definite values of the kinetic energy, which suggests the term "kinetic Thomas-Fermi" (KTF) solutions. Despite their formal simplicity, KTF-solutions can possess complex current density fields with unconventional topology. We point out that a large class of light-shift potentials gives rise to KTF-solutions. As elementary examples, we consider one-dimensional and two-dimensional optical lattice scenarios, obtained by means of the superposition of two, three and four laser beams, and discuss the stability properties of the corresponding KTF-solutions. A general method is proposed to excite two-dimensional KTF-solutions in experiments by means of time-modulated light-shift potentials.

  1. Convective kinetic energy equation under the mass-flux subgrid-scale parameterization

    NASA Astrophysics Data System (ADS)

    Yano, Jun-Ichi

    2015-03-01

    The present paper originally derives the convective kinetic energy equation under mass-flux subgrid-scale parameterization in a formal manner based on the segmentally-constant approximation (SCA). Though this equation is long since presented by Arakawa and Schubert (1974), a formal derivation is not known in the literature. The derivation of this formulation is of increasing interests in recent years due to the fact that it can explain basic aspects of the convective dynamics such as discharge-recharge and transition from shallow to deep convection. The derivation is presented in two manners: (i) for the case that only the vertical component of the velocity is considered and (ii) the case that both the horizontal and vertical components are considered. The equation reduces to the same form as originally presented by Arakwa and Schubert in both cases, but with the energy dissipation term defined differently. In both cases, nevertheless, the energy "dissipation" (loss) term consists of the three principal contributions: (i) entrainment-detrainment, (ii) outflow from top of convection, and (iii) pressure effects. Additionally, inflow from the bottom of convection contributing to a growth of convection is also formally counted as a part of the dissipation term. The eddy dissipation is also included for a completeness. The order-of-magnitude analysis shows that the convective kinetic energy "dissipation" is dominated by the pressure effects, and it may be approximately described by Rayleigh damping with a constant time scale of the order of 102-103 s. The conclusion is also supported by a supplementary analysis of a cloud-resolving model (CRM) simulation. The Appendix discusses how the loss term ("dissipation") of the convective kinetic energy is qualitatively different from the conventional eddy-dissipation process found in turbulent flows.

  2. An asymptotic preserving unified gas kinetic scheme for gray radiative transfer equations

    SciTech Connect

    Sun, Wenjun; Jiang, Song; Xu, Kun

    2015-03-15

    The solutions of radiative transport equations can cover both optical thin and optical thick regimes due to the large variation of photon's mean-free path and its interaction with the material. In the small mean free path limit, the nonlinear time-dependent radiative transfer equations can converge to an equilibrium diffusion equation due to the intensive interaction between radiation and material. In the optical thin limit, the photon free transport mechanism will emerge. In this paper, we are going to develop an accurate and robust asymptotic preserving unified gas kinetic scheme (AP-UGKS) for the gray radiative transfer equations, where the radiation transport equation is coupled with the material thermal energy equation. The current work is based on the UGKS framework for the rarefied gas dynamics [14], and is an extension of a recent work [12] from a one-dimensional linear radiation transport equation to a nonlinear two-dimensional gray radiative system. The newly developed scheme has the asymptotic preserving (AP) property in the optically thick regime in the capturing of diffusive solution without using a cell size being smaller than the photon's mean free path and time step being less than the photon collision time. Besides the diffusion limit, the scheme can capture the exact solution in the optical thin regime as well. The current scheme is a finite volume method. Due to the direct modeling for the time evolution solution of the interface radiative intensity, a smooth transition of the transport physics from optical thin to optical thick can be accurately recovered. Many numerical examples are included to validate the current approach.

  3. From Lévy flights to the fractional kinetic equation for dynamical chaos

    NASA Astrophysics Data System (ADS)

    Zaslavsky, G. M.

    Chaotic dynamics of Hamiltonian systems can be described by the random process which resembles the Lévy-type flights and trappings in the phase space of a system. The probability distribution function satisfies the fractional in space and time generalization of the Fokker-Planck-Kolmogorov equation. Orders of the fractional derivatives in space and time can be connected to the Pesin's dimensions of the trajectories. A new look on the problem of Maxwell's Demon is discussed in the context of the anomalous ("strange") kinetics.

  4. Simulation of light generation in cholesteric liquid crystals using kinetic equations: Time-independent solution

    SciTech Connect

    Shtykov, N. M. Palto, S. P.; Umanskii, B. A.

    2013-08-15

    We report on the results of calculating the conditions for light generation in cholesteric liquid crystals doped with fluorescent dyes using kinetic equations. Specific features of spectral properties of the chiral cholesteric medium as a photonic structure and spatially distributed type of the feedback in the active medium are taken into account. The expression is derived for the threshold pump radiation intensity as a function of the dye concentration and sample thickness. The importance of taking into account the distributed loss level in the active medium for calculating the optimal parameters of the medium and for matching the calculated values with the results of experiments is demonstrated.

  5. Theory of warm ionized gases: equation of state and kinetic Schottky anomaly.

    PubMed

    Capolupo, A; Giampaolo, S M; Illuminati, F

    2013-10-01

    Based on accurate Lennard-Jones-type interaction potentials, we derive a closed set of state equations for the description of warm atomic gases in the presence of ionization processes. The specific heat is predicted to exhibit peaks in correspondence to single and multiple ionizations. Such kinetic analog in atomic gases of the Schottky anomaly in solids is enhanced at intermediate and low atomic densities. The case of adiabatic compression of noble gases is analyzed in detail and the implications on sonoluminescence are discussed. In particular, the predicted plasma electron density in a sonoluminescent bubble turns out to be in good agreement with the value measured in recent experiments. PMID:24229140

  6. Application of convergence acceleration to the reactor kinetic equations: A comparative study

    SciTech Connect

    Picca, P.; Furfaro, R.; Ganapol, B. D.

    2013-07-01

    This presentation provides a comparison of two methodologies for the solution of reactor kinetic equations, namely for a standard finite difference and a semi-analytical approach. The above-mentioned methods are implemented in a convergence acceleration framework to enhance their efficiency and a comparative study is reported to verify whether it is more convenient to use a rudimentary but fast algorithm (finite difference) with respect to the more refined but computationally intense approach of the semi-analytical method. Performance on several test cases from the literature are compared. (authors)

  7. Berry curvature and four-dimensional monopoles in the relativistic chiral kinetic equation.

    PubMed

    Chen, Jiunn-Wei; Pu, Shi; Wang, Qun; Wang, Xin-Nian

    2013-06-28

    We derive a relativistic chiral kinetic equation with manifest Lorentz covariance from Wigner functions of spin-1/2 massless fermions in a constant background electromagnetic field. It contains vorticity terms and a four-dimensional Euclidean Berry monopole which gives an axial anomaly. By integrating out the zeroth component of the 4-momentum p, we reproduce the previous three-dimensional results derived from the Hamiltonian approach, together with the newly derived vorticity terms. The phase space continuity equation has an anomalous source term proportional to the product of electric and magnetic fields (FσρF[over ˜]σρ∼EσBσ). This provides a unified interpretation of the chiral magnetic and vortical effects, chiral anomaly, Berry curvature, and the Berry monopole in the framework of Wigner functions. PMID:23848865

  8. On Self-Similar Solutions to a Kinetic Equation Arising in Weak Turbulence Theory for the Nonlinear Schrödinger Equation

    NASA Astrophysics Data System (ADS)

    Kierkels, A. H. M.; Velázquez, J. J. L.

    2016-06-01

    We construct a family of self-similar solutions with fat tails to a quadratic kinetic equation. This equation describes the long time behaviour of weak solutions with finite mass to the weak turbulence equation associated to the nonlinear Schrödinger equation. The solutions that we construct have finite mass, but infinite energy. In Kierkels and Velázquez (J Stat Phys 159:668-712, 2015) self-similar solutions with finite mass and energy were constructed. Here we prove upper and lower exponential bounds on the tails of these solutions.

  9. On Self-Similar Solutions to a Kinetic Equation Arising in Weak Turbulence Theory for the Nonlinear Schrödinger Equation

    NASA Astrophysics Data System (ADS)

    Kierkels, A. H. M.; Velázquez, J. J. L.

    2016-04-01

    We construct a family of self-similar solutions with fat tails to a quadratic kinetic equation. This equation describes the long time behaviour of weak solutions with finite mass to the weak turbulence equation associated to the nonlinear Schrödinger equation. The solutions that we construct have finite mass, but infinite energy. In Kierkels and Velázquez (J Stat Phys 159:668-712, 2015) self-similar solutions with finite mass and energy were constructed. Here we prove upper and lower exponential bounds on the tails of these solutions.

  10. Improved Experimental and Computational Methodology for Determining the Kinetic Equation and the Extant Kinetic Constants of Fe(II) Oxidation by Acidithiobacillus ferrooxidans▿

    PubMed Central

    Molchanov, Sharon; Gendel, Yuri; Ioslvich, Ilya; Lahav, Ori

    2007-01-01

    The variety of kinetics expressions encountered in the literature and the unreasonably broad range of values reported for the kinetics constants of Acidithiobacillus ferrooxidans underscore the need for a unifying experimental procedure and for the development of a reliable kinetics equation. Following an extensive and critical review of reported experimental techniques, a method based on batch pH-controlled kinetics experiments lasting less than one doubling time was developed for the determination of extant kinetics constants. The Fe(II) concentration in the experiments was measured by a method insensitive to Fe(III) interference. Kinetics parameters were determined by nonlinear fitting of the integrated form of the Monod equation to yield a KS of 31 ± 4 mg Fe2+ liter−1 (mean ± standard deviation), a KP of 139 ± 20 mg Fe3+ liter−1, and a μmax of 0.082 ± 0.002 h−1. The corresponding kinetics equation was as follows: \\documentclass[10pt]{article} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{pmc} \\pagestyle{empty} \\oddsidemargin -1.0in \\begin{document} \\begin{equation*}\\frac{dS}{dt}= \\left \\left(-\\frac{0.082}{2.3{\\cdot}10^{7}}\\right) \\right \\frac{S{\\cdot}X}{31(1+\\displaystyle\\frac{P_{0}+S_{0}-S}{139})+S}\\end{equation*}\\end{document} where S represents the Fe(II) concentration in mg liter−1, P0 represents the initial Fe(III) concentration in mg liter−1, X represents the suspended bacterial cell concentration in cells ml−1, and t represents time in hours. The measured data fit this equation exceptionally well, with an R2 of >0.99. Fe(III) inhibition was found to be of a competitive nature. Contrary to previous reports, the results show that the concentration of Acidithiobacillus ferrooxidans cells has no affect on the kinetics constants. The kinetics equation can be considered applicable only to A. ferrooxidans cells grown under

  11. The solution of the point kinetics equations via converged accelerated Taylor series (CATS)

    SciTech Connect

    Ganapol, B.; Picca, P.; Previti, A.; Mostacci, D.

    2012-07-01

    This paper deals with finding accurate solutions of the point kinetics equations including non-linear feedback, in a fast, efficient and straightforward way. A truncated Taylor series is coupled to continuous analytical continuation to provide the recurrence relations to solve the ordinary differential equations of point kinetics. Non-linear (Wynn-epsilon) and linear (Romberg) convergence accelerations are employed to provide highly accurate results for the evaluation of Taylor series expansions and extrapolated values of neutron and precursor densities at desired edits. The proposed Converged Accelerated Taylor Series, or CATS, algorithm automatically performs successive mesh refinements until the desired accuracy is obtained, making use of the intermediate results for converged initial values at each interval. Numerical performance is evaluated using case studies available from the literature. Nearly perfect agreement is found with the literature results generally considered most accurate. Benchmark quality results are reported for several cases of interest including step, ramp, zigzag and sinusoidal prescribed insertions and insertions with adiabatic Doppler feedback. A larger than usual (9) number of digits is included to encourage honest benchmarking. The benchmark is then applied to the enhanced piecewise constant algorithm (EPCA) currently being developed by the second author. (authors)

  12. Analysis of the methods for the derivation of binary kinetic equations in the theory of fluorescence concentration quenching

    SciTech Connect

    Doktorov, A. B.

    2014-09-14

    In the framework of unified many-particle approach the familiar problem of fluorescence concentration quenching in the presence of pumping (light pulse) of arbitrary intensity is considered. This process is a vivid and the simplest example of multistage bulk reaction including bimolecular irreversible quenching reaction and reversible monomolecular transformation as elementary stages. General relation between the kinetics of multistage bulk reaction and that of the elementary stage of quenching has been established. This allows one to derive general kinetic equations (of two types) for the multistage reaction in question on the basis of general kinetic equations (differential and integro-differential) of elementary stage of quenching. Relying on the same unified many-particle approach we have developed binary approximations with the use of two (frequently employed in the literature) many-particle methods (such as simple superposition approximation and the method of extracting pair channels in three-particle correlation evolution) to the derivation of non-Markovian binary kinetic equations. The possibility of reducing the obtained binary equations to the Markovian equations of formal chemical kinetics has been considered. As an example the exact solution of the problem (for the specific case) is examined, and the applicability of two many particle methods of derivation of binary equations is analyzed.

  13. Kinetic flux-vector splitting for the Navier-Stokes equations

    SciTech Connect

    Chou, S.Y.; Baganoff, D.

    1997-01-15

    Before a hybrid scheme can be developed combining the direct simulation Monte Carlo (DSMC) method and a Navier-Stokes (NS) representation, one must have access to compatible kinetic-split fluxes from the NS portion of the hybrid scheme. The kinetic theory basis is given for the development of the required fluxes from the Chapman-Enskog velocity distribution function for a simple gas; and these are then extended to a polyatomic gas by use of the Eucken approximation. The derived fluxes are then used to implement boundary conditions at solid surfaces that are based on concepts associated with kinetic theory and the DSMC method. This approach is shown to lead to temperature slip and velocity slip as a natural outcome of the new formulation, a requirement for use in the near-continuum regime where DSMC and NS must be joined. Several different flows, for which solid boundaries are not present, are computed using the derived fluxes, together with a second-order finite-volume scheme, and the results are shown to agree well with several established numerical schemes for the NS equations. 22 refs., 12 figs.

  14. Kinetic boundary layers in gas mixtures: Systems described by nonlinearly coupled kinetic and hydrodynamic equations and applications to droplet condensation and evaporation

    SciTech Connect

    Widder, M.E.; Titulaer, U.M. )

    1993-03-01

    The authors consider a mixture of heavy vapor molecules and a light carrier gas surrounding a liquid droplet. The vapor is described by a variant of the Klein-Kramers equation; the gas is described by the Navier-Stokes equations; the droplet acts as a heat source due to the released heat of condensation. The exchange of momentum and energy between the constituents of the mixture is taken into account by force terms in the kinetic equation and source terms in the Navier-Stokes equations. These are chosen to obtain maximal agreement with the irreversible thermodynamics of a gas mixture. The structure of the kinetic boundary layer around the sphere is determined from the self-consistent solution of this set of coupled equations with appropriate boundary conditions at the surface of the sphere. The kinetic equation is rewritten as a set of coupled moment equations. A complete set of solutions of these moment equations is constructed by numerical integration inward from the region far away from the droplet, where the background inhomogeneities are small. A technique developed earlier is used to deal with the numerical instability of the moment equations. The solutions obtained for given temperature and pressure profiles in the gas are then combined linearly such that they obey the boundary conditions at the droplet surface; from this solution source terms for the Navier-Stokes equation of the gas are constructed and used to determine improved temperature and pressure profiles for the background gas. For not too large temperature differneces between the droplet and the gas at infinity, self-consistency is reached after a few iterations. The method is applied to the condensation of droplets from a supersaturated vapor as well as to strong evaporation of droplets under the influence of an external heat source, where corrections of up to 40% are obtained.

  15. Coalescence kinetics in surfactant stabilized emulsions: evolution equations from direct numerical simulations.

    PubMed

    Skartlien, R; Grimes, B; Meakin, P; Sjöblom, J; Sollum, E

    2012-12-01

    Lattice Boltzmann simulations were used to study the coalescence kinetics in emulsions with amphiphilic surfactant, under neutrally buoyant conditions, and with a significant kinematic viscosity contrast between the phases (emulating water in oil emulsions). The 3D simulation domain was large enough (256(3) ~ 10(7) grid points) to obtain good statistics with droplet numbers ranging from a few thousand at early times to a few hundred near equilibrium. Increased surfactant contents slowed down the coalescence rate between droplets due to the Gibbs-Marangoni effect, and the coalescence was driven by a quasi-turbulent velocity field. The kinetic energy decayed at a relatively slow rate at early times, due to conversion of interfacial energy to kinetic energy in the flow during coalescence. Phenomenological, coupled differential equations for the mean droplet diameter D(t) and the number density n(d)(t) were obtained from the simulation data and from film draining theories. Local (in time) power law exponents for the growth of the mean diameter (and for the concomitant decrease of n(d)) were established in terms of the instantaneous values of the kinetic energy, coalescence probability, Gibbs elasticity, and interfacial area. The model studies indicated that true power laws for the growth of the droplet size and decrease of the number of droplets with time may not be justified, since the exponents derived using the phenomenological model were time dependent. In contrast to earlier simulation results for symmetric blends with surfactant, we found no evidence for stretched logarithmic scaling of the form D ~ [ln (ct)](α) for the morphology length, or exponential scalings associated with arrested growth, on the basis of the phenomenological model. PMID:23231250

  16. Coalescence kinetics in surfactant stabilized emulsions: Evolution equations from direct numerical simulations

    SciTech Connect

    R. Skartlien; E. Sollum; A. Akselsen; P. Meakin; B. Grimes; J. Sjoblom

    2012-12-01

    Lattice Boltzmann simulations were used to study the coalescence kinetics in emulsions with amphiphilic surfactant, under neutrally buoyant conditions, and with a significant kinematic viscosity contrast between the phases (emulating water in oil emulsions). The 3D simulation domain was large enough (256 3rd power -- 10 7th power grid points) to obtain good statistics with droplet numbers ranging from a few thousand at early times to a few hundred near equilibrium. Increased surfactant contents slowed down the coalescence rate between droplets due to the Gibbs-Marangoni effect, and the coalescence was driven by a quasi-turbulent velocity field. The kinetic energy decayed at a relatively slow rate at early times, due to conversion of interfacial energy to kinetic energy in the flow during coalescence. Phenomenological, coupled differential equations for the mean droplet diameter D(t) and the number density nd(t) were obtained from the simulation data and from film draining theories. Local (in time) power law exponents for the growth of the mean diameter (and for the concomitant decrease of nd) were established in terms of the instantaneous values of the kinetic energy, coalescence probability, Gibbs elasticity, and interfacial area. The model studies indicated that true power laws for the growth of the droplet size and decrease of the number of droplets with time may not be justified, since the exponents derived using the phenomenological model were time dependent. In contrast to earlier simulation results for symmetric blends with surfactant, we found no evidence for stretched logarithmic scaling of the formD -- [ln (ct)]a for the morphology length, or exponential scalings associated with arrested growth, on the basis of the phenomenological model.

  17. Elementary solutions of coupled model equations in the kinetic theory of gases

    NASA Technical Reports Server (NTRS)

    Kriese, J. T.; Siewert, C. E.; Chang, T. S.

    1974-01-01

    The method of elementary solutions is employed to solve two coupled integrodifferential equations sufficient for determining temperature-density effects in a linearized BGK model in the kinetic theory of gases. Full-range completeness and orthogonality theorems are proved for the developed normal modes and the infinite-medium Green's function is constructed as an illustration of the full-range formalism. The appropriate homogeneous matrix Riemann problem is discussed, and half-range completeness and orthogonality theorems are proved for a certain subset of the normal modes. The required existence and uniqueness theorems relevant to the H matrix, basic to the half-range analysis, are proved, and an accurate and efficient computational method is discussed. The half-space temperature-slip problem is solved analytically, and a highly accurate value of the temperature-slip coefficient is reported.

  18. A robust low diffusive kinetic scheme for the Navier-Stokes/Euler equations

    SciTech Connect

    Moschetta, J.M.; Pullin, D.I.

    1997-05-15

    A new kinetic scheme based on the equilibrium flux method (EFM) and modified using Osher intermediate states is proposed. This new scheme called EFMO combines the robustness of the equilibrium flux method and the accuracy of flux-difference splitting schemes. The original EFM scheme is expressed in terms of simple wave decomposition in which only the linearly degenerate subpath is calculated from Osher numerical flux while nonlinear waves are still evaluated from the regular EFM splitting. Owing to its capability of withstanding intense nonlinear waves and yet exactly resolving contact discontinuities, EFMO is particularly well suited for the resolution of the Navier-Stokes equations as demonstrated by a series of severe test cases including the high-speed viscous flow around a cone, a shock-boundary layer interaction problem, a vacuum apparition problem, the hypersonic flow around a circular cylinder at Mach 100, and the forward-facing step at Mach 3. 36 refs., 7 figs.

  19. Approximate method for stochastic chemical kinetics with two-time scales by chemical Langevin equations

    NASA Astrophysics Data System (ADS)

    Wu, Fuke; Tian, Tianhai; Rawlings, James B.; Yin, George

    2016-05-01

    The frequently used reduction technique is based on the chemical master equation for stochastic chemical kinetics with two-time scales, which yields the modified stochastic simulation algorithm (SSA). For the chemical reaction processes involving a large number of molecular species and reactions, the collection of slow reactions may still include a large number of molecular species and reactions. Consequently, the SSA is still computationally expensive. Because the chemical Langevin equations (CLEs) can effectively work for a large number of molecular species and reactions, this paper develops a reduction method based on the CLE by the stochastic averaging principle developed in the work of Khasminskii and Yin [SIAM J. Appl. Math. 56, 1766-1793 (1996); ibid. 56, 1794-1819 (1996)] to average out the fast-reacting variables. This reduction method leads to a limit averaging system, which is an approximation of the slow reactions. Because in the stochastic chemical kinetics, the CLE is seen as the approximation of the SSA, the limit averaging system can be treated as the approximation of the slow reactions. As an application, we examine the reduction of computation complexity for the gene regulatory networks with two-time scales driven by intrinsic noise. For linear and nonlinear protein production functions, the simulations show that the sample average (expectation) of the limit averaging system is close to that of the slow-reaction process based on the SSA. It demonstrates that the limit averaging system is an efficient approximation of the slow-reaction process in the sense of the weak convergence.

  20. Approximate method for stochastic chemical kinetics with two-time scales by chemical Langevin equations.

    PubMed

    Wu, Fuke; Tian, Tianhai; Rawlings, James B; Yin, George

    2016-05-01

    The frequently used reduction technique is based on the chemical master equation for stochastic chemical kinetics with two-time scales, which yields the modified stochastic simulation algorithm (SSA). For the chemical reaction processes involving a large number of molecular species and reactions, the collection of slow reactions may still include a large number of molecular species and reactions. Consequently, the SSA is still computationally expensive. Because the chemical Langevin equations (CLEs) can effectively work for a large number of molecular species and reactions, this paper develops a reduction method based on the CLE by the stochastic averaging principle developed in the work of Khasminskii and Yin [SIAM J. Appl. Math. 56, 1766-1793 (1996); ibid. 56, 1794-1819 (1996)] to average out the fast-reacting variables. This reduction method leads to a limit averaging system, which is an approximation of the slow reactions. Because in the stochastic chemical kinetics, the CLE is seen as the approximation of the SSA, the limit averaging system can be treated as the approximation of the slow reactions. As an application, we examine the reduction of computation complexity for the gene regulatory networks with two-time scales driven by intrinsic noise. For linear and nonlinear protein production functions, the simulations show that the sample average (expectation) of the limit averaging system is close to that of the slow-reaction process based on the SSA. It demonstrates that the limit averaging system is an efficient approximation of the slow-reaction process in the sense of the weak convergence. PMID:27155630

  1. Conditions for critical effects in the mass action kinetics equations for water radiolysis

    SciTech Connect

    Wittman, Richard S.; Buck, Edgar C.; Mausolf, Edward J.; McNamara, Bruce K.; Smith, Frances N.; Soderquist, Chuck Z.

    2014-11-25

    We report on a subtle global feature of the mass action kinetics equations for water radiolysis that results in predictions of a critical behavior in H2O2 and associated radical concentrations. While radiolysis kinetics has been studied extensively in the past, it is only in recent years that high speed computing has allowed the rapid exploration of the solution over widely varying dose and compositional conditions. We explore the radiolytic production of H2O2 under various externally fixed conditions of molecular H2 and O2 that have been regarded as problematic in the literature – specifically, “jumps” in predicted concentrations, and inconsistencies between predictions and experiments have been reported for alpha radiolysis. We computationally map-out a critical concentration behavior for alpha radiolysis kinetics using a comprehensive set of reactions. We then show that all features of interest are accurately reproduced with 15 reactions. An analytical solution for steady-state concentrations of the 15 reactions reveals regions in [H2] and [O2] where the H2O2 concentration is not unique – both stable and unstable concentrations exist. The boundary of this region can be characterized analytically as a function of G-values and rate constants independent of dose rate. Physically, the boundary can be understood as separating a region where a steady-state H2O2 concentration exists, from one where it does not exist without a direct decomposition reaction. We show that this behavior is consistent with reported alpha radiolysis data and that no such behavior should occur for gamma radiolysis. We suggest experiments that could verify or discredit a critical concentration behavior for alpha radiolysis and could place more restrictive ranges on G-values from derived relationships between them.

  2. Conditions for critical effects in the mass action kinetics equations for water radiolysis

    SciTech Connect

    Wittman, Richard S.; Buck, Edgar C.; Mausolf, Edward J.; McNamara, Bruce K.; Smith, Frances N.; Soderquist, Chuck Z.

    2014-12-26

    We report on a subtle global feature of the mass action kinetics equations for water radiolysis that results in predictions of a critical behavior in H2O2 and associated radical concentrations. While radiolysis kinetics has been studied extensively in the past, it is only in recent years that high speed computing has allowed the rapid exploration of the solution over widely varying dose and compositional conditions. We explore the radiolytic production of H2O2 under various externally fixed conditions of molecular H2 and O2 that have been regarded as problematic in the literature – specifically, “jumps” in predicted concentrations, and inconsistencies between predictions and experiments have been reported for alpha radiolysis. We computationally map-out a critical concentration behavior for alpha radiolysis kinetics using a comprehensive set of reactions. We then show that all features of interest are accurately reproduced with 15 reactions. An analytical solution for steady-state concentrations of the 15 reactions reveals regions in [H2] and [O2] where the H2O2 concentration is not unique – both stable and unstable concentrations exist. The boundary of this region can be characterized analytically as a function of G-values and rate constants independent of dose rate. Physically, the boundary can be understood as separating a region where a steady-state H2O2 concentration exists, from one where it does not exist without a direct decomposition reaction. We show that this behavior is consistent with reported alpha radiolysis data and that no such behavior should occur for gamma radiolysis. We suggest experiments that could verify or discredit a critical concentration behavior for alpha radiolysis and could place more restrictive ranges on G-values from derived relationships between them.

  3. Verification of particle-in-cell simulations against exact solutions of kinetic equations

    NASA Astrophysics Data System (ADS)

    Turner, Miles

    2015-09-01

    Demonstrating correctness of computer simulations (or verification) has become a matter of increasing concern in recent years. The strongest type of verification is a demonstration that the simulation converges to an exact solution of the mathematical model that is supposed to be solved. Of course, this is possible only if such an exact solution is available. In this paper, we are interested in kinetic simulation using the particle-in-cell method, and consequently a relevant exact solution must be a solution of a kinetic equation. While we know of no such solutions that exercise all the features of a typical particle-in-cell simulation, in this paper we show that the mathematical literature contains several such solutions that involve a large fraction of the functionality of such a code, and which collectively exercise essentially all of the code functionality. These solutions include the plane diode, the neutron criticality problem, and the calculation of ion energy distribution functions in oscillating fields. In each of theses cases, we can show the the particle-in-cell simulation converges to the exact solution in the expected way. These demonstrations are strong evidence of correct implementation. Work supported by Science Foundation Ireland under grant 08/SRC/I1411.

  4. Padé Approximants for the Equation of State for Relativistic Hydrodynamics by Kinetic Theory

    NASA Astrophysics Data System (ADS)

    Tsai, Shang-Hsi; Yang, Jaw-Yen

    2015-07-01

    A two-point Padé approximant (TPPA) algorithm is developed for the equation of state (EOS) for relativistic hydrodynamic systems, which are described by the classical Maxwell-Boltzmann statistics and the semiclassical Fermi-Dirac statistics with complete degeneracy. The underlying rational function is determined by the ratios of the macroscopic state variables with various orders of accuracy taken at the extreme relativistic limits. The nonunique TPPAs are validated by Taub's inequality for the consistency of the kinetic theory and the special theory of relativity. The proposed TPPA is utilized in deriving the EOS of the dilute gas and in calculating the specific heat capacity, the adiabatic index function, and the isentropic sound speed of the ideal gas. Some general guidelines are provided for the application of an arbitrary accuracy requirement. The superiority of the proposed TPPA is manifested in manipulating the constituent polynomials of the approximants, which avoids the arithmetic complexity of struggling with the modified Bessel functions and the hyperbolic trigonometric functions arising from the relativistic kinetic theory.

  5. High order asymptotic preserving DG-IMEX schemes for discrete-velocity kinetic equations in a diffusive scaling

    NASA Astrophysics Data System (ADS)

    Jang, Juhi; Li, Fengyan; Qiu, Jing-Mei; Xiong, Tao

    2015-01-01

    In this paper, we develop a family of high order asymptotic preserving schemes for some discrete-velocity kinetic equations under a diffusive scaling, that in the asymptotic limit lead to macroscopic models such as the heat equation, the porous media equation, the advection-diffusion equation, and the viscous Burgers' equation. Our approach is based on the micro-macro reformulation of the kinetic equation which involves a natural decomposition of the equation to the equilibrium and non-equilibrium parts. To achieve high order accuracy and uniform stability as well as to capture the correct asymptotic limit, two new ingredients are employed in the proposed methods: discontinuous Galerkin (DG) spatial discretization of arbitrary order of accuracy with suitable numerical fluxes; high order globally stiffly accurate implicit-explicit (IMEX) Runge-Kutta scheme in time equipped with a properly chosen implicit-explicit strategy. Formal asymptotic analysis shows that the proposed scheme in the limit of ε → 0 is a consistent high order discretization for the limiting equation. Numerical results are presented to demonstrate the stability and high order accuracy of the proposed schemes together with their performance in the limit. Our methods are also tested for the continuous-velocity one-group transport equation in slab geometry and for several examples with spatially varying parameters.

  6. Thermal electron heating rate - A derivation. [from kinetic equation for earth ionosphere

    NASA Technical Reports Server (NTRS)

    Hoegy, W. R.

    1984-01-01

    The thermal electron heating rate is an important heat source term in the ionospheric electron energy balance equation, representing heating by photoelectrons or by precipitating higher energy electrons. A formula for the thermal electron heating rate is derived from the kinetic equation using the electron-electron collision operator as given by the unified theory of Kihara and Aono. This collision operator includes collective interactions to produce a finite collision operator with an exact Coulomb logarithm term. The derived heating rate O(e) is the sum of three terms, O(e) = O(P) + S + O(int), which are respectively: (1) primary electron production term giving the heating from newly created electrons that have not yet suffered collisions with the ambient electrons; (2) a heating term evaluated on the energy surface m(e)/2 = E(T) at the transition between Maxwellian and tail electrons at E(T); and (3) the integral term representing heating of Maxwellian electrons by eneegetic tail electrons at energies ET. Published ionospheric electron temperature studies used only the integral term O(int) with differing lower integration limits. Use of the incomplete heating rate could lead to erroneous conclusions regarding electron heat balance, since O(e) is greater than O(int) by as much as a factor of two. Previously announced in STAR as N84-15941

  7. Gas-kinetic unified algorithm for hypersonic flows covering various flow regimes solving Boltzmann model equation in nonequilibrium effect

    SciTech Connect

    Li, Zhihui; Ma, Qiang; Wu, Junlin; Jiang, Xinyu; Zhang, Hanxin

    2014-12-09

    Based on the Gas-Kinetic Unified Algorithm (GKUA) directly solving the Boltzmann model equation, the effect of rotational non-equilibrium is investigated recurring to the kinetic Rykov model with relaxation property of rotational degrees of freedom. The spin movement of diatomic molecule is described by moment of inertia, and the conservation of total angle momentum is taken as a new Boltzmann collision invariant. The molecular velocity distribution function is integrated by the weight factor on the internal energy, and the closed system of two kinetic controlling equations is obtained with inelastic and elastic collisions. The optimization selection technique of discrete velocity ordinate points and numerical quadrature rules for macroscopic flow variables with dynamic updating evolvement are developed to simulate hypersonic flows, and the gas-kinetic numerical scheme is constructed to capture the time evolution of the discretized velocity distribution functions. The gas-kinetic boundary conditions in thermodynamic non-equilibrium and numerical procedures are studied and implemented by directly acting on the velocity distribution function, and then the unified algorithm of Boltzmann model equation involving non-equilibrium effect is presented for the whole range of flow regimes. The hypersonic flows involving non-equilibrium effect are numerically simulated including the inner flows of shock wave structures in nitrogen with different Mach numbers of 1.5-Ma-25, the planar ramp flow with the whole range of Knudsen numbers of 0.0009-Kn-10 and the three-dimensional re-entering flows around tine double-cone body.

  8. Formal analogy between the Dirac equation in its Majorana form and the discrete-velocity version of the Boltzmann kinetic equation.

    PubMed

    Fillion-Gourdeau, F; Herrmann, H J; Mendoza, M; Palpacelli, S; Succi, S

    2013-10-18

    We point out a formal analogy between the Dirac equation in Majorana form and the discrete-velocity version of the Boltzmann kinetic equation. By a systematic analysis based on the theory of operator splitting, this analogy is shown to turn into a concrete and efficient computational method, providing a unified treatment of relativistic and nonrelativistic quantum mechanics. This might have potentially far-reaching implications for both classical and quantum computing, because it shows that, by splitting time along the three spatial directions, quantum information (Dirac-Majorana wave function) propagates in space-time as a classical statistical process (Boltzmann distribution). PMID:24182245

  9. The generalized kinetic equation as a model for the nonlinear transfer in third-generation wave models

    NASA Astrophysics Data System (ADS)

    Gramstad, Odin; Babanin, Alexander

    2016-04-01

    An alternative model for the nonlinear interaction term S n l in spectral wave models, the so called generalized kinetic equation (Janssen J Phys Oceanogr 33(4):863-884, 2003; Annenkov and Shrira J Fluid Mech 561:181-207, 2006b; Gramstad and Stiassnie J Fluid Mech 718:280-303, 2013), is discussed and implemented in the third generation wave model WAVEWATCH-III. The generalized kinetic equation includes the effects of near-resonant nonlinear interactions, and is therefore able, in theory, to describe faster nonlinear evolution than the existing forms of S n l which are based on the standard Hasselmann kinetic equation (Hasselmann J Fluid Mech 12:481-500, 1962). Numerical simulations with WAVEWATCH have been carried out to thoroughly test the performance of the new form of S n l , and to compare it to the existing models for S n l in WAVEWATCH; the DIA and WRT. Some differences between the different models for S n l are observed. As expected, the DIA is shown to perform less well compared to the exact terms in certain situations, in particular for narrow wave spectra. Also for the case of turning wind significant differences between the different models are observed. Nevertheless, different from the case of unidirectional waves where the generalized kinetic equation represents a obvious improvement to the standard forms of S n l (Gramstad and Stiassnie 2013), the differences seems to be less pronounced for the more realistic cases considered in this paper.

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

  11. Numerical solutions of the semiclassical Boltzmann ellipsoidal-statistical kinetic model equation

    PubMed Central

    Yang, Jaw-Yen; Yan, Chin-Yuan; Huang, Juan-Chen; Li, Zhihui

    2014-01-01

    Computations of rarefied gas dynamical flows governed by the semiclassical Boltzmann ellipsoidal-statistical (ES) kinetic model equation using an accurate numerical method are presented. The semiclassical ES model was derived through the maximum entropy principle and conserves not only the mass, momentum and energy, but also contains additional higher order moments that differ from the standard quantum distributions. A different decoding procedure to obtain the necessary parameters for determining the ES distribution is also devised. The numerical method in phase space combines the discrete-ordinate method in momentum space and the high-resolution shock capturing method in physical space. Numerical solutions of two-dimensional Riemann problems for two configurations covering various degrees of rarefaction are presented and various contours of the quantities unique to this new model are illustrated. When the relaxation time becomes very small, the main flow features a display similar to that of ideal quantum gas dynamics, and the present solutions are found to be consistent with existing calculations for classical gas. The effect of a parameter that permits an adjustable Prandtl number in the flow is also studied. PMID:25104904

  12. Rate Equations and Kinetic Parameters of the Reactions Involved in Pyrite Oxidation by Thiobacillus ferrooxidans.

    PubMed

    Lizama, H M; Suzuki, I

    1989-11-01

    Rate equations and kinetic parameters were obtained for various reactions involved in the bacterial oxidation of pyrite. The rate constants were 3.5 muM Fe per min per FeS(2) percent pulp density for the spontaneous pyrite dissolution, 10 muM Fe per min per mM Fe for the indirect leaching with Fe, 90 muM O(2) per min per mg of wet cells per ml for the Thiobacillus ferrooxidans oxidation of washed pyrite, and 250 muM O(2) per min per mg of wet cells per ml for the T. ferrooxidans oxidation of unwashed pyrite. The K(m) values for pyrite concentration were similar and were 1.9, 2.5, and 2.75% pulp density for indirect leaching, washed pyrite oxidation by T. ferrooxidans, and unwashed pyrite oxidation by T. ferrooxidans, respectively. The last reaction was competitively inhibited by increasing concentrations of cells, with a K(i) value of 0.13 mg of wet cells per ml. T. ferrooxidans cells also increased the rate of Fe production from Fe plus pyrite. PMID:16348054

  13. A kinetic scheme for the Navier-Stokes equations and high-order methods for hyperbolic conservation laws

    NASA Astrophysics Data System (ADS)

    May, Georg

    This dissertation revolves around algorithm development in the context of numerical methods for hyperbolic conservation laws and the compressible Navier-Stokes equations, with particular emphasis on unstructured meshes. Three distinct topics may be identified: Firstly, a new kinetic scheme for the compressible Navier-Stokes equations is developed. Kinetic numerical schemes are based on the discretization of a probability density function. In the context of fluid flow such schemes have a natural basis rooted in the kinetic theory of gases. A significant advantage of kinetic schemes is that they allow a compact, completely mesh-independent discretization of the Navier-Stokes equations, which makes them well suited for next-generation solvers on general unstructured meshes. The new kinetic scheme is based on the Xu-Prenderaast BGK scheme, and achieves a dramatic reduction in computational cost, while also improving and clarifying the formulation with respect to the underlying kinetic gas theory. The second topic addresses high-order numerical methods for conservation laws on unstructured meshes. High-order methods potentially produce higher accuracy with fewer degrees of freedom, compared to standard first or second order accurate schemes, while formulation for unstructured meshes makes complex computational domains amenable. The Spectral Difference Method offers a remarkably simple alternative to such high-order schemes for unstructured meshes as the Discontinuous Galerkin and Spectral Volume Methods. Significant contributions to the development of the Spectral Difference Method are presented, including stability analysis, viscous formulation, and h/p-multigrid convergence acceleration. Finally, the theory of Gibbs-complementary reconstruction is utilized in the context of high-order numerical methods for hyperbolic equations. Gibbs-complementary reconstruction makes it possible to extract pointwise high-order convergence in the spectral approximation of non

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

    NASA Astrophysics Data System (ADS)

    Sun, Wenjun; Jiang, Song; Xu, Kun; Li, Shu

    2015-12-01

    This paper presents an extension of previous work (Sun et al., 2015 [22]) of the unified gas kinetic scheme (UGKS) for the gray radiative transfer equations to the frequency-dependent (multi-group) radiative transfer system. Different from the gray radiative transfer equations, where the optical opacity is only a function of local material temperature, the simulation of frequency-dependent radiative transfer is associated with additional difficulties from the frequency-dependent opacity. For the multiple frequency radiation, the opacity depends on both the spatial location and the frequency. For example, the opacity is typically a decreasing function of frequency. At the same spatial region the transport physics can be optically thick for the low frequency photons, and optically thin for high frequency ones. Therefore, the optical thickness is not a simple function of space location. In this paper, the UGKS for frequency-dependent radiative system is developed. The UGKS is a finite volume method and the transport physics is modeled according to the ratio of the cell size to the photon's frequency-dependent mean free path. When the cell size is much larger than the photon's mean free path, a diffusion solution for such a frequency radiation will be obtained. On the other hand, when the cell size is much smaller than the photon's mean free path, a free transport mechanism will be recovered. In the regime between the above two limits, with the variation of the ratio between the local cell size and photon's mean free path, the UGKS provides a smooth transition in the physical and frequency space to capture the corresponding transport physics accurately. The seemingly straightforward extension of the UGKS from the gray to multiple frequency radiation system is due to its intrinsic consistent multiple scale transport modeling, but it still involves lots of work to properly discretize the multiple groups in order to design an asymptotic preserving (AP) scheme in all

  15. Gas-kinetic theory and Boltzmann equation of share price within an equilibrium market hypothesis and ad hoc strategy

    NASA Astrophysics Data System (ADS)

    Ausloos, M.

    2000-09-01

    Recent observations have indicated that the traditional equilibrium market hypothesis (EMH; also known as Efficient Market Hypothesis) is unrealistic. It is shown here that it is the analog of a Boltzmann equation in physics, thus having some bad properties of mean-field approximations like a Gaussian distribution of price fluctuations. A kinetic theory for prices can be simply derived, considering in a first approach that market actors have all identical relaxation times, and solved within a Chapman-Enskog like formalism. In closing the set of equations, (i) an equation of state with a pressure and (ii) the equilibrium (isothermal) equation for the price (taken as the order parameter) of a stock as a function of the volume of money available are obtained.

  16. A multilevel method for coupling the neutron kinetics and heat transfer equations

    SciTech Connect

    Tamang, A.; Anistratov, D. Y.

    2013-07-01

    We present a computational method for adequate and efficient coupling the neutron transport equation with the precursor and heat transfer equations. It is based on the multilevel nonlinear quasidiffusion (QD) method for solving the multigroup transport equation. The system of equations includes the time-dependent high-order transport equation and time-dependent multigroup and effective one-group low-order QD equations. We also formulate a method applying the {alpha}-approximation for the time-dependent high-order transport equation. This approach enables one to avoid storing the angular flux from the previous time step. Numerical results for a model transient problem are presented. (authors)

  17. Collision integral in the kinetic equation for a rarefied electron gas with allowance for its spin polarization

    SciTech Connect

    Sasorov, P. V.; Fomin, I. V.

    2015-06-15

    The collision integral in the kinetic equation for a rarefied spin-polarized gas of fermions (electrons) is derived. The collisions between these fermions and the collisions with much heavier particles (ions) forming a randomly located stationary background (gas) are taken into account. An important new circumstance is that the particle-particle scattering amplitude is not assumed to be small, which could be obtained, for example, in the first Born approximation. The derived collision integral can be used in the kinetic equation, including that for a relatively cold rarefied spin-polarized plasma with a characteristic electron energy below α{sup 2}m{sub e}c{sup 2}, where α is the fine-structure constant.

  18. Analysis of atmospheric flow over a surface protrusion using the turbulence kinetic energy equation with reference to aeronautical operating systems

    NASA Technical Reports Server (NTRS)

    Frost, W.; Harper, W. L.

    1975-01-01

    Flow over surface obstructions can produce significantly large wind shears such that adverse flying conditions can occur for aeronautical systems (helicopters, STOL vehicles, etc.). Atmospheric flow fields resulting from a semi-elliptical surface obstruction in an otherwise horizontally homogeneous statistically stationary flow are modelled with the boundary-layer/Boussinesq-approximation of the governing equation of fluid mechanics. The turbulence kinetic energy equation is used to determine the dissipative effects of turbulent shear on the mean flow. Iso-lines of turbulence kinetic energy and turbulence intensity are plotted in the plane of the flow and highlight regions of high turbulence intensity in the stagnation zone and sharp gradients in intensity along the transition from adverse to favourable pressure gradient. Discussion of the effects of the disturbed wind field in CTOL and STOL aircraft flight path and obstruction clearance standards is given. The results indicate that closer inspection of these presently recommended standards as influenced by wind over irregular terrains is required.

  19. Simulation of rarefied gas flows on the basis of the Boltzmann kinetic equation solved by applying a conservative projection method

    NASA Astrophysics Data System (ADS)

    Dodulad, O. I.; Kloss, Yu. Yu.; Potapov, A. P.; Tcheremissine, F. G.; Shuvalov, P. V.

    2016-06-01

    Flows of a simple rarefied gas and gas mixtures are computed on the basis of the Boltzmann kinetic equation, which is solved by applying various versions of the conservative projection method, namely, a two-point method for a simple gas and gas mixtures with a small difference between the molecular masses and a multipoint method in the case of a large mass difference. Examples of steady and unsteady flows are computed in a wide range of Mach and Knudsen numbers.

  20. Extracellular microbial polysaccharides. I. Substrate, biomass, and product kinetic equations for batch xanthan gum fermentation

    SciTech Connect

    Weiss, R.M.; Ollis, D.F.

    1980-04-01

    Various biomass (X), product (P), and substrate (S) rate equations are investigated in order to synthesize a general xanthan fermentation model from literature data. Analytical forms that provide reasonable descriptions for the X, P, and S behaviors reported by Moraine and Rogovin are shown to be the logistic equation, the Luedeking-Piret equation, and a modified Leudeking-Piret equation, respectively. The autonomous logistic equation allows the serial evaluation of parameters for all three equations, rather than a simultaneous determination required by nonautonomous models. 22 references.

  1. A unified gas-kinetic scheme for continuum and rarefied flows IV: Full Boltzmann and model equations

    NASA Astrophysics Data System (ADS)

    Liu, Chang; Xu, Kun; Sun, Quanhua; Cai, Qingdong

    2016-06-01

    Fluid dynamic equations are valid in their respective modeling scales, such as the particle mean free path scale of the Boltzmann equation and the hydrodynamic scale of the Navier-Stokes (NS) equations. With a variation of the modeling scales, theoretically there should have a continuous spectrum of fluid dynamic equations. Even though the Boltzmann equation is claimed to be valid in all scales, many Boltzmann solvers, including direct simulation Monte Carlo method, require the cell resolution to the order of particle mean free path scale. Therefore, they are still single scale methods. In order to study multiscale flow evolution efficiently, the dynamics in the computational fluid has to be changed with the scales. A direct modeling of flow physics with a changeable scale may become an appropriate approach. The unified gas-kinetic scheme (UGKS) is a direct modeling method in the mesh size scale, and its underlying flow physics depends on the resolution of the cell size relative to the particle mean free path. The cell size of UGKS is not limited by the particle mean free path. With the variation of the ratio between the numerical cell size and local particle mean free path, the UGKS recovers the flow dynamics from the particle transport and collision in the kinetic scale to the wave propagation in the hydrodynamic scale. The previous UGKS is mostly constructed from the evolution solution of kinetic model equations. Even though the UGKS is very accurate and effective in the low transition and continuum flow regimes with the time step being much larger than the particle mean free time, it still has space to develop more accurate flow solver in the region, where the time step is comparable with the local particle mean free time. In such a scale, there is dynamic difference from the full Boltzmann collision term and the model equations. This work is about the further development of the UGKS with the implementation of the full Boltzmann collision term in the region

  2. Utilization of Integrated Michaelis-Menten Equation to Determine Kinetic Constants

    ERIC Educational Resources Information Center

    Bezerra, Rui M. F.; Dias, Albino A.

    2007-01-01

    Students of biochemistry and related biosciences are urged to solve problems where kinetic parameters are calculated from initial rates obtained at different substrate concentrations. Troubles begin when they go to the laboratory to perform kinetic experiments and realize that usual laboratory instruments do not measure initial rates but only…

  3. Kinetic theory of transport processes in partially ionized reactive plasma, I: General transport equations

    NASA Astrophysics Data System (ADS)

    Zhdanov, V. M.; Stepanenko, A. A.

    2016-03-01

    In this paper we derive the set of general transport equations for multicomponent partially ionized reactive plasma in the presence of electric and magnetic fields taking into account the internal degrees of freedom and electronic excitation of plasma particles. Our starting point is a generalized Boltzmann equation with the collision integral in the Wang-Chang and Uhlenbeck form and a reactive collision integral. We obtain a set of conservation equations for such plasma and employ a linearized variant of Grad's moment method to derive the system of moment (or transport) equations for the plasma species nonequilibrium parameters. Full and reduced transport equations, resulting from the linearized system of moment equations, are presented, which can be used to obtain transport relations and expressions for transport coefficients of electrons and heavy plasma particles (molecules, atoms and ions) in partially ionized reactive plasma.

  4. Rarefied gas flow simulations using high-order gas-kinetic unified algorithms for Boltzmann model equations

    NASA Astrophysics Data System (ADS)

    Li, Zhi-Hui; Peng, Ao-Ping; Zhang, Han-Xin; Yang, Jaw-Yen

    2015-04-01

    This article reviews rarefied gas flow computations based on nonlinear model Boltzmann equations using deterministic high-order gas-kinetic unified algorithms (GKUA) in phase space. The nonlinear Boltzmann model equations considered include the BGK model, the Shakhov model, the Ellipsoidal Statistical model and the Morse model. Several high-order gas-kinetic unified algorithms, which combine the discrete velocity ordinate method in velocity space and the compact high-order finite-difference schemes in physical space, are developed. The parallel strategies implemented with the accompanying algorithms are of equal importance. Accurate computations of rarefied gas flow problems using various kinetic models over wide ranges of Mach numbers 1.2-20 and Knudsen numbers 0.0001-5 are reported. The effects of different high resolution schemes on the flow resolution under the same discrete velocity ordinate method are studied. A conservative discrete velocity ordinate method to ensure the kinetic compatibility condition is also implemented. The present algorithms are tested for the one-dimensional unsteady shock-tube problems with various Knudsen numbers, the steady normal shock wave structures for different Mach numbers, the two-dimensional flows past a circular cylinder and a NACA 0012 airfoil to verify the present methodology and to simulate gas transport phenomena covering various flow regimes. Illustrations of large scale parallel computations of three-dimensional hypersonic rarefied flows over the reusable sphere-cone satellite and the re-entry spacecraft using almost the largest computer systems available in China are also reported. The present computed results are compared with the theoretical prediction from gas dynamics, related DSMC results, slip N-S solutions and experimental data, and good agreement can be found. The numerical experience indicates that although the direct model Boltzmann equation solver in phase space can be computationally expensive

  5. Steady-state benchmarks of DK4D: A time-dependent, axisymmetric drift-kinetic equation solver

    SciTech Connect

    Lyons, B. C.; Jardin, S. C.; Ramos, J. J.

    2015-05-15

    The DK4D code has been written to solve a set of time-dependent, axisymmetric, finite-Larmor-radius drift-kinetic equations (DKEs) for the non-Maxwellian part of the electron and ion distribution functions using the full, linearized Fokker–Planck–Landau collision operator. The plasma is assumed to be in the low- to finite-collisionality regime, as is found in the cores of modern and future magnetic confinement fusion experiments. Each DKE is formulated such that the perturbed distribution function carries no net density, parallel momentum, or kinetic energy. Rather, these quantities are contained within the background Maxwellians and would be evolved by an appropriate set of extended magnetohydrodynamic (MHD) equations. This formulation allows for straight-forward coupling of DK4D to existing extended MHD time evolution codes. DK4D uses a mix of implicit and explicit temporal representations and finite element and spectral spatial representations. These, along with other computational methods used, are discussed extensively. Steady-state benchmarks are then presented comparing the results of DK4D to expected analytic results at low collisionality, qualitatively, and to the Sauter analytic fits for the neoclassical conductivity and bootstrap current, quantitatively. These benchmarks confirm that DK4D is capable of solving for the correct, gyroaveraged distribution function in stationary magnetic equilibria. Furthermore, the results presented demonstrate how the exact drift-kinetic solution varies with collisionality as a function of the magnetic moment and the poloidal angle.

  6. Physical and numerical sources of computational inefficiency in integration of chemical kinetic rate equations: Etiology, treatment and prognosis

    NASA Technical Reports Server (NTRS)

    Pratt, D. T.; Radhakrishnan, K.

    1986-01-01

    The design of a very fast, automatic black-box code for homogeneous, gas-phase chemical kinetics problems requires an understanding of the physical and numerical sources of computational inefficiency. Some major sources reviewed in this report are stiffness of the governing ordinary differential equations (ODE's) and its detection, choice of appropriate method (i.e., integration algorithm plus step-size control strategy), nonphysical initial conditions, and too frequent evaluation of thermochemical and kinetic properties. Specific techniques are recommended (and some advised against) for improving or overcoming the identified problem areas. It is argued that, because reactive species increase exponentially with time during induction, and all species exhibit asymptotic, exponential decay with time during equilibration, exponential-fitted integration algorithms are inherently more accurate for kinetics modeling than classical, polynomial-interpolant methods for the same computational work. But current codes using the exponential-fitted method lack the sophisticated stepsize-control logic of existing black-box ODE solver codes, such as EPISODE and LSODE. The ultimate chemical kinetics code does not exist yet, but the general characteristics of such a code are becoming apparent.

  7. AdS/CFT connection between Boltzmann and Einstein equations: Kinetic theory and pure gravity in AdS space

    SciTech Connect

    Iyer, Ramakrishnan; Mukhopadhyay, Ayan

    2010-04-15

    The AdS/CFT correspondence defines a sector with universal strongly coupled dynamics in the field theory as the dual of pure gravity in AdS described by Einstein's equation with a negative cosmological constant. We explain here, from the field-theoretic viewpoint how the dynamics in this sector gets determined by the expectation value of the energy-momentum tensor alone. We first show that the Boltzmann equation has very special solutions which could be functionally completely determined in terms of the energy-momentum tensor alone. We call these solutions conservative solutions. We indicate why conservative solutions should also exist when we refine this kinetic description to go closer to the exact microscopic theory or even move away from the regime of weak coupling so that no kinetic description could be employed. We argue that these conservative solutions form the universal sector dual to pure gravity at strong coupling and large N. Based on this observation, we propose a regularity condition on the energy-momentum tensor so that the dual solution in pure gravity has a smooth future horizon. We also study if irreversibility emerges only at long time scales of observation, unlike the case of the Boltzmann equation.

  8. Permeability and kinetic coefficients for mesoscale BCF surface step dynamics: Discrete two-dimensional deposition-diffusion equation analysis

    DOE PAGESBeta

    Zhao, Renjie; Evans, James W.; Oliveira, Tiago J.

    2016-04-08

    Here, a discrete version of deposition-diffusion equations appropriate for description of step flow on a vicinal surface is analyzed for a two-dimensional grid of adsorption sites representing the stepped surface and explicitly incorporating kinks along the step edges. Model energetics and kinetics appropriately account for binding of adatoms at steps and kinks, distinct terrace and edge diffusion rates, and possible additional barriers for attachment to steps. Analysis of adatom attachment fluxes as well as limiting values of adatom densities at step edges for nonuniform deposition scenarios allows determination of both permeability and kinetic coefficients. Behavior of these quantities is assessedmore » as a function of key system parameters including kink density, step attachment barriers, and the step edge diffusion rate.« less

  9. Using Lamm-Equation Modeling of Sedimentation Velocity Data to Determine the Kinetic and Thermodynamic Properties of Macromolecular Interactions

    PubMed Central

    Brautigam, Chad A.

    2011-01-01

    The interaction of macromolecules with themselves and with other macromolecules is fundamental to the functioning of living systems. Recent advances in the analysis of sedimentation velocity (SV) data obtained by analytical ultracentrifugation allow the experimenter to determine important features of such interactions, including the equilibrium association constant and information about the kinetic off-rate of the interaction. The determination of these parameters is made possible by the ability of modern software to fit numerical solutions of the Lamm Equation with kinetic considerations directly to SV data. Herein, the SV analytical advances implemented in the software package SEDPHAT are summarized. Detailed analyses of SV data using these strategies are presented. Finally, a few highlights of recent literature reports that feature this type of SV data analysis are surveyed. PMID:21187153

  10. Collision frequencies in density-matrix kinetic equations describing nonlinear effects in the wings of spectral lines

    SciTech Connect

    Parkhomenko, A I; Shalagin, Anatolii M

    2011-11-30

    Using the eikonal approximation, we have calculated effective collision frequencies in density-matrix kinetic equations describing nonlinear effects in the wings of spectral lines. We have established the relation between the probabilities of absorption and stimulated emission and the characteristics of the radiation and elementary scattering event. The example of the power interaction potential shows that quantum mechanical calculation of the collision frequencies in the eikonal approximation and previously known spectral line wing theory give similar results for the probability of radiation absorption.

  11. Analytical and semianalytical solutions to the kinetic equation with Coulomb collision term and a monoenergetic source function

    SciTech Connect

    Goncharov, P. R.; Kuteev, B. V.; Ozaki, T.; Sudo, S.

    2010-11-15

    Analytical and semianalytical solutions have been obtained using a practical dimensionless form of Boltzmann kinetic equation assuming spatial homogeneity, azimuthal symmetry, and Maxwellian distributions of target plasma species. In contrast with formerly considered simplified equations with truncated collision terms, the exact Landau-Boltzmann collision operator is used, which conserves the number of particles, nullifies the collision term at statistical equilibrium, and describes the Maxwellization process naturally observed in correct solutions. A comparison with previous stationary and time-dependent analytical solutions is given. The new semianalytical results can be used in numerical modeling, for verification of solutions in more complex models, and in experimental data analysis, especially concerning nuclear processes and advanced localized, angle-resolved suprathermal particle diagnostics.

  12. A Gas-Kinetic Method for Hyperbolic-Elliptic Equations and Its Application in Two-Phase Fluid Flow

    NASA Technical Reports Server (NTRS)

    Xu, Kun

    1999-01-01

    A gas-kinetic method for the hyperbolic-elliptic equations is presented in this paper. In the mixed type system, the co-existence and the phase transition between liquid and gas are described by the van der Waals-type equation of state (EOS). Due to the unstable mechanism for a fluid in the elliptic region, interface between the liquid and gas can be kept sharp through the condensation and evaporation process to remove the "averaged" numerical fluid away from the elliptic region, and the interface thickness depends on the numerical diffusion and stiffness of the phase change. A few examples are presented in this paper for both phase transition and multifluid interface problems.

  13. A covariant Fokker-Planck equation for a simple gas from relativistic kinetic theory

    SciTech Connect

    Chacon-Acosta, Guillermo; Dagdug, Leonardo; Morales-Tecotl, Hugo A.

    2010-12-14

    A manifestly covariant Fokker-Planck differential equation is derived for the case of a relativistic simple gas by taking a small momentum transfer approximation within the collision integral of the relativistic Boltzmann equation. We follow closely previous work, with the main difference that we keep manifest covariance at every stage of the analysis. In addition, we use the covariant Juettner distribution function to find a relativistic generalization of the Einstein's fluctuation-dissipation relation.

  14. A Reynolds-averaged turbulence modeling approach using three transport equations for the turbulent viscosity, kinetic energy, and dissipation rate

    NASA Astrophysics Data System (ADS)

    Yoshizawa, Akira; Abe, Hiroyuki; Matsuo, Yuichi; Fujiwara, Hitoshi; Mizobuchi, Yasuhiro

    2012-07-01

    A Reynolds-averaged approach to turbulent shear flows is sought with resort to a three-equation method. Its novelty is the introduction of a turbulent-viscosity transport equation through the transport equation for the Reynolds stress in addition to those for the turbulent kinetic energy and the dissipation rate. The latter two equations are used for evaluating the dimensional coefficients in the former. The aim of this model is to enhance the capability to cope with nonstationary and advection effects in various turbulent flows. The adaptability to them is confirmed through the application to homogeneous-shear and supersonic free-shear flows. In particular, the reasonable prediction is obtained in the latter where the growth rate of the shear layer is suppressed with the increase in the convective Mach number. The present model is also applied to a three-dimensional flow past a wing tip as an instance of complex aeronautical flows, and the excessive diffusion of the trailing vortices is shown to be suppressed. The turbulent-viscosity representation for the Reynolds stress is systematically supplemented with nonlinear effects of mean-velocity gradient tensors, and its adequacy is verified in a channel flow.

  15. Approaches to the solution of coupled multiexponential transient-state rate kinetic equations: A critical review.

    PubMed

    Fisher, Harvey F

    2016-08-01

    The transient-state kinetic approach has failed to reach its full potential despite its advantage over the steady-state approach in its ability to observe mechanistic events directly and in real time. This failure has been due in part to the lack of any rigorously derived and readily applicable body of theory corresponding to that which currently characterizes the steady-state approach. In order to clarify the causes of this discrepancy and to suggest a route to its solution we examine the capabilities and limitations of the various forms of transient-state kinetic approaches to the mathematical resolution of enzymatic reaction mechanisms currently available. We document a lack of validity inherent in their basic assumptions and suggest the need for a potentially more rigorous analytic approach. PMID:27173736

  16. Reduced fluid-kinetic equations for low-frequency dynamics, magnetic reconnection, and electron heating in low-beta plasmas

    SciTech Connect

    Zocco, Alessandro; Schekochihin, Alexander A.

    2011-10-15

    A minimal model for magnetic reconnection and, generally, low-frequency dynamics in low-beta plasmas is proposed. The model combines analytical and computational simplicity with physical realizability: it is a rigorous limit of gyrokinetics for plasma beta of order the electron-ion mass ratio. The model contains collisions and can be used both in the collisional and collisionless reconnection regimes. It includes gyrokinetic ions (not assumed cold) and allows for the topological rearrangement of the magnetic field lines by either resistivity or electron inertia, whichever predominates. The two-fluid dynamics are coupled to electron kinetics--electrons are not assumed isothermal and are described by a reduced drift-kinetic equation. The model, therefore allows for irreversibility and conversion of magnetic energy into electron heat via parallel phase mixing in velocity space. An analysis of the exchanges between various forms of free energy and its conversion into electron heat is provided. It is shown how all relevant linear waves and regimes of the tearing instability (collisionless, semicollisional, and fully resistive) are recovered in various limits of our model. An efficient way to simulate our equations numerically is proposed, via the Hermite representation of the velocity space. It is shown that small scales in velocity space will form, giving rise to a shallow Hermite-space spectrum, whence it is inferred that, for steady-state or sufficiently slow dynamics, the electron heating rate will remain finite in the limit of vanishing collisionality.

  17. Exact solution of the one- and three-dimensional quantum kinetic equations with velocity-dependent collision rates: Comparative analysis

    NASA Astrophysics Data System (ADS)

    Privalov, T.; Shalagin, A.

    1999-06-01

    The interaction of a plane monochromatic traveling wave with two-level particles suffering collisions with buffer-gas particles is considered. Collision rates are assumed to be velocity dependent. The collision integral is obtained on the basis of the strong-collision model, generalized to the case of velocity-dependent collision rates (the so-called ``kangaroo'' model). We obtained the exact analytical solution of the problem for arbitrary intensity of radiation, arbitrary ratio of homogeneous and Doppler widths of the absorption line, and arbitrary mass ratio between absorbing- and buffer-gas particles. The obtained analytical solutions of the quantum kinetic equations allowed us to analyze the spectral shape of the strong-field absorption line as well as the probe-field absorption line (the nonlinear part of the work done by the probe field) and the frequency dependence of the light-induced drift (LID) velocity. A comprehensive comparative analysis for the three- and one-dimensional versions of the model is given. On the basis of this analysis, we reach the conclusion that the one-dimensional quantum kinetic equation has quite a wide range of application. We also reveal the conditions for the strongest manifestation of the velocity dependence of the collision rates, which affects most strongly the anomalous LID.

  18. Evaluation of the recrystallization kinetics of hot-melt extruded polymeric solid dispersions using an improved Avrami equation.

    PubMed

    Feng, Xin; Ye, Xingyou; Park, Jun-Bom; Lu, Wenli; Morott, Joe; Beissner, Brad; Lian, Zhuoyang John; Pinto, Elanor; Bi, Vivian; Porter, Stu; Durig, Tom; Majumdar, Soumyajit; Repka, Michael A

    2015-01-01

    The recrystallization of an amorphous drug in a solid dispersion system could lead to a loss in the drug solubility and bioavailability. The primary objective of the current research was to use an improved kinetic model to evaluate the recrystallization kinetics of amorphous structures and to further understand the factors influencing the physical stability of amorphous solid dispersions. Amorphous solid dispersions of fenofibrate with different molecular weights of hydroxypropylcellulose, HPC (Klucel™ LF, EF, ELF) were prepared utilizing hot-melt extrusion technology. Differential scanning calorimetry was utilized to quantitatively analyze the extent of recrystallization in the samples stored at different temperatures and relative humidity (RH) conditions. The experimental data were fitted into the improved kinetics model of a modified Avrami equation to calculate the recrystallization rate constants. Klucel LF, the largest molecular weight among the HPCs used, demonstrated the greatest inhibition of fenofibrate recrystallization. Additionally, the recrystallization rate (k) decreased with increasing polymer content, however exponentially increased with higher temperature. Also k increased linearly rather than exponentially over the range of RH studied. PMID:25224341

  19. A numerical method for integrating the kinetic equations of droplet spectra evolution by condensation/evaporation and by coalescence/breakup processes

    NASA Technical Reports Server (NTRS)

    Emukashvily, I. M.

    1982-01-01

    An extension of the method of moments is developed for the numerical integration of the kinetic equations of droplet spectra evolution by condensation/evaporation and by coalescence/breakup processes. The number density function n sub k (x,t) in each separate droplet packet between droplet mass grid points (x sub k, x sub k+1) is represented by an expansion in orthogonal polynomials with a given weighting function. In this way droplet number concentrations, liquid water contents and other moments in each droplet packet are conserved and the problem of solving the kinetic equations is replaced by one of solving a set of coupled differential equations for the number density function moments. The method is tested against analytic solutions of the corresponding kinetic equations. Numerical results are obtained for different coalescence/breakup and condensation/evaporation kernels and for different initial droplet spectra. Also droplet mass grid intervals, weighting functions, and time steps are varied.

  20. Thermostatted kinetic equations as models for complex systems in physics and life sciences

    NASA Astrophysics Data System (ADS)

    Bianca, Carlo

    2012-12-01

    Statistical mechanics is a powerful method for understanding equilibrium thermodynamics. An equivalent theoretical framework for nonequilibrium systems has remained elusive. The thermodynamic forces driving the system away from equilibrium introduce energy that must be dissipated if nonequilibrium steady states are to be obtained. Historically, further terms were introduced, collectively called a thermostat, whose original application was to generate constant-temperature equilibrium ensembles. This review surveys kinetic models coupled with time-reversible deterministic thermostats for the modeling of large systems composed both by inert matter particles and living entities. The introduction of deterministic thermostats allows to model the onset of nonequilibrium stationary states that are typical of most real-world complex systems. The first part of the paper is focused on a general presentation of the main physical and mathematical definitions and tools: nonequilibrium phenomena, Gauss least constraint principle and Gaussian thermostats. The second part provides a review of a variety of thermostatted mathematical models in physics and life sciences, including Kac, Boltzmann, Jager-Segel and the thermostatted (continuous and discrete) kinetic for active particles models. Applications refer to semiconductor devices, nanosciences, biological phenomena, vehicular traffic, social and economics systems, crowds and swarms dynamics.

  1. Treatment of Severe Metabolic Alkalosis with Continuous Renal Replacement Therapy: Bicarbonate Kinetic Equations of Clinical Value.

    PubMed

    Yessayan, Lenar; Yee, Jerry; Frinak, Stan; Kwon, David; Szamosfalvi, Balazs

    2015-01-01

    Concomitant severe metabolic alkalosis, hypernatremia, and kidney failure pose a therapeutic challenge. Hemodialysis to correct azotemia and abnormal electrolytes results in rapid correction of serum sodium, bicarbonate, and urea but presents a risk for dialysis disequilibrium and brain edema. We describe a patient with Zollinger-Ellison syndrome with persistent encephalopathy, severe metabolic alkalosis (highest bicarbonate 81 mEq/L), hypernatremia (sodium 157 mEq/L), and kidney failure despite 30 hours of intravenous crystalloids and proton pump inhibitor. We used continuous renal replacement therapy (RRT) with delivered hourly urea clearance of ~3 L/hour (24 hour sustained low efficiency dialysis with regional citrate anticoagulation protocol at blood flow rate 60 ml/min and dialysate flow rate 400 ml/min). To mitigate a pronounced decrease in plasma osmolality while removing urea from this hypernatremic patient, dialysate sodium was set to start at 155 mEq/L then at 150 mEq/L after 6 hours. Serum bicarbonate, urea, and sodium were slowly corrected over 26 hours. This case demonstrates how to regulate and predict the systemic bicarbonate level using single pool kinetic modeling during convective or diffusive RRT. Kinetic modeling provides a valuable tool for systemic blood pH control in future combined use of extracorporeal CO2 removal and continuous RRT systems. PMID:25794247

  2. Master equation for a kinetic model of a trading market and its analytic solution

    NASA Astrophysics Data System (ADS)

    Chatterjee, Arnab; Chakrabarti, Bikas K.; Stinchcombe, Robin B.

    2005-08-01

    We analyze an ideal-gas-like model of a trading market with quenched random saving factors for its agents and show that the steady state income (m) distribution P(m) in the model has a power law tail with Pareto index ν exactly equal to unity, confirming the earlier numerical studies on this model. The analysis starts with the development of a master equation for the time development of P(m) . Precise solutions are then obtained in some special cases.

  3. Geochemical kinetics via the Swift-Connick equations and solution NMR

    NASA Astrophysics Data System (ADS)

    Harley, Steven J.; Ohlin, C. André; Casey, William H.

    2011-07-01

    Signal analysis in Nuclear Magnetic Resonance spectroscopy is among the most powerful methods to quantify reaction rates in aqueous solutions. To this end, the Swift-Connick approximations to the Bloch-McConnell equations have been used extensively to estimate rate parameters for elementary reactions. The method is primarily used for 17O NMR in aqueous solutions, but the list of geochemically relevant nuclei that can be used is long, and includes 29Si, 27Al, 19F, 13C and many others of particular interest to geochemists. Here we review the derivation of both the Swift-Connick and Bloch-McConnell equations and emphasize assumptions and quirks. For example, the equations were derived for CW-NMR, but are used with modern pulse FT-NMR and can be applied to systems that have exchange rates that are shorter than the lifetime of a typical pulse. The method requires a dilute solution where the minor reacting species contributes a negligible amount of total magnetization. We evaluate the sensitivity of results to this dilute-solution requirement and also highlight the need for chemically well-defined systems if reliable data are to be obtained. The limitations in using longitudinal relaxation to estimate reaction rate parameters are discussed. Finally, we provide examples of the application of the method, including ligand exchanges from aqua ions and hydrolysis complexes, that emphasize its flexibility. Once the basic requirements of the Swift-Connick method are met, it allows geochemists to establish rates of elementary reactions. Reactions at this scale lend themselves well to methods of computational simulation and could provide key tests of accuracy.

  4. Characterization of Weak Convergence of Probability-Valued Solutions of General One-Dimensional Kinetic Equations

    NASA Astrophysics Data System (ADS)

    Perversi, Eleonora; Regazzini, Eugenio

    2015-05-01

    For a general inelastic Kac-like equation recently proposed, this paper studies the long-time behaviour of its probability-valued solution. In particular, the paper provides necessary and sufficient conditions for the initial datum in order that the corresponding solution converges to equilibrium. The proofs rest on the general CLT for independent summands applied to a suitable Skorokhod representation of the original solution evaluated at an increasing and divergent sequence of times. It turns out that, roughly speaking, the initial datum must belong to the standard domain of attraction of a stable law, while the equilibrium is presentable as a mixture of stable laws.

  5. Analysis of Electromagnetic Wave Propagation in a Magnetized Re-Entry Plasma Sheath Via the Kinetic Equation

    NASA Technical Reports Server (NTRS)

    Manning, Robert M.

    2009-01-01

    Based on a theoretical model of the propagation of electromagnetic waves through a hypersonically induced plasma, it has been demonstrated that the classical radiofrequency communications blackout that is experienced during atmospheric reentry can be mitigated through the appropriate control of an external magnetic field of nominal magnitude. The model is based on the kinetic equation treatment of Vlasov and involves an analytical solution for the electric and magnetic fields within the plasma allowing for a description of the attendant transmission, reflection and absorption coefficients. The ability to transmit through the magnetized plasma is due to the magnetic windows that are created within the plasma via the well-known whistler modes of propagation. The case of 2 GHz transmission through a re-entry plasma is considered. The coefficients are found to be highly sensitive to the prevailing electron density and will thus require a dynamic control mechanism to vary the magnetic field as the plasma evolves through the re-entry phase.

  6. A Bhatnagar-Gross-Krook-like Model Kinetic Equation for a Granular Gas of Inelastic Rough Hard Spheres

    SciTech Connect

    Santos, Andres

    2011-05-20

    The Boltzmann collision operator for a dilute granular gas of inelastic rough hard spheres is much more intricate than its counterpart for inelastic smooth spheres. Now the one-body distribution function depends not only on the translational velocity v of the center of mass but also on the angular velocity {omega} of the particle. Moreover, the collision rules couple v and {omega}, involving not only the coefficient of normal restitution {alpha} but also the coefficient of tangential restitution {beta}. The aim of this paper is to propose an extension to inelastic rough particles of a Bhatnagar-Gross-Krook-like kinetic model previously proposed for inelastic smooth particles. The Boltzmann collision operator is replaced by the sum of three terms representing: (i) the relaxation to a two-temperature local equilibrium distribution, (ii) the action of a nonconservative drag force F proportional to v-u(u being the flow velocity), and (iii) the action of a nonconservative torque M equal to a linear combination of {omega} and {Omega}({Omega} being the mean angular velocity). The three coefficients in F and M are fixed to reproduce the Boltzmann collisional rates of change of {Omega} and of the two granular temperatures (translational and rotational). A simpler version of the model is also constructed in the form of two coupled kinetic equations for the translational and rotational velocity distributions. The kinetic model is applied to the simple shear flow steady state and the combined influence of {alpha} and {beta} on the shear and normal stresses and on the translational velocity distribution function is analyzed.

  7. Pharmaceutical solid-state kinetic stability investigation by using moisture-modified Arrhenius equation and JMP statistical software.

    PubMed

    Fu, Mingkun; Perlman, Michael; Lu, Qing; Varga, Csanad

    2015-03-25

    An accelerated stress approach utilizing the moisture-modified Arrhenius equation and JMP statistical software was utilized to quantitatively assess the solid state stability of an investigational oncology drug MLNA under the influence of temperature (1/T) and humidity (%RH). Physical stability of MLNA under stress conditions was evaluated by using XRPD, DSC, TGA, and DVS, while chemical stability was evaluated by using HPLC. The major chemical degradation product was identified as a hydrolysis product of MLNA drug substance, and was subsequently subjected to an investigation of kinetics based on the isoconversion concept. A mathematical model (ln k=-11,991×(1/T)+0.0298×(%RH)+29.8823) based on the initial linear kinetics observed for the formation of this degradant at all seven stress conditions was built by using the moisture-modified Arrhenius equation and JMP statistical software. Comparison of the predicted versus experimental lnk values gave a mean deviation value of 5.8%, an R(2) value of 0.94, a p-value of 0.0038, and a coefficient of variation of the root mean square error CV(RMSE) of 7.9%. These statistics all indicated a good fit to the model for the stress data of MLNA. Both temperature and humidity were shown to have a statistically significant impact on stability by using effect leverage plots (p-value<0.05 for both 1/T and %RH). Inclusion of a term representing the interaction of relative humidity and temperature (%RH×1/T) was shown not to be justified by using Analysis of Covariance (ANCOVA), which supported the use of the moisture-corrected Arrhenius equation modeling theory. The model was found to be of value to aid setting of specifications and retest period, and storage condition selection. A model was also generated using only four conditions, as an example from a resource saving perspective, which was found to provide a good fit to the entire set of data. PMID:25656488

  8. Inference for Stochastic Chemical Kinetics Using Moment Equations and System Size Expansion

    PubMed Central

    Thomas, Philipp; Kazeroonian, Atefeh; Theis, Fabian J.; Grima, Ramon; Hasenauer, Jan

    2016-01-01

    Quantitative mechanistic models are valuable tools for disentangling biochemical pathways and for achieving a comprehensive understanding of biological systems. However, to be quantitative the parameters of these models have to be estimated from experimental data. In the presence of significant stochastic fluctuations this is a challenging task as stochastic simulations are usually too time-consuming and a macroscopic description using reaction rate equations (RREs) is no longer accurate. In this manuscript, we therefore consider moment-closure approximation (MA) and the system size expansion (SSE), which approximate the statistical moments of stochastic processes and tend to be more precise than macroscopic descriptions. We introduce gradient-based parameter optimization methods and uncertainty analysis methods for MA and SSE. Efficiency and reliability of the methods are assessed using simulation examples as well as by an application to data for Epo-induced JAK/STAT signaling. The application revealed that even if merely population-average data are available, MA and SSE improve parameter identifiability in comparison to RRE. Furthermore, the simulation examples revealed that the resulting estimates are more reliable for an intermediate volume regime. In this regime the estimation error is reduced and we propose methods to determine the regime boundaries. These results illustrate that inference using MA and SSE is feasible and possesses a high sensitivity. PMID:27447730

  9. Inference for Stochastic Chemical Kinetics Using Moment Equations and System Size Expansion.

    PubMed

    Fröhlich, Fabian; Thomas, Philipp; Kazeroonian, Atefeh; Theis, Fabian J; Grima, Ramon; Hasenauer, Jan

    2016-07-01

    Quantitative mechanistic models are valuable tools for disentangling biochemical pathways and for achieving a comprehensive understanding of biological systems. However, to be quantitative the parameters of these models have to be estimated from experimental data. In the presence of significant stochastic fluctuations this is a challenging task as stochastic simulations are usually too time-consuming and a macroscopic description using reaction rate equations (RREs) is no longer accurate. In this manuscript, we therefore consider moment-closure approximation (MA) and the system size expansion (SSE), which approximate the statistical moments of stochastic processes and tend to be more precise than macroscopic descriptions. We introduce gradient-based parameter optimization methods and uncertainty analysis methods for MA and SSE. Efficiency and reliability of the methods are assessed using simulation examples as well as by an application to data for Epo-induced JAK/STAT signaling. The application revealed that even if merely population-average data are available, MA and SSE improve parameter identifiability in comparison to RRE. Furthermore, the simulation examples revealed that the resulting estimates are more reliable for an intermediate volume regime. In this regime the estimation error is reduced and we propose methods to determine the regime boundaries. These results illustrate that inference using MA and SSE is feasible and possesses a high sensitivity. PMID:27447730

  10. Kinetic theory of spin-polarized systems in electric and magnetic fields with spin-orbit coupling. I. Kinetic equation and anomalous Hall and spin-Hall effects

    NASA Astrophysics Data System (ADS)

    Morawetz, K.

    2015-12-01

    The coupled kinetic equation for density and spin Wigner functions is derived including spin-orbit coupling, electric and magnetic fields, and self-consistent Hartree mean fields suited for SU(2) transport. The interactions are assumed to be with scalar and magnetic impurities as well as scalar and spin-flip potentials among the particles. The spin-orbit interaction is used in a form suitable for solid state physics with Rashba or Dresselhaus coupling, graphene, extrinsic spin-orbit coupling, and effective nuclear matter coupling. The deficiencies of the two-fluid model are worked out consisting of the appearance of an effective in-medium spin precession. The stationary solution of all these systems shows a band splitting controlled by an effective medium-dependent Zeeman field. The self-consistent precession direction is discussed and a cancellation of linear spin-orbit coupling at zero temperature is reported. The precession of spin around this effective direction caused by spin-orbit coupling leads to anomalous charge and spin currents in an electric field. Anomalous Hall conductivity is shown to consist of the known results obtained from the Kubo formula or Berry phases and a symmetric part interpreted as an inverse Hall effect. Analogously the spin-Hall and inverse spin-Hall effects of spin currents are discussed which are present even without magnetic fields showing a spin accumulation triggered by currents. The analytical dynamical expressions for zero temperature are derived and discussed in dependence on the magnetic field and effective magnetizations. The anomalous Hall and spin-Hall effect changes sign at higher than a critical frequency dependent on the relaxation time.

  11. Self-Consistent System of Equations for a Kinetic Description of the Low-Pressure Discharges Accounting for the Nonlocal and Collisionless Electron Dynamics

    SciTech Connect

    Igor D. Kaganovich; Oleg Polomarov

    2003-05-19

    In low-pressure discharges, when the electron mean free path is larger or comparable with the discharge length, the electron dynamics is essentially non-local. Moreover, the electron energy distribution function (EEDF) deviates considerably from a Maxwellian. Therefore, an accurate kinetic description of the low-pressure discharges requires knowledge of the non-local conductivity operator and calculation of the non-Maxwellian EEDF. The previous treatments made use of simplifying assumptions: a uniform density profile and a Maxwellian EEDF. In the present study a self-consistent system of equations for the kinetic description of nonlocal, non-uniform, nearly collisionless plasmas of low-pressure discharges is derived. It consists of the nonlocal conductivity operator and the averaged kinetic equation for calculation of the non-Maxwellian EEDF. The importance of accounting for the non-uniform plasma density profile on both the current density profile and the EEDF is demonstrated.

  12. Exponential-fitted methods for integrating stiff systems of ordinary differential equations: Applications to homogeneous gas-phase chemical kinetics

    NASA Technical Reports Server (NTRS)

    Pratt, D. T.

    1984-01-01

    Conventional algorithms for the numerical integration of ordinary differential equations (ODEs) are based on the use of polynomial functions as interpolants. However, the exact solutions of stiff ODEs behave like decaying exponential functions, which are poorly approximated by polynomials. An obvious choice of interpolant are the exponential functions themselves, or their low-order diagonal Pade (rational function) approximants. A number of explicit, A-stable, integration algorithms were derived from the use of a three-parameter exponential function as interpolant, and their relationship to low-order, polynomial-based and rational-function-based implicit and explicit methods were shown by examining their low-order diagonal Pade approximants. A robust implicit formula was derived by exponential fitting the trapezoidal rule. Application of these algorithms to integration of the ODEs governing homogenous, gas-phase chemical kinetics was demonstrated in a developmental code CREK1D, which compares favorably with the Gear-Hindmarsh code LSODE in spite of the use of a primitive stepsize control strategy.

  13. The solution to the nuclear reactor kinetic equations with a continuous slowing down model using numerical Laplace transform inversion

    SciTech Connect

    Ravetto, P.; Sumini, M.; Ganapol, B.D.

    1988-01-01

    In an attempt to better understand the influence of prompt and delayed neutrons on nuclear reactor dynamics, a continuous slowing down model based on Fermi age theory was developed several years ago. This model was easily incorporated into the one-group diffusion equation and provided a realistic physical picture of how delayed and prompt neutrons slow down and simultaneously diffuse throughout a medium. The model allows for different slowing down times for each delayed neutron group as well as for prompt neutrons and for spectral differences between the two typed of neutrons. Because of its generality, this model serves not only a a useful predictive tool to anticipate reactor transients, but also as an excellent educational tool to demonstrate the effect of delayed neutrons in reactor kinetics. However, because of numerical complications, the slowing down model could not be developed to its full potential. In particular, the major limitation was the inversion of the Laplace transform, which relied on a knowledge of the poles associated with the resulting transformed flux. For this reason, only one group of delayed neutrons and times longer than the slowing down times could be considered. As is shown, the new inversion procedure removes the short time limitation as well as allows for any number of delayed neutron groups. The inversion technique is versatile and is useful in teaching numerical methods in nuclear science.

  14. Numerical simulation of flood inundation using a well-balanced kinetic scheme for the shallow water equations with bulk recharge and discharge

    NASA Astrophysics Data System (ADS)

    Ersoy, Mehmet; Lakkis, Omar; Townsend, Philip

    2016-04-01

    The flow of water in rivers and oceans can, under general assumptions, be efficiently modelled using Saint-Venant's shallow water system of equations (SWE). SWE is a hyperbolic system of conservation laws (HSCL) which can be derived from a starting point of incompressible Navier-Stokes. A common difficulty in the numerical simulation of HSCLs is the conservation of physical entropy. Work by Audusse, Bristeau, Perthame (2000) and Perthame, Simeoni (2001), proposed numerical SWE solvers known as kinetic schemes (KSs), which can be shown to have desirable entropy-consistent properties, and are thus called well-balanced schemes. A KS is derived from kinetic equations that can be integrated into the SWE. In flood risk assessment models the SWE must be coupled with other equations describing interacting meteorological and hydrogeological phenomena such as rain and groundwater flows. The SWE must therefore be appropriately modified to accommodate source and sink terms, so kinetic schemes are no longer valid. While modifications of SWE in this direction have been recently proposed, e.g., Delestre (2010), we depart from the extant literature by proposing a novel model that is "entropy-consistent" and naturally extends the SWE by respecting its kinetic formulation connections. This allows us to derive a system of partial differential equations modelling flow of a one-dimensional river with both a precipitation term and a groundwater flow model to account for potential infiltration and recharge. We exhibit numerical simulations of the corresponding kinetic schemes. These simulations can be applied to both real world flood prediction and the tackling of wider issues on how climate and societal change are affecting flood risk.

  15. Kinlsq: a program for fitting kinetics data with numerically integrated rate equations and its application to the analysis of slow, tight-binding inhibition data.

    PubMed

    Gutheil, W G; Kettner, C A; Bachovchin, W W

    1994-11-15

    Kinlsq, a Matlab-based computer program for the least-squares fitting of parameters to kinetics data described by numerically integrated rate equations, is described, and three applications to the analysis of enzyme kinetics data are given. The first application was to the analysis of a simple bimolecular enzyme plus inhibitor binding curve. The kinlsq fit to these data was essentially identical to that obtained with the corresponding analytically integrated rate equation, validating kinlsq. The second application was to the fit of a numerically integrated Michaelis-Menten model to the progress curve for dipeptidyl peptidase IV-catalyzed hydrolysis of Ala-Pro-p-nitroanilide as a demonstration of the analysis of steady-state enzyme kinetics data. The results obtained with kinlsq were compared with the results obtained by fitting this time course with the integrated Michaelis-Menten equation, and with the results obtained by fitting the (S,dP/dt) transform of the data with the Michaelis-Menten equation. The third application was to the analysis of the inhibition of chymotrypsin by the slow, tight-binding inhibitor MeOSuc-Ala-Ala-Pro-boroPhe, data not readily amenable to other methods of analysis. These applications demonstrate how kinlsq can be used to fit rate constants, equilibrium constants, steady-state constants, and the stoichiometric relationships between components. PMID:7695087

  16. Numerical modeling of the sound propagation through a rarefied gas in a semi-infinite space on the basis of linearized kinetic equation.

    PubMed

    Sharipov, Felix; Kalempa, Denize

    2008-10-01

    The sound propagation through a rarefied gas is investigated on the basis of the linearized kinetic equation. A plate oscillating in the direction normal to its own plane is considered as a sound wave source. It is assumed a fully established oscillation so that the solution of the kinetic equation depends on time harmonically, while its dependence on the spatial coordinates is obtained numerically. The problem is solved over a wide range of the oscillation speed parameter defined as a ratio of the intermolecular collision frequency to the sound frequency. In order to evaluate the influence of the momentum and energy accommodation coefficients on the solution of the problem, the Cercignani-Lampis scattering kernel is applied as the boundary condition. An analysis of wave characteristics near the source surface shows that they are significantly different from those far from the surface even if the oscillation is slow, i.e., the solution is not harmonic in the space. PMID:19062839

  17. Equations describing the effects of the transient stability of subpopulations of nascent hnRNAs on the kinetics of turnover of hnRNA.

    PubMed Central

    Kleene, K C

    1986-01-01

    The equations that have been used previously to analyse the rate of decay of hnRNA implicitly assume that nascent hnRNAs are degraded stochastically. This assumption is inconsistent with electron-microscopic studies of transcription cited here, which show that nascent hnRNAs are not degraded during transcription, implying that hnRNA degradation occurs only after termination of transcription and release of the hnRNA from chromatin. Equations are derived describing the accumulation of radioactivity hnRNA during continuous labelling assuming that nascent hnRNAs are stable and that hnRNAs decay with first-order kinetics only after completion of transcription. The effects of the transient stability of nascent hnRNAs on the kinetics of hnRNA turnover can become important when the half-life of the hnRNA is shorter than the time to transcribe an hnRNA from the point of initiation to the point of termination. These equations should prove useful in studies of hnRNA turnover that require a precise description of the labelling kinetics of nascent and completed subpopulations of hnRNA. PMID:3707533

  18. Equations describing the effects of the transient stability of subpopulations of nascent hnRNAs on the kinetics of turnover of hnRNA.

    PubMed

    Kleene, K C

    1986-02-01

    The equations that have been used previously to analyse the rate of decay of hnRNA implicitly assume that nascent hnRNAs are degraded stochastically. This assumption is inconsistent with electron-microscopic studies of transcription cited here, which show that nascent hnRNAs are not degraded during transcription, implying that hnRNA degradation occurs only after termination of transcription and release of the hnRNA from chromatin. Equations are derived describing the accumulation of radioactivity hnRNA during continuous labelling assuming that nascent hnRNAs are stable and that hnRNAs decay with first-order kinetics only after completion of transcription. The effects of the transient stability of nascent hnRNAs on the kinetics of hnRNA turnover can become important when the half-life of the hnRNA is shorter than the time to transcribe an hnRNA from the point of initiation to the point of termination. These equations should prove useful in studies of hnRNA turnover that require a precise description of the labelling kinetics of nascent and completed subpopulations of hnRNA. PMID:3707533

  19. Boltzmann equation modelling of Learning Dynamics. Comment on "Collective learning modeling based on the kinetic theory of active particles" by D. Burini et al.

    NASA Astrophysics Data System (ADS)

    Shizgal, Bernie

    2016-03-01

    The paper by Burini et al. [7] presents an interesting use of the Boltzmann equation of kinetic theory to model real learning processes. The authors provide a comprehensive discussion of the basic concepts involved in their modelling work. The Boltzmann equation as used by physicists and chemists to model a variety of transport processes in many diverse fields is based on the notion of the binary collisions between identifiable particles in the defined system [9]. The particles exchange energy on collision and the distribution function, which depends on the three velocity components and the three spatial coordinates, varies with time. The classical or quantum collision dynamics between particles play a central role in the definition of the kernels in the integral operators that define the Boltzmann equation [8].

  20. On the accuracy of a one-dimensional approach to the solution of kinetic equations with velocity-dependent collision frequencies

    SciTech Connect

    Parkhomenko, A. I. Shalagin, A. M.

    2014-11-15

    The solution of many problems in light-induced gas kinetics can be simplified significantly using quantum kinetic equations in the context of the so-called one-dimensional approximation, in which the initial equations are averaged over transverse (relative to the direction of radiation) velocities. The errors introduced in such an approach are usually assumed to be small; however, this has been confirmed quantitatively only on the basis of the simplest (two- and three-level) particle models. We analyze the accuracy of the one-dimensional approximation for multilevel particles quantitatively for the light-induced drift (LID) effect in cesium atoms in the atmosphere of inert buffer gases. It is shown that in the case of the so-called “normal” LID, one-dimensional kinetic equations can always be used instead of three-dimensional equations without a risk of losing some important fine details in the dependence of the drift velocity on the radiation frequency. In the case of anomalous LID, the error of the one-dimensional approximation is also insignificant, but it can be disregarded only in the case of light buffer particles. For comparable masses of resonant and buffer particles, the one-dimensional approximation may give a noticeable error in determination of drift velocity amplitudes; however, the positions of drift velocity zeros and extrema depending on radiation-frequency detuning can be described successfully. Results show that the error introduced by using the one-dimensional approximation for multilevel particles turns out to be more significant than for the simplest particle models.

  1. Solving Kinetic Equations for the Laser Flash Photolysis Experiment on Nitric Oxide Synthases: Effect of Conformational Dynamics on the Interdomain Electron Transfer.

    PubMed

    Astashkin, Andrei V; Feng, Changjian

    2015-11-12

    The production of nitric oxide by the nitric oxide synthase (NOS) enzyme depends on the interdomain electron transfer (IET) between the flavin mononucleotide (FMN) and heme domains. Although the rate of this IET has been measured by laser flash photolysis (LFP) for various NOS proteins, no rigorous analysis of the relevant kinetic equations was performed so far. In this work, we provide an analytical solution of the kinetic equations underlying the LFP approach. The derived expressions reveal that the bulk IET rate is significantly affected by the conformational dynamics that determines the formation and dissociation rates of the docking complex between the FMN and heme domains. We show that in order to informatively study the electron transfer across the NOS enzyme, LFP should be used in combination with other spectroscopic methods that could directly probe the docking equilibrium and the conformational change rate constants. The implications of the obtained analytical expressions for the interpretation of the LFP results from various native and modified NOS proteins are discussed. The mathematical formulas derived in this work should also be applicable for interpreting the IET kinetics in other modular redox enzymes. PMID:26477677

  2. Inverse Kinetics

    Energy Science and Technology Software Center (ESTSC)

    2000-03-20

    Given the space-independent, one energy group reactor kinetics equations and the initial conditions, this prgram determines the time variation of reactivity required to produce the given input of flux-time data.

  3. Quantitative determination of the steady-state kinetics of multienzyme reactions using the algebraic rate equations for the component single-enzyme reactions.

    PubMed Central

    Stoner, C D

    1993-01-01

    Methods are given whereby the steady-state kinetic characteristics of multienzyme reactions consisting of individual single-enzyme reactions linked by freely diffusible intermediates can be determined quantitatively from the experimentally determined complete algebraic rate equations for the individual reactions. The approach is based on the fact that a valid steady-state rate equation for such a multienzyme reaction, in terms of the rate equations for the individual reactions, can be obtained simply from knowledge of the relative rates of the individual reactions when the multienzyme reaction is in the steady state. A number of model multienzyme reactions, which differ as to structural arrangement of the individual reactions, are examined by this approach. Simple mathematical methods which are applicable to most of these models are given for direct calculation of dependent variables. It is either pointed out or demonstrated with Mathematica that the rate equations for all of these models can be handled very easily with the aid of a personal computer equipped with appropriate equation-solving software. Since the approach permits evaluation of all dependent variables for any specific combination of values for the kinetic parameters and independent variables, numerical values for the flux control coefficients of the individual enzymes can be obtained by direct calculation for a wide variety of conditions and can be compared with those obtained according to the methods of Metabolic Control Analysis. Several such comparisons have been made and in all cases identical results were obtained. The intuitive notion that the individual enzymes of a multienzyme reaction would be equally rate limiting if the total amount of enzyme were being used with maximum efficiency is tested and shown to be incorrect. In the course of this test the flux control coefficient for the individual enzymes were found to be appropriate indicators of relative rate limitation or control by the

  4. Use of Coupled Rate Equations to Model NIR-to-Visible Upconversion Kinetics in Er3+, Yb3+: NaYF4 Nanocrystals

    NASA Astrophysics Data System (ADS)

    Yao, Ge

    The Er3+, Yb3+: NaYF4 system is the most efficient upconversion (UC) phosphor known, and yet the kinetic details of the mechanism responsible for upconversion are poorly understood. In this work, the dynamics of the photo-physical processes leading to NIR-to-visible upconversion luminescence in Er3+, Yb3+: NaYF 4 nanocrystals are investigated using a coupled-rate-equation model. The rate equations used in the simulation contain parameters representing the microscopic rate constants for individual mechanistic steps. Following initial population of the excited Yb3+ 2F5/2 state in the NIR, the population density of the excited states of the Er 3+ and Yb3+ ions are propagated as a function of time. Experimental spectroscopic characterization performed in our research group provides the data with which to study the kinetics of the UC process. The data include the time evolution of green, red, 1.0microm, and 1.5microm emission following pulsed excitation (decay curves), and the relative integrated intensity ratios of green, red, 1.0microm, and 1.5microm emission as a function of excitation power. Simulated and experimental decay curves and intensity ratios are compared to determine the optimum set of parameterized kinetic rate constants. The curve fits and integrated-relative-intensity-ratio fits are optimized using a Nelder-Meade simplex search method, in which the parameters are chosen to minimize the chi2 value calculated from the difference of the simulated and measure intensity ratios and the simulated and measured decay curves. The optimized parameter values obtained from the fits are physically reasonable, which suggests that the rate-equation modeling is an appropriate approach to describe the UC process in Er 3+, Yb3+: NaYF4 nanocrystals.

  5. Gaussian short-time propagators and electron kinetics: Numerical evaluation of path-sum solutions to Fokker{endash}Planck equations for rf heating and current drive

    SciTech Connect

    Bizarro, J.P.; Belo, J.H.; Figueiredo, A.C.

    1997-06-01

    Knowing that short-time propagators for Fokker{endash}Planck equations are Gaussian, and based on a path-sum formulation, an efficient and simple numerical method is presented to solve the initial-value problem for electron kinetics during rf heating and current drive. The formulation is thoroughly presented and discussed, its advantages are stressed, and general, practical criteria for its implementation are derived regarding the time step and grid spacing. The new approach is illustrated and validated by solving the one-dimensional model for lower-hybrid current drive, which has a well-known steady-state analytical solution. {copyright} {ital 1997 American Institute of Physics.}

  6. Model of 2,3-bisphosphoglycerate metabolism in the human erythrocyte based on detailed enzyme kinetic equations: equations and parameter refinement.

    PubMed

    Mulquiney, P J; Kuchel, P W

    1999-09-15

    Over the last 25 years, several mathematical models of erythrocyte metabolism have been developed. Although these models have identified the key features in the regulation and control of erythrocyte metabolism, many important aspects remain unexplained. In particular, none of these models have satisfactorily accounted for 2,3-bisphosphoglycerate (2,3-BPG) metabolism. 2,3-BPG is an important modulator of haemoglobin oxygen affinity, and hence an understanding of the regulation of 2,3-BPG concentration is important for understanding blood oxygen transport. A detailed, comprehensive, and hence realistic mathematical model of erythrocyte metabolism is presented that can explain the regulation and control of 2,3-BPG concentration and turnover. The model is restricted to the core metabolic pathways, namely glycolysis, the 2,3-BPG shunt and the pentose phosphate pathway (PPP), and includes membrane transport of metabolites, the binding of metabolites to haemoglobin and Mg(2+), as well as pH effects on key enzymic reactions and binding processes. The model is necessarily complex, since it is intended to describe the regulation and control of 2,3-BPG metabolism under a wide variety of physiological and experimental conditions. In addition, since H(+) and blood oxygen tension are important external effectors of 2,3-BPG concentration, it was important that the model take into account the large array of kinetic and binding phenomena that result from changes in these effectors. Through an iterative loop of experimental and simulation analysis many values of enzyme-kinetic parameters of the model were refined to yield close conformity between model simulations and 'real' experimental data. This iterative process enabled a single set of parameters to be found which described well the metabolic behaviour of the erythrocyte under a wide variety of conditions. PMID:10477269

  7. Model of 2,3-bisphosphoglycerate metabolism in the human erythrocyte based on detailed enzyme kinetic equations: equations and parameter refinement.

    PubMed Central

    Mulquiney, P J; Kuchel, P W

    1999-01-01

    Over the last 25 years, several mathematical models of erythrocyte metabolism have been developed. Although these models have identified the key features in the regulation and control of erythrocyte metabolism, many important aspects remain unexplained. In particular, none of these models have satisfactorily accounted for 2,3-bisphosphoglycerate (2,3-BPG) metabolism. 2,3-BPG is an important modulator of haemoglobin oxygen affinity, and hence an understanding of the regulation of 2,3-BPG concentration is important for understanding blood oxygen transport. A detailed, comprehensive, and hence realistic mathematical model of erythrocyte metabolism is presented that can explain the regulation and control of 2,3-BPG concentration and turnover. The model is restricted to the core metabolic pathways, namely glycolysis, the 2,3-BPG shunt and the pentose phosphate pathway (PPP), and includes membrane transport of metabolites, the binding of metabolites to haemoglobin and Mg(2+), as well as pH effects on key enzymic reactions and binding processes. The model is necessarily complex, since it is intended to describe the regulation and control of 2,3-BPG metabolism under a wide variety of physiological and experimental conditions. In addition, since H(+) and blood oxygen tension are important external effectors of 2,3-BPG concentration, it was important that the model take into account the large array of kinetic and binding phenomena that result from changes in these effectors. Through an iterative loop of experimental and simulation analysis many values of enzyme-kinetic parameters of the model were refined to yield close conformity between model simulations and 'real' experimental data. This iterative process enabled a single set of parameters to be found which described well the metabolic behaviour of the erythrocyte under a wide variety of conditions. PMID:10477269

  8. Development and analysis of prognostic equations for mesoscale kinetic energy and mesoscale (subgrid scale) fluxes for large-scale atmospheric models

    NASA Technical Reports Server (NTRS)

    Avissar, Roni; Chen, Fei

    1993-01-01

    Generated by landscape discontinuities (e.g., sea breezes) mesoscale circulation processes are not represented in large-scale atmospheric models (e.g., general circulation models), which have an inappropiate grid-scale resolution. With the assumption that atmospheric variables can be separated into large scale, mesoscale, and turbulent scale, a set of prognostic equations applicable in large-scale atmospheric models for momentum, temperature, moisture, and any other gaseous or aerosol material, which includes both mesoscale and turbulent fluxes is developed. Prognostic equations are also developed for these mesoscale fluxes, which indicate a closure problem and, therefore, require a parameterization. For this purpose, the mean mesoscale kinetic energy (MKE) per unit of mass is used, defined as E-tilde = 0.5 (the mean value of u'(sub i exp 2), where u'(sub i) represents the three Cartesian components of a mesoscale circulation (the angle bracket symbol is the grid-scale, horizontal averaging operator in the large-scale model, and a tilde indicates a corresponding large-scale mean value). A prognostic equation is developed for E-tilde, and an analysis of the different terms of this equation indicates that the mesoscale vertical heat flux, the mesoscale pressure correlation, and the interaction between turbulence and mesoscale perturbations are the major terms that affect the time tendency of E-tilde. A-state-of-the-art mesoscale atmospheric model is used to investigate the relationship between MKE, landscape discontinuities (as characterized by the spatial distribution of heat fluxes at the earth's surface), and mesoscale sensible and latent heat fluxes in the atmosphere. MKE is compared with turbulence kinetic energy to illustrate the importance of mesoscale processes as compared to turbulent processes. This analysis emphasizes the potential use of MKE to bridge between landscape discontinuities and mesoscale fluxes and, therefore, to parameterize mesoscale fluxes

  9. On the interpretation of the inverted kinetics equation and space-time calculations of the effectiveness of the VVER-1000 reactor scram system

    NASA Astrophysics Data System (ADS)

    Zizin, M. N.; Ivanov, L. D.

    2013-12-01

    In the present paper, an attempt is made to analyze the accuracy of calculating the effectiveness of the VVER-1000 reactor scram system by means of the inverted solution of the kinetics equation (ISKE). In the numerical studies in the intellectual ShIPR software system, the actuation of the reactor scram system with the possible jamming of one of the two most effective rods is simulated. First, the connection of functionals calculated in the space-time computation in different approximations with the kinetics equation is considered on the theoretical level. The formulas are presented in a manner facilitating their coding. Then, the results of processing of several such functions by the ISKE are presented. For estimating the effectiveness of the VVER-1000 reactor scram system, it is proposed to use the measured currents of ionization chambers (IC) jointly with calculated readings of IC imitators. In addition, the integral of the delayed neutron (DN) generation rate multiplied by the adjoint DN source over the volume of the reactor, calculated for the instant of time when insertion of safety rods ends, is used. This integral is necessary for taking into account the spatial reactivity effects. Reasonable agreement was attained for the considered example between the effectiveness of the scram system evaluated by this method and the values obtained by steady-state calculations as the difference of the reciprocal effective multiplication factors with withdrawn and inserted control rods. This agreement was attained with the use of eight-group DN parameters.

  10. Solution of the Krook kinetic equation model and non-equilibrium thermodynamics of a rarefied gas affected by a non-linear thermal radiation field

    NASA Astrophysics Data System (ADS)

    Abourabia, Aly Maher; Wahid, Taha Zakaraia Abdel

    2011-05-01

    A new approach for studying the influence of a thermal radiation field upon a rarefied neutral gas is introduced. We insert the radiation field effect in the force term of the Boltzmann equation. In a frame co-moving with the fluid, the BGK (Bhatnager-Gross-Krook) model kinetic equation is applied analytically. The one-dimensional steady problem is studied using the Liu-Lees model. We apply the moment method to follow the behavior of the macroscopic properties of the gas, such as the temperature and concentration. They are substituted into the corresponding two-stream Maxwellian distribution functions, permitting the investigation of the non-equilibrium thermodynamic properties of the system (gas + heated plate). The entropy, entropy flux, entropy production, thermodynamic forces and the kinetic coefficients are obtained. We verify the celebrated Onsager reciprocity relations for the system. The ratios between the different contributions of the internal energy changes based upon the total derivatives of the extensive parameters are estimated via the Gibbs formula. The results are applied to the Helium gas for various radiation field intensities due to different plate temperatures. Figures illustrating the calculated variables are drawn to predict their behavior and the results are discussed.

  11. Modeling and simulation of Streptomyces peucetius var. caesius N47 cultivation and epsilon-rhodomycinone production with kinetic equations and neural networks.

    PubMed

    Kiviharju, Kristiina; Salonen, Kalle; Leisola, Matti; Eerikäinen, Tero

    2006-11-10

    This study focuses on comparing different kinetic growth models and the use of neural networks in the batch cultivation of Streptomyces peucetius var. caesius producing epsilon-rhodomycinone. Contois, Monod and Teissier microbial growth models were used as well as the logistic growth modeling approach, which was found best in the simulations of growth and glucose consumption in the batch growth phase. The lag phase was included in the kinetic model with a CO2 trigger and a delay factor. Substrate consumption and product formation were included as Luedeking-Piret and logistic type equations, respectively. Biomass formation was modeled successfully with a 6-8-2 network, and the network was capable of biomass prediction with an R2-value of 0.983. Epsilon-rhodomycinone production was successfully modeled with a recursive 8-3-1 network capable of epsilon-rhodomycinone prediction with an R2-value of 0.903. The predictive power of the neural networks was superior to the kinetic models, which could not be used in predictive modeling of arbitrary batch cultivations. PMID:16797766

  12. Extending the range of validity of Fourier's law into the kinetic transport regime via asymptotic solution of the phonon Boltzmann transport equation

    NASA Astrophysics Data System (ADS)

    Péraud, Jean-Philippe M.; Hadjiconstantinou, Nicolas G.

    2016-01-01

    We derive the continuum equations and boundary conditions governing phonon-mediated heat transfer in the limit of a small but finite mean-free path from the asymptotic solution of the linearized Boltzmann equation in the relaxation time approximation. Our approach uses the ratio of the mean-free path to the characteristic system length scale, also known as the Knudsen number, as the expansion parameter to study the effects of boundaries on the breakdown of the Fourier description. We show that, in the bulk, the traditional heat conduction equation using Fourier's law as a constitutive relation is valid at least up to second order in the Knudsen number for steady problems and first order for time-dependent problems. However, this description does not hold within distances on the order of a few mean-free paths from the boundary; this breakdown is a result of kinetic effects that are always present in the boundary vicinity and require solution of a Boltzmann boundary layer problem to be determined. Matching the inner, boundary layer solution to the outer, bulk solution yields boundary conditions for the Fourier description as well as additive corrections in the form of universal kinetic boundary layers; both are found to be proportional to the bulk-solution gradients at the boundary and parametrized by the material model and the phonon-boundary interaction model (Boltzmann boundary condition). Our derivation shows that the traditional no-jump boundary condition for prescribed temperature boundaries and the no-flux boundary condition for diffusely reflecting boundaries are appropriate only to zeroth order in the Knudsen number; at higher order, boundary conditions are of the jump type. We illustrate the utility of the asymptotic solution procedure by demonstrating that it can be used to predict the Kapitza resistance (and temperature jump) associated with an interface between two materials. All results are validated via comparisons with low-variance deviational Monte

  13. Energy approach to kinetics equations for dislocations and twins and its application for high strain rate collision problems

    NASA Astrophysics Data System (ADS)

    Borodin, E. N.; Mayer, A. E.

    2015-11-01

    A computational plasticity model with accounting of coupled evolution of the dislocations and twins in metals under the dynamic loading is presented. The model is based on our previous results for the dislocation plasticity, but generalizes them and accounts mechanical twinning in addition. It includes equations of the mechanics of continua for elastic-plastic medium, where the plastic deformation tensor is determined through the structural defects evolution in the material. The model is self-consistent and allows determining of mechanical properties in wide range of strain rates and thermodynamic conditions as well as modification of the defect subsystems. The equations and parameters, its numerical implementation and some of obtained results are presented.

  14. Direct discrete simulation of the kinetic Landau-Fokker-Planck equation with an alternating external electromagnetic field

    NASA Astrophysics Data System (ADS)

    Karpov, S. A.; Potapenko, I. F.

    2015-10-01

    A stochastic method of simulation of Coulomb interaction is considered. The main idea of the method is to approximate the nonlinear Landau kinetic collision integral by the Boltzmann integral. In its realization, the method can be attributed to a wide class of Monte Carlo-type methods. It is easily combined with the existing particle methods used to simulate collisionless plasmas. This is important for simulation of the dynamics of both laboratory and space plasmas when the mean free path of plasma particles is comparable with the plasma inhomogeneity scale length. Illustrative examples of relaxation of two-temperature plasma being subject to a high-frequency alternating electric field are given, and differences from their classical description are considered. The method satisfies the conservation laws for the number of particles, momentum, and energy and is simple and efficient in implementation.

  15. Direct discrete simulation of the kinetic Landau-Fokker-Planck equation with an alternating external electromagnetic field

    SciTech Connect

    Karpov, S. A.; Potapenko, I. F.

    2015-10-15

    A stochastic method of simulation of Coulomb interaction is considered. The main idea of the method is to approximate the nonlinear Landau kinetic collision integral by the Boltzmann integral. In its realization, the method can be attributed to a wide class of Monte Carlo-type methods. It is easily combined with the existing particle methods used to simulate collisionless plasmas. This is important for simulation of the dynamics of both laboratory and space plasmas when the mean free path of plasma particles is comparable with the plasma inhomogeneity scale length. Illustrative examples of relaxation of two-temperature plasma being subject to a high-frequency alternating electric field are given, and differences from their classical description are considered. The method satisfies the conservation laws for the number of particles, momentum, and energy and is simple and efficient in implementation.

  16. Solution of a kinetic equation for diatomic gas with the use of differential scattering cross sections computed by the method of classical trajectories

    NASA Astrophysics Data System (ADS)

    Anikin, Yu. A.; Dodulad, O. I.

    2013-07-01

    A collision integral is constructed taking into account the rotational degrees of freedom of the gas molecules. Its truncation error is shown to be second order in the rotational velocity mesh size. In the solution of the kinetic equation, the resulting collision integral is directly computed using a projection method. Preliminarily, the differential scattering cross sections of nitrogen molecules are computed by applying the method of classical trajectories. The resulting cross section values are tabulated in multimillion data arrays. The one-dimensional problems of shock wave structure and heat transfer between two plates are computed as tests, and the results are compared with experimental data. The convergence of the results with decreasing rotational velocity mesh size is analyzed.

  17. Kinetic description of intense nonneutral beam propagation through a periodic solenoidal focusing field based on the nonlinear Vlasov-Maxwell equations

    SciTech Connect

    Davidson, R.C.; Chen, C.

    1997-08-01

    A kinetic description of intense nonneutral beam propagation through a periodic solenoidal focusing field B{sup sol}({rvec x}) is developed. The analysis is carried out for a thin beam with characteristic beam radius r{sub b} {much_lt} S, and directed axial momentum {gamma}{sub b}m{beta}{sub b}c (in the z-direction) large compared with the transverse momentum and axial momentum spread of the beam particles. Making use of the nonlinear Vlasov-Maxwell equations for general distribution function f{sub b}({rvec x},{rvec p},t) and self-consistent electrostatic field consistent with the thin-beam approximation, the kinetic model is used to investigate detailed beam equilibrium properties for a variety of distribution functions. Examples are presented both for the case of a uniform solenoidal focusing field B{sub z}(z) = B{sub 0} = const. and for the case of a periodic solenoidal focusing field B{sub z}(z + S) = B{sub z}(z). The nonlinear Vlasov-Maxwell equations are simplified in the thin-beam approximation, and an alternative Hamiltonian formulation is developed that is particularly well-suited to intense beam propagation in periodic focusing systems. Based on the present analysis, the Vlasov-Maxwell description of intense nonneutral beam propagation through a periodic solenoidal focusing field {rvec B}{sup sol}({rvec x}) is found to be remarkably tractable and rich in physics content. The Vlasov-Maxwell formalism developed here can be extended in a straightforward manner to investigate detailed stability behavior for perturbations about specific choices of beam equilibria.

  18. Effect of the adsorbate kinetic diameter on the accuracy of the Dubinin-Radushkevich equation for modeling adsorption of organic vapors on activated carbon.

    PubMed

    Jahandar Lashaki, Masoud; Fayaz, Mohammadreza; Niknaddaf, Saeid; Hashisho, Zaher

    2012-11-30

    This paper investigates the effect of the kinetic diameter (KD) of the reference adsorbate on the accuracy of the Dubinin-Radushkevich (D-R) equation for predicting the adsorption isotherms of organic vapors on microporous activated carbon. Adsorption isotherms for 13 organic compounds on microporous beaded activated carbon were experimentally measured, and predicted using the D-R model and affinity coefficients. The affinity coefficients calculated based on molar volumes, molecular polarizabilities, and molecular parachors were used to predict the isotherms based on four reference compounds (4.3≤KD≤6.8 Å). The results show that the affinity coefficients are independent of the calculation method if the reference and test adsorbates are from the same organic group. Choosing a reference adsorbate with a KD similar to that of the test adsorbate results in better prediction of the adsorption isotherm. The relative error between the predicted and the measured adsorption isotherms increases as the absolute difference in the kinetic diameters of the reference and test adsorbates increases. Finally, the proposed hypothesis was used to explain reports of inconsistent findings among published articles. The results from this study are important because they allow a more accurate prediction of adsorption capacities of adsorbents which allow for better design of adsorption systems. PMID:23044198

  19. Stability of Planar Fronts for a Non-Local Phase Kinetics Equation with a Conservation Law in D ≤ 3

    NASA Astrophysics Data System (ADS)

    Carlen, Eric A.; Orlandi, Enza

    2012-05-01

    We consider, in a D-dimensional cylinder, a non-local evolution equation that describes the evolution of the local magnetization in a continuum limit of an Ising spin system with Kawasaki dynamics and Kac potentials. We consider sub-critical temperatures, for which there are two local spatially homogeneous equilibria, and show a local nonlinear stability result for the minimum free energy profiles for the magnetization at the interface between regions of these two different local equilibrium; i.e. the planar fronts: We show that an initial perturbation of a front that is sufficiently small in L2 norm, and sufficiently localized yields a solution that relaxes to another front, selected by a conservation law, in the L1 norm at an algebraic rate that we explicitly estimate. We also obtain rates for the relaxation in the L2 norm and the rate of decrease of the excess free energy.

  20. Kinetic model framework for aerosol and cloud surface chemistry and gas-particle interactions: Part 1 - general equations, parameters, and terminology

    NASA Astrophysics Data System (ADS)

    Pöschl, U.; Rudich, Y.; Ammann, M.

    2005-04-01

    Aerosols and clouds play central roles in atmospheric chemistry and physics, climate, air pollution, and public health. The mechanistic understanding and predictability of aerosol and cloud properties, interactions, transformations, and effects are, however, still very limited. This is due not only to the limited availability of measurement data, but also to the limited applicability and compatibility of model formalisms used for the analysis, interpretation, and description of heterogeneous and multiphase processes. To support the investigation and elucidation of atmospheric aerosol and cloud surface chemistry and gas-particle interactions, we present a comprehensive kinetic model framework with consistent and unambiguous terminology and universally applicable rate equations and parameters. It allows to describe mass transport and chemical reactions at the gas-particle interface and to link aerosol and cloud surface processes with gas phase and particle bulk processes in systems with multiple chemical components and competing physicochemical processes. The key elements and essential aspects of the presented framework are: a simple and descriptive double-layer surface model (sorption layer and quasi-static layer); straightforward flux-based mass balance and rate equations; clear separation of mass transport and chemical reactions; well-defined rate parameters (uptake and accommodation coefficients, reaction and transport rate coefficients); clear distinction between gas phase, gas-surface, and surface-bulk transport (gas phase diffusion correction, surface and bulk accommodation); clear distinction between gas-surface, surface layer, and surface-bulk reactions (Langmuir-Hinshelwood and Eley-Rideal mechanisms); mechanistic description of concentration and time dependencies; flexible inclusion/omission of chemical species and physicochemical processes; flexible convolution/deconvolution of species and processes; and full compatibility with traditional resistor model

  1. Departures of the electron energy distribution from a Maxwellian in hydrogen. I - Formulation and solution of the electron kinetic equation. II - Consequences

    NASA Technical Reports Server (NTRS)

    Shoub, E. C.

    1977-01-01

    The problem of calculating the steady-state free-electron energy distribution in a hydrogen gas is considered in order to study departures of that distribution from a Maxwellian at sufficiently low degrees of ionization. A model kinetic equation is formulated and solved analytically for the one-particle electron distribution function in a steady-state partially ionized hydrogen gas, and it is shown that the formal solution can be accurately approximated by using the WKB method. The solutions obtained indicate that the high-energy tail of the distribution is susceptible to distortion by imbalanced inelastic collisions for ionization fractions not exceeding about 0.1 and that such departures from a Maxwellian can lead to significant changes in the collisional excitation and ionization rates of ground-state hydrogen atoms. Expressions for the electron-hydrogen collision rates are derived which explicitly display their dependence on the hydrogen departure coefficients. The results are applied in order to compare self-consistent predictions with those based on the a priori assumption of a Maxwellian distribution for models of the thermal ionization equilibrium of hydrogen in the optically thin limit, spectral-line formation by a gas consisting of two-level atoms, and radiative transfer in finite slabs by a gas of four-level hydrogen atoms.

  2. The Chemical Master Equation Approach to Nonequilibrium Steady-State of Open Biochemical Systems: Linear Single-Molecule Enzyme Kinetics and Nonlinear Biochemical Reaction Networks

    PubMed Central

    Qian, Hong; Bishop, Lisa M.

    2010-01-01

    We develop the stochastic, chemical master equation as a unifying approach to the dynamics of biochemical reaction systems in a mesoscopic volume under a living environment. A living environment provides a continuous chemical energy input that sustains the reaction system in a nonequilibrium steady state with concentration fluctuations. We discuss the linear, unimolecular single-molecule enzyme kinetics, phosphorylation-dephosphorylation cycle (PdPC) with bistability, and network exhibiting oscillations. Emphasis is paid to the comparison between the stochastic dynamics and the prediction based on the traditional approach based on the Law of Mass Action. We introduce the difference between nonlinear bistability and stochastic bistability, the latter has no deterministic counterpart. For systems with nonlinear bistability, there are three different time scales: (a) individual biochemical reactions, (b) nonlinear network dynamics approaching to attractors, and (c) cellular evolution. For mesoscopic systems with size of a living cell, dynamics in (a) and (c) are stochastic while that with (b) is dominantly deterministic. Both (b) and (c) are emergent properties of a dynamic biochemical network; We suggest that the (c) is most relevant to major cellular biochemical processes such as epi-genetic regulation, apoptosis, and cancer immunoediting. The cellular evolution proceeds with transitions among the attractors of (b) in a “punctuated equilibrium” manner. PMID:20957107

  3. The chemical master equation approach to nonequilibrium steady-state of open biochemical systems: linear single-molecule enzyme kinetics and nonlinear biochemical reaction networks.

    PubMed

    Qian, Hong; Bishop, Lisa M

    2010-01-01

    We develop the stochastic, chemical master equation as a unifying approach to the dynamics of biochemical reaction systems in a mesoscopic volume under a living environment. A living environment provides a continuous chemical energy input that sustains the reaction system in a nonequilibrium steady state with concentration fluctuations. We discuss the linear, unimolecular single-molecule enzyme kinetics, phosphorylation-dephosphorylation cycle (PdPC) with bistability, and network exhibiting oscillations. Emphasis is paid to the comparison between the stochastic dynamics and the prediction based on the traditional approach based on the Law of Mass Action. We introduce the difference between nonlinear bistability and stochastic bistability, the latter has no deterministic counterpart. For systems with nonlinear bistability, there are three different time scales: (a) individual biochemical reactions, (b) nonlinear network dynamics approaching to attractors, and (c) cellular evolution. For mesoscopic systems with size of a living cell, dynamics in (a) and (c) are stochastic while that with (b) is dominantly deterministic. Both (b) and (c) are emergent properties of a dynamic biochemical network; We suggest that the (c) is most relevant to major cellular biochemical processes such as epi-genetic regulation, apoptosis, and cancer immunoediting. The cellular evolution proceeds with transitions among the attractors of (b) in a "punctuated equilibrium" manner. PMID:20957107

  4. A realizability-preserving discontinuous Galerkin scheme for entropy-based moment closures for linear kinetic equations in one space dimension

    NASA Astrophysics Data System (ADS)

    Alldredge, Graham; Schneider, Florian

    2015-08-01

    We implement a high-order numerical scheme for the entropy-based moment closure, the so-called MN model, for linear kinetic equations in slab geometry. A discontinuous Galerkin (DG) scheme in space along with a strong-stability preserving Runge-Kutta time integrator is a natural choice to achieve a third-order scheme, but so far, the challenge for such a scheme in this context is the implementation of a linear scaling limiter when the numerical solution leaves the set of realizable moments (that is, those moments associated with a positive underlying distribution). The difficulty for such a limiter lies in the computation of the intersection of a ray with the set of realizable moments. We avoid this computation by using quadrature to generate a convex polytope which approximates this set. The half-space representation of this polytope is used to compute an approximation of the required intersection straightforwardly, and with this limiter in hand, the rest of the DG scheme is constructed using standard techniques. We consider the resulting numerical scheme on a new manufactured solution and standard benchmark problems for both traditional MN models and the so-called mixed-moment models. The manufactured solution allows us to observe the expected convergence rates and explore the effects of the regularization in the optimization.

  5. Kinetic model framework for aerosol and cloud surface chemistry and gas-particle interactions - Part 1: General equations, parameters, and terminology

    NASA Astrophysics Data System (ADS)

    Pöschl, U.; Rudich, Y.; Ammann, M.

    2007-12-01

    Aerosols and clouds play central roles in atmospheric chemistry and physics, climate, air pollution, and public health. The mechanistic understanding and predictability of aerosol and cloud properties, interactions, transformations, and effects are, however, still very limited. This is due not only to the limited availability of measurement data, but also to the limited applicability and compatibility of model formalisms used for the analysis, interpretation, and description of heterogeneous and multiphase processes. To support the investigation and elucidation of atmospheric aerosol and cloud surface chemistry and gas-particle interactions, we present a comprehensive kinetic model framework with consistent and unambiguous terminology and universally applicable rate equations and parameters. It enables a detailed description of mass transport and chemical reactions at the gas-particle interface, and it allows linking aerosol and cloud surface processes with gas phase and particle bulk processes in systems with multiple chemical components and competing physicochemical processes. The key elements and essential aspects of the presented framework are: a simple and descriptive double-layer surface model (sorption layer and quasi-static layer); straightforward flux-based mass balance and rate equations; clear separation of mass transport and chemical reactions; well-defined and consistent rate parameters (uptake and accommodation coefficients, reaction and transport rate coefficients); clear distinction between gas phase, gas-surface, and surface-bulk transport (gas phase diffusion, surface and bulk accommodation); clear distinction between gas-surface, surface layer, and surface-bulk reactions (Langmuir-Hinshelwood and Eley-Rideal mechanisms); mechanistic description of concentration and time dependences (transient and steady-state conditions); flexible addition of unlimited numbers of chemical species and physicochemical processes; optional aggregation or resolution

  6. STOMP Subsurface Transport Over Multiple Phases Version 1.0 Addendum: ECKEChem Equilibrium-Conservation-Kinetic Equation Chemistry and Reactive Transport

    SciTech Connect

    White, Mark D.; McGrail, B. Peter

    2005-12-01

    flow and transport simulator, STOMP (Subsurface Transport Over Multiple Phases). Prior to these code development activities, the STOMP simulator included sequential and scalable implementations for numerically simulating the injection of supercritical CO2 into deep saline aquifers. Additionally, the sequential implementations included operational modes that considered nonisothermal conditions and kinetic dissolution of CO2 into the saline aqueous phase. This addendum documents the advancement of these numerical simulation capabilities to include reactive transport in the STOMP simulator through the inclusion of the recently PNNL developed batch geochemistry solution module ECKEChem (Equilibrium-Conservation-Kinetic Equation Chemistry). Potential geologic reservoirs for sequestering CO2 include deep saline aquifers, hydrate-bearing formations, depleted or partially depleted natural gas and petroleum reservoirs, and coal beds. The mechanisms for sequestering carbon dioxide in geologic reservoirs include physical trapping, dissolution in the reservoir fluids, hydraulic trapping (hysteretic entrapment of nonwetting fluids), and chemical reaction. This document and the associated code development and verification work are concerned with the chemistry of injecting CO2 into geologic reservoirs. As geologic sequestration of CO2 via chemical reaction, namely precipitation reactions, are most dominate in deep saline aquifers, the principal focus of this document is the numerical simulation of CO2 injection, migration, and geochemical reaction in deep saline aquifers. The ECKEChem batch chemistry module was developed in a fashion that would allow its implementation into all operational modes of the STOMP simulator, making it a more versatile chemistry component. Additionally, this approach allows for verification of the ECKEChem module against more classical reactive transport problems involving aqueous systems.

  7. Brownian motion from Boltzmann's equation.

    NASA Technical Reports Server (NTRS)

    Montgomery, D.

    1971-01-01

    Two apparently disparate lines of inquiry in kinetic theory are shown to be equivalent: (1) Brownian motion as treated by the (stochastic) Langevin equation and Fokker-Planck equation; and (2) Boltzmann's equation. The method is to derive the kinetic equation for Brownian motion from the Boltzmann equation for a two-component neutral gas by a simultaneous expansion in the density and mass ratios.

  8. O2(a1Δg) + Mg, Fe, and Ca: Experimental kinetics and formulation of a weak collision, multiwell master equation with spin-hopping

    NASA Astrophysics Data System (ADS)

    Plane, J. M. C.; Whalley, C. L.; Frances-Soriano, L.; Goddard, A.; Harvey, J. N.; Glowacki, D. R.; Viggiano, A. A.

    2012-07-01

    utilized density functional theory along with multireference and explicitly correlated CCSD(T)-F12 electronic structure calculations to examine the lowest lying singlet and triplet surfaces. In addition to mapping stationary points, we used a genetic algorithm to locate minimum energy crossing points between the two surfaces. Simulations of the Ca + O2(a) kinetics were then carried out using a combination of both standard and non-adiabatic Rice-Ramsperger-Kassel-Marcus (RRKM) theory implemented within a weak collision, multiwell master equation model. In terms of atmospheric significance, only in the case of Ca does reaction with O2(a) compete with O3 during the daytime between 85 and 110 km.

  9. A comparative examination of the adsorption mechanism of an anionic textile dye (RBY 3GL) onto the powdered activated carbon (PAC) using various the isotherm models and kinetics equations with linear and non-linear methods

    NASA Astrophysics Data System (ADS)

    Açıkyıldız, Metin; Gürses, Ahmet; Güneş, Kübra; Yalvaç, Duygu

    2015-11-01

    The present study was designed to compare the linear and non-linear methods used to check the compliance of the experimental data corresponding to the isotherm models (Langmuir, Freundlich, and Redlich-Peterson) and kinetics equations (pseudo-first order and pseudo-second order). In this context, adsorption experiments were carried out to remove an anionic dye, Remazol Brillant Yellow 3GL (RBY), from its aqueous solutions using a commercial activated carbon as a sorbent. The effects of contact time, initial RBY concentration, and temperature onto adsorbed amount were investigated. The amount of dye adsorbed increased with increased adsorption time and the adsorption equilibrium was attained after 240 min. The amount of dye adsorbed enhanced with increased temperature, suggesting that the adsorption process is endothermic. The experimental data was analyzed using the Langmuir, Freundlich, and Redlich-Peterson isotherm equations in order to predict adsorption isotherm. It was determined that the isotherm data were fitted to the Langmuir and Redlich-Peterson isotherms. The adsorption process was also found to follow a pseudo second-order kinetic model. According to the kinetic and isotherm data, it was found that the determination coefficients obtained from linear method were higher than those obtained from non-linear method.

  10. Two-dimensional equations of the surface harmonics method for solving problems of spatial neutron kinetics in square-lattice reactors

    SciTech Connect

    Boyarinov, V. F. Kondrushin, A. E. Fomichenko, P. A.

    2014-12-15

    Two-dimensional time-dependent finite-difference equations of the surface harmonics method (SHM) for the description of the neutron transport are derived for square-lattice reactors. These equations are implemented in the SUHAM-TD code. Verification of the derived equations and the developed code are performed by the example of known test problems, and the potential and efficiency of the SHM as applied to the solution of the time-dependent neutron transport equation in the diffusion approximation in two-dimensional geometry are demonstrated. These results show the substantial advantage of SHM over direct finite-difference modeling in computational costs.

  11. Two-dimensional equations of the surface harmonics method for solving problems of spatial neutron kinetics in square-lattice reactors

    NASA Astrophysics Data System (ADS)

    Boyarinov, V. F.; Kondrushin, A. E.; Fomichenko, P. A.

    2014-12-01

    Two-dimensional time-dependent finite-difference equations of the surface harmonics method (SHM) for the description of the neutron transport are derived for square-lattice reactors. These equations are implemented in the SUHAM-TD code. Verification of the derived equations and the developed code are performed by the example of known test problems, and the potential and efficiency of the SHM as applied to the solution of the time-dependent neutron transport equation in the diffusion approximation in two-dimensional geometry are demonstrated. These results show the substantial advantage of SHM over direct finite-difference modeling in computational costs.

  12. Solving Simple Kinetics without Integrals

    ERIC Educational Resources Information Center

    de la Pen~a, Lisandro Herna´ndez

    2016-01-01

    The solution of simple kinetic equations is analyzed without referencing any topic from differential equations or integral calculus. Guided by the physical meaning of the rate equation, a systematic procedure is used to generate an approximate solution that converges uniformly to the exact solution in the case of zero, first, and second order…

  13. Spatial kinetics in fast reactors

    NASA Astrophysics Data System (ADS)

    Seleznev, E. F.; Belov, A. A.; Panova, I. S.; Matvienko, I. P.; Zhukov, A. M.

    2013-12-01

    The analysis of the solution to the spatial nonstationary equation of neutron transport is presented by the example of a fast reactor. Experiments in spatial kinetics conducted recently at the complex of critical assemblies (fast physical stand) and computations of their data using the TIMER code (for solving the nonstationary equation in multidimensional diffusion approximation for direct and inverse problems of reactor kinetics) have shown that kinetics of fast reactors substantially differs from kinetics of thermal reactors. The difference is connected with influence of the delayed neutron spectrum on rates of the process in a fast reactor.

  14. Nuclear Reactor Kinetics and Control.

    SciTech Connect

    JEFFERY,; LEWINS, D.

    2009-07-27

    Version 00 Dr. J.D. Lewins has now released the following legacy book for free distribution: Nuclear Reactor Kinetics and Control, Pergamon Press, London, 275 pages, 1978. 1. Introductory Review 2. Neutron and Precursor Equations 3. Elementary Solutions of the Kinetics Equations at Low Power 4. Linear Reactor Process Dynamics with Feedback 5. Power Reactor Control Systems 6. Fluctuations and Reactor Noise 7. Safety and Reliability 8. Non Linear Systems; Stability and Control 9. Analogue Computing Addendum: Jay Basken and Jeffery D. Lewins: Power Series Solution of the Reactor Kinetics Equations, Nuclear Science and Engineering: 122, 407-436 (1996) (authorized for distribution with the book: courtesy of the American Nuclear Society)

  15. Algebraic Rate of Decay for the Excess Free Energy and Stability of Fronts for a Nonlocal Phase Kinetics Equation with a Conservation Law. I

    NASA Astrophysics Data System (ADS)

    Carlen, E. A.; Carvalho, M. C.; Orlandi, E.

    1999-06-01

    This is the first of two papers devoted to the study of a nonlocal evolution equation that describes the evolution of the local magnetization in a continuum limit of an Ising spin system with Kawasaki dynamics and Kac potentials. We consider subcritical temperatures, for which there are two local equilibria, and begin the proof of a local nonlinear stability result for the minimum free energy profiles for the magnetization at the interface between regions of these two different local equilibria; i.e., the fronts. We shall show in the second paper that an initial perturbation v 0 of a front that is sufficiently small in L 2 norm, and sufficiently localized that ∫ x 2 v 0( x)2 dx<∞, yields a solution that relaxes to another front, selected by a conservation law, in the L 1 norm at an algebraic rate that we explicitly estimate. There we also obtain rates for the relaxation in the L 2 norm and the rate of decrease of the excess free energy. Here we prove a number of estimates essential for this result. Moreover, the estimates proved here suffice to establish the main result in an important special case.

  16. Enskog-like kinetic models for vehicular traffic

    SciTech Connect

    Klar, A.; Wegener, R.

    1997-04-01

    In the present paper a general criticism of kinetic equations for vehicular traffic is given. The necessity of introducing an Enskog-type correction into these equations is shown. An Enskog-like kinetic traffic flow equation is presented and fluid dynamic equations are derived. This derivation yields new coefficients for the standard fluid dynamic equations of vehicular traffic. Numerical simulations for inhomogeneous traffic flow situations are shown together with a comparison between kinetic and fluid dynamics models.

  17. Kinetic Demonstration.

    ERIC Educational Resources Information Center

    Burgardt, Erik D.; Ryan, Hank

    1996-01-01

    Presents a unit on chemical reaction kinetics that consists of a predemonstration activity, the demonstration, and a set of postdemonstration activities that help students transfer the concepts to actual chemical reactions. Simulates various aspects of chemical reaction kinetics. (JRH)

  18. Kinetic Atom.

    ERIC Educational Resources Information Center

    Wilson, David B.

    1981-01-01

    Surveys the research of scientists like Joule, Kelvin, Maxwell, Clausius, and Boltzmann as it comments on the basic conceptual issues involved in the development of a more precise kinetic theory and the idea of a kinetic atom. (Author/SK)

  19. Enzyme Kinetics.

    ERIC Educational Resources Information Center

    Moe, Owen; Cornelius, Richard

    1988-01-01

    Conveys an appreciation of enzyme kinetic analysis by using a practical and intuitive approach. Discusses enzyme assays, kinetic models and rate laws, the kinetic constants (V, velocity, and Km, Michaels constant), evaluation of V and Km from experimental data, and enzyme inhibition. (CW)

  20. A Kinetic-fluid Model

    SciTech Connect

    First Author = C.Z. Cheng; Jay R. Johnson

    1998-07-10

    A nonlinear kinetic-fluid model for high-beta plasmas with multiple ion species which can be applied to multiscale phenomena is presented. The model embeds important kinetic effects due to finite ion Larmor radius (FLR), wave-particle resonances, magnetic particle trapping, etc. in the framework of simple fluid descriptions. When further restricting to low frequency phenomena with frequencies less than the ion cyclotron frequency the kinetic-fluid model takes a simpler form in which the fluid equations of multiple ion species collapse into single-fluid density and momentum equations and a low frequency generalized Ohm's law. The kinetic effects are introduced via plasma pressure tensors for ions and electrons which are computed from particle distribution functions that are governed by the Vlasov equation or simplified plasma dynamics equations such as the gyrokinetic equation. The ion FLR effects provide a finite parallel electric field, a perpendicular velocity that modifies the ExB drift, and a gyroviscosity tensor, all of which are neglected in the usual one-fluid MHD description. Eigenmode equations are derived which include magnetosphere-ionosphere coupling effects for low frequency waves (e.g., kinetic/inertial Alfven waves and ballooning-mirror instabilities).

  1. Nonlinear gyrokinetic equations

    SciTech Connect

    Dubin, D.H.E.; Krommes, J.A.; Oberman, C.; Lee, W.W.

    1983-03-01

    Nonlinear gyrokinetic equations are derived from a systematic Hamiltonian theory. The derivation employs Lie transforms and a noncanonical perturbation theory first used by Littlejohn for the simpler problem of asymptotically small gyroradius. For definiteness, we emphasize the limit of electrostatic fluctuations in slab geometry; however, there is a straight-forward generalization to arbitrary field geometry and electromagnetic perturbations. An energy invariant for the nonlinear system is derived, and various of its limits are considered. The weak turbulence theory of the equations is examined. In particular, the wave kinetic equation of Galeev and Sagdeev is derived from an asystematic truncation of the equations, implying that this equation fails to consider all gyrokinetic effects. The equations are simplified for the case of small but finite gyroradius and put in a form suitable for efficient computer simulation. Although it is possible to derive the Terry-Horton and Hasegawa-Mima equations as limiting cases of our theory, several new nonlinear terms absent from conventional theories appear and are discussed.

  2. Nuclear Reactor Kinetics and Control.

    Energy Science and Technology Software Center (ESTSC)

    2009-07-27

    Version 00 Dr. J.D. Lewins has now released the following legacy book for free distribution: Nuclear Reactor Kinetics and Control, Pergamon Press, London, 275 pages, 1978. 1. Introductory Review 2. Neutron and Precursor Equations 3. Elementary Solutions of the Kinetics Equations at Low Power 4. Linear Reactor Process Dynamics with Feedback 5. Power Reactor Control Systems 6. Fluctuations and Reactor Noise 7. Safety and Reliability 8. Non Linear Systems; Stability and Control 9. Analogue Computingmore » Addendum: Jay Basken and Jeffery D. Lewins: Power Series Solution of the Reactor Kinetics Equations, Nuclear Science and Engineering: 122, 407-436 (1996) (authorized for distribution with the book: courtesy of the American Nuclear Society)« less

  3. Kinetic theory of vehicular traffic

    NASA Astrophysics Data System (ADS)

    Iannini, M. L. L.; Dickman, Ronald

    2016-02-01

    We review the kinetic theory of traffic proposed by Prigogine and Herman in which the Boltzmann equation is adapted to vehicular traffic. The kinetic equation and its solution are discussed, and a novel distribution of desired velocities that is more suitable for describing real traffic conditions is analyzed. We also study the stationary velocity distribution at the transition between individual and collective flow patterns. At this transition, the distribution splits into a smoothly varying regular part, in which vehicles have nonzero velocities, and a singular one, corresponding to stopped vehicles. Computational methods for obtaining the stationary velocity distribution and the full space-time evolution of the vehicular distribution are explained.

  4. Relativistic kinetic theory of magnetoplasmas

    SciTech Connect

    Beklemishev, Alexei; Nicolini, Piero; Tessarotto, Massimo

    2005-05-16

    Recently, an increasing interest in astrophysical as well as laboratory plasmas has been manifested in reference to the existence of relativistic flows, related in turn to the production of intense electric fields in magnetized systems. Such phenomena require their description in the framework of a consistent relativistic kinetic theory, rather than on relativistic MHD equations, subject to specific closure conditions. The purpose of this work is to apply the relativistic single-particle guiding-center theory developed by Beklemishev and Tessarotto, including the nonlinear treatment of small-wavelength EM perturbations which may naturally arise in such systems. As a result, a closed set of relativistic gyrokinetic equations, consisting of the collisionless relativistic kinetic equation, expressed in hybrid gyrokinetic variables, and the averaged Maxwell's equations, is derived for an arbitrary four-dimensional coordinate system.

  5. Kinetics of reactive wetting

    SciTech Connect

    Yost, F.G.

    2000-04-14

    The importance of interfacial processes in materials joining has a long history. A significant amount of work has suggested that processes collateral to wetting can affect the extent of wetting and moderate or retard wetting rate. Even very small additions of a constituent, known to react with the substrate, cause pronounced improvement in wetting and are exploited in braze alloys, especially those used for joining to ceramics. In the following a model will be constructed for the wetting kinetics of a small droplet of metal containing a constituent that diffuses to the wetting line and chemically reacts with a flat, smooth substrate. The model is similar to that of Voitovitch et al. and Mortensen et al. but incorporates chemical reaction kinetics such that the result contains both diffusion and reaction kinetics. The model is constructed in the circular cylinder coordinate system, satisfies the diffusion equation under conditions of slow flow, and considers diffusion and reaction at the wetting line to be processes in series. This is done by solving the diffusion equation with proper initial and boundary conditions, computing the diffusive flux at the wetting line, and equating this to both the convective flux and reaction flux. This procedure is similar to equating the current flowing in components of a series circuit. The wetting rate will be computed versus time for a variety of diffusion and reaction conditions. A transition is observed from nonlinear (diffusive) to linear (reactive) behavior as the control parameters (such as the diffusion coefficient) are modified. This is in agreement with experimental observations. The adequacy of the slow flow condition, used in this type of analysis, is discussed and an amended procedure is suggested.

  6. Multiflow approach to plasma kinetics

    SciTech Connect

    Ignatov, A. M.

    2015-10-15

    Instead of the commonly used Vlasov equation, one is able to treat kinetic phenomena in collisionless plasma with the help of the infinite set of hydrodynamic equations. The present paper deals with the linear approximation of multiflow hydrodynamics. It is shown that single-particle and collective excitations analogous to Van Kampen waves are explicitly separated. Expressions for the energy of all eigenmodes are obtained.

  7. Analytic solutions of the relativistic Boltzmann equation

    NASA Astrophysics Data System (ADS)

    Hatta, Yoshitaka; Martinez, Mauricio; Xiao, Bo-Wen

    2015-04-01

    We present new analytic solutions to the relativistic Boltzmann equation within the relaxation time approximation. We first obtain spherically expanding solutions which are the kinetic counterparts of the exact solutions of the Israel-Stewart equation in the literature. This allows us to compare the solutions of the kinetic and hydrodynamic equations at an analytical level. We then derive a novel boost-invariant solution of the Boltzmann equation which has an unconventional dependence on the proper time. The existence of such a solution is also suggested in second-order hydrodynamics and fluid-gravity correspondence.

  8. THE BERNOULLI EQUATION AND COMPRESSIBLE FLOW THEORIES

    EPA Science Inventory

    The incompressible Bernoulli equation is an analytical relationship between pressure, kinetic energy, and potential energy. As perhaps the simplest and most useful statement for describing laminar flow, it buttresses numerous incompressible flow models that have been developed ...

  9. Kinetics of Reactive Wetting

    SciTech Connect

    YOST, FREDERICK G.

    1999-09-09

    The importance of interfacial processes in materials joining has a long history. A significant amount of work has suggested that processes collateral to wetting can affect the extent of wetting and moderate or retard wetting rate. Even very small additions of a constituent, known to react with the substrate, cause pronounced improvement in wetting and are exploited in braze alloys, especially those used for joining to ceramics. The wide diversity of processes, such as diffusion, chemical reaction, and fluxing, and their possible combinations suggest that various rate laws should be expected for wetting kinetics depending on the controlling processes. These rate laws are expected to differ crucially from the standard fluid controlled wetting models found in the literature. Voitovitch et al. and Mortensen et al. have shown data that suggests diffusion control for some systems and reaction control for others. They also presented a model of wetting kinetics controlled by the diffusion of a constituent contained by the wetting fluid. In the following a model will be constructed for the wetting kinetics of a small droplet of metal containing a constituent that diffuses to the wetting line and chemically reacts with a flat, smooth substrate. The model is similar to that of Voitovitch et al. and Mortensen et al. but incorporates chemical reaction kinetics such that the result contains both diffusion and reaction kinetics. The model is constructed in the circular cylinder coordinate system, satisfies the diffusion equation under conditions of slow flow, and considers diffusion and reaction at the wetting line to be processes in series. This is done by solving the diffusion equation with proper initial and boundary conditions, computing the diffusive flux at the wetting line and equating this to both the convective flux and reaction flux. This procedure is similar to equating the current flowing in components of a series circuit. The wetting rate will be computed versus time

  10. O{sub 2}(a{sup 1}{Delta}{sub g}) + Mg, Fe, and Ca: Experimental kinetics and formulation of a weak collision, multiwell master equation with spin-hopping

    SciTech Connect

    Plane, J. M. C.; Whalley, C. L.; Goddard, A.; Frances-Soriano, L.; Harvey, J. N.; Glowacki, D. R.; Viggiano, A. A.

    2012-07-07

    complex, suggesting that this reaction proceeds mostly by near-resonant electronic energy transfer to Fe(a{sup 5}F) + O{sub 2}(X). The reaction of Ca + O{sub 2}(a) occurs in an intermediate regime with two competing pressure dependent channels: (1) a recombination to produce CaO{sub 2}({sup 1}A{sub 1}), and (2) a singlet/triplet non-adiabatic hopping channel leading to CaO + O({sup 3}P). In order to interpret the Ca + O{sub 2}(a) results, we utilized density functional theory along with multireference and explicitly correlated CCSD(T)-F12 electronic structure calculations to examine the lowest lying singlet and triplet surfaces. In addition to mapping stationary points, we used a genetic algorithm to locate minimum energy crossing points between the two surfaces. Simulations of the Ca + O{sub 2}(a) kinetics were then carried out using a combination of both standard and non-adiabatic Rice-Ramsperger-Kassel-Marcus (RRKM) theory implemented within a weak collision, multiwell master equation model. In terms of atmospheric significance, only in the case of Ca does reaction with O{sub 2}(a) compete with O{sub 3} during the daytime between 85 and 110 km.

  11. Mass Conservation and Chemical Kinetics.

    ERIC Educational Resources Information Center

    Barbara, Thomas M.; Corio, P. L.

    1980-01-01

    Presents a method for obtaining all mass conservation conditions implied by a given mechanism in which the conditions are used to simplify integration of the rate equations and to derive stoichiometric relations. Discusses possibilities of faulty inference of kinetic information from a given stoichiometry. (CS)

  12. Teaching Equations.

    ERIC Educational Resources Information Center

    Nibbelink, William H.

    1990-01-01

    Proposed is a gradual transition from arithmetic to the idea of an equation with variables in the elementary grades. Vertical and horizontal formats of open sentences, the instructional sequence, vocabulary, and levels of understanding are discussed in this article. (KR)

  13. Nucleation and growth transformation kinetics

    NASA Astrophysics Data System (ADS)

    Erukhimovitch, V.; Baram, J.

    1995-03-01

    As a result of the reassessment of the Kolmogorov-Johnson-Mehl-Avrami (KJMA) theory for the kinetics of nucleation and growth transformations, an integral-equation formulation has been developed instead of the well-known and widely used Avrami equation. The presented formulation considers interfacial and diffusional growths, in one, two, and three dimensions, with both time-dependent and time-invariant nucleation and growth rates. The integral-equation model corrects reported inadequacies of the KJMA theory when applied in numerous experiments and various solid-state transformations. It is shown that in the example cases examined in this paper, crystallization from the amorphous state in melt-spun ribbons, isothermal aging of CuAlZn, pearlitic transition in an eutectoid steel, and crystallization in a PEKK polymer, the thermodynamic and kinetic interpretation and parameters extracted from best fits of the Avrami equations to the experimental data are erroneous. The KJMA formulation is a simplification of the real physical conditions. The main limitation of the new model is that almost all the integral equations representing the kinetics of solid-state transformations have no analytical solutions.

  14. Chemical kinetics of geminal recombination

    SciTech Connect

    Levin, P.P.; Khudyakov, I.V.; Brin, E.F.; Kuz'min, V.A.

    1988-09-01

    The kinetics of geminal recombination of triplet radical pairs formed in photoreduction of benzophenone by p-cresol in glycerin solution was studied by pulsed laser photolysis. The experiments were conducted at several temperatures and in a constant magnetic field of H = 0.34 T. The parameters in six kinetic equations describing geminal recombination were determined with a computer. The values of the sums of the squares of the residual deviations of the approximation were obtained. It was found that the kinetics are best described by the functions proposed by Noyes and Shushin. It was shown that it is necessary to use the mutual diffusion coefficient of the radicals, which is significantly smaller than the sum of the estimations of the experimental values of the radical diffusion coefficients, for describing the kinetics due to the correlations of the molecular motions of the radicals in the cage.

  15. Population kinetics in dense plasmas

    SciTech Connect

    Schlanges, M.; Bornath, T.; Prenzel, R.; Kremp, D.

    1996-07-01

    Starting from quantum kinetic equations, rate equations for the number densities of the different atomic states and equations for the energy density are derived which are valid for dense nonideal plasmas. Statistical expressions are presented for the rate coefficients taking into account many-body effects as dynamical screening, lowering of the ionization energy and Pauli-blocking. Based on these generalized expressions, the coefficients of impact ionization, three-body recombination, excitation and deexcitation are calculated for nonideal hydrogen and carbon plasmas. As a result, higher ionization and recombination rates are obtained in the dense plasma region. The influence of the many-body effects on the population kinetics, including density and temperature relaxation, is shown then for a dense hydrogen plasma. {copyright} {ital 1996 American Institute of Physics.}

  16. Aerosol kinetic code "AERFORM": Model, validation and simulation results

    NASA Astrophysics Data System (ADS)

    Gainullin, K. G.; Golubev, A. I.; Petrov, A. M.; Piskunov, V. N.

    2016-06-01

    The aerosol kinetic code "AERFORM" is modified to simulate droplet and ice particle formation in mixed clouds. The splitting method is used to calculate condensation and coagulation simultaneously. The method is calibrated with analytic solutions of kinetic equations. Condensation kinetic model is based on cloud particle growth equation, mass and heat balance equations. The coagulation kinetic model includes Brownian, turbulent and precipitation effects. The real values are used for condensation and coagulation growth of water droplets and ice particles. The model and the simulation results for two full-scale cloud experiments are presented. The simulation model and code may be used autonomously or as an element of another code.

  17. A Gas-Kinetic Scheme for Reactive Flows

    NASA Technical Reports Server (NTRS)

    Lian,Youg-Sheng; Xu, Kun

    1998-01-01

    In this paper, the gas-kinetic BGK scheme for the compressible flow equations is extended to chemical reactive flow. The mass fraction of the unburnt gas is implemented into the gas kinetic equation by assigning a new internal degree of freedom to the particle distribution function. The new variable can be also used to describe fluid trajectory for the nonreactive flows. Due to the gas-kinetic BGK model, the current scheme basically solves the Navier-Stokes chemical reactive flow equations. Numerical tests validate the accuracy and robustness of the current kinetic method.

  18. Beautiful equations

    NASA Astrophysics Data System (ADS)

    Viljamaa, Panu; Jacobs, J. Richard; Chris; JamesHyman; Halma, Matthew; EricNolan; Coxon, Paul

    2014-07-01

    In reply to a Physics World infographic (part of which is given above) about a study showing that Euler's equation was deemed most beautiful by a group of mathematicians who had been hooked up to a functional magnetic-resonance image (fMRI) machine while viewing mathematical expressions (14 May, http://ow.ly/xHUFi).

  19. Production of a sterile species: Quantum kinetics

    SciTech Connect

    Ho, Chiu Man; Boyanovsky, D.; Ho, C.M.

    2007-04-23

    Production of a sterile species is studied within an effective model of active-sterile neutrino mixing in a medium in thermal equilibrium. The quantum kinetic equations for the distribution functions and coherences are obtained from two independent methods: the effective action and the quantum master equation. The decoherence time scale for active-sterile oscillations is tau(dec)=2/Gamma(aa), but the evolution of the distribution functions is determined by the two different time scales associated with the damping rates of the quasiparticle modes in the medium: Gamma(1)=Gamma(aa)cos^2theta(m); Gamma(2)=Gamma(aa)sin^2theta(m) where Gamma(aa) is the interaction rate of the active species in the absence of mixing and theta(m) the mixing angle in the medium. These two time scales are widely different away from Mikheyev-Smirnov-Wolfenstein resonances and preclude the kinetic description of active-sterile production in terms of a simple rate equation. We give the complete set of quantum kinetic equations for the active and sterile populations and coherences and discuss in detail the various approximations. A generalization of the active-sterile transition probability in a medium is provided via the quantum master equation. We derive explicitly the usual quantum kinetic equations in terms of the"polarization vector" and show their equivalence to those obtained from the quantum master equation and effective action.

  20. Marcus equation

    SciTech Connect

    1998-11-01

    In the late 1950s to early 1960s Rudolph A. Marcus developed a theory for treating the rates of outer-sphere electron-transfer reactions. Outer-sphere reactions are reactions in which an electron is transferred from a donor to an acceptor without any chemical bonds being made or broken. (Electron-transfer reactions in which bonds are made or broken are referred to as inner-sphere reactions.) Marcus derived several very useful expressions, one of which has come to be known as the Marcus cross-relation or, more simply, as the Marcus equation. It is widely used for correlating and predicting electron-transfer rates. For his contributions to the understanding of electron-transfer reactions, Marcus received the 1992 Nobel Prize in Chemistry. This paper discusses the development and use of the Marcus equation. Topics include self-exchange reactions; net electron-transfer reactions; Marcus cross-relation; and proton, hydride, atom and group transfers.

  1. Marcus equation

    DOE R&D Accomplishments Database

    1998-09-21

    In the late 1950s to early 1960s Rudolph A. Marcus developed a theory for treating the rates of outer-sphere electron-transfer reactions. Outer-sphere reactions are reactions in which an electron is transferred from a donor to an acceptor without any chemical bonds being made or broken. (Electron-transfer reactions in which bonds are made or broken are referred to as inner-sphere reactions.) Marcus derived several very useful expressions, one of which has come to be known as the Marcus cross-relation or, more simply, as the Marcus equation. It is widely used for correlating and predicting electron-transfer rates. For his contributions to the understanding of electron-transfer reactions, Marcus received the 1992 Nobel Prize in Chemistry. This paper discusses the development and use of the Marcus equation. Topics include self-exchange reactions; net electron-transfer reactions; Marcus cross-relation; and proton, hydride, atom and group transfers.

  2. Busulfan kinetics.

    PubMed

    Ehrsson, H; Hassan, M; Ehrnebo, M; Beran, M

    1983-07-01

    Busulfan kinetics were studied in patients with chronic myelocytic leukemia after oral doses of 2, 4, and 6 mg. The plasma concentration-time data could be fitted to a zero-order absorption one-compartment open model. The elimination rate constant averaged 0.27 +/- 0.05 hr-1 (SD). The plasma AUC was linearly related to the dose. The lag time for the start of absorption, the time absorption ends, and the absorption rate constant showed some interindividual variations. About 1% of busulfan is excreted unchanged in urine over 24 hr. PMID:6574831

  3. Group-kinetic theory of turbulence

    NASA Technical Reports Server (NTRS)

    Tchen, C. M.

    1986-01-01

    The two phases are governed by two coupled systems of Navier-Stokes equations. The couplings are nonlinear. These equations describe the microdynamical state of turbulence, and are transformed into a master equation. By scaling, a kinetic hierarchy is generated in the form of groups, representing the spectral evolution, the diffusivity and the relaxation. The loss of memory in formulating the relaxation yields the closure. The network of sub-distributions that participates in the relaxation is simulated by a self-consistent porous medium, so that the average effect on the diffusivity is to make it approach equilibrium. The kinetic equation of turbulence is derived. The method of moments reverts it to the continuum. The equation of spectral evolution is obtained and the transport properties are calculated. In inertia turbulence, the Kolmogoroff law for weak coupling and the spectrum for the strong coupling are found. As the fluid analog, the nonlinear Schrodinger equation has a driving force in the form of emission of solitons by velocity fluctuations, and is used to describe the microdynamical state of turbulence. In order for the emission together with the modulation to participate in the transport processes, the non-homogeneous Schrodinger equation is transformed into a homogeneous master equation. By group-scaling, the master equation is decomposed into a system of transport equations, replacing the Bogoliubov system of equations of many-particle distributions. It is in the relaxation that the memory is lost when the ensemble of higher-order distributions is simulated by an effective porous medium. The closure is thus found. The kinetic equation is derived and transformed into the equation of spectral flow.

  4. Tolrestat kinetics

    SciTech Connect

    Hicks, D.R.; Kraml, M.; Cayen, M.N.; Dubuc, J.; Ryder, S.; Dvornik, D.

    1984-10-01

    The kinetics of tolrestat, a potent inhibitor of aldose reductase, were examined. Serum concentrations of tolrestat and of total /sup 14/C were measured after dosing normal subjects and subjects with diabetes with /sup 14/C-labeled tolrestat. In normal subjects, tolrestat was rapidly absorbed and disappearance from serum was biphasic. Distribution and elimination t 1/2s were approximately 2 and 10 to 12 hr, respectively, after single and multiple doses. Unchanged tolrestat accounted for the major portion of /sup 14/C in serum. Radioactivity was rapidly and completely excreted in urine and feces in an approximate ratio of 2:1. Findings were much the same in subjects with diabetes. In normal subjects, the kinetics of oral tolrestat were independent of dose in the 10 to 800 mg range. Repetitive dosing did not result in unexpected cumulation. Tolrestat was more than 99% bound to serum protein; it did not compete with warfarin for binding sites but was displaced to some extent by high concentrations of tolbutamide or salicylate.

  5. On fast reactor kinetics studies

    SciTech Connect

    Seleznev, E. F.; Belov, A. A.; Matveenko, I. P.; Zhukov, A. M.; Raskach, K. F.

    2012-07-01

    The results and the program of fast reactor core time and space kinetics experiments performed and planned to be performed at the IPPE critical facility is presented. The TIMER code was taken as computation support of the experimental work, which allows transient equations to be solved in 3-D geometry with multi-group diffusion approximation. The number of delayed neutron groups varies from 6 to 8. The code implements the solution of both transient neutron transfer problems: a direct one, where neutron flux density and its derivatives, such as reactor power, etc, are determined at each time step, and an inverse one for the point kinetics equation form, where such a parameter as reactivity is determined with a well-known reactor power time variation function. (authors)

  6. Consistent lattice Boltzmann equations for phase transitions.

    PubMed

    Siebert, D N; Philippi, P C; Mattila, K K

    2014-11-01

    Unlike conventional computational fluid dynamics methods, the lattice Boltzmann method (LBM) describes the dynamic behavior of fluids in a mesoscopic scale based on discrete forms of kinetic equations. In this scale, complex macroscopic phenomena like the formation and collapse of interfaces can be naturally described as related to source terms incorporated into the kinetic equations. In this context, a novel athermal lattice Boltzmann scheme for the simulation of phase transition is proposed. The continuous kinetic model obtained from the Liouville equation using the mean-field interaction force approach is shown to be consistent with diffuse interface model using the Helmholtz free energy. Density profiles, interface thickness, and surface tension are analytically derived for a plane liquid-vapor interface. A discrete form of the kinetic equation is then obtained by applying the quadrature method based on prescribed abscissas together with a third-order scheme for the discretization of the streaming or advection term in the Boltzmann equation. Spatial derivatives in the source terms are approximated with high-order schemes. The numerical validation of the method is performed by measuring the speed of sound as well as by retrieving the coexistence curve and the interface density profiles. The appearance of spurious currents near the interface is investigated. The simulations are performed with the equations of state of Van der Waals, Redlich-Kwong, Redlich-Kwong-Soave, Peng-Robinson, and Carnahan-Starling. PMID:25493907

  7. A kinetic-theory approach to turbulent chemically reacting flows

    NASA Technical Reports Server (NTRS)

    Chung, P. M.

    1976-01-01

    The paper examines the mathematical and physical foundations for the kinetic theory of reactive turbulent flows, discussing the differences and relation between the kinetic and averaged equations, and comparing some solutions of the kinetic equations obtained by the Green's function method with those obtained by the approximate bimodal method. The kinetic method described consists essentially in constructing the probability density functions of the chemical species on the basis of solutions of the Langevin stochastic equation for the influence of eddies on the behavior of fluid elements. When the kinetic equations are solved for the structure of the diffusion flame established in a shear layer by the bimodal method, discontinuities in gradients of the mean concentrations at the two flame edges appear. This is a consequence of the bimodal approximation of all distribution functions by two dissimilar half-Maxwellian functions, which is a very crude approximation. These discontinuities do not appear when the solutions are constructed by the Green's function method described here.

  8. Dual Diagonalization of Reactive Transport Equations

    NASA Astrophysics Data System (ADS)

    Yeh, G.; Tsai, C.

    2013-12-01

    One solves a system of species transport equations in the primitive approach to reactive transport modeling. This approach is not able to decouple equilibrium reaction rates from species concentrations. This problem has been overcome with the approach to diagonalizing the reaction matrix since mid 1990's, which yields the same number of transport equations for reaction-extents. In the diagonalization approach, first, a subset of reaction- extent transport equations is solved for concentrations of components and kinetic-variables. Then, the component, kinetic-variable, and mass action equations are solved for all species concentrations. Finally, the equilibrium reaction rates are posterior computed. The difficulty in this approach is that the solution of species concentrations in the second step is a stiff problem when the concentrations of master species are small compared to those of equilibrium species. To overcome the problem of stiffness, we propose a dual diagonalization approach. Here, a second diagonalization is performed on the decomposed unit matrix to yield species concentrations, each as a linear function of reaction extents. In this dual diagonalization approach, four steps are needed to complete the modeling. First, component and kinetic-variable transport equations are solved for the concentrations of components (a subset of reaction-extents) and kinetic-variables (another subset of reaction-extents). Second, the set of mass action equations written in terms of reaction extents are solved for equilibrium-variables (yet another subset of reaction-extents). Third, species concentrations are posterior obtained by solving the set of linear equations defining reaction-extents. Fourth, equilibrium rates are posterior calculated with transport equations for equilibrium-variables. Several example problems will be used to demonstrate the efficiency of this approach. Keywords: Reactive Transport, Reaction-Extent, Component, Kinetic-Variable, Equilibrium

  9. On the relationships between Michaelis-Menten kinetics, reverse Michaelis-Menten kinetics, Equilibrium Chemistry Approximation kinetics and quadratic kinetics

    NASA Astrophysics Data System (ADS)

    Tang, J. Y.

    2015-09-01

    The Michaelis-Menten kinetics and the reverse Michaelis-Menten kinetics are two popular mathematical formulations used in many land biogeochemical models to describe how microbes and plants would respond to changes in substrate abundance. However, the criteria of when to use which of the two are often ambiguous. Here I show that these two kinetics are special approximations to the Equilibrium Chemistry Approximation kinetics, which is the first order approximation to the quadratic kinetics that solves the equation of enzyme-substrate complex exactly for a single enzyme single substrate biogeochemical reaction with the law of mass action and the assumption of quasi-steady-state for the enzyme-substrate complex and that the product genesis from enzyme-substrate complex is much slower than the equilibration between enzyme-substrate complexes, substrates and enzymes. In particular, I showed that the derivation of the Michaelis-Menten kinetics does not consider the mass balance constraint of the substrate, and the reverse Michaelis-Menten kinetics does not consider the mass balance constraint of the enzyme, whereas both of these constraints are taken into account in the Equilibrium Chemistry Approximation kinetics. By benchmarking against predictions from the quadratic kinetics for a wide range of substrate and enzyme concentrations, the Michaelis-Menten kinetics was found to persistently under-predict the normalized sensitivity ∂ ln v / ∂ ln k2+ of the reaction velocity v with respect to the maximum product genesis rate k2+, persistently over-predict the normalized sensitivity ∂ ln v / ∂ ln k1+ of v with respect to the intrinsic substrate affinity k1+, persistently over-predict the normalized sensitivity ∂ ln v / ∂ ln [ E ]T of v with respect the total enzyme concentration [ E ]T and persistently under-predict the normalized sensitivity ∂ ln v / ∂ ln [ S ]T of v with respect to the total substrate concentration [ S ]T. Meanwhile

  10. On the relationships between Michaelis–Menten kinetics, reverse Michaelis–Menten kinetics, Equilibrium Chemistry Approximation kinetics and quadratic kinetics

    DOE PAGESBeta

    Tang, J. Y.

    2015-09-03

    The Michaelis–Menten kinetics and the reverse Michaelis–Menten kinetics are two popular mathematical formulations used in many land biogeochemical models to describe how microbes and plants would respond to changes in substrate abundance. However, the criteria of when to use which of the two are often ambiguous. Here I show that these two kinetics are special approximations to the Equilibrium Chemistry Approximation kinetics, which is the first order approximation to the quadratic kinetics that solves the equation of enzyme-substrate complex exactly for a single enzyme single substrate biogeochemical reaction with the law of mass action and the assumption of quasi-steady-state formore » the enzyme-substrate complex and that the product genesis from enzyme-substrate complex is much slower than the equilibration between enzyme-substrate complexes, substrates and enzymes. In particular, I showed that the derivation of the Michaelis–Menten kinetics does not consider the mass balance constraint of the substrate, and the reverse Michaelis–Menten kinetics does not consider the mass balance constraint of the enzyme, whereas both of these constraints are taken into account in the Equilibrium Chemistry Approximation kinetics. By benchmarking against predictions from the quadratic kinetics for a wide range of substrate and enzyme concentrations, the Michaelis–Menten kinetics was found to persistently under-predict the normalized sensitivity ∂ ln v / ∂ ln k2+ of the reaction velocity v with respect to the maximum product genesis rate k2+, persistently over-predict the normalized sensitivity ∂ ln v / ∂ ln k1+ of v with respect to the intrinsic substrate affinity k1+, persistently over-predict the normalized sensitivity ∂ ln v / ∂ ln [ E ]T of v with respect the total enzyme concentration [ E ]T and persistently under-predict the normalized sensitivity ∂ ln v / ∂ ln [ S ]T of v with respect to the total substrate concentration [ S ]T. Meanwhile, the

  11. Algebraic operator approach to gas kinetic models

    NASA Astrophysics Data System (ADS)

    Il'ichov, L. V.

    1997-02-01

    Some general properties of the linear Boltzmann kinetic equation are used to present it in the form ∂ tϕ = - †Âϕ with the operators Âand† possessing some nontrivial algebraic properties. When applied to the Keilson-Storer kinetic model, this method gives an example of quantum ( q-deformed) Lie algebra. This approach provides also a natural generalization of the “kangaroo model”.

  12. Perceptions of the Schrodinger equation

    NASA Astrophysics Data System (ADS)

    Efthimiades, Spyros

    2014-03-01

    The Schrodinger equation has been considered to be a postulate of quantum physics, but it is also perceived as the quantum equivalent of the non-relativistic classical energy relation. We argue that the Schrodinger equation cannot be a physical postulate, and we show explicitly that its second space derivative term is wrongly associated with the kinetic energy of the particle. The kinetic energy of a particle at a point is proportional to the square of the momentum, that is, to the square of the first space derivative of the wavefunction. Analyzing particle interactions, we realize that particles have multiple virtual motions and that each motion is accompanied by a wave that has constant amplitude. Accordingly, we define the wavefunction as the superposition of the virtual waves of the particle. In simple interaction settings we can tell what particle motions arise and can explain the outcomes in direct and tangible terms. Most importantly, the mathematical foundation of quantum mechanics becomes clear and justified, and we derive the Schrodinger, Dirac, etc. equations as the conditions the wavefunction must satisfy at each space-time point in order to fulfill the respective total energy equation.

  13. Hybrid fluid/kinetic model for parallel heat conduction

    SciTech Connect

    Callen, J.D.; Hegna, C.C.; Held, E.D.

    1998-12-31

    It is argued that in order to use fluid-like equations to model low frequency ({omega} < {nu}) phenomena such as neoclassical tearing modes in low collisionality ({nu} < {omega}{sub b}) tokamak plasmas, a Chapman-Enskog-like approach is most appropriate for developing an equation for the kinetic distortion (F) of the distribution function whose velocity-space moments lead to the needed fluid moment closure relations. Further, parallel heat conduction in a long collision mean free path regime can be described through a combination of a reduced phase space Chapman-Enskog-like approach for the kinetics and a multiple-time-scale analysis for the fluid and kinetic equations.

  14. Eulerian action principles for linearized reduced dynamical equations

    NASA Astrophysics Data System (ADS)

    Brizard, Alain

    1994-08-01

    New Eulerian action principles for the linearized gyrokinetic Maxwell-Vlasov equations and the linearized kinetic-magnetohydrodynamic (kinetic-MHD) equations are presented. The variational fields for the linearized gyrokinetic Vlasov-Maxwell equations are the perturbed electromagnetic potentials (φ1,A1) and the gyroangle-independent gyrocenter (gy) function Sgy, while the variational fields for the linearized kinetic-MHD equations are the ideal MHD fluid displacement ξ and the gyroangle-independent drift-kinetic (dk) function Sdk (defined as the drift-kinetic limit of Sgy). According to the Lie-transform approach to Vlasov perturbation theory, Sgy generates first-order perturbations in the gyrocenter distribution F1≡{Sgy, F0}gc, where F1 satisfies the linearized gyrokinetic Vlasov equation and {, }gc denotes the unperturbed guiding-center (gc) Poisson bracket. Previous quadratic variational forms were constructed ad hoc from the linearized equations, and required the linearized gyrokinetic (or drift-kinetic) Vlasov equation to be solved a priori (e.g., by integration along an unperturbed guiding-center orbit) through the use of the normal-mode and ballooning-mode representations. The presented action principles ignore these requirements and, thus, apply to more general perturbations.

  15. Kinetic theory of spatially inhomogeneous stellar systems without collective effects

    NASA Astrophysics Data System (ADS)

    Chavanis, P.-H.

    2013-08-01

    We review and complete the kinetic theory of spatially inhomogeneous stellar systems when collective effects (dressing of the stars by their polarization cloud) are neglected. We start from the BBGKY hierarchy issued from the Liouville equation and consider an expansion in powers of 1/N in a proper thermodynamic limit. For N → +∞, we obtain the Vlasov equation describing the evolution of collisionless stellar systems like elliptical galaxies. This corresponds to the mean field approximation. At the order 1/N, we obtain a kinetic equation describing the evolution of collisional stellar systems like globular clusters. This corresponds to the weak coupling approximation. This equation coincides with the generalized Landau equation derived from a more abstract projection operator formalism. This equation does not suffer logarithmic divergences at large scales since spatial inhomogeneity is explicitly taken into account. Making a local approximation, and introducing an upper cut-off at the Jeans length, it reduces to the Vlasov-Landau equation which is the standard kinetic equation of stellar systems. Our approach provides a simple and pedagogical derivation of these important equations from the BBGKY hierarchy which is more rigorous for systems with long-range interactions than the two-body encounters theory. Making an adiabatic approximation, we write the generalized Landau equation in angle-action variables and obtain a Landau-type kinetic equation that is valid for fully inhomogeneous stellar systems and is free of divergences at large scales. This equation is less general than the recently derived Lenard-Balescu-type kinetic equation since it neglects collective effects, but it is substantially simpler and could be useful as a first step. We discuss the evolution of the system as a whole and the relaxation of a test star in a bath of field stars. We derive the corresponding Fokker-Planck equation in angle-action variables and provide expressions for the

  16. Extended rate equations

    SciTech Connect

    Shore, B.W.

    1981-01-30

    The equations of motion are discussed which describe time dependent population flows in an N-level system, reviewing the relationship between incoherent (rate) equations, coherent (Schrodinger) equations, and more general partially coherent (Bloch) equations. Approximations are discussed which replace the elaborate Bloch equations by simpler rate equations whose coefficients incorporate long-time consequences of coherence.

  17. Zakharov equations in quantum dusty plasmas

    SciTech Connect

    Sayed, F.; Vladimirov, S. V.; Ishihara, O.

    2015-08-15

    By generalizing the formalism of modulational interactions in quantum dusty plasmas, we derive the kinetic quantum Zakharov equations in dusty plasmas that describe nonlinear coupling of high frequency Langmuir waves to low frequency plasma density variations, for cases of non-degenerate and degenerate plasma electrons.

  18. Reciprocal relations between kinetic curves

    NASA Astrophysics Data System (ADS)

    Yablonsky, G. S.; Gorban, A. N.; Constales, D.; Galvita, V. V.; Marin, G. B.

    2011-01-01

    We study coupled irreversible processes. For linear or linearized kinetics with microreversibility, \\dot{x}=Kx , the kinetic operator K is symmetric in the entropic inner product. This form of Onsager's reciprocal relations implies that the shift in time, exp(Kt), is also a symmetric operator. This generates the reciprocity relations between the kinetic curves. For example, for the Master equation, if we start the process from the i-th pure state and measure the probability pj(t) of the j-th state (j≠i), and, similarly, measure pi(t) for the process, which starts at the j-th pure state, then the ratio of these two probabilities pj(t)/pi(t) is constant in time and coincides with the ratio of the equilibrium probabilities. We study similar and more general reciprocal relations between the kinetic curves. The experimental evidence provided as an example is from the reversible water gas shift reaction over iron oxide catalyst. The experimental data are obtained using Temporal Analysis of Products (TAP) pulse-response studies. These offer excellent confirmation within the experimental error.

  19. Kinetics of the humid aging of magnetic recording tape

    NASA Technical Reports Server (NTRS)

    Bertram, H. N.; Cuddihy, E. F.

    1982-01-01

    The kinetics of the hydrolysis of polyester urethane binders of magnetic recording tape is described. Kinetic data were generated from measurements of acetone-extractable hydrolyzed binder products versus time for various humidity-temperature environments. These data can be described by a linear, single product, reversible rate equation. This equation, coupled with measurements on the effect of hydrolysis on recorded tape performance, is used to predict proper environmental storage conditions for magnetic tape.

  20. Kinetic theory of electromagnetic ion waves in relativistic plasmas

    SciTech Connect

    Marklund, Mattias; Shukla, Padma K.

    2006-09-15

    A kinetic theory for electromagnetic ion waves in a cold relativistic plasma is derived. The kinetic equation for the broadband electromagnetic ion waves is coupled to the slow density response via an acoustic equation driven by a ponderomotive force-like term linear in the electromagnetic field amplitude. The modulational instability growth rate is derived for an arbitrary spectrum of waves. The monochromatic and random phase cases are studied.

  1. Simplification of the unified gas kinetic scheme.

    PubMed

    Chen, Songze; Guo, Zhaoli; Xu, Kun

    2016-08-01

    The unified gas kinetic scheme (UGKS) is an asymptotic preserving (AP) scheme for kinetic equations. It is superior for transition flow simulation and has been validated in the past years. However, compared to the well-known discrete ordinate method (DOM), which is a classical numerical method solving the kinetic equations, the UGKS needs more computational resources. In this study, we propose a simplification of the unified gas kinetic scheme. It allows almost identical numerical cost as the DOM, but predicts numerical results as accurate as the UGKS. In the simplified scheme, the numerical flux for the velocity distribution function and the numerical flux for the macroscopic conservative quantities are evaluated separately. The equilibrium part of the UGKS flux is calculated by analytical solution instead of the numerical quadrature in velocity space. The simplification is equivalent to a flux hybridization of the gas kinetic scheme for the Navier-Stokes (NS) equations and the conventional discrete ordinate method. Several simplification strategies are tested, through which we can identify the key ingredient of the Navier-Stokes asymptotic preserving property. Numerical tests show that, as long as the collision effect is built into the macroscopic numerical flux, the numerical scheme is Navier-Stokes asymptotic preserving, regardless the accuracy of the microscopic numerical flux for the velocity distribution function. PMID:27627418

  2. Kinetic theory of relativistic plasmas

    NASA Technical Reports Server (NTRS)

    Gould, R. J.

    1981-01-01

    The thermalization of particle kinetic motion by binary collisions is considered for a plasma with a Boltzmann constant-temperature product approximately equal to 10 to 100 times the product of the electron mass with the square of the speed of light. At this temperature, the principal mechanism for relaxation of electron motion is via radiationless electron-electron collisions (Moller scattering). Ions are nonrelativistic, but are energetic enough so that their Coulomb scattering can be treated in the Born approximation. Relaxation times are computed and Boltzmann-equation Fokker-Planck operators are derived for the various binary-collision processes. The expression for the rate of kinetic energy exchange between electron and ion gases is derived for the case where the gases are at different temperatures.

  3. Fulvenallene decomposition kinetics.

    PubMed

    Polino, Daniela; Cavallotti, Carlo

    2011-09-22

    While the decomposition kinetics of the benzyl radical has been studied in depth both from the experimental and the theoretical standpoint, much less is known about the reactivity of what is likely to be its main decomposition product, fulvenallene. In this work the high temperature reactivity of fulvenallene was investigated on a Potential Energy Surface (PES) consisting of 10 wells interconnected through 11 transition states using a 1 D Master Equation (ME). Rate constants were calculated using RRKM theory and the ME was integrated using a stochastic kinetic Monte Carlo code. It was found that two main decomposition channels are possible, the first is active on the singlet PES and leads to the formation of the fulvenallenyl radical and atomic hydrogen. The second requires intersystem crossing to the triplet PES and leads to acetylene and cyclopentadienylidene. ME simulations were performed calculating the microcanonical intersystem crossing frequency using Landau-Zener theory convoluting the crossing probability with RRKM rates evaluated at the conical intersection. It was found that the reaction channel leading to the cyclopentadienylidene diradical is only slightly faster than that leading to the fulvenallenyl radical, so that it can be concluded that both reactions are likely to be active in the investigated temperature (1500-2000 K) and pressure (0.05-50 bar) ranges. However, the simulations show that intersystem crossing is rate limiting for the first reaction channel, as the removal of this barrier leads to an increase of the rate constant by a factor of 2-3. Channel specific rate constants are reported as a function of temperature and pressure. PMID:21819060

  4. Lattice kinetic simulation of nonisothermal magnetohydrodynamics.

    PubMed

    Chatterjee, Dipankar; Amiroudine, Sakir

    2010-06-01

    In this paper, a lattice kinetic algorithm is presented to simulate nonisothermal magnetohydrodynamics in the low-Mach number incompressible limit. The flow and thermal fields are described by two separate distribution functions through respective scalar kinetic equations and the magnetic field is governed by a vector distribution function through a vector kinetic equation. The distribution functions are only coupled via the macroscopic density, momentum, magnetic field, and temperature computed at the lattice points. The novelty of the work is the computation of the thermal field in conjunction with the hydromagnetic fields in the lattice Boltzmann framework. A 9-bit two-dimensional (2D) lattice scheme is used for the numerical computation of the hydrodynamic and thermal fields, whereas the magnetic field is simulated in a 5-bit 2D lattice. Simulation of Hartmann flow in a channel provides excellent agreement with corresponding analytical results. PMID:20866540

  5. Classical non-Markovian Boltzmann equation

    SciTech Connect

    Alexanian, Moorad

    2014-08-01

    The modeling of particle transport involves anomalous diffusion, (x²(t) ) ∝ t{sup α} with α ≠ 1, with subdiffusive transport corresponding to 0 < α < 1 and superdiffusive transport to α > 1. These anomalies give rise to fractional advection-dispersion equations with memory in space and time. The usual Boltzmann equation, with only isolated binary collisions, is Markovian and, in particular, the contributions of the three-particle distribution function are neglected. We show that the inclusion of higher-order distribution functions give rise to an exact, non-Markovian Boltzmann equation with resulting transport equations for mass, momentum, and kinetic energy with memory in both time and space. The two- and the three-particle distribution functions are considered under the assumption that the two- and the three-particle correlation functions are translationally invariant that allows us to obtain advection-dispersion equations for modeling transport in terms of spatial and temporal fractional derivatives.

  6. Quantum logics and chemical kinetics

    NASA Astrophysics Data System (ADS)

    Ivanov, C. I.

    1981-06-01

    A statistical theory of chemical kinetics is presented based on the quantum logical concept of chemical observables. The apparatus of Boolean algebra B is applied for the construction of appropriate composition polynomials referring to any stipulated arrangement of the atomic constituents. A physically motivated probability measure μ( F) is introduced on the field B of chemical observables, which considers the occurrence of the yes response of a given F ɛ B. The equations for the time evolution of the species density operators and the master equations for the corresponding number densities are derived. The general treatment is applied to a superposition of elementary substitution reactions (AB) α + C ⇄ (AC) β + B. The expressions for the reaction rate coefficients are established.

  7. Solid State Kinetic Parameters and Chemical Mechanism of the Dehydration of CoCl2.6H2O.

    ERIC Educational Resources Information Center

    Ribas, Joan; And Others

    1988-01-01

    Presents an experimental example illustrating the most common methods for the determination of kinetic parameters. Discusses the different theories and equations to be applied and the mechanism derived from the kinetic results. (CW)

  8. Representing Rate Equations for Enzyme-Catalyzed Reactions

    ERIC Educational Resources Information Center

    Ault, Addison

    2011-01-01

    Rate equations for enzyme-catalyzed reactions are derived and presented in a way that makes it easier for the nonspecialist to see how the rate of an enzyme-catalyzed reaction depends upon kinetic constants and concentrations. This is done with distribution equations that show how the rate of the reaction depends upon the relative quantities of…

  9. Transport equations in tokamak plasmas

    SciTech Connect

    Callen, J. D.; Hegna, C. C.; Cole, A. J.

    2010-05-15

    Tokamak plasma transport equations are usually obtained by flux surface averaging the collisional Braginskii equations. However, tokamak plasmas are not in collisional regimes. Also, ad hoc terms are added for neoclassical effects on the parallel Ohm's law, fluctuation-induced transport, heating, current-drive and flow sources and sinks, small magnetic field nonaxisymmetries, magnetic field transients, etc. A set of self-consistent second order in gyroradius fluid-moment-based transport equations for nearly axisymmetric tokamak plasmas has been developed using a kinetic-based approach. The derivation uses neoclassical-based parallel viscous force closures, and includes all the effects noted above. Plasma processes on successive time scales and constraints they impose are considered sequentially: compressional Alfven waves (Grad-Shafranov equilibrium, ion radial force balance), sound waves (pressure constant along field lines, incompressible flows within a flux surface), and collisions (electrons, parallel Ohm's law; ions, damping of poloidal flow). Radial particle fluxes are driven by the many second order in gyroradius toroidal angular torques on a plasma species: seven ambipolar collision-based ones (classical, neoclassical, etc.) and eight nonambipolar ones (fluctuation-induced, polarization flows from toroidal rotation transients, etc.). The plasma toroidal rotation equation results from setting to zero the net radial current induced by the nonambipolar fluxes. The radial particle flux consists of the collision-based intrinsically ambipolar fluxes plus the nonambipolar fluxes evaluated at the ambipolarity-enforcing toroidal plasma rotation (radial electric field). The energy transport equations do not involve an ambipolar constraint and hence are more directly obtained. The 'mean field' effects of microturbulence on the parallel Ohm's law, poloidal ion flow, particle fluxes, and toroidal momentum and energy transport are all included self-consistently. The

  10. Spin field equations and Heun's equations

    NASA Astrophysics Data System (ADS)

    Jiang, Min; Wang, Xuejing; Li, Zhongheng

    2015-06-01

    The Kerr-Newman-(anti) de Sitter metric is the most general stationary black hole solution to the Einstein-Maxwell equation with a cosmological constant. We study the separability of the equations of the massless scalar (spin s=0), neutrino ( s=1/2), electromagnetic ( s=1), Rarita-Schwinger ( s=3/2), and gravitational ( s=2) fields propagating on this background. We obtain the angular and radial master equations, and show that the master equations are transformed to Heun's equation. Meanwhile, we give the condition of existence of event horizons for Kerr-Newman-(anti) de Sitter spacetime by using Sturm theorem.