NASA Astrophysics Data System (ADS)
Blanco, R.; Pesquera, L.; Santos, E.
1983-03-01
The motion of a charged pointlike relativistic particle under the action of a given force field plus a random electromagnetic radiation is studied. It is assumed that the given force field alone should produce a multiply periodic motion, which is perturbed by the action of both the random radiation and the reaction damping. The random radiation is represented by a stochastic process and an equation is obtained for the equilibrium probability density of the particle in phase space. In the particular case of a random radiation with Rayleigh-Jeans spectrum, it is shown that the stationary solution, corresponding to radiation-matter equilibrium, is given by the Maxwell-Boltzmann distribution.
Boltzmann equation solver adapted to emergent chemical non-equilibrium
Birrell, Jeremiah; Wilkening, Jon; Rafelski, Johann
2015-01-15
We present a novel method to solve the spatially homogeneous and isotropic relativistic Boltzmann equation. We employ a basis set of orthogonal polynomials dynamically adapted to allow for emergence of chemical non-equilibrium. Two time dependent parameters characterize the set of orthogonal polynomials, the effective temperature T(t) and phase space occupation factor ϒ(t). In this first paper we address (effectively) massless fermions and derive dynamical equations for T(t) and ϒ(t) such that the zeroth order term of the basis alone captures the particle number density and energy density of each particle distribution. We validate our method and illustrate the reduced computational cost and the ability to easily represent final state chemical non-equilibrium by studying a model problem that is motivated by the physics of the neutrino freeze-out processes in the early Universe, where the essential physical characteristics include reheating from another disappearing particle component (e{sup ±}-annihilation)
Revised lattice Boltzmann model for traffic flow with equilibrium traffic pressure
NASA Astrophysics Data System (ADS)
Shi, Wei; Lu, Wei-Zhen; Xue, Yu; He, Hong-Di
2016-02-01
A revised lattice Boltzmann model concerning the equilibrium traffic pressure is proposed in this study to tackle the phase transition phenomena of traffic flow system. The traditional lattice Boltzmann model has limitation to investigate the complex traffic phase transitions due to its difficulty for modeling the equilibrium velocity distribution. Concerning this drawback, the equilibrium traffic pressure is taken into account to derive the equilibrium velocity distribution in the revised lattice Boltzmann model. In the proposed model, a three-dimensional velocity-space is assumed to determine the equilibrium velocity distribution functions and an alternative, new derivative approach is introduced to deduct the macroscopic equations with the first-order accuracy level from the lattice Boltzmann model. Based on the linear stability theory, the stability conditions of the corresponding macroscopic equations can be obtained. The outputs indicate that the stability curve is divided into three regions, i.e., the stable region, the neutral stability region, and the unstable region. In the stable region, small disturbance appears in the initial uniform flow and will vanish after long term evolution, while in the unstable region, the disturbance will be enlarged and finally leads to the traffic system entering the congested state. In the neutral stability region, small disturbance does not vanish with time and maintains its amplitude in the traffic system. Conclusively, the stability of traffic system is found to be enhanced as the equilibrium traffic pressure increases. Finally, the numerical outputs of the proposed model are found to be consistent with the recognized, theoretical results.
Operational derivation of Boltzmann distribution with Maxwell's demon model.
Hosoya, Akio; Maruyama, Koji; Shikano, Yutaka
2015-01-01
The resolution of the Maxwell's demon paradox linked thermodynamics with information theory through information erasure principle. By considering a demon endowed with a Turing-machine consisting of a memory tape and a processor, we attempt to explore the link towards the foundations of statistical mechanics and to derive results therein in an operational manner. Here, we present a derivation of the Boltzmann distribution in equilibrium as an example, without hypothesizing the principle of maximum entropy. Further, since the model can be applied to non-equilibrium processes, in principle, we demonstrate the dissipation-fluctuation relation to show the possibility in this direction. PMID:26598363
Exponential trend to equilibrium for the inelastic Boltzmann equation driven by a particle bath
NASA Astrophysics Data System (ADS)
Cañizo, José A.; Lods, Bertrand
2016-05-01
We consider the spatially homogeneous Boltzmann equation for inelastic hard spheres (with constant restitution coefficient α \\in (0,1) ) under the thermalization induced by a host medium with a fixed Maxwellian distribution. We prove that the solution to the associated initial-value problem converges exponentially fast towards the unique equilibrium solution. The proof combines a careful spectral analysis of the linearised semigroup as well as entropy estimates. The trend towards equilibrium holds in the weakly inelastic regime in which α is close to 1, and the rate of convergence is explicit and depends solely on the spectral gap of the elastic linear collision operator.
NASA Astrophysics Data System (ADS)
Asinari, P.
2011-03-01
Boltzmann equation is one the most powerful paradigms for explaining transport phenomena in fluids. Since early fifties, it received a lot of attention due to aerodynamic requirements for high altitude vehicles, vacuum technology requirements and nowadays, micro-electro-mechanical systems (MEMs). Because of the intrinsic mathematical complexity of the problem, Boltzmann himself started his work by considering first the case when the distribution function does not depend on space (homogeneous case), but only on time and the magnitude of the molecular velocity (isotropic collisional integral). The interest with regards to the homogeneous isotropic Boltzmann equation goes beyond simple dilute gases. In the so-called econophysics, a Boltzmann type model is sometimes introduced for studying the distribution of wealth in a simple market. Another recent application of the homogeneous isotropic Boltzmann equation is given by opinion formation modeling in quantitative sociology, also called socio-dynamics or sociophysics. The present work [1] aims to improve the deterministic method for solving homogenous isotropic Boltzmann equation proposed by Aristov [2] by two ideas: (a) the homogeneous isotropic problem is reformulated first in terms of particle kinetic energy (this allows one to ensure exact particle number and energy conservation during microscopic collisions) and (b) a DVM-like correction (where DVM stands for Discrete Velocity Model) is adopted for improving the relaxation rates (this allows one to satisfy exactly the conservation laws at macroscopic level, which is particularly important for describing the late dynamics in the relaxation towards the equilibrium).
An alternative approach to the Boltzmann distribution through the chemical potential
NASA Astrophysics Data System (ADS)
D'Anna, Michele; Job, Georg
2016-05-01
The Boltzmann distribution is one of the most significant results of classical physics. Despite its importance and its wide range of application, at high school level it is mostly presented without any derivation or link to some basic ideas. In this contribution we present an approach based on the chemical potential that allows to derive it directly from the basic idea of thermodynamical equilibrium.
Permit allocation in emissions trading using the Boltzmann distribution
NASA Astrophysics Data System (ADS)
Park, Ji-Won; Kim, Chae Un; Isard, Walter
2012-10-01
In emissions trading, the initial allocation of permits is an intractable issue because it needs to be essentially fair to the participating countries. There are many ways to distribute a given total amount of emissions permits among countries, but the existing distribution methods, such as auctioning and grandfathering, have been debated. In this paper we describe a new method for allocating permits in emissions trading using the Boltzmann distribution. We introduce the Boltzmann distribution to permit allocation by combining it with concepts in emissions trading. We then demonstrate through empirical data analysis how emissions permits can be allocated in practice among participating countries. The new allocation method using the Boltzmann distribution describes the most probable, natural, and unbiased distribution of emissions permits among multiple countries. Simple and versatile, this new method holds potential for many economic and environmental applications.
Liñán-García, Ernesto; Sánchez-Hernández, Juan Paulo; González-Barbosa, J. Javier; González-Flores, Carlos
2016-01-01
A new hybrid Multiphase Simulated Annealing Algorithm using Boltzmann and Bose-Einstein distributions (MPSABBE) is proposed. MPSABBE was designed for solving the Protein Folding Problem (PFP) instances. This new approach has four phases: (i) Multiquenching Phase (MQP), (ii) Boltzmann Annealing Phase (BAP), (iii) Bose-Einstein Annealing Phase (BEAP), and (iv) Dynamical Equilibrium Phase (DEP). BAP and BEAP are simulated annealing searching procedures based on Boltzmann and Bose-Einstein distributions, respectively. DEP is also a simulated annealing search procedure, which is applied at the final temperature of the fourth phase, which can be seen as a second Bose-Einstein phase. MQP is a search process that ranges from extremely high to high temperatures, applying a very fast cooling process, and is not very restrictive to accept new solutions. However, BAP and BEAP range from high to low and from low to very low temperatures, respectively. They are more restrictive for accepting new solutions. DEP uses a particular heuristic to detect the stochastic equilibrium by applying a least squares method during its execution. MPSABBE parameters are tuned with an analytical method, which considers the maximal and minimal deterioration of problem instances. MPSABBE was tested with several instances of PFP, showing that the use of both distributions is better than using only the Boltzmann distribution on the classical SA. PMID:27413369
Lin, X.
1991-01-01
This paper reports the development of an object-oriented programming methodology for particle simulations. It is established on the (m reductionist) view that many physical phenomena cana be reduced to many-body problems. By doing the reduction, many seemly unrelated physical phenomena can be simulated in a systematic way and a high-level programming system can be constructed to facilitate the programming and the solution of the simulations. In the object-oriented particle simulation methodology, a hierarchy of abstract particles is defined to represent a variety of characteristics in physical system simulations. A simulation program is constructed from particles derived from the abstract particles. The object- oriented particle simulation methodology provides a unifying modeling and simulation framework for a variety of simulation applications with the use of particle methods. It allows easy composition of simulation programs from predefined software modules and facilitates software reusability. It greatly increase the productivity of simulation program constructions. Boltzmann (after Ludwig Boltzmann, 1844-1906) is a prototype programming system in the object-oriented particle simulation methodology. Boltzmann is implemented in C++ and the X Window System. It contains a library of data types and functions that support simulations in particle methods. Moreover, it provides a visualization window to support friendly user-computer interaction. Examples of the application of the Boltzmann programming system are presented. The effectiveness of the object-oriented particle simulation methodology is demonstrated. A user's manual is included in the appendix.
Operational derivation of Boltzmann distribution with Maxwell’s demon model
NASA Astrophysics Data System (ADS)
Hosoya, Akio; Maruyama, Koji; Shikano, Yutaka
2015-11-01
The resolution of the Maxwell’s demon paradox linked thermodynamics with information theory through information erasure principle. By considering a demon endowed with a Turing-machine consisting of a memory tape and a processor, we attempt to explore the link towards the foundations of statistical mechanics and to derive results therein in an operational manner. Here, we present a derivation of the Boltzmann distribution in equilibrium as an example, without hypothesizing the principle of maximum entropy. Further, since the model can be applied to non-equilibrium processes, in principle, we demonstrate the dissipation-fluctuation relation to show the possibility in this direction.
Operational derivation of Boltzmann distribution with Maxwell’s demon model
Hosoya, Akio; Maruyama, Koji; Shikano, Yutaka
2015-01-01
The resolution of the Maxwell’s demon paradox linked thermodynamics with information theory through information erasure principle. By considering a demon endowed with a Turing-machine consisting of a memory tape and a processor, we attempt to explore the link towards the foundations of statistical mechanics and to derive results therein in an operational manner. Here, we present a derivation of the Boltzmann distribution in equilibrium as an example, without hypothesizing the principle of maximum entropy. Further, since the model can be applied to non-equilibrium processes, in principle, we demonstrate the dissipation-fluctuation relation to show the possibility in this direction. PMID:26598363
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.
The Boltzmann Equation for a Multi-species Mixture Close to Global Equilibrium
NASA Astrophysics Data System (ADS)
Briant, Marc; Daus, Esther S.
2016-07-01
We study the Cauchy theory for a multi-species mixture, where the different species can have different masses, in a perturbative setting on the three dimensional torus. The ultimate aim of this work is to obtain the existence, uniqueness and exponential trend to equilibrium of solutions to the multi-species Boltzmann equation in {L^1_vL^∞_x(m)} , where {m˜ (1+ |v|^k)} is a polynomial weight. We prove the existence of a spectral gap for the linear multi-species Boltzmann operator allowing different masses, and then we establish a semigroup property thanks to a new explicit coercive estimate for the Boltzmann operator. Then we develop an {L^2-L^E&infty}; theory à la Guo for the linear perturbed equation. Finally, we combine the latter results with a decomposition of the multi-species Boltzmann equation in order to deal with the full equation. We emphasize that dealing with different masses induces a loss of symmetry in the Boltzmann operator which prevents the direct adaptation of standard mono-species methods (for example Carleman representation, Povzner inequality). Of important note is the fact that all methods used and developed in this work are constructive. Moreover, they do not require any Sobolev regularity and the {L^1_vL^∞_x} framework is dealt with for any {k > k_0} , recovering the optimal physical threshold of finite energy {k_0=2} in the particular case of a multi-species hard spheres mixture with the same masses.
Beyond Boltzmann-Gibbs statistics: Maximum entropy hyperensemblesout-of-equilibrium
Crooks, Gavin E.
2006-02-23
What is the best description that we can construct of athermodynamic system that is not in equilibrium, given only one, or afew, extra parameters over and above those needed for a description ofthe same system at equilibrium? Here, we argue the most appropriateadditional parameter is the non-equilibrium entropy of the system, andthat we should not attempt to estimate the probability distribution ofthe system, but rather the metaprobability (or hyperensemble) that thesystem is described by a particular probability distribution. The resultis an entropic distribution with two parameters, one a non-equilibriumtemperature, and the other a measure of distance from equilibrium. Thisdispersion parameter smoothly interpolates between certainty of acanonical distribution at equilibrium and great uncertainty as to theprobability distribution as we move away from equilibrium. We deducethat, in general, large, rare fluctuations become far more common as wemove away from equilibrium.
Boltzmann Gibbs distribution of fortune and broken time reversible symmetry in econodynamics
NASA Astrophysics Data System (ADS)
Ao, P.
2007-08-01
Within the framework of stochastic differential equations it is demonstrated that the existence of Boltzmann-Gibbs type distribution in economy is independent of the time reversal symmetry in econodynamics. Both power law and exponential distributions can be accommodated naturally. The demonstration is based on a mathematical structure discovered during a study in gene regulatory network dynamics. Further possible analogy between equilibrium economy and thermodynamics is explored, suggesting that statistical physics methods can indeed play an important role in the study of complex systems.
Polyelectrolyte Microcapsules: Ion Distributions from a Poisson-Boltzmann Model
NASA Astrophysics Data System (ADS)
Tang, Qiyun; Denton, Alan R.; Rozairo, Damith; Croll, Andrew B.
2014-03-01
Recent experiments have shown that polystyrene-polyacrylic-acid-polystyrene (PS-PAA-PS) triblock copolymers in a solvent mixture of water and toluene can self-assemble into spherical microcapsules. Suspended in water, the microcapsules have a toluene core surrounded by an elastomer triblock shell. The longer, hydrophilic PAA blocks remain near the outer surface of the shell, becoming charged through dissociation of OH functional groups in water, while the shorter, hydrophobic PS blocks form a networked (glass or gel) structure. Within a mean-field Poisson-Boltzmann theory, we model these polyelectrolyte microcapsules as spherical charged shells, assuming different dielectric constants inside and outside the capsule. By numerically solving the nonlinear Poisson-Boltzmann equation, we calculate the radial distribution of anions and cations and the osmotic pressure within the shell as a function of salt concentration. Our predictions, which can be tested by comparison with experiments, may guide the design of microcapsules for practical applications, such as drug delivery. This work was supported by the National Science Foundation under Grant No. DMR-1106331.
Derivation of the Second Law of Thermodynamics from Boltzmann's Distribution Law.
ERIC Educational Resources Information Center
Nelson, P. G.
1988-01-01
Shows how the thermodynamic condition for equilibrium in an isolated system can be derived by the application of Boltzmann's law to a simple physical system. States that this derivation could be included in an introductory course on chemical equilibrium to help prepare students for a statistical mechanical treatment presented in the curriculum.…
Liu, Bo; Liu, Pei; Xu, Zhenli; Zhou, Shenggao
2013-01-01
Near a charged surface, counterions of different valences and sizes cluster; and their concentration profiles stratify. At a distance from such a surface larger than the Debye length, the electric field is screened by counterions. Recent studies by a variational mean-field approach that includes ionic size effects and by Monte Carlo simulations both suggest that the counterion stratification is determined by the ionic valence-to-volume ratios. Central in the mean-field approach is a free-energy functional of ionic concentrations in which the ionic size effects are included through the entropic effect of solvent molecules. The corresponding equilibrium conditions define the generalized Boltzmann distributions relating the ionic concentrations to the electrostatic potential. This paper presents a detailed analysis and numerical calculations of such a free-energy functional to understand the dependence of the ionic charge density on the electrostatic potential through the generalized Boltzmann distributions, the role of ionic valence-to-volume ratios in the counterion stratification, and the modification of Debye length due to the effect of ionic sizes. PMID:24465094
Equilibrium Tail Distribution Due to Touschek Scattering
Nash,B.; Krinsky, S.
2009-05-04
Single large angle Coulomb scattering is referred to as Touschek scattering. In addition to causing particle loss when the scattered particles are outside the momentum aperture, the process also results in a non-Gaussian tail, which is an equilibrium between the Touschek scattering and radiation damping. Here we present an analytical calculation for this equilibrium distribution.
NASA Astrophysics Data System (ADS)
Qin, Feng; Zhao, Hua; Cai, Wei; Zhang, Zhiguo; Cao, Wenwu
2016-06-01
Noncontact monitoring temperature is very important in modern medicine, science, and technologies. The fluorescence intensity ratio (FIR) technique based on the Boltzmann distribution law exhibits excellent application potential, but the observed FIR deviates from the Boltzmann distribution law in the low temperature range. We propose a fluorescence intensity ratio relation FIR* = ηFIR by introducing a quantity η representing thermal population degree, which can be obtained from measured fluorescence decay curves of the upper emitting level. Using Eu3+ as an example, the method is confirmed that the deviated FIR is able to be corrected and return to follow the Boltzmann law.
NASA Astrophysics Data System (ADS)
Liang, Jun; Liu, Yan-Chun; Zhu, Qiao
2014-02-01
In order to further explore the effects of non-Gaussian smeared mass distribution on the thermodynamical properties of noncommutative black holes, we consider noncommutative black holes based on Maxwell-Boltzmann smeared mass distribution in (2+1)-dimensional spacetime. The thermodynamical properties of the black holes are investigated, including Hawking temperature, heat capacity, entropy and free energy. We find that multiple black holes with the same temperature do not exist, while there exists a possible decay of the noncommutative black hole based on Maxwell-Boltzmann smeared mass distribution into the rotating (commutative) BTZ black hole.
[Equilibrium surface charge distribution in phospholipid vesicles. I. Method of calculation].
Tenchov, B G; Raĭchev, B D
1977-01-01
This paper presents a method of calculation of the surface charge equilibrium distribution between the two surfaces of a spherically closed phospholipid bilayer suspended in aqueous electrolyte solution. The net surface charge is supposed to be provided by the ionized polar groups of the phospholipid molecules. Its equilibrium distribution is found by minimization of the free electrostatic energy. The procedure of minimization utilizes the solution of the Poisson-Boltzmann equation which describes the double electric layers of the membrane and an expression for the membrane potential derived under the assumption of absence of charges in the membrane phase. An analytical solution of the problem in the range of validity of the linearized Poisson-Boltzman equation is obtained. It is shown that in this case an equilibrium transmembrane potential exists, and the surface charge density is greater at the outer surface of the vesicle. PMID:588604
NASA Astrophysics Data System (ADS)
Lucia, Umberto
2016-02-01
The balance of forces and processes between the system and the environment and the processes inside the system are the result of the flows of the quanta. Moreover, the transition between two thermodynamic states is the consequence of absorption or emission of quanta, but, during the transition, the entropy variation due to the irreversibility occurs and it breaks any symmetry of time. Consequently, the irreversibility is the result of a transition, a process, an interaction between the system and its environment. This interaction results completely time-irreversible for any real process because of irreversibility. As a consequence, a proof of the third law is obtained proving that the zero temperature state can be achieved only for an infinite work lost for dissipation or in an infinite time. The fundamental role of time both in equilibrium and in non equilibrium analysis is pointed out. Moreover, the non equilibrium temperature is related to the entropy generation and its fluctuation rate; indeed, non-stationary temperature means that the system has not yet attained free energy minimum state, i.e., the maximum entropy state; the consequence is that the zero temperature state can be achieved only for an infinite work lost for dissipation or in an infinite time. In engineering thermodynamics the efficiency is always obtained without any consideration on time, while, here, just the time is introduced as a fundamental quantity of the analysis of non equilibrium states.
Consistent Application of the Boltzmann Distribution to Residual Entropy in Crystals
ERIC Educational Resources Information Center
Kozliak, Evguenii I.
2007-01-01
Four different approaches to residual entropy (the entropy remaining in crystals comprised of nonsymmetric molecules like CO, N[subscript 2]O, FClO[subscript 3], and H[subscript 2]O as temperatures approach 0 K) are analyzed and a new method of its calculation is developed based on application of the Boltzmann distribution. The inherent connection…
ERIC Educational Resources Information Center
Castle, Karen J.
2007-01-01
In this undergraduate physical chemistry laboratory experiment, students acquire a high-resolution infrared absorption spectrum of carbon dioxide and use their data to show that the rotational-vibrational state populations follow a Boltzmann distribution. Data are acquired with a mid-infrared laser source and infrared detector. Appropriate…
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.
Far-from-equilibrium distribution from near-steady-state work fluctuations.
Marsland, Robert; England, Jeremy
2015-11-01
A long-standing goal of nonequilibrium statistical mechanics has been to extend the conceptual power of the Boltzmann distribution to driven systems. We report some new progress towards this goal. Instead of writing the nonequilibrium steady-state distribution in terms of perturbations around thermal equilibrium, we start from the linearized driven dynamics of observables about their stable fixed point, and expand in the strength of the nonlinearities encountered during typical fluctuations away from the fixed point. The first terms in this expansion retain the simplicity of known expansions about equilibrium, but can correctly describe the statistics of a certain class of systems even under strong driving. We illustrate this approach by comparison with a numerical simulation of a sheared Brownian colloid, where we find that the first two terms in our expansion are sufficient to account for the shear thinning behavior at high shear rates. PMID:26651660
Energy Distributions in Small Populations: Pascal versus Boltzmann
ERIC Educational Resources Information Center
Kugel, Roger W.; Weiner, Paul A.
2010-01-01
The theoretical distributions of a limited amount of energy among small numbers of particles with discrete, evenly-spaced quantum levels are examined systematically. The average populations of energy states reveal the pattern of Pascal's triangle. An exact formula for the probability that a particle will be in any given energy state is derived.…
NASA Astrophysics Data System (ADS)
Feng, Yue
Plasma is currently a hot topic and it has many significant applications due to its composition of both positively and negatively charged particles. The energy distribution function is important in plasma science since it characterizes the ability of the plasma to affect chemical reactions, affect physical outcomes, and drive various applications. The Boltzmann Transport Equation is an important kinetic equation that provides an accurate basis for characterizing the distribution function---both in energy and space. This dissertation research proposes a multi-term approximation to solve the Boltzmann Transport Equation by treating the relaxation process using an expansion of the electron distribution function in Legendre polynomials. The elastic and 29 inelastic cross sections for electron collisions with nitrogen molecules (N2) and singly ionized nitrogen molecules ( N+2 ) have been used in this application of the Boltzmann Transport Equation. Different numerical methods have been considered to compare the results. The numerical methods discussed in this thesis are the implicit time-independent method, the time-dependent Euler method, the time-dependent Runge-Kutta method, and finally the implicit time-dependent relaxation method by generating the 4-way grid with a matrix solver. The results show that the implicit time-dependent relaxation method is the most accurate and stable method for obtaining reliable results. The results were observed to match with the published experimental data rather well.
NASA Astrophysics Data System (ADS)
Bergeron, H.; Curado, E. M. F.; Gazeau, J. P.; Rodrigues, Ligia M. C. S.
2016-02-01
Asymptotic behavior (with respect to the number of trials) of symmetric generalizations of binomial distributions and their related entropies is studied through three examples. The first one has the q-exponential as the generating function, the second one involves the modified Abel polynomials, and the third one has Hermite polynomials. We prove analytically that the Rényi entropy is extensive for these three cases, i.e., it is proportional (asymptotically) to the number n of events and that q-exponential and Hermite cases have also extensive Boltzmann-Gibbs. The Abel case is exceptional in the sense that its Boltzmann-Gibbs entropy is not extensive and behaves asymptotically as the square root of n. This result is obtained numerically and also confirmed analytically, under reasonable assumptions, by using a regularization of the beta function and its derivative. Probabilistic urn and genetic models are presented for illustrating this remarkable case.
Yuan Jiankui; Jette, David; Chen Weimin
2008-09-15
A photon transport algorithm for fully three-dimensional radiotherapy treatment planning has been developed based on the discrete ordinates (S{sub N}) solution of the Boltzmann equation. The algorithm is characterized by orthogonal adaptive meshes, which place additional points where large gradients occur and a procedure to evaluate the collided flux using the representation of spherical harmonic expansion instead of the summation of the volume-weighted contribution from discrete angles. The Boltzmann equation was solved in the form of S{sub N} spatial, energy, and angular discretization with mitigation of ray effects by the first-collision source method. Unlike existing S{sub N} codes, which were designed for general purpose for multiparticle transport in areas such as nuclear engineering, our code is optimized for medical radiation transport. To validate the algorithm, several examples were employed to calculate the photon flux distribution. Numerical results show good agreement with the Monte Carlo calculations using EGSnrc.
Fully relativistic lattice Boltzmann algorithm
Romatschke, P.; Mendoza, M.; Succi, S.
2011-09-15
Starting from the Maxwell-Juettner equilibrium distribution, we develop a relativistic lattice Boltzmann (LB) algorithm capable of handling ultrarelativistic systems with flat, but expanding, spacetimes. The algorithm is validated through simulations of a quark-gluon plasma, yielding excellent agreement with hydrodynamic simulations. The present scheme opens the possibility of transferring the recognized computational advantages of lattice kinetic theory to the context of both weakly and ultrarelativistic systems.
Verheest, Frank; Hellberg, Manfred A.; Kourakis, Ioannis
2013-08-15
Acoustic supersolitons arise when a plasma model is able to support three consecutive local extrema of the Sagdeev pseudopotential between the undisturbed conditions and an accessible root. This leads to a characteristic electric field signature, where a simple bipolar shape is enriched by subsidiary maxima. Large-amplitude nonlinear acoustic modes are investigated, using a pseudopotential approach, for plasmas containing two-temperature electrons having Boltzmann or kappa distributions, in the presence of cold fluid ions. The existence domains for positive supersolitons are derived in a methodological way, both for structure velocities and amplitudes, in terms of plasma compositional parameters. In addition, typical pseudopotentials, soliton, and electric field profiles have been given to illustrate that positive supersolitons can be found in the whole range of electron distributions from Maxwellian to a very hard nonthermal spectrum in kappa. However, it is found that the parameter ranges that support supersolitons vary significantly over the wide range of kappa considered.
Technology Transfer Automated Retrieval System (TEKTRAN)
The distribution coefficient (KD) for the human drug carbamazepine was measured using a non-equilibrium technique. Repacked soil columns were prepared using an Airport silt loam (Typic Natrustalf) with an average organic matter content of 2.45%. Carbamazepine solutions were then leached through th...
NASA Astrophysics Data System (ADS)
García-García, J.; Martín, F.
2000-11-01
From a coupling model between the Boltzmann transport equation and the quantum Liouville equation, we have developed a simulator based on the Wigner distribution function (WDF) approach that can be applied to resonant tunneling diodes (RTDs) and other vertical transport quantum devices. In comparison to previous WDF simulators, the tool allows one to extend the simulation domains up to hundreds of nanometers, which are the typical dimensions required for the study of actual multilayer structures. With these improvements, a level of agreement between theory and experiment comparable to that obtained by using other simulators based on Green functions has been achieved. The results of this work reveal that the WDF formalism can be alternatively used to study the behavior of actual multilayered RTDs.
ERIC Educational Resources Information Center
Angeli, Celestino; Cimiraglia, Renzo; Dallo, Federico; Guareschi, Riccardo; Tenti, Lorenzo
2013-01-01
The dependence on the temperature of the population of the "i"th state, "P"[subscript "i"], in the Boltzmann distribution is analyzed by studying its derivative with respect to the temperature, "T." A simple expression is found, involving "P"[subscript "i"], the energy of the state,…
The Equilibrium Distribution of Income and the Market for Status.
ERIC Educational Resources Information Center
Becker, Gary S.; Murphy, Kevin M.; Werning, Ivan
2005-01-01
This paper explores the implications for risk-taking behavior and the equilibrium distribution of income of assuming that the desire for status positions is a powerful motive and that it raises the marginal utility of consumption. In contrast to previous analyses, we consider the case in which status positions are sold in a hedonic market. We show…
NASA Astrophysics Data System (ADS)
Norouzi, Ali; Esfahani, Javad Abolfazli
2016-02-01
In this study, gaseous flow through a micro/nano-channel is investigated via a novel two relaxation time lattice Boltzmann method. In this method, the slip velocity at the fluid-solid interface is realized by defining the free relaxation parameter. Furthermore, in order to capture the non-linear phenomena associated with the Knudsen layer, the wall function correction is employed. To this respect, different available wall functions are implemented. The objective of the study is to provide a comparative study on the accuracy, range of applicability and computational efficiency of these wall functions in a wide range of Knudsen numbers. The results of the present study are compared against direct simulation Mont Carlo and information preservation data. It is found that only a few of the implemented wall functions are capable of predicting the flow behavior with reasonable accuracy, particularly when the Knudsen number lies in the transition flow regime.
Lattice Boltzmann model for compressible fluids
NASA Technical Reports Server (NTRS)
Alexander, F. J.; Chen, H.; Chen, S.; Doolen, G. D.
1992-01-01
A lattice Boltzmann model is derived which simulates compressible fluids. By choosing the parameters of the equilibrium distribution appropriately, the sound speed (which may be set arbitrarily low), bulk viscosity, and kinematic viscosity can be selected. This model simulates compressible flows and can include shocks. With a proper rescaling and zero-sound speed, this model simulates Burgers's equation. The viscosity determined by a Chapman-Enskog expansion compares well with that measured form simulations. The exact solutions of Burgers's equation on the unit circle are compared to solutions of lattice Boltzmann model finding reasonable agreement.
Unsteady electroosmosis in a microchannel with Poisson-Boltzmann charge distribution.
Chang, Chien C; Kuo, Chih-Yu; Wang, Chang-Yi
2011-11-01
The present study is concerned with unsteady electroosmotic flow (EOF) in a microchannel with the electric charge distribution described by the Poisson-Boltzmann (PB) equation. The nonlinear PB equation is solved by a systematic perturbation with respect to the parameter λ which measures the strength of the wall zeta potential relative to the thermal potential. In the small λ limits (λ<1), we recover the linearized PB equation - the Debye-Hückel approximation. The solutions obtained by using only three terms in the perturbation series are shown to be accurate with errors <1% for λ up to 2. The accurate solution to the PB equation is then used to solve the electrokinetic fluid transport equation for two types of unsteady flow: transient flow driven by a suddenly applied voltage and oscillatory flow driven by a time-harmonic voltage. The solution for the transient flow has important implications on EOF as an effective means for transporting electrolytes in microchannels with various electrokinetic widths. On the other hand, the solution for the oscillatory flow is shown to have important physical implications on EOF in mixing electrolytes in terms of the amplitude and phase of the resulting time-harmonic EOF rate, which depends on the applied frequency and the electrokinetic width of the microchannel as well as on the parameter λ. PMID:22072500
Application of the multi distribution function lattice Boltzmann approach to thermal flows
NASA Astrophysics Data System (ADS)
Parmigiani, A.; Huber, C.; Chopard, B.; Latt, J.; Bachmann, O.
2009-04-01
Numerical methods able to model high Rayleigh ( Ra) and high Prandtl ( Pr) number thermal convection are important to study large-scale geophysical phenomena occuring in very viscous fluids such as magma chamber dynamics (104 < Pr < 107 and 107 < Ra < 1011). The important variable to quantify the thermal state of a convective fluid is a generalized dimensionless heat transfer coefficient (the Nusselt number) whose measure indicates the relative efficiency of the thermal convection. In this paper we test the ability of Multi-distribution Function approach (MDF) Thermal Lattice Boltzmann method to study the well-established scaling result for the Nusselt number ( Nu ∝ Ra 1/3) in Rayleigh Bénard convection for 104 ≤ Ra ≤ 109 and 101 ≤ Pr ≤ 104. We explore its main drawbacks in the range of Pr and Ra number under investigation: (1) high computational time N c required for the algorithm to converge and (2) high spatial accuracy needed to resolve the thickness of thermal plumes and both thermal and velocity boundary layer. We try to decrease the computational demands of the method using a multiscale approach based on the implicit dependence of the Pr number on the relaxation time, the spatial and temporal resolution characteristic of the MDF thermal model.
Spatial distributions at equilibrium under heterogeneous transient subdiffusion
Berry, Hugues; Soula, Hédi A.
2014-01-01
Experimental measurements of the mobility of macromolecules, especially proteins, in cells and their membranes consistently report transient subdiffusion with possibly position-dependent—non-homogeneous—properties. However, the spatiotemporal dynamics of protein mobility when transient subdiffusion is restricted to a subregion of space is still unclear. Here, we investigated the spatial distribution at equilibrium of proteins undergoing transient subdiffusion due to continuous-time random walks (CTRW) in a restricted subregion of a two-dimensional space. Our Monte-Carlo simulations suggest that this process leads to a non-homogeneous spatial distribution of the proteins at equilibrium, where proteins increasingly accumulate in the CTRW subregion as its anomalous properties are increasingly marked. In the case of transient CTRW, we show that this accumulation is dictated by the asymptotic Brownian regime and not by the initial anomalous transient dynamics. Moreover, our results also show that this dominance of the asymptotic Brownian regime cannot be simply generalized to other scenarios of transient subdiffusion. In particular, non-homogeneous transient subdiffusion due to hindrance by randomly-located immobile obstacles does not lead to such a strong local accumulation. These results suggest that, even though they exhibit the same time-dependence of the mean-squared displacement, the different scenarios proposed to account for subdiffusion in the cell lead to different protein distribution in space, even at equilibrium and without coupling with reaction. PMID:25429273
Yousfi, M.; Himoudi, A.; Gaouar, A. )
1992-12-15
Electron distribution functions in nonthermal cold plasmas generated by classical electrical discharges have been calculated from a powerful Boltzmann equation solution and an original Monte Carlo simulation. In these two methods both classical (i.e., elastic, inelastic, and superelastic) electron-atom (or molecule) collisions and electron-electron interactions are taken into account. The approximations considered to include long-range (electron-electron) and short-range (electron-atom) interactions in the same Monte Carlo algorithm are first validated by comparing with Boltzmann equation results. Then, the influence of electron-electron interactions on electron distribution functions, swarm parameters, and reaction rates under nonthermal cold plasma conditions are analyzed and discussed as a function of reduced electric field [ital E]/[ital N] and ionization degree [ital n][sub [ital e
Equilibrium Distributions and the Nanostructure Diagram for Epitaxial Quantum Dots
Rudd, R E; Briggs, G D; Sutton, A P; Medeiros-Ribeiro, G; Williams, R S
2006-05-01
We present in detail a thermodynamic equilibrium model for the growth of nanostructures on semiconductor substrates in heteroepitaxy and its application to germanium deposition on silicon. Some results of this model have been published previously, but the details of the formulation of the model are given here for the first time. The model allows the computation of the shape and size distributions of the surface nanostructures, as well as other properties of the system. We discuss the results of the model, and their incorporation into a nanostructure diagram that summarizes the relative stability of domes and pyramids in the bimodal size distributions.
Dechant, Andreas; Shafier, Shalom Tzvi; Kessler, David A; Barkai, Eli
2016-08-01
The Boltzmann-Gibbs density, a central result of equilibrium statistical mechanics, relates the energy of a system in contact with a thermal bath to its equilibrium statistics. This relation is lost for nonthermal systems such as cold atoms in optical lattices, where the heat bath is replaced with the laser beams of the lattice. We investigate in detail the stationary phase-space probability for Sisyphus cooling under harmonic confinement. In particular, we elucidate whether the total energy of the system still describes its stationary state statistics. We find that this is true for the center part of the phase-space density for deep lattices, where the Boltzmann-Gibbs density provides an approximate description. The relation between energy and statistics also persists for strong confinement and in the limit of high energies, where the system becomes underdamped. However, the phase-space density now exhibits heavy power-law tails. In all three cases we find expressions for the leading-order phase-space density and corrections which break the equivalence of probability and energy and violate energy equipartition. The nonequilibrium nature of the steady state is corroborated by explicit violations of detailed balance. We complement these analytical results with numerical simulations to map out the intricate structure of the phase-space density. PMID:27627290
NASA Astrophysics Data System (ADS)
Blanco, R.; Pesquera, L.
1986-01-01
We analyze in detail the validity of the compatibility between the Maxwell-Boltzmann (MB) distribution and the Rayleigh-Jeans (RJ) spectrum obtained in recent works by the authors for classical relativistic systems. We show that the MB distribution and the RJ spectrum are not compatible if we do not remove high enough frequencies. By analyzing the applicability of the approximation methods used in previous works to obtain the MB distribution from the RJ spectrum, we conclude that these methods are valid only if we introduce a high-frequency cutoff in the RJ spectrum. A short discussion is made on the meaning of this cutoff.
NASA Technical Reports Server (NTRS)
Lanyi, Gabor E.
2003-01-01
This viewgraph presentation reviews the 1901 work in Planck's constant and blackbody radiation law and the 1916 Einstein rederivation of the blackbody radiation law. It also reviews Wien's law. It also presents equations that demonstrate the thermal balance between radiation and matter.
Williams, C F; Watson, J E; Nelson, S D
2014-01-01
The distribution coefficient (KD) for the human drug carbamazepine was measured using a non-equilibrium technique. Repacked soil columns were prepared using an Airport silt loam (Typic Natrustalf) with an average organic matter content of 2.45%. Carbamazepine solutions were then leached through the columns at 0.5, 1.0 and 1.5 mL min(-1) representing average linear velocities of 1.8, 3.5 and 5.3 cm h(-1) respectively. Each flow rate was replicated three times and three carbamazepine pulses were applied to each column resulting in a total of 9 columns with 27 total carbamazepine pulses. Breakthrough curves were used to determine KD using the parameter fitting software CXTFIT. Results indicate that as flow rate decreased from 5.3 to 1.8 cm h(-1), KD increased an average of 21%. Additionally, KD determined by column leaching (14.7-22.7 L kg(-1)) was greater than KD determined by a 2h batch equilibrium adsorption (12.6 L kg(-1)). Based on these KD's carbamazepine would be generally characterized as non-mobile in the soil investigated. However, repeated carbamazepine applications resulted in an average 22% decrease in KD between the first and third applications. Decreasing KD is attributed to differences in sorption site kinetics and carbamazepine residence time in contact with the soil. This would indicate that the repeated use of reclaimed wastewater at high application rates for long-term irrigation or groundwater recharge has the potential to lead to greater transport of carbamazepine than KD determined by batch equilibrium would predict. PMID:24050717
Truncated Thermal Equilibrium Distribution for Intense Beam Propagation
Ronald C. Davidson; Hong Qin; Steven M. Lund
2003-02-26
An intense charged-particle beam with directed kinetic energy ({lambda}{sub b}-1)m{sub b}c{sup 2} propagates in the z-direction through an applied focusing field with transverse focusing force modeled by F{sub foc} = -{lambda}{sub b}m{sub b}{omega}{sub beta}{sup 2} {perpendicular} x {perpendicular} in the smooth focusing approximation. This paper examines properties of the axisymmetric, truncated thermal equilibrium distribution F(sub)b(r,p perpendicular) = A exp (-H Perpendicular/T perpendicular (sub)b) = (H perpendicular-E(sub)b), where A, T perpendicular (sub)b, and E (sub)b are positive constants, and H perpendicular is the Hamiltonian for transverse particle motion. The equilibrium profiles for beam number density, n(sub)b(r) = * d{sup 2}pF(sub)b(r,p perpendicular), and transverse temperature, T perpendicular (sub)b(r) = * d{sup 2}p(p{sup 2} perpendicular/2 lambda (sbu)bm (sub)b)F(sub)b(r,p perpendicular), are calculated self-consistently including space-charge effects. Several properties of the equilibrium profiles are noteworthy. For example, the beam has a sharp outer edge radius r(sub)b with n(sub)b(r greater than or equal to rb) = 0, where r(sub)b depends on the value of E(sub)b/T (sub)perpendicular(sub)b. In addition, unlike the choice of a semi-Gaussian distribution, F{sup SG}(sub)b = A exp (-p{sup 2}(sub)perpendicular/2lambda(sub)bm(sub)bTperpendicular(sub)b) = (r-r(sub)b), the truncated thermal equilibrium distribution F(sub)b(r,p) depends on (r,p) only through the single-particle constant of the motion Hperpendiuclar and is therefore a true steady-state solution (*/*t = 0) of the nonlinear Vlasov-Maxwell equations.
Extended Tonks-Langmuir-type model with non-Boltzmann-distributed electrons and cold ion sources
NASA Astrophysics Data System (ADS)
Kamran, M.; Kuhn, S.; Tskhakaya, D. D.; Khan, M.; Khan
2013-04-01
kinetic Tonks-Langmuir model. Phys. Plasmas 13, 063508) or bi-Maxwellian (Godyak, V. A. et al. 1995 Tonks-Langmuir problem for a bi-Maxwellian plasma. IEEE Trans. Plasma Sci. 23, 728) electron velocity distribution functions (VDFs), which satisfy the zero-CSS-term (Vlasov) kinetic equation and imply zero electron currents, we here propose a more general class of electron VDFs allowing, in an approximate manner, for non-zero CSS terms and finite electron currents inside the plasma region. The sheath-edge and floating-wall potentials are calculated by balancing the ion and electron current densities at sheath-edge singularities. In a first detailed application, the type-t and type-p electron VDFs are assumed to be `inner' and `outer' cut-off Maxwellians respectively, with different amplitudes and `formal' temperatures, implying the perfectly CSS-free limit. For the special case of equal type-t and type-p electron VDF amplitudes and formal temperatures, the classical Boltzmann distribution for electrons is formally retrieved. Special cases with other amplitude and formal-temperature ratios show significant deviations from the classical case.
Modeling flue pipes: Subsonic flow, lattice Boltzmann, and parallel distributed computers
NASA Astrophysics Data System (ADS)
Skordos, Panayotis A.
1995-01-01
The problem of simulating the hydrodynamics and the acoustic waves inside wind musical instruments such as the recorder the organ, and the flute is considered. The problem is attacked by developing suitable local-interaction algorithms and a parallel simulation system on a cluster of non-dedicated workstations. Physical measurements of the acoustic signal of various flue pipes show good agreement with the simulations. Previous attempts at this problem have been frustrated because the modeling of acoustic waves requires small integration time steps which make the simulation very compute-intensive. In addition, the simulation of subsonic viscous compressible flow at high Reynolds numbers is susceptible to slow-growing numerical instabilities which are triggered by high-frequency acoustic modes. The numerical instabilities are mitigated by employing suitable explicit algorithms: lattice Boltzmann method, compressible finite differences, and fourth-order artificial-viscosity filter. Further, a technique for accurate initial and boundary conditions for the lattice Boltzmann method is developed, and the second-order accuracy of the lattice Boltzmann method is demonstrated. The compute-intensive requirements are handled by developing a parallel simulation system on a cluster of non-dedicated workstations. The system achieves 80 percent parallel efficiency (speedup/processors) using 20 HP-Apollo workstations. The system is built on UNIX and TCP/IP communication routines, and includes automatic process migration from busy hosts to free hosts.
Modeling Flue Pipes: Subsonic Flow, Lattice Boltzmann, and Parallel Distributed Computers.
NASA Astrophysics Data System (ADS)
Skordos, Panayotis A.
1995-01-01
The problem of simulating the hydrodynamics and the acoustic waves inside wind musical instruments such as the recorder, the organ, and the flute is considered. The problem is attacked by developing suitable local-interaction algorithms and a parallel simulation system on a cluster of non-dedicated workstations. Physical measurements of the acoustic signal of various flue pipes show good agreement with the simulations. Previous attempts at this problem have been frustrated because the modeling of acoustic waves requires small integration time steps which make the simulation very compute-intensive. In addition, the simulation of subsonic viscous compressible flow at high Reynolds numbers is susceptible to slow-growing numerical instabilities which are triggered by high-frequency acoustic modes. The numerical instabilities are mitigated by employing suitable explicit algorithms: lattice Boltzmann method, compressible finite differences, and fourth-order artificial -viscosity filter. Further, a technique for accurate initial and boundary conditions for the lattice Boltzmann method is developed, and the second-order accuracy of the lattice Boltzmann method is demonstrated. The compute-intensive requirements are handled by developing a parallel simulation system on a cluster of non-dedicated workstations. The system achieves 80% parallel efficiency (speedup/processors) using 20 HP-Apollo workstations. The system is built on UNIX and TCP/IP communication routines, and includes automatic process migration from busy hosts to free hosts. (Copies available exclusively from MIT Libraries, Rm. 14-0551, Cambridge, MA 02139-4307. Ph. 617 -253-5668; Fax 617-253-1690.).
On the electron equilibrium distribution function in the kinetic theory of electron cyclotron maser
NASA Astrophysics Data System (ADS)
Shenggang, Liu
1981-11-01
The problems concerning the specification of electron equilibrium distribution function for the kinetic theory of ECRM are investigated in this paper. After detailed analysis of the published equilibium distribution functions, several conclusion have been achieved.
Equilibrium distribution of shapes for linear and star macromolecules
NASA Astrophysics Data System (ADS)
Cannon, Joel W.; Aronovitz, Joseph A.; Goldbart, Paul
1991-05-01
We investigate the equilibrium distribution of three-dimensional shapes adopted by isolated linear and star-shaped macromolecules, both with and without intramolecular interactions, using an implementation of the Monte Carlo method suitable for macromolecules with branches. We compute the joint probability distribution function for two quantities which together characterise invariant features of the normalised radius of gyration tensor associated with the shape (rather than size) of macromolecular configurations. Amongst other things, knowledge of this distribution function allows us to compute the expectation values <~ngle Δ0rangle and <~ngle S0rangle introduced by Aronovitz and Nelson (J. Phys. France 47 (1986) 1445) to characterise the extent and nature of anisotropy of typical shapes drawn from the ensemble of macromolecular configurations. We also compute a third expectation value <~ngle Σ0rangle which isolates the nature of the anisotropy from its extent. Furthermore, our simulation permits a comparison of <~ngle Δ0rangle and <~ngle S0rangle with the less natural (but analytically tractable) alternative quantities Δ, the asphericity examined by Rudnick and Gaspari (J. Phys. A 19 (1986) L191) and by Aronovitz and Nelson, and S, examined by Aronovitz and Nelson, which have enhanced sensititvity to larger configurations and therefore convolve shape information with size information. It is found that although Δ and S do provide some characterisation of anisotropy, they differ considerably from the natural measures <~ngle Δ0rangle and <~ngle S0rangle. In particular, if Δ and S are regarded as approximations to Δ0rangle and <~ngle S0rangle then, for both linear and branched macromolecules, they severely underestimate the increase (or overestimate the decrease) in extent and prolateness of anisotropy due to intramolecular interactions. Nous examinons la distribution à l'équilibre des formes tri-dimensionnelles prises par des macromolécules isolées de
Modified lattice Boltzmann method for compressible fluid simulations.
Hinton, F L; Rosenbluth, M N; Wong, S K; Lin-Liu, Y R; Miller, R L
2001-06-01
A modified lattice Boltzmann algorithm is shown to have much better stability to growing temperature perturbations, when compared with the standard lattice Boltzmann algorithm. The damping rates of long-wavelength waves, which determine stability, are derived using a collisional equilibrium distribution function which has the property that the Euler equations are obtained exactly in the limit of zero time step. Using this equilibrium distribution function, we show that our algorithm has inherent positive hyperviscosity and hyperdiffusivity, for very small values of viscosity and thermal diffusivity, which are lacking in the standard algorithm. Short-wavelength modes are shown to be stable for temperatures greater than a lower limit. Results from a computer code are used to compare these algorithms, and to confirm the damping rate predictions made analytically. Finite amplitude sound waves in the simulated fluid steepen, as expected from gas dynamic theory. PMID:11415085
Kikuchi, Takashi; Horioka, Kazuhiko
2009-05-15
Possible emittance growths of intense, nonuniform beams during a transport in a focusing channel are derived as a function of nonlinear field energy and space charge tune depression factors. The nonlinear field energy of the beam with thermal equilibrium distribution is estimated by considering the particle distribution across the cross section of the beam. The results show that the possible emittance growth can be suppressed by keeping the beam particle in thermal equilibrium distribution during the beam transport.
IS THE SIZE DISTRIBUTION OF URBAN AEROSOLS DETERMINED BY THERMODYNAMIC EQUILIBRIUM? (R826371C005)
A size-resolved equilibrium model, SELIQUID, is presented and used to simulate the size–composition distribution of semi-volatile inorganic aerosol in an urban environment. The model uses the efflorescence branch of aerosol behavior to predict the equilibrium partitioni...
NASA Astrophysics Data System (ADS)
Imai, M.; Sataka, M.; Matsuda, M.; Okayasu, S.; Kawatsura, K.; Takahiro, K.; Komaki, K.; Shibata, H.; Nishio, K.
2015-07-01
Both equilibrium and non-equilibrium charge-state distributions were studied experimentally for 2.0 MeV/u carbon ions after passing through carbon foils. Measured charge-state distribution established the equilibrium at a target thickness of 10 μg/cm2 and this remained unchanged until a maximum target thickness of 98 μg/cm2. The equilibrium charge-state distribution, the equilibrium mean charge-state, and the width and skewness of the equilibrium distribution were compared with predictions using existing semi-empirical formulae as well as simulation results, including the ETACHA code. It was found that charge-state distributions, mean charge states, and distribution widths for C2+, C3+, and C4+ incident ions merged into quasi-equilibrium values at a target thickness of 5.7 μg/cm2 in the pre-equilibrium region and evolved simultaneously to the 'real equilibrium' values for all of the initial charge states, including C5+ and C6+ ions, as previously demonstrated for sulfur projectile ions at the same velocity (Imai et al., 2009). Two kinds of simulation, ETACHA and solution of rate equations taking only single electron transfers into account, were used, and both of them reproduced the measured charge evolution qualitatively. The quasi-equilibrium behavior could be reproduced with the ETACHA code, but not with solution of elementary rate equations.
NASA Astrophysics Data System (ADS)
Torrisi, L.
2016-02-01
The charge production from laser-generated plasmas generates not isotropically ion acceleration in vacuum and with mean kinetic energy proportional to the ion charge state. The ion velocity depends on many factors of which the most important are the plasma temperature, the adiabatic gas expansion in vacuum and the Coulomb acceleration. The ion energy distributions of the emitted ions from the plasma can be well explained by the Coulomb-Boltzmann-Shifted function, with a cut-off limitation at high energy for a wide range of laser intensities. It can be applied for intensities of 1010 W/cm2, when plasma is produced only in the backward direction from thick targets (backward plasma acceleration regime), as well as at intensities of the order of 1019 W/cm2, when plasma is produced in the forward direction from thin targets in target-normal sheath acceleration regime. It loses of validity in radiation pressure acceleration regime, at which ions are emitted near mono-energetically.
Sonnad, Kiran G.; Cary, John R.
2015-04-15
A procedure to obtain a near equilibrium phase space distribution function has been derived for beams with space charge effects in a generalized periodic focusing transport channel. The method utilizes the Lie transform perturbation theory to canonically transform to slowly oscillating phase space coordinates. The procedure results in transforming the periodic focusing system to a constant focusing one, where equilibrium distributions can be found. Transforming back to the original phase space coordinates yields an equilibrium distribution function corresponding to a constant focusing system along with perturbations resulting from the periodicity in the focusing. Examples used here include linear and nonlinear alternating gradient focusing systems. It is shown that the nonlinear focusing components can be chosen such that the system is close to integrability. The equilibrium distribution functions are numerically calculated, and their properties associated with the corresponding focusing system are discussed.
Polar Coordinate Lattice Boltzmann Kinetic Modeling of Detonation Phenomena
NASA Astrophysics Data System (ADS)
Lin, Chuan-Dong; Xu, Ai-Guo; Zhang, Guang-Cai; Li, Ying-Jun
2014-11-01
A novel polar coordinate lattice Boltzmann kinetic model for detonation phenomena is presented and applied to investigate typical implosion and explosion processes. In this model, the change of discrete distribution function due to local chemical reaction is dynamically coupled into the modified lattice Boltzmann equation which could recover the Navier—Stokes equations, including contribution of chemical reaction, via the Chapman—Enskog expansion. For the numerical investigations, the main focuses are the nonequilibrium behaviors in these processes. The system at the disc center is always in its thermodynamic equilibrium in the highly symmetric case. The internal kinetic energies in different degrees of freedom around the detonation front do not coincide. The dependence of the reaction rate on the pressure, influences of the shock strength and reaction rate on the departure amplitude of the system from its local thermodynamic equilibrium are probed.
Size-frequency distribution of crater populations in equilibrium on the Moon
NASA Astrophysics Data System (ADS)
Xiao, Zhiyong; Werner, Stephanie C.
2015-12-01
Overprinting of craters by subsequent impacts and topographic degradation complicates crater statistics, especially for old surfaces and small-diameter crater populations. A crater population is regarded as in equilibrium at a particular crater size when smaller craters are being produced at the same rate at which they are being destroyed. Evaluating the equilibrium state of crater populations is challenging, and empirical equilibrium densities are frequently inferred. By performing careful crater counts and cross comparisons on several lunar surfaces, we study the size-frequency distributions (SFD) for the crater populations, which have portions in equilibrium. The results are one of the few observational constraints on the SFD of crater populations in equilibrium, showing that referring to empirical equilibrium densities is not safe for evaluating the equilibrium states of crater populations. Equilibrium densities are not positively correlated with the ages of crater populations, and some populations in equilibrium have crater densities less than those previously believed to represent equilibrium conditions. Besides the SFD of the production population, different crater removal rates at different diameters also affect the SFD of crater populations in equilibrium. The equilibrium onset diameter (Deq) of a crater population can be translated to model ages because older populations have larger Deq, and those for same-aged surfaces are comparable. We show that the crater populations studied here are in equilibrium at much smaller diameters than those predicted for same-aged surfaces by crater degradation models, thus indicating lower crater degradation rates on the Moon, and/or younger ages of the counting areas.
Transition in the Equilibrium Distribution Function of Relativistic Particles
Mendoza, M.; Araújo, N. A. M.; Succi, S.; Herrmann, H. J.
2012-01-01
We analyze a transition from single peaked to bimodal velocity distribution in a relativistic fluid under increasing temperature, in contrast with a non-relativistic gas, where only a monotonic broadening of the bell-shaped distribution is observed. Such transition results from the interplay between the raise in thermal energy and the constraint of maximum velocity imposed by the speed of light. We study the Bose-Einstein, the Fermi-Dirac, and the Maxwell-Jüttner distributions, and show that they all exhibit the same qualitative behavior. We characterize the nature of the transition in the framework of critical phenomena and show that it is either continuous or discontinuous, depending on the group velocity. We analyze the transition in one, two, and three dimensions, with special emphasis on twodimensions, for which a possible experiment in graphene, based on the measurement of the Johnson-Nyquist noise, is proposed. PMID:22937220
Transition in the Equilibrium Distribution Function of Relativistic Particles
NASA Astrophysics Data System (ADS)
Mendoza, M.; Araújo, N. A. M.; Succi, S.; Herrmann, H. J.
2012-08-01
We analyze a transition from single peaked to bimodal velocity distribution in a relativistic fluid under increasing temperature, in contrast with a non-relativistic gas, where only a monotonic broadening of the bell-shaped distribution is observed. Such transition results from the interplay between the raise in thermal energy and the constraint of maximum velocity imposed by the speed of light. We study the Bose-Einstein, the Fermi-Dirac, and the Maxwell-Jüttner distributions, and show that they all exhibit the same qualitative behavior. We characterize the nature of the transition in the framework of critical phenomena and show that it is either continuous or discontinuous, depending on the group velocity. We analyze the transition in one, two, and three dimensions, with special emphasis on twodimensions, for which a possible experiment in graphene, based on the measurement of the Johnson-Nyquist noise, is proposed.
NASA Astrophysics Data System (ADS)
Imai, M.; Sataka, M.; Kawatsura, K.; Takahiro, K.; Komaki, K.; Shibata, H.; Sugai, H.; Nishio, K.
2009-08-01
Both equilibrium and non-equilibrium charge-state distributions for 2.0 MeV/u sulfur ions after passing through carbon foils were studied experimentally. For the equilibrium charge-state distribution, incident ions of S 7+, S 12+, S 14+ and S 16+ were injected into carbon foils 54, 98, 150 and 200 μg/cm 2 in thickness, whereas for the non-equilibrium distributions, new measurements for S 15+ and S 16+ incidences were made through carbon foils of 0.9-10 μg/cm 2 to supplement our previous experiments regarding S 6+-S 14+ incidences [M. Imai, M. Sataka, K. Kawatsura, K. Takahiro, K. Komaki, H. Shibata, H. Sugai, K. Nishio, Nucl. Instr. and Meth. B 230 (2005) 63; M. Imai, M. Sataka, K. Kawatsura, K. Takahiro, K. Komaki, H. Shibata, H. Sugai, K. Nishio, Nucl. Instr. and Meth. B 256 (2007) 11]. Mean charge states for S 6+-S 14+ incidences as functions of the penetration thickness merged at 6.9 μg/cm 2 and changed together until reaching equilibrium at around 100 μg/cm 2, while those for S 15+ and S 16+ incidences took different paths to equilibrium, which was also the case for distribution widths for S 6+-S 14+, S 15+ and S 16+ incidences. An equilibrium mean charge state of 12.68 and distribution width of 1.11 were attained with equilibrium charge distributions between 6+ and 16+.
NASA Astrophysics Data System (ADS)
Salimi, M. R.; Taeibi-Rahni, M.
2015-12-01
Due to superior accuracy and stability of multiple relaxation time (MRT) collision operator over its single relaxation time (SRT) counterpart, new lifting relations are proposed here to construct single particle distribution functions for MRT-LBM from macroscopic variables. Using these lifting relations, a new hybrid FVM-LB method is presented (called Finite type-LB hybrid method), which is consistent with MRT-LBM. In this new hybrid method, single-particle distribution functions in MRT-LBM sub-domain boundaries are computed, using equilibrium and non-equilibrium moments. These moments are computed in Navier-Stokes/FVM sub-domain boundaries, using macroscopic variables and their derivatives. The new method is validated by solving three benchmark problems, i.e., two- and three-dimensional lid driven cavity flows and two-dimensional unsteady flow around a squared section cylinder. These problems are analyzed with pure FVM, pure LBM, and Finite type-LB hybrid method (FTLBHM) and the related results are compared with each other and with benchmark data. These comparisons clearly demonstrate the accuracy of the present novel methodology for simulating steady/unsteady flow fields in two and three dimensions.
Silvestre-Alcantara, Whasington; Bhuiyan, Lutful B.; Outhwaite, Christopher W.; Henderson, Douglas
2010-01-01
The properties of the singlet ion distributions at and around contact in a restricted primitive model double layer are characterized in the modified Poisson–Boltzmann theory. Comparisons are made with the corresponding exact Monte Carlo simulation data, the results from the Gouy–Chapman–Stern theory coupled to an exclusion volume term, and the mean spherical approximation. Particular emphasis is given to the behaviour of the theoretical predictions in relation to the contact value theorem involving the charge profile. The simultaneous behaviour of the coion and counterion contact values is also examined. The performance of the modified Poisson–Boltzmann theory in regard to the contact value theorems is very reasonable with the contact characteristics showing semi-quantitative or better agreement overall with the simulation results. The exclusion-volume-treated Gouy–Chapman–Stern theory reveals a fortuitous cancellation of errors, while the mean spherical approximation is poor. PMID:20664814
A Lattice-Boltzmann model for suspensions of self-propelling colloidal particles.
Ramachandran, S; Sunil Kumar, P B; Pagonabarraga, I
2006-06-01
We present a Lattice-Boltzmann method for simulating self-propelling (active) colloidal particles in two dimensions. Active particles with symmetric and asymmetric force distribution on their surface are considered. The velocity field generated by a single active particle, changing its orientation randomly, and the different time scales involved are characterized in detail. The steady-state speed distribution in the fluid, resulting from the activity, is shown to deviate considerably from the equilibrium distribution. PMID:16779527
A Lattice-Boltzmann model for suspensions of self-propelling colloidal particles
NASA Astrophysics Data System (ADS)
Ramachandran, S.; Kumar, P. B. Sunil; Pagonabarraga, I.
2006-06-01
We present a Lattice-Boltzmann method for simulating self-propelling (active) colloidal particles in two dimensions. Active particles with symmetric and asymmetric force distribution on their surface are considered. The velocity field generated by a single active particle, changing its orientation randomly, and the different time scales involved are characterized in detail. The steady-state speed distribution in the fluid, resulting from the activity, is shown to deviate considerably from the equilibrium distribution.
NASA Astrophysics Data System (ADS)
Blanco, R.; Pesquera, L.
1986-08-01
The problem of radiation-matter equilibrium for classical nonrelativistic rigid extended charges is studied by using the dipolar approximation for the fields. An exact solution is obtained for a large class of spherically symmetric distribution charges. It is shown that equilibrium holds with the Rayleigh-Jeans radiation. As concerns the equilibrium distribution for the matter, the Maxwell-Boltzmann law is obtained only when the radius is large enough for the electromagnetic mass to be negligible.
Equilibrium Distribution of Subgrid Convection: A Grand Canonic Ensemble Approach
NASA Astrophysics Data System (ADS)
Bao, J.; Penland, M. C.
2011-12-01
Moist convection on scales smaller than the horizontal grid spacing that is commonly used in operational numerical weather and climate prediction models is turbulent and therefore its interaction with the environment is stochastic. Traditionally in operational weather and climate prediction models, the effect of unresolved subgrid convection on the prediction of resolved scales is parameterized deterministically as an ensemble mean, and the stochastic fluctuations about this ensemble mean are ignored. It has recently been advocated that the stochastic fluctuations should be properly accounted for in the subgrid parameterization in order to address a persistent issue in operational ensemble prediction: the spread of ensemble members tends to be underestimated. In this study, the probability of requiring n mutually independently convective plumes and a total cloud-base mass flux M for subgrid convection to occur in a given grid box is derived based on the concept of the grand canonical ensemble, which is well known in classic statistical mechanics. The probability distribution functions of the cloud-base mass flux and the number of subgrid convective plumes are dependent on the average of each of the two quantities. For a large number of such grid boxes in a given area, the concept can be extended to a homogenous stochastic situation. In this situation, the probability of finding exact k subgrid convective plumes in one of the grid boxes is given by the binomial distribution, which converges to the Poisson distribution when the number of the boxes approaches to infinity. The latter result provides an alternative way to derive and interpret the previous theoretical results obtained by Craig and Cohen (2006, JAS, Vol. 63, p. 1996-2015).
Equilibrium distribution of the wave energy in a carbyne chain
NASA Astrophysics Data System (ADS)
Kovriguine, D. A.; Nikitenkova, S. P.
2016-03-01
The steady-state energy distribution of thermal vibrations at a given ambient temperature has been investigated based on a simple mathematical model that takes into account central and noncentral interactions between carbon atoms in a one-dimensional carbyne chain. The investigation has been performed using standard asymptotic methods of nonlinear dynamics in terms of the classical mechanics. In the first-order nonlinear approximation, there have been revealed resonant wave triads that are formed at a typical nonlinearity of the system under phase matching conditions. Each resonant triad consists of one longitudinal and two transverse vibration modes. In the general case, the chain is characterized by a superposition of similar resonant triplets of different spectral scales. It has been found that the energy equipartition of nonlinear stationary waves in the carbyne chain at a given temperature completely obeys the standard Rayleigh-Jeans law due to the proportional amplitude dispersion. The possibility of spontaneous formation of three-frequency envelope solitons in carbyne has been demonstrated. Heat in the form of such solitons can propagate in a chain of carbon atoms without diffusion, like localized waves.
Modelling spectral properties of non-equilibrium atomic hydrogen plasma
NASA Astrophysics Data System (ADS)
D'Ammando, G.; Pietanza, L. D.; Colonna, G.; Longo, S.; Capitelli, M.
2010-02-01
A model to predict the emissivity and absorption coefficient of atomic hydrogen plasma is presented in detail. Non-equilibrium plasma is studied through coupling of the model with a collisional-radiative code for the excited states population as well as with the Boltzmann equation for the electron energy distribution function.
The Approach to Equilibrium: Detailed Balance and the Master Equation
ERIC Educational Resources Information Center
Alexander, Millard H.; Hall, Gregory E.; Dagdigian, Paul J.
2011-01-01
The approach to the equilibrium (Boltzmann) distribution of populations of internal states of a molecule is governed by inelastic collisions in the gas phase and with surfaces. The set of differential equations governing the time evolution of the internal state populations is commonly called the master equation. An analytic solution to the master…
NASA Astrophysics Data System (ADS)
Gao, Haixia; Li, Ting; Xiao, Changming
2016-05-01
When a simple system is in its nonequilibrium state, it will shift to its equilibrium state. Obviously, in this process, there are a series of nonequilibrium states. With the assistance of Bayesian statistics and hyperensemble, a probable probability distribution of these nonequilibrium states can be determined by maximizing the hyperensemble entropy. It is known that the largest probability is the equilibrium state, and the far a nonequilibrium state is away from the equilibrium one, the smaller the probability will be, and the same conclusion can also be obtained in the multi-state space. Furthermore, if the probability stands for the relative time the corresponding nonequilibrium state can stay, then the velocity of a nonequilibrium state returning back to its equilibrium can also be determined through the reciprocal of the derivative of this probability. It tells us that the far away the state from the equilibrium is, the faster the returning velocity will be; if the system is near to its equilibrium state, the velocity will tend to be smaller and smaller, and finally tends to 0 when it gets the equilibrium state.
High-frequency spectral distribution of the equilibrium radiation energy in a plasma
NASA Astrophysics Data System (ADS)
Bobrov, V. B.; Trigger, S. A.
2016-04-01
We establish that the difference of the spectral distribution of the equilibrium radiation energy in matter from the Planck formula in the high-frequency range is determined by the imaginary part of the transverse dielectric permittivity of the matter. Based on this, we show that in a rarified high-temperature fully ionized nonrelativistic plasma, the high-frequency spectral distribution of the equilibrium radiation energy differs essentially from the Planck formula because of the power-law character of the decrease in the frequency, which is due to the presence of matter.
Relativistic distribution function for particles with spin at local thermodynamical equilibrium
Becattini, F.; Chandra, V.; Del Zanna, L.; Grossi, E.
2013-11-15
We present an extension of relativistic single-particle distribution function for weakly interacting particles at local thermodynamical equilibrium including spin degrees of freedom, for massive spin 1/2 particles. We infer, on the basis of the global equilibrium case, that at local thermodynamical equilibrium particles acquire a net polarization proportional to the vorticity of the inverse temperature four-vector field. The obtained formula for polarization also implies that a steady gradient of temperature entails a polarization orthogonal to particle momentum. The single-particle distribution function in momentum space extends the so-called Cooper–Frye formula to particles with spin 1/2 and allows us to predict their polarization in relativistic heavy ion collisions at the freeze-out. -- Highlights: •Single-particle distribution function in local thermodynamical equilibrium with spin. •Polarization of spin 1/2 particles in a fluid at local thermodynamical equilibrium. •Prediction of a new effect: a steady gradient of temperature induces a polarization. •Application to the calculation of polarization in relativistic heavy ion collisions.
Numerical scheme for a spatially inhomogeneous matrix-valued quantum Boltzmann equation
NASA Astrophysics Data System (ADS)
Lu, Jianfeng; Mendl, Christian B.
2015-06-01
We develop an efficient algorithm for a spatially inhomogeneous matrix-valued quantum Boltzmann equation derived from the Hubbard model. The distribution functions are 2 × 2 matrix-valued to accommodate the spin degree of freedom, and the scalar quantum Boltzmann equation is recovered as a special case when all matrices are proportional to the identity. We use Fourier discretization and fast Fourier transform to efficiently evaluate the collision kernel with spectral accuracy, and numerically investigate periodic, Dirichlet and Maxwell boundary conditions. Model simulations quantify the convergence to local and global thermal equilibrium.
A multispeed Discrete Boltzmann Model for transcritical 2D shallow water flows
NASA Astrophysics Data System (ADS)
La Rocca, Michele; Montessori, Andrea; Prestininzi, Pietro; Succi, Sauro
2015-03-01
In this work a Discrete Boltzmann Model for the solution of transcritical 2D shallow water flows is presented and validated. In order to provide the model with transcritical capabilities, a particular multispeed velocity set has been employed for the discretization of the Boltzmann equation. It is shown that this particular set naturally yields a simple and closed procedure to determine higher order equilibrium distribution functions needed to simulate transcritical flow. The model is validated through several classical benchmarks and is proven to correctly and accurately simulate both 1D and 2D transitions between the two flow regimes.
Stability analysis of lattice Boltzmann methods
Sterling, J.D.; Chen, Shiyi
1996-01-01
The lattice Boltzmann equation describes the evolution of the velocity distribution function on a lattice in a manner that macroscopic fluid dynamical behavior is recovered. Although the equation is a derivative of lattice gas automata, it may be interpreted as a Lagrangian finite-difference method for the numerical simulation of the discrete-velocity Boltzmann equation that makes use of a BGK collision operator. As a result, it is not surprising that numericaI instability of lattice Boltzmann methods have been frequently encountered by researchers. We present an analysis of the stability of perturbations of the particle populations linearized about equilibrium values corresponding to a constant-density uniform mean flow. The linear stability depends on the following parameters: the distribution of the mass at a site between the different discrete speeds, the BGK relaxation time, the mean velocity, and the wave-number of the perturbations. This parameter space is too large to compute the complete stability characteristics. We report some stability results for a subset of the parameter space for a 7-velocity hexagonal lattice, a 9-velocity square lattice, and a 15-velocity cubic lattice. Results common to all three lattices are (1) the BGK relaxation time {tau} must be greater than 1/2 corresponding to positive shear viscosity, (2) there exists a maximum stable mean velocity for fixed values of theother parameters, and (3) as {tau} is increased from 1/2 the maximum stable velocity increases monotonically until some fixed velocity is reached which does not change for larger {tau}.
Relaxation of non-equilibrium hydrogen distributions in a-Si:H
Adriaenssens, G.J.; Zhang, Q.
1998-12-31
It is shown that a reduction of the three-level energy diagram proposed by Van de Walle (Phys. Rev. B 53, 11292, 1996) to describe the relaxation of non-equilibrium hydrogen distributions, to just the interstitial transport level and a distribution of traps, allows an essentially equivalent formulation of the hydrogen relaxation kinetics. The modified formulation offers the possibility of accounting for dispersive diffusion while preserving the essential multiple retrapping aspect of the original proposal.
On return to thermal equilibrium for the excess charge distribution in semiconductors
NASA Astrophysics Data System (ADS)
Paranjape, V. V.
1985-04-01
The relaxation time for the return to thermal equilibrium for the excess charge distribution in metals and in semiconductors is shown to be longer than the relaxation time derived in many textbooks. In semiconductors the relaxation of charge is accompanied by an increase in temperature of the charged carriers. The time of decay τT for the temperature is considerably longer than the time of decay for the excess charge. Hence the overall relaxation time for attaining thermal equilibrium is determined by τT.
Equilibrium distribution of hard-sphere systems and revised Enskog theory
NASA Astrophysics Data System (ADS)
van Beijeren, H.
1983-10-01
A revised Enskog theory (RET) is shown to lead to a correct equilibrium distribution in hard-sphere systems in a stationary external potential, while the standard Enskog theory (SET) does not. Attention is given to the s-component hard-sphere mixture with constant external potential acting on particles of a particular species. The different definition of the pair correlation function at a contact point of two different species used in RET is demonstrated to be consistent with equilibrium statistical mechanics, whereas the definition chosen for the contact point in SET is not.
The Equilibrium Allele Frequency Distribution for a Population with Reproductive Skew
Der, Ricky; Plotkin, Joshua B.
2014-01-01
We study the population genetics of two neutral alleles under reversible mutation in a model that features a skewed offspring distribution, called the Λ-Fleming–Viot process. We describe the shape of the equilibrium allele frequency distribution as a function of the model parameters. We show that the mutation rates can be uniquely identified from this equilibrium distribution, but the form of the offspring distribution cannot itself always be so identified. We introduce an estimator for the mutation rate that is consistent, independent of the form of reproductive skew. We also introduce a two-allele infinite-sites version of the Λ-Fleming–Viot process, and we use it to study how reproductive skew influences standing genetic diversity in a population. We derive asymptotic formulas for the expected number of segregating sites as a function of sample size and offspring distribution. We find that the Wright–Fisher model minimizes the equilibrium genetic diversity, for a given mutation rate and variance effective population size, compared to all other Λ-processes. PMID:24473932
Radtke, Gregg A; Hadjiconstantinou, Nicolas G
2009-05-01
We present an efficient variance-reduced particle simulation technique for solving the linearized Boltzmann transport equation in the relaxation-time approximation used for phonon, electron, and radiative transport, as well as for kinetic gas flows. The variance reduction is achieved by simulating only the deviation from equilibrium. We show that in the limit of small deviation from equilibrium of interest here, the proposed formulation achieves low relative statistical uncertainty that is also independent of the magnitude of the deviation from equilibrium, in stark contrast to standard particle simulation methods. Our results demonstrate that a space-dependent equilibrium distribution improves the variance reduction achieved, especially in the collision-dominated regime where local equilibrium conditions prevail. We also show that by exploiting the physics of relaxation to equilibrium inherent in the relaxation-time approximation, a very simple collision algorithm with a clear physical interpretation can be formulated. PMID:19518597
Temperature based Restricted Boltzmann Machines
NASA Astrophysics Data System (ADS)
Li, Guoqi; Deng, Lei; Xu, Yi; Wen, Changyun; Wang, Wei; Pei, Jing; Shi, Luping
2016-01-01
Restricted Boltzmann machines (RBMs), which apply graphical models to learning probability distribution over a set of inputs, have attracted much attention recently since being proposed as building blocks of multi-layer learning systems called deep belief networks (DBNs). Note that temperature is a key factor of the Boltzmann distribution that RBMs originate from. However, none of existing schemes have considered the impact of temperature in the graphical model of DBNs. In this work, we propose temperature based restricted Boltzmann machines (TRBMs) which reveals that temperature is an essential parameter controlling the selectivity of the firing neurons in the hidden layers. We theoretically prove that the effect of temperature can be adjusted by setting the parameter of the sharpness of the logistic function in the proposed TRBMs. The performance of RBMs can be improved by adjusting the temperature parameter of TRBMs. This work provides a comprehensive insights into the deep belief networks and deep learning architectures from a physical point of view.
NASA Astrophysics Data System (ADS)
Frank, T. D.
The virial theorem and the concept of canonical-statistical distributions represent two fundamental elements of statistical physics. We apply these concepts to hand tremor oscillations recorded from six Parkinson patients. We find that the virial theorem holds for Parkinson tremor oscillations. In contrast, we find that the concept of canonical distributions fails to a certain extent and needs to be replaced by the notion of non-canonical (i.e., canonical-dissipative) distributions. In doing so, our analysis reveals both general statistical aspects and non-equilibrium aspects of Parkinson hand tremor.
Three-dimensional lattice Boltzmann model for compressible flows.
Sun, Chenghai; Hsu, Andrew T
2003-07-01
A three-dimensional compressible lattice Boltzmann model is formulated on a cubic lattice. A very large particle-velocity set is incorporated in order to enable a greater variation in the mean velocity. Meanwhile, the support set of the equilibrium distribution has only six directions. Therefore, this model can efficiently handle flows over a wide range of Mach numbers and capture shock waves. Due to the simple form of the equilibrium distribution, the fourth-order velocity tensors are not involved in the formulation. Unlike the standard lattice Boltzmann model, no special treatment is required for the homogeneity of fourth-order velocity tensors on square lattices. The Navier-Stokes equations were recovered, using the Chapman-Enskog method from the Bhatnagar-Gross-Krook (BGK) lattice Boltzmann equation. The second-order discretization error of the fluctuation velocity in the macroscopic conservation equation was eliminated by means of a modified collision invariant. The model is suitable for both viscous and inviscid compressible flows with or without shocks. Since the present scheme deals only with the equilibrium distribution that depends only on fluid density, velocity, and internal energy, boundary conditions on curved wall are easily implemented by an extrapolation of macroscopic variables. To verify the scheme for inviscid flows, we have successfully simulated a three-dimensional shock-wave propagation in a box and a normal shock of Mach number 10 over a wedge. As an application to viscous flows, we have simulated a flat plate boundary layer flow, flow over a cylinder, and a transonic flow over a NACA0012 airfoil cascade. PMID:12935242
Importance of Pore Size Distribution of Fine-grained Sediments on Gas Hydrate Equilibrium
NASA Astrophysics Data System (ADS)
Kwon, T. H.; Kim, H. S.; Cho, G. C.; Park, T. H.
2015-12-01
Gas hydrates have been considered as a new source of natural gases. For the gas hydrate production, the gas hydrate reservoir should be depressurized below the equilibrium pressure of gas hydrates. Therefore, it is important to predict the equilibrium of gas hydrates in the reservoir conditions because it can be affected by the pore size of the host sediments due to the capillary effect. In this study, gas hydrates were synthesized in fine-grained sediment samples including a pure silt sample and a natural clayey silt sample cored from a hydrate occurrence region in Ulleung Basin, East Sea, offshore Korea. Pore size distributions of the samples were obtained by the nitrogen adsorption and desorption test and the mercury intrusion porosimetry. The equilibrium curve of gas hydrates in the fine-grained sediments were found to be significantly influenced by the clay fraction and the corresponding small pores (>50 nm in diameter). For the clayey silt sample, the equilibrium pressure was higher by ~1.4 MPa than the bulk equilibrium pressure. In most cases of oceanic gas hydrate reservoirs, sandy layers are found interbedded with fine-grained sediment layers while gas hydrates are intensively accumulated in the sandy layers. Our experiment results reveal the inhibition effect of fine-grained sediments against gas hydrate formation, in which greater driving forces (e.g., higher pressure or lower temperature) are required during natural gas migration. Therefore, gas hydrate distribution in interbedded layers of sandy and fine-grained sediments can be explained by such capillary effect induced by the pore size distribution of host sediments.
High-precision work distributions for extreme non-equilibrium processes in large systems
NASA Astrophysics Data System (ADS)
Hartmann, Alexander
2014-03-01
The distributions of work for strongly non-equilibrium processes are studied using a very general form of a large-deviation approach, which allows one to study distributions down to extremely small probabilities of almost arbitrary quantities of interest for equilibrium, non-equilibrium stationary and even non-stationary processes. The method is applied to varying quickly the external field in a wide range B = 3 <--> 0 for critical (T = 2 . 269) two-dimensional Ising system of size L × L = 128 × 128 . To obtain free energy differences from the work distributions, they must be studied in ranges where the probabilities are as small as 10-240, which is not possible using direct simulation approaches. By comparison with the exact free energies, one sees that the present approach allows one to obtain the free energy with a very high relative precision of 10-4. This works well also for non-zero field, i.e., for a case where standard umbrella-sampling methods seem to be not so efficient to calculate free energies. Furthermore, for the present case it is verified that the resulting distributions of work fulfill Crooks theorem with high precision. Finally, the free energy for the Ising magnet as a function of the field strength is obtained.
The lattice Boltzmann model for the second-order Benjamin-Ono equations
NASA Astrophysics Data System (ADS)
Lai, Huilin; Ma, Changfeng
2010-04-01
In this paper, in order to extend the lattice Boltzmann method to deal with more complicated nonlinear equations, we propose a 1D lattice Boltzmann scheme with an amending function for the second-order (1 + 1)-dimensional Benjamin-Ono equation. With the Taylor expansion and the Chapman-Enskog expansion, the governing evolution equation is recovered correctly from the continuous Boltzmann equation. The equilibrium distribution function and the amending function are obtained. Numerical simulations are carried out for the 'good' Boussinesq equation and the 'bad' one to validate the proposed model. It is found that the numerical results agree well with the analytical solutions. The present model can be used to solve more kinds of nonlinear partial differential equations.
The standard map: From Boltzmann-Gibbs statistics to Tsallis statistics
NASA Astrophysics Data System (ADS)
Tirnakli, Ugur; Borges, Ernesto P.
2016-03-01
As well known, Boltzmann-Gibbs statistics is the correct way of thermostatistically approaching ergodic systems. On the other hand, nontrivial ergodicity breakdown and strong correlations typically drag the system into out-of-equilibrium states where Boltzmann-Gibbs statistics fails. For a wide class of such systems, it has been shown in recent years that the correct approach is to use Tsallis statistics instead. Here we show how the dynamics of the paradigmatic conservative (area-preserving) stan-dard map exhibits, in an exceptionally clear manner, the crossing from one statistics to the other. Our results unambiguously illustrate the domains of validity of both Boltzmann-Gibbs and Tsallis statistical distributions. Since various important physical systems from particle confinement in magnetic traps to autoionization of molecular Rydberg states, through particle dynamics in accelerators and comet dynamics, can be reduced to the standard map, our results are expected to enlighten and enable an improved interpretation of diverse experimental and observational results.
Lattice Boltzmann method for bosons and fermions and the fourth-order Hermite polynomial expansion.
Coelho, Rodrigo C V; Ilha, Anderson; Doria, Mauro M; Pereira, R M; Aibe, Valter Yoshihiko
2014-04-01
The Boltzmann equation with the Bhatnagar-Gross-Krook collision operator is considered for the Bose-Einstein and Fermi-Dirac equilibrium distribution functions. We show that the expansion of the microscopic velocity in terms of Hermite polynomials must be carried to the fourth order to correctly describe the energy equation. The viscosity and thermal coefficients, previously obtained by Yang et al. [Shi and Yang, J. Comput. Phys. 227, 9389 (2008); Yang and Hung, Phys. Rev. E 79, 056708 (2009)] through the Uehling-Uhlenbeck approach, are also derived here. Thus the construction of a lattice Boltzmann method for the quantum fluid is possible provided that the Bose-Einstein and Fermi-Dirac equilibrium distribution functions are expanded to fourth order in the Hermite polynomials. PMID:24827360
Equilibrium distribution of samarium and europium between fluoride salt melts and liquid bismuth
NASA Astrophysics Data System (ADS)
Zagnit'ko, A. V.; Ignat'ev, V. V.
2016-01-01
The extraction of samarium and europium from a melt of a molar composition 73LiF-27BeF2 into liquid bismuth with additions of lithium as a reducing agent at a temperature of 600-610°C was studied. The equilibrium distribution coefficients of samarium and europium were measured. In the metal fluoride salt melt under study, the valence of samarium and europium was shown to be equal to two.
Mélykúti, Bence; Hespanha, João P.; Khammash, Mustafa
2014-01-01
Many biochemical reaction networks are inherently multiscale in time and in the counts of participating molecular species. A standard technique to treat different time scales in the stochastic kinetics framework is averaging or quasi-steady-state analysis: it is assumed that the fast dynamics reaches its equilibrium (stationary) distribution on a time scale where the slowly varying molecular counts are unlikely to have changed. We derive analytic equilibrium distributions for various simple biochemical systems, such as enzymatic reactions and gene regulation models. These can be directly inserted into simulations of the slow time-scale dynamics. They also provide insight into the stimulus–response of these systems. An important model for which we derive the analytic equilibrium distribution is the binding of dimer transcription factors (TFs) that first have to form from monomers. This gene regulation mechanism is compared to the cases of the binding of simple monomer TFs to one gene or to multiple copies of a gene, and to the cases of the cooperative binding of two or multiple TFs to a gene. The results apply equally to ligands binding to enzyme molecules. PMID:24920118
Mélykúti, Bence; Hespanha, João P; Khammash, Mustafa
2014-08-01
Many biochemical reaction networks are inherently multiscale in time and in the counts of participating molecular species. A standard technique to treat different time scales in the stochastic kinetics framework is averaging or quasi-steady-state analysis: it is assumed that the fast dynamics reaches its equilibrium (stationary) distribution on a time scale where the slowly varying molecular counts are unlikely to have changed. We derive analytic equilibrium distributions for various simple biochemical systems, such as enzymatic reactions and gene regulation models. These can be directly inserted into simulations of the slow time-scale dynamics. They also provide insight into the stimulus-response of these systems. An important model for which we derive the analytic equilibrium distribution is the binding of dimer transcription factors (TFs) that first have to form from monomers. This gene regulation mechanism is compared to the cases of the binding of simple monomer TFs to one gene or to multiple copies of a gene, and to the cases of the cooperative binding of two or multiple TFs to a gene. The results apply equally to ligands binding to enzyme molecules. PMID:24920118
NASA Astrophysics Data System (ADS)
Zagnit'ko, A. V.; Ignat'ev, V. V.
2013-04-01
The distribution of lanthanum, neodymium, and thorium between a lithium chloride melt and liquid bismuth with additions of lithium as a reducing agent are investigated at 650°C. Equilibrium values of their distribution constants are measured. It is shown that in contrast to neodymium and lanthanum, thorium cannot be extracted from bismuth into lithium chloride. This allows us to propose an efficient scheme for separating lanthanides and thorium in a system for the extraction of fuel salts in molten-salt nuclear reactors.
Boltzmann's Approach to Statistical Mechanics
NASA Astrophysics Data System (ADS)
Goldstein, Sheldon
In the last quarter of the nineteenth century, Ludwig Boltzmann explained how irreversible macroscopic laws, in particular the second law of thermodynamics, originate in the time-reversible laws of microscopic physics. Boltzmann's analysis, the essence of which I shall review here, is basically correct. The most famous criticisms of Boltzmann's later work on the subject have little merit. Most twentieth century innovations - such as the identification of the state of a physical system with a probability distribution \\varrho on its phase space, of its thermodynamic entropy with the Gibbs entropy of \\varrho, and the invocation of the notions of ergodicity and mixing for the justification of the foundations of statistical mechanics - are thoroughly misguided.
Zheng, Xiliang; Wang, Jin
2015-04-01
We uncovered the universal statistical laws for the biomolecular recognition/binding process. We quantified the statistical energy landscapes for binding, from which we can characterize the distributions of the binding free energy (affinity), the equilibrium constants, the kinetics and the specificity by exploring the different ligands binding with a particular receptor. The results of the analytical studies are confirmed by the microscopic flexible docking simulations. The distribution of binding affinity is Gaussian around the mean and becomes exponential near the tail. The equilibrium constants of the binding follow a log-normal distribution around the mean and a power law distribution in the tail. The intrinsic specificity for biomolecular recognition measures the degree of discrimination of native versus non-native binding and the optimization of which becomes the maximization of the ratio of the free energy gap between the native state and the average of non-native states versus the roughness measured by the variance of the free energy landscape around its mean. The intrinsic specificity obeys a Gaussian distribution near the mean and an exponential distribution near the tail. Furthermore, the kinetics of binding follows a log-normal distribution near the mean and a power law distribution at the tail. Our study provides new insights into the statistical nature of thermodynamics, kinetics and function from different ligands binding with a specific receptor or equivalently specific ligand binding with different receptors. The elucidation of distributions of the kinetics and free energy has guiding roles in studying biomolecular recognition and function through small-molecule evolution and chemical genetics. PMID:25885453
Zheng, Xiliang; Wang, Jin
2015-01-01
We uncovered the universal statistical laws for the biomolecular recognition/binding process. We quantified the statistical energy landscapes for binding, from which we can characterize the distributions of the binding free energy (affinity), the equilibrium constants, the kinetics and the specificity by exploring the different ligands binding with a particular receptor. The results of the analytical studies are confirmed by the microscopic flexible docking simulations. The distribution of binding affinity is Gaussian around the mean and becomes exponential near the tail. The equilibrium constants of the binding follow a log-normal distribution around the mean and a power law distribution in the tail. The intrinsic specificity for biomolecular recognition measures the degree of discrimination of native versus non-native binding and the optimization of which becomes the maximization of the ratio of the free energy gap between the native state and the average of non-native states versus the roughness measured by the variance of the free energy landscape around its mean. The intrinsic specificity obeys a Gaussian distribution near the mean and an exponential distribution near the tail. Furthermore, the kinetics of binding follows a log-normal distribution near the mean and a power law distribution at the tail. Our study provides new insights into the statistical nature of thermodynamics, kinetics and function from different ligands binding with a specific receptor or equivalently specific ligand binding with different receptors. The elucidation of distributions of the kinetics and free energy has guiding roles in studying biomolecular recognition and function through small-molecule evolution and chemical genetics. PMID:25885453
NASA Astrophysics Data System (ADS)
Dzifčáková, Elena; Dudík, Jaroslav
2015-12-01
We use the latest available atomic data to calculate the ionisation and recombination rates for the non-Maxwellian n-distributions, which were shown previously to provide a good fit to the enhanced intensities of dielectronic satellite lines during solar flares. The ionisation and recombination coefficients are subsequently used to derive the ionisation equilibrium. To do so, we consider odd values of n ranging from 1 to 19, i.e., from Maxwellian to strongly non-Maxwellian cases. These calculations involve all elements with proton number up to 30, i.e., H to Zn. The n-distributions modify both the ionisation and the recombination rates. The ionisation rates decrease more steeply at lower pseudo-temperatures, while the radiative recombination rate is reduced due to a lower number of low-energy electrons. The peaks of the dielectronic recombination rates become narrower. These changes are reflected in the ionisation equilibrium. Ion abundance peaks become narrower and can also be shifted, mostly towards higher temperatures. The He-like ions are an important exception, as they are formed in a larger temperature range than that for the Maxwellian distribution. The ions Si xiii - xiv used previously for the diagnostics of the n-distributions are affected only weakly, confirming the determination of n. The ionisation equilibria are available as the electronic supplementary material in a format compatible with the CHIANTI database.
Lattice Boltzmann modeling of phonon transport
NASA Astrophysics Data System (ADS)
Guo, Yangyu; Wang, Moran
2016-06-01
A novel lattice Boltzmann scheme is proposed for phonon transport based on the phonon Boltzmann equation. Through the Chapman-Enskog expansion, the phonon lattice Boltzmann equation under the gray relaxation time approximation recovers the classical Fourier's law in the diffusive limit. The numerical parameters in the lattice Boltzmann model are therefore rigorously correlated to the bulk material properties. The new scheme does not only eliminate the fictitious phonon speed in the diagonal direction of a square lattice system in the previous lattice Boltzmann models, but also displays very robust performances in predicting both temperature and heat flux distributions consistent with analytical solutions for diverse numerical cases, including steady-state and transient, macroscale and microscale, one-dimensional and multi-dimensional phonon heat transport. This method may provide a powerful numerical tool for deep studies of nonlinear and nonlocal heat transports in nanosystems.
U.S. stock market interaction network as learned by the Boltzmann machine
NASA Astrophysics Data System (ADS)
Borysov, Stanislav S.; Roudi, Yasser; Balatsky, Alexander V.
2015-12-01
We study historical dynamics of joint equilibrium distribution of stock returns in the U.S. stock market using the Boltzmann distribution model being parametrized by external fields and pairwise couplings. Within Boltzmann learning framework for statistical inference, we analyze historical behavior of the parameters inferred using exact and approximate learning algorithms. Since the model and inference methods require use of binary variables, effect of this mapping of continuous returns to the discrete domain is studied. The presented results show that binarization preserves the correlation structure of the market. Properties of distributions of external fields and couplings as well as the market interaction network and industry sector clustering structure are studied for different historical dates and moving window sizes. We demonstrate that the observed positive heavy tail in distribution of couplings is related to the sparse clustering structure of the market. We also show that discrepancies between the model's parameters might be used as a precursor of financial instabilities.
U.S. stock market interaction network as learned by the Boltzmann machine
Borysov, Stanislav S.; Roudi, Yasser; Balatsky, Alexander V.
2015-12-07
Here, we study historical dynamics of joint equilibrium distribution of stock returns in the U.S. stock market using the Boltzmann distribution model being parametrized by external fields and pairwise couplings. Within Boltzmann learning framework for statistical inference, we analyze historical behavior of the parameters inferred using exact and approximate learning algorithms. Since the model and inference methods require use of binary variables, effect of this mapping of continuous returns to the discrete domain is studied. The presented results show that binarization preserves the correlation structure of the market. Properties of distributions of external fields and couplings as well as the market interaction network and industry sector clustering structure are studied for different historical dates and moving window sizes. We demonstrate that the observed positive heavy tail in distribution of couplings is related to the sparse clustering structure of the market. We also show that discrepancies between the model’s parameters might be used as a precursor of financial instabilities.
Effects of non-equilibrium particle distributions in deuterium-tritium burning
Michta, D; Graziani, F; Pruet, J; Luu, T
2009-08-18
We investigate the effects of non-equilibrium particle distributions resulting from rapid deuterium-tritium burning in plasmas using a Fokker-Planck code that incorporates small-angle Coulomb scattering, Brehmsstrahlung, Compton scattering, and thermal-nuclear burning. We find that in inertial confinement fusion environments, deviations away from Maxwellian distributions for either deuterium or tritium ions are small and result in 1% changes in the energy production rates. The deuterium and tritium effective temperatures are not equal, but differ by only about 2.5% near the time of peak burn rate. Simulations with high Z (Xe) dopants show that the dopant temperature closely tracks that of the fuel. On the other hand, fusion product ion distributions are highly non-Maxwellian, and careful treatments of energy-exchange between these ions and other particles is important for determining burn rates.
H to Zn Ionization Equilibrium for the Non-Maxwellian Electron κ-distributions: Updated Calculations
NASA Astrophysics Data System (ADS)
Dzifčáková, E.; Dudík, J.
2013-05-01
New data for the calculation of ionization and recombination rates have been published in the past few years, most of which are included in the CHIANTI database. We used these data to calculate collisional ionization and recombination rates for the non-Maxwellian κ-distributions with an enhanced number of particles in the high-energy tail, which have been detected in the solar transition region and the solar wind. Ionization equilibria for elements H to Zn are derived. The κ-distributions significantly influence both the ionization and recombination rates and widen the ion abundance peaks. In comparison with the Maxwellian distribution, the ion abundance peaks can also be shifted to lower or higher temperatures. The updated ionization equilibrium calculations result in large changes for several ions, notably Fe VIII-Fe XIV. The results are supplied in electronic form compatible with the CHIANTI database.
NASA Astrophysics Data System (ADS)
Di Troia, C.
2015-11-01
A class of parametric distribution functions was proposed in (Di Troia 2012 Plasma Phys. Control. Fusion 54 105017) as equilibrium distribution functions (EDFs) for charged particles in fusion plasmas, representing supra-thermal particles in anisotropic equilibria for Neutral Beam Injection and Ion Cyclotron Heating scenarios. Moreover, those EDFs can be used to represent also nearly isotropic equilibria for Slowing-Down alpha particles and core thermal plasma populations. Such EDFs depend on constants of motion (COMs). In axisymmetric system with no equilibrium electric field, they depend on toroidal canonical momentum {{P}φ} , kinetic energy w and magnetic moment μ. In the present work, the same EDFs are obtained from first principles and general hypothesis. The derivation is probabilistic and makes use of the Bayes’ Theorem. The bayesian argument is used to describe how the plasma is far from the prior probability distribution function (pdf), e.g. Maxwellian, based on the information obtained from magnetic moment and guiding center velocity pdf. Once the general functional form of the EDF has been settled, it is shown how to associate a modified Landau collision operator in the Fokker-Planck equation, to describe the system relaxation towards the proposed EDF.
Diffusion and near-equilibrium distribution of MRI and CT contrast agents in articular cartilage.
Silvast, Tuomo S; Kokkonen, Harri T; Jurvelin, Jukka S; Quinn, Thomas M; Nieminen, Miika T; Töyräs, Juha
2009-11-21
Charged contrast agents have been used both in vitro and in vivo for estimation of the fixed charge density (FCD) in articular cartilage. In the present study, the effects of molecular size and charge on the diffusion and equilibrium distribution of several magnetic resonance imaging (MRI) and computed tomography (CT) contrast agents were investigated. Full thickness cartilage disks (Ø = 4.0 mm, n = 64) were prepared from fresh bovine patellae. Contrast agent (gadopentetate: Magnevist((R)), gadodiamide: Omniscan, ioxaglate: Hexabrix or sodium iodide: NaI) diffusion was allowed either through the articular surface or through the deep cartilage. CT imaging of the samples was conducted before contrast agent administration and after 1, 5, 9, 16, 25 and 29 h (and with three samples after 2, 3, 4 and 5 days) diffusion using a clinical peripheral quantitative computed tomography (pQCT) instrument. With all contrast agents, the diffusion through the deep cartilage was slower when compared to the diffusion through the articular surface. With ioxaglate, gadopentetate and gadodiamide it took over 29 h for diffusion to reach the near-equilibrium state. The slow diffusion of the contrast agents raise concerns regarding the validity of techniques for FCD estimation, as these contrast agents may not reach the equilibrium state that is assumed. However, since cartilage composition, i.e. deep versus superficial, had a significant effect on diffusion, imaging of the nonequilibrium diffusion process might enable more accurate assessment of cartilage integrity. PMID:19864699
Equilibrium poloidal field distributions in reversed-field-pinch toroidal discharges
Baker, D.A.; Mann, L.W.; Schoenberg, K.F.
1982-04-01
A comparison between the analytic formulae of Shafranov for equilibrium in axisymmetric toroidal reversed field pinch (RFP) systems and fully toroidal numerical solutions of the Grad-Shafranov equation is presented as a function of poloidal beta, internal plasma inductance, and aspect ratio. The Shafranov formula for the equilibrium poloidal field distribution is accurate to within 5% for aspect ratios greater than 2, poloidal betas less than 50%, and for plasma current channels that exceed one-third of the minor toroidal radius. The analytic description for the center shift of the innermost flux surface that encloses the plasma current (the Shafranov shift) is accurate to within 15% for aspect ratios greater than 2 and poloidal betas below 50%, provided the shift does not exceed one-tenth of the minor conducting boundary radius. The behavior of the magnetic axis shift as a function of plasma parameters is included. The Shafranov formulae provide a convenient method for describing the equilibrium behavior of an RFP discharge. Examples illustrating the application of the analytic formulae to the Los Alamos ZT-40M RFP experiment are given.
Non-equilibrium steady-state distributions of colloids in a tilted periodic potential
NASA Astrophysics Data System (ADS)
Ma, Xiaoguang; Lai, Pik-Yin; Ackerson, Bruce; Tong, Penger
A two-layer colloidal system is constructed to study the effects of the external force F on the non-equilibrium steady-state (NESS) dynamics of the diffusing particles over a tilted periodic potential, in which detailed balance is broken due to the presence of a steady particle flux. The periodic potential is provided by the bottom layer colloidal spheres forming a fixed crystalline pattern on a glass substrate. The corrugated surface of the bottom colloidal crystal provides a gravitational potential field for the top layer diffusing particles. By tilting the sample with respect to gravity, a tangential component F is applied to the diffusing particles. The measured NESS probability density function Pss (x , y) of the particles is found to deviate from the equilibrium distribution depending on the driving or distance from equilibrium. The experimental results are compared with the exact solution of the 1D Smoluchowski equation and the numerical results of the 2D Smoluchowski equation. Moreover, from the obtained exact 1D solution, we develop an analytical method to accurately extract the 1D potential U0 (x) from the measured Pss (x) . Work supported in part by the Research Grants Council of Hong Kong SAR.
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.
Robbins, Timothy J.; Ziebarth, Jesse D.; Wang, Yongmei
2014-01-01
The ion atmosphere created by monovalent (Na+) or divalent (Mg2+) cations surrounding a B-form DNA duplex were examined using atomistic molecular dynamics (MD) simulations and the nonlinear Poisson-Boltzmann (PB) equation. The ion distributions predicted by the two methods were compared using plots of radial and two-dimensional cation concentrations and by calculating the total number of cations and net solution charge surrounding the DNA. Na+ ion distributions near the DNA were more diffuse in PB calculations than in corresponding MD simulations, with PB calculations predicting lower concentrations near DNA groove sites and phosphate groups and a higher concentration in the region between these locations. Other than this difference, the Na+ distributions generated by the two methods largely agreed, as both predicted similar locations of high Na+ concentration and nearly identical values of the number of cations and the net solution charge at all distances from the DNA. In contrast, there was greater disagreement between the two methods for Mg2+ cation concentration profiles, as both the locations and magnitudes of peaks in Mg2+ concentration were different. Despite experimental and simulation observations that Mg2+ typically maintains its first solvation shell when interacting with nucleic acids, modeling Mg2+ as an unsolvated ion during PB calculations improved the agreement of the Mg2+ ion atmosphere predicted by the two methods and allowed for values of the number of bound ions and net solution charge surrounding the DNA from PB calculations that approached the values observed in MD simulations. PMID:24443090
Equilibrium and stability in a heliotron with anisotropic hot particle slowing-down distribution
Cooper, W. A.; Asahi, Y.; Narushima, Y.; Suzuki, Y.; Watanabe, K. Y.; Graves, J. P.; Isaev, M. Yu.
2012-10-15
The equilibrium and linear fluid Magnetohydrodynamic (MHD) stability in an inward-shifted large helical device heliotron configuration are investigated with the 3D ANIMEC and TERPSICHORE codes, respectively. A modified slowing-down distribution function is invoked to study anisotropic pressure conditions. An appropriate choice of coefficients and exponents allows the simulation of neutral beam injection in which the angle of injection is varied from parallel to perpendicular. The fluid stability analysis concentrates on the application of the Johnson-Kulsrud-Weimer energy principle. The growth rates are maximum at <{beta}>{approx}2%, decrease significantly at <{beta}>{approx}4.5%, do not vary significantly with variations of the injection angle and are similar to those predicted with a bi-Maxwellian hot particle distribution function model. Stability is predicted at <{beta}>{approx}2.5% with a sufficiently peaked energetic particle pressure profile. Electrostatic potential forms from the MHD instability necessary for guiding centre orbit following are calculated.
[Equilibrium surface charge distribution in phospholipid vesicles. II. Results of calculations].
Tenchov, B G; Raĭchev, B D
1977-01-01
The results of the calculations of the equilibrium potential and surface charge distribution in a vesicle with radii r1=55 A and r2=100 A are presented. The calculations were carried out for the electrolyte concentrations 100 mM, 10mM, 1mM and 0.1 mM. The digitally obtained solutions of the Poisson-Boltzman equation for the case of spherical symmetry were utilized in the exact solution of the problem. The predictions of the exact solution differ significantly from these of the linear approximation in the range of low electrolyte concentrations. The influence of the membrane dielectric permeability on the charge distribution is negligible. PMID:588605
Non-equilibrium work distribution for interacting colloidal particles under friction
NASA Astrophysics Data System (ADS)
Ruben Gomez-Solano, Juan; July, Christoph; Mehl, Jakob; Bechinger, Clemens
2015-04-01
We experimentally investigate the non-equilibrium steady-state distribution of the work done by an external force on a mesoscopic system with many coupled degrees of freedom: a colloidal crystal mechanically driven across a commensurate periodic light field. Since this system mimics the spatiotemporal dynamics of a crystalline surface moving on a corrugated substrate, our results show general properties of the work distribution for atomically flat surfaces undergoing friction. We address the role of several parameters which can influence the shape of the work distribution, e.g. the number of particles used to locally probe the properties of the system and the time interval to measure the work. We find that, when tuning the control parameters to induce particle depinning from the substrate, there is an abrupt change of the shape of the work distribution. While in the completely static and sliding friction regimes the work distribution is Gaussian, non-Gaussian tails show up due to the spatiotemporal heterogeneity of the particle dynamics during the transition between these two regimes.
Suzuki, Hideyuki; Imura, Jun-ichi; Horio, Yoshihiko; Aihara, Kazuyuki
2013-01-01
The chaotic Boltzmann machine proposed in this paper is a chaotic pseudo-billiard system that works as a Boltzmann machine. Chaotic Boltzmann machines are shown numerically to have computing abilities comparable to conventional (stochastic) Boltzmann machines. Since no randomness is required, efficient hardware implementation is expected. Moreover, the ferromagnetic phase transition of the Ising model is shown to be characterised by the largest Lyapunov exponent of the proposed system. In general, a method to relate probabilistic models to nonlinear dynamics by derandomising Gibbs sampling is presented. PMID:23558425
Entropic lattice Boltzmann model for compressible flows.
Frapolli, N; Chikatamarla, S S; Karlin, I V
2015-12-01
We present a lattice Boltzmann model (LBM) that covers the entire range of fluid flows, from low Mach weakly compressible to transonic and supersonic flows. One of the most restrictive limitations of the lattice Boltzmann method, the low Mach number limit, is overcome here by three fundamental changes to the LBM scheme: use of an appropriately chosen multispeed lattice, accurate evaluation of the equilibrium, and the entropic relaxation for the collision. The range of applications is demonstrated through the simulation of a bow shock in front of an airfoil and the simulation of decaying compressible turbulence with shocklets. PMID:26764625
Entropic lattice Boltzmann model for compressible flows
NASA Astrophysics Data System (ADS)
Frapolli, N.; Chikatamarla, S. S.; Karlin, I. V.
2015-12-01
We present a lattice Boltzmann model (LBM) that covers the entire range of fluid flows, from low Mach weakly compressible to transonic and supersonic flows. One of the most restrictive limitations of the lattice Boltzmann method, the low Mach number limit, is overcome here by three fundamental changes to the LBM scheme: use of an appropriately chosen multispeed lattice, accurate evaluation of the equilibrium, and the entropic relaxation for the collision. The range of applications is demonstrated through the simulation of a bow shock in front of an airfoil and the simulation of decaying compressible turbulence with shocklets.
Derivation of anisotropic dissipative fluid dynamics from the Boltzmann equation
NASA Astrophysics Data System (ADS)
Molnár, Etele; Niemi, Harri; Rischke, Dirk H.
2016-06-01
Fluid-dynamical equations of motion can be derived from the Boltzmann equation in terms of an expansion around a single-particle distribution function which is in local thermodynamical equilibrium, i.e., isotropic in momentum space in the rest frame of a fluid element. However, in situations where the single-particle distribution function is highly anisotropic in momentum space, such as the initial stage of heavy-ion collisions at relativistic energies, such an expansion is bound to break down. Nevertheless, one can still derive a fluid-dynamical theory, called anisotropic dissipative fluid dynamics, in terms of an expansion around a single-particle distribution function, f^0 k, which incorporates (at least parts of) the momentum anisotropy via a suitable parametrization. We construct such an expansion in terms of polynomials in energy and momentum in the direction of the anisotropy and of irreducible tensors in the two-dimensional momentum subspace orthogonal to both the fluid velocity and the direction of the anisotropy. From the Boltzmann equation we then derive the set of equations of motion for the irreducible moments of the deviation of the single-particle distribution function from f^0 k. Truncating this set via the 14-moment approximation, we obtain the equations of motion of anisotropic dissipative fluid dynamics.
Lattice Boltzmann formulation for Braginskii magnetohydrodynamics
NASA Astrophysics Data System (ADS)
Dellar, Paul
2012-03-01
We present a lattice Boltzmann formulation of the Braginskii magnetohydrodynamic equations that describe large-scale motions in strongly magnetised plasmas. Fluid quantities, density, velocity and stress, are represented by a finite set of distribution functions associated with particles moving on a square or cubic lattice. Equilibrium distributions are constructed from Hermite moment expansions, so slowly varying solutions of the discrete kinetic equation exactly satisfy the Navier--Stokes or MHD momentum equations. Electromagnetic quantities are represented by a second kinetic equation for a set of vector-valued distribution functions. Maxwell's equations and the resistive MHD induction equation may be recovered from slowly varying solutions using different scalings. The resulting algorithm, comprising only local operations at grid points and data copying between adjacent points, readily lends itself to large-scale parallel computations. We modify the collision operator to apply different relaxation times to components of the stress parallel and perpendicular to the local magnetic field, simulating a form of the Braginskii MHD equations encountered in astrophysics. Large shears develop in simulations where the fluid velocity perpendicular to the field lines reverses.
Boltzmann kinetic equation for filtered fluid turbulence.
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
NASA Astrophysics Data System (ADS)
Bihani, A. D.; Daigle, H.; Cook, A.; Glosser, D.; Shushtarian, A.
2015-12-01
Coexistence of three methane phases (liquid (L), gas (G), hydrate (H)) in marine gas hydrate systems may occur according to in-situ pressure, temperature, salinity and pore size. In sediments with salinity close to seawater, a discrete zone of three-phase (3P) equilibrium may occur near the base of the regional hydrate stability zone (RHSZ) due to capillary effects. The existence of a 3P zone influences the location of the bottom-simulating reflection (BSR) and has implications for methane fluxes at the base of the RHSZ. We studied hydrate stability conditions in two wells, WR313-G and WR313-H, at Walker Ridge Block 313 in the northern Gulf of Mexico. We determined pore size distributions (PSD) by constructing a synthetic nuclear magnetic resonance (NMR) relaxation time distribution. Correlations were obtained by non-linear regression on NMR, gamma ray, and bulk density logs from well KC-151 at Keathley Canyon. The correlations enabled construction of relaxation time distributions for WR313-G and WR313-H, which were used to predict PSD through comparison with mercury injection capillary pressure measurements. With the computed PSD, L+H and L+G methane solubility was determined from in-situ pressure and temperature. The intersection of the L+G and L+H curves for various pore sizes allowed calculation of the depth range of the 3P equilibrium zone. As in previous studies at Blake Ridge and Hydrate Ridge, the top of the 3P zone moves upwards with increasing water depth and overlies the bulk 3P equilibrium depth. In clays at Walker Ridge, the predicted thickness of the 3P zone is approximately 35 m, but in coarse sands it is only a few meters due to the difference in absolute pore sizes and the width of the PSD. The thick 3P zone in the clays may explain in part why the BSR is only observed in the sand layers at Walker Ridge, although other factors may influence the presence or absence of a BSR.
Bihani, Abhishek; Daigle, Hugh; Cook, Ann; Glosser, Deborah; Shushtarian, Arash
2015-12-15
Coexistence of three methane phases (liquid (L), gas (G), hydrate (H)) in marine gas hydrate systems may occur according to in-situ pressure, temperature, salinity and pore size. In sediments with salinity close to seawater, a discrete zone of three-phase (3P) equilibrium may occur near the base of the regional hydrate stability zone (RHSZ) due to capillary effects. The existence of a 3P zone influences the location of the bottom-simulating reflection (BSR) and has implications for methane fluxes at the base of the RHSZ. We studied hydrate stability conditions in two wells, WR313-G and WR313-H, at Walker Ridge Block 313 in the northern Gulf of Mexico. We determined pore size distributions (PSD) by constructing a synthetic nuclear magnetic resonance (NMR) relaxation time distribution. Correlations were obtained by non-linear regression on NMR, gamma ray, and bulk density logs from well KC-151 at Keathley Canyon. The correlations enabled construction of relaxation time distributions for WR313-G and WR313-H, which were used to predict PSD through comparison with mercury injection capillary pressure measurements. With the computed PSD, L+H and L+G methane solubility was determined from in-situ pressure and temperature. The intersection of the L+G and L+H curves for various pore sizes allowed calculation of the depth range of the 3P equilibrium zone. As in previous studies at Blake Ridge and Hydrate Ridge, the top of the 3P zone moves upwards with increasing water depth and overlies the bulk 3P equilibrium depth. In clays at Walker Ridge, the predicted thickness of the 3P zone is approximately 35 m, but in coarse sands it is only a few meters due to the difference in absolute pore sizes and the width of the PSD. The thick 3P zone in the clays may explain in part why the BSR is only observed in the sand layers at Walker Ridge, although other factors may influence the presence or absence of a BSR.
Raut, L K
1991-01-01
A study is conducted in attempts to increase the understanding of the links between macroeconomic effects and causes of population growth in formulating policy. An overlapping generations general equilibrium model is employed aggregating household decisions about fertility, savings, and investment in the human capital of children with the objective of studying intertemporal relationships among population growth, income distribution, inter-generation social mobility, skill composition of the labor force, and household income. As a result of endogenous fertility, the equilibrium path attains steady state from the second generation. Income tax transfer, child taxation, and social security taxation policies are also examined in the paper. A structural explanation is given for the inverse household income-child quantity and negative child quality-quantity relationships seen in developing countries. In a Cobb-Douglas economy, these relationships hold in the short-run, potentially working over the long-run in other economies. Overall, the model shows that group interests may hinder emergence of perfect capital markets with private initiatives. Where developing countries are concerned, these results have strong implications for population policy. A policy mix of building good quality schools, or subsidizing rural education, introducing a formal social security program, and providing high-yield, risk-free investments, banking, and insurance services to the poor is recommended. PMID:12284076
Hu, Wei; Zhang, Chao; Hu, Wenjin; Fang, Yun; Hou, Wenjie
2012-10-01
Aerosol delivery is a route which is advantageous to the therapy of pulmonary diseases, such as lung cancer. The pharmacokinetics and tissue distribution after aerosol delivery of carboxylate form of hydroxycamptothecin (C-HCPT) were investigated. The concentrations of the three different types (lactone, carboxylate and the total of both forms) of HCPT were measured by HPLC analysis. The initial experiment showed no evident difference between lactone and carboxylate in the lungs during the aerosol treatment, compared with the HCPT content in plasma. The AUC(inf) value of lactone in the lungs was higher than that of carboxylate, which was 138,176.00 min ng g⁻¹ and 128,460.00 min ng g⁻¹, respectively. Meanwhile, AUC(inf) in the plasma during the entire treatment indicated that the lactone content was always at a lower level, and the carboxylate form tended to predominate, as shown by the lactone/carboxylate (L/C) equilibrium. The tissue distribution results showed that the lactone proportion in the liver increased up to the maximum value of 69.69% after aerosol administration, whereas the mean L/C equilibrium index for the liver was 2.07±1.06, and the C(max) and AUC(0-∞) values of the total HCPT were highest in the tissues. Based on these results we speculated that the initial wholly carboxylate form of the HCPT atomized liquid did not influence the transformation to lactone form. Moreover, the deposition of the total HCPT and lactone was higher in the lungs and other tissues than in the plasma after the aerosol treatment. This study will be beneficial to the therapy of pulmonary carcinoma. PMID:22858157
NASA Astrophysics Data System (ADS)
Barthes, Laurent; Mallet, Cécile
2010-05-01
Keywords: Rain Drop Size Distribution, Breakup, coalescence, disdrometer The study of the vertical evolution of raindrop size distributions (DSDs) during rainfall, from the freezing level isotherm to ground level, is a key to improving our understanding of the microphysics of rain. In numerous domains such as remote sensing, telecommunications, soil erosion, and the study of the rain's efficiency in 'washing' the atmosphere, the DSD plays an important role. Among the different processes affecting the evolution of DSD, breakup and coalescence are two of the most significant. Models of coalescence and breakup lead to equilibrium of the raindrop size distribution (DSD) after a fall through sufficient vertical height. At equilibrium, the DSD no longer evolves, and its shape is unique whatever the rain rate or LWC. This implies that the DSD is known, to within a multiplication constant. These models based on experimental measurements have been developed over the past 40 years. The Low and List (1982a,b) parameterization (hereinafter LL82) and the Greg M. McFarquhar (2004) model are both based on the same laboratory experiments, which lead to an equilibrium drop size distribution (EDSD) with two or three peaks, and an exponential tail with a slope of approximately Λ=65 cm-1. Numerous measurements using disdrometer collected in different climatic areas: Paris, France (Mars to October 2000), Iowa-City (April to October 2002), and Djougou (Benin June to September 2006) corresponding to 537 hours of rain period have shown that for high rain rates, close to a state of equilibrium, this slope lies between Λ=20 - 22 cm-1. This latter value is corroborated by others measurements found in the literature (Hu & Srivastava, 1995). Hu & Srivastava suggested that the Low and List parameterization may overestimate the effects of the breakup process. This hypothesis is in adequation with recent laboratory experiments (A.P. Barros 2008) in which the authors conclude that the number of
He, Xiaoyi; Lou, Li-Shi Lou, Li-Shi
1997-12-01
In this paper, the lattice Boltzmann equation is directly derived from the Boltzmann equation. It is shown that the lattice Boltzmann equation is a special discretized form of the Boltzmann equation. Various approximations for the discretization of the Boltzmann equation in both time and phase space are discussed in detail. A general procedure to derive the lattice Boltzmann model from the continuous Boltzmann equation is demonstrated explicitly. The lattice Boltzmann models derived include the two-dimensional 6-bit, 7-bit, and 9-bit, and three-dimensional 27-bit models. {copyright} {ital 1997} {ital The American Physical Society}
Temperature based Restricted Boltzmann Machines
Li, Guoqi; Deng, Lei; Xu, Yi; Wen, Changyun; Wang, Wei; Pei, Jing; Shi, Luping
2016-01-01
Restricted Boltzmann machines (RBMs), which apply graphical models to learning probability distribution over a set of inputs, have attracted much attention recently since being proposed as building blocks of multi-layer learning systems called deep belief networks (DBNs). Note that temperature is a key factor of the Boltzmann distribution that RBMs originate from. However, none of existing schemes have considered the impact of temperature in the graphical model of DBNs. In this work, we propose temperature based restricted Boltzmann machines (TRBMs) which reveals that temperature is an essential parameter controlling the selectivity of the firing neurons in the hidden layers. We theoretically prove that the effect of temperature can be adjusted by setting the parameter of the sharpness of the logistic function in the proposed TRBMs. The performance of RBMs can be improved by adjusting the temperature parameter of TRBMs. This work provides a comprehensive insights into the deep belief networks and deep learning architectures from a physical point of view. PMID:26758235
A lattice Boltzmann model for the Burgers-Fisher equation.
Zhang, Jianying; Yan, Guangwu
2010-06-01
A lattice Boltzmann model is developed for the one- and two-dimensional Burgers-Fisher equation based on the method of the higher-order moment of equilibrium distribution functions and a series of partial differential equations in different time scales. In order to obtain the two-dimensional Burgers-Fisher equation, vector sigma(j) has been used. And in order to overcome the drawbacks of "error rebound," a new assumption of additional distribution is presented, where two additional terms, in first order and second order separately, are used. Comparisons with the results obtained by other methods reveal that the numerical solutions obtained by the proposed method converge to exact solutions. The model under new assumption gives better results than that with second order assumption. PMID:20590325
Energy distribution and heat transfer mechanisms in atmospheric pressure non-equilibrium plasmas
NASA Astrophysics Data System (ADS)
Nozaki, Tomohiro; Miyazaki, Yu; Unno, Yasuko; Okazaki, Ken
2001-12-01
Energy distribution and heat transfer mechanisms in atmospheric pressure non-equilibrium plasmas were investigated extensively through energy balance analysis, emission spectroscopy of the rotational band of CH (A2Δ→X2Π), and gas chromatographic analysis. Two plasma sources were examined: methane-fed dielectric barrier discharge (DBD) and atmospheric pressure glow-discharge (APG). The DBD features filamentary microdischarges accompanied by surface discharge along a dielectric barrier. As a result, 60% of the input power was measured as heat transfer to the dielectric electrode, whereas 20% was to the metallic electrode. Consequently, feed gas average temperature was increased only by 20-40 K. On the other hand, rotational temperature of the corresponding emission region exceeded average gas temperature by 100 K. In APG, heat transfer to electrodes was dominated by formation of negative glow regardless of whether the electrode was covered by a dielectric. However, negative glow tended to be thinner and more intense when it formed on a metallic electrode, leading to slightly higher metallic heating. Rotational temperature in APG was close to average gas temperature since APG does not show radial localization of plasma. Energy efficiency for methane decomposition process to produce ethane, ethylene, and hydrogen was about 1% regardless of the plasma source. Energy distribution and heat transfer mechanisms depend strongly on the plasma spatial structure rather than flow fields or feed gas physical properties.
NASA Astrophysics Data System (ADS)
Hahn, M.; Savin, D. W.
2015-02-01
We describe the influence of electron-impact multiple ionization (EIMI) on the ionization balance of collisionally ionized plasmas. Previous ionization balance calculations have largely neglected EIMI. Here, EIMI cross-section data are incorporated into calculations of both equilibrium and non-equilibrium charge-state distributions (CSDs). For equilibrium CSDs, we find that EIMI has only a small effect and can usually be ignored. However, for non-equilibrium plasmas the influence of EIMI can be important. In particular, we find that for plasmas in which the temperature oscillates there are significant differences in the CSD when including versus neglecting EIMI. These results have implications for modeling and spectroscopy of impulsively heated plasmas, such as nanoflare heating of the solar corona.
Hahn, M.; Savin, D. W.
2015-02-10
We describe the influence of electron-impact multiple ionization (EIMI) on the ionization balance of collisionally ionized plasmas. Previous ionization balance calculations have largely neglected EIMI. Here, EIMI cross-section data are incorporated into calculations of both equilibrium and non-equilibrium charge-state distributions (CSDs). For equilibrium CSDs, we find that EIMI has only a small effect and can usually be ignored. However, for non-equilibrium plasmas the influence of EIMI can be important. In particular, we find that for plasmas in which the temperature oscillates there are significant differences in the CSD when including versus neglecting EIMI. These results have implications for modeling and spectroscopy of impulsively heated plasmas, such as nanoflare heating of the solar corona.
Multiple anisotropic collisions for advection-diffusion Lattice Boltzmann schemes
NASA Astrophysics Data System (ADS)
Ginzburg, Irina
2013-01-01
This paper develops a symmetrized framework for the analysis of the anisotropic advection-diffusion Lattice Boltzmann schemes. Two main approaches build the anisotropic diffusion coefficients either from the anisotropic anti-symmetric collision matrix or from the anisotropic symmetric equilibrium distribution. We combine and extend existing approaches for all commonly used velocity sets, prescribe most general equilibrium and build the diffusion and numerical-diffusion forms, then derive and compare solvability conditions, examine available anisotropy and stable velocity magnitudes in the presence of advection. Besides the deterioration of accuracy, the numerical diffusion dictates the stable velocity range. Three techniques are proposed for its elimination: (i) velocity-dependent relaxation entries; (ii) their combination with the coordinate-link equilibrium correction; and (iii) equilibrium correction for all links. Two first techniques are also available for the minimal (coordinate) velocity sets. Even then, the two-relaxation-times model with the isotropic rates often gains in effective stability and accuracy. The key point is that the symmetric collision mode does not modify the modeled diffusion tensor but it controls the effective accuracy and stability, via eigenvalue combinations of the opposite parity eigenmodes. We propose to reduce the eigenvalue spectrum by properly combining different anisotropic collision elements. The stability role of the symmetric, multiple-relaxation-times component, is further investigated with the exact von Neumann stability analysis developed in diffusion-dominant limit.
Nicholls, David C.; Dopita, Michael A.; Sutherland, Ralph S.; Kewley, Lisa J.; Palay, Ethan
2013-08-15
In this paper we develop tools for observers to use when analyzing nebular spectra for temperatures and metallicities, with two goals: to present a new, simple method to calculate equilibrium electron temperatures for collisionally excited line flux ratios, using the latest atomic data; and to adapt current methods to include the effects of possible non-equilibrium ''{kappa}'' electron energy distributions. Adopting recent collision strength data for [O III], [S III], [O II], [S II], and [N II], we find that existing methods based on older atomic data seriously overestimate the electron temperatures, even when considering purely Maxwellian statistics. If {kappa} distributions exist in H II regions and planetary nebulae as they do in solar system plasmas, it is important to investigate the observational consequences. This paper continues our previous work on the {kappa} distribution. We present simple formulaic methods that allow observers to (1) measure equilibrium electron temperatures and atomic abundances using the latest atomic data, and (2) to apply simple corrections to existing equilibrium analysis techniques to allow for possible non-equilibrium effects. These tools should lead to better consistency in temperature and abundance measurements, and a clearer understanding of the physics of H II regions and planetary nebulae.
Tong, Weijun; Song, Haiqing; Gao, Changyou; Möhwald, Helmuth
2006-07-01
The effects of ionic strength and solvent polarity on the equilibrium distribution of fluorescein (FL) and FITC-dextran between the interior of polyelectrolyte multilayer microcapsules filled with negatively charged strong polyelectrolyte and the bulk solution were systematically investigated. A negatively charged strong polyelectrolyte, poly(styrene sulfonate) (PSS), used for CaCO3 core fabrication, was entrapped inside the capsules. Due to the semipermeability of the capsule wall, a Donnan equilibrium between the inner solution within the capsules and the bulk solution was created. The equilibrium distribution of the negatively charged permeants was investigated by means of confocal laser scanning microscopy as a function of ionic strength and solvent polarity. The equilibrium distribution of the negatively charged permeants could be tuned by increasing the bulk ionic strength to decrease the Donnan potential. Decreasing the solvent polarity also could enhance the permeation of FL, which induces a sudden increase of permeation when the ethanol volume fraction was higher than 0.7. This is mainly attributed to the precipitation of PSS. A theoretical model combining the Donnan equilibrium and Manning counterion condensation was employed to discuss the results. PMID:16805590
Theory for non-equilibrium statistical mechanics.
Attard, Phil
2006-08-21
This paper reviews a new theory for non-equilibrium statistical mechanics. This gives the non-equilibrium analogue of the Boltzmann probability distribution, and the generalization of entropy to dynamic states. It is shown that this so-called second entropy is maximized in the steady state, in contrast to the rate of production of the conventional entropy, which is not an extremum. The relationships of the new theory to Onsager's regression hypothesis, Prigogine's minimal entropy production theorem, the Langevin equation, the formula of Green and Kubo, the Kawasaki distribution, and the non-equilibrium fluctuation and work theorems, are discussed. The theory is worked through in full detail for the case of steady heat flow down an imposed temperature gradient. A Monte Carlo algorithm based upon the steady state probability density is summarized, and results for the thermal conductivity of a Lennard-Jones fluid are shown to be in agreement with known values. Also discussed is the generalization to non-equilibrium mechanical work, and to non-equilibrium quantum statistical mechanics. As examples of the new theory two general applications are briefly explored: a non-equilibrium version of the second law of thermodynamics, and the origin and evolution of life. PMID:16883388
Student understanding of the Boltzmann factor
NASA Astrophysics Data System (ADS)
Smith, Trevor I.; Mountcastle, Donald B.; Thompson, John R.
2015-12-01
[This paper is part of the Focused Collection on Upper Division Physics Courses.] We present results of our investigation into student understanding of the physical significance and utility of the Boltzmann factor in several simple models. We identify various justifications, both correct and incorrect, that students use when answering written questions that require application of the Boltzmann factor. Results from written data as well as teaching interviews suggest that many students can neither recognize situations in which the Boltzmann factor is applicable nor articulate the physical significance of the Boltzmann factor as an expression for multiplicity, a fundamental quantity of statistical mechanics. The specific student difficulties seen in the written data led us to develop a guided-inquiry tutorial activity, centered around the derivation of the Boltzmann factor, for use in undergraduate statistical mechanics courses. We report on the development process of our tutorial, including data from teaching interviews and classroom observations of student discussions about the Boltzmann factor and its derivation during the tutorial development process. This additional information informed modifications that improved students' abilities to complete the tutorial during the allowed class time without sacrificing the effectiveness as we have measured it. These data also show an increase in students' appreciation of the origin and significance of the Boltzmann factor during the student discussions. Our findings provide evidence that working in groups to better understand the physical origins of the canonical probability distribution helps students gain a better understanding of when the Boltzmann factor is applicable and how to use it appropriately in answering relevant questions.
Force-field functor theory: classical force-fields which reproduce equilibrium quantum distributions
Babbush, Ryan; Parkhill, John; Aspuru-Guzik, Alán
2013-01-01
Feynman and Hibbs were the first to variationally determine an effective potential whose associated classical canonical ensemble approximates the exact quantum partition function. We examine the existence of a map between the local potential and an effective classical potential which matches the exact quantum equilibrium density and partition function. The usefulness of such a mapping rests in its ability to readily improve Born-Oppenheimer potentials for use with classical sampling. We show that such a map is unique and must exist. To explore the feasibility of using this result to improve classical molecular mechanics, we numerically produce a map from a library of randomly generated one-dimensional potential/effective potential pairs then evaluate its performance on independent test problems. We also apply the map to simulate liquid para-hydrogen, finding that the resulting radial pair distribution functions agree well with path integral Monte Carlo simulations. The surprising accessibility and transferability of the technique suggest a quantitative route to adapting Born-Oppenheimer potentials, with a motivation similar in spirit to the powerful ideas and approximations of density functional theory. PMID:24790954
Non-equilibrium Statistical Mechanics and the Sea Ice Thickness Distribution
NASA Astrophysics Data System (ADS)
Wettlaufer, John; Toppaladoddi, Srikanth
We use concepts from non-equilibrium statistical physics to transform the original evolution equation for the sea ice thickness distribution g (h) due to Thorndike et al., (1975) into a Fokker-Planck like conservation law. The steady solution is g (h) = calN (q) hqe - h / H , where q and H are expressible in terms of moments over the transition probabilities between thickness categories. The solution exhibits the functional form used in observational fits and shows that for h << 1 , g (h) is controlled by both thermodynamics and mechanics, whereas for h >> 1 only mechanics controls g (h) . Finally, we derive the underlying Langevin equation governing the dynamics of the ice thickness h, from which we predict the observed g (h) . This allows us to demonstrate that the ice thickness field is ergodic. The genericity of our approach provides a framework for studying the geophysical scale structure of the ice pack using methods of broad relevance in statistical mechanics. Swedish Research Council Grant No. 638-2013-9243, NASA Grant NNH13ZDA001N-CRYO and the National Science Foundation and the Office of Naval Research under OCE-1332750 for support.
Force-Field Functor Theory: Classical Force-Fields which Reproduce Equilibrium Quantum Distributions
NASA Astrophysics Data System (ADS)
Babbush, Ryan; Parkhill, John; Aspuru-Guzik, Alan
2013-10-01
Feynman and Hibbs were the first to variationally determine an effective potential whose associated classical canonical ensemble approximates the exact quantum partition function. We examine the existence of a map between the local potential and an effective classical potential which matches the exact quantum equilibrium density and partition function. The usefulness of such a mapping rests in its ability to readily improve Born-Oppenheimer potentials for use with classical sampling. We show that such a map is unique and must exist. To explore the feasibility of using this result to improve classical molecular mechanics, we numerically produce a map from a library of randomly generated one-dimensional potential/effective potential pairs then evaluate its performance on independent test problems. We also apply the map to simulate liquid para-hydrogen, finding that the resulting radial pair distribution functions agree well with path integral Monte Carlo simulations. The surprising accessibility and transferability of the technique suggest a quantitative route to adapting Born-Oppenheimer potentials, with a motivation similar in spirit to the powerful ideas and approximations of density functional theory.
Moment Closure Approximations of the Boltzmann Equation Based on \\varphi -Divergences
NASA Astrophysics Data System (ADS)
Abdelmalik, M. R. A.; van Brummelen, E. H.
2016-07-01
This paper is concerned with approximations of the Boltzmann equation based on the method of moments. We propose a generalization of the setting of the moment-closure problem from relative entropy to \\varphi -divergences and a corresponding closure procedure based on minimization of \\varphi -divergences. The proposed description encapsulates as special cases Grad's classical closure based on expansion in Hermite polynomials and Levermore's entropy-based closure. We establish that the generalization to divergence-based closures enables the construction of extended thermodynamic theories that avoid essential limitations of the standard moment-closure formulations such as inadmissibility of the approximate phase-space distribution, potential loss of hyperbolicity and singularity of flux functions at local equilibrium. The divergence-based closure leads to a hierarchy of tractable symmetric hyperbolic systems that retain the fundamental structural properties of the Boltzmann equation.
Wang, Chen; Goss, Kai-Uwe; Lei, Ying Duan; Abbatt, Jonathan P D; Wania, Frank
2015-07-21
Challenges in the parametrization of compound distribution between the gas and particle phase contribute significantly to the uncertainty in the prediction of secondary organic aerosol (SOA) formation and are rooted in the complexity and variability of atmospheric condensed matter, which includes water, salts, and a multitude of organic oxidation products, often in two separated phases. Here, we explore the use of the commercial quantum-chemistry-based software COSMOtherm to predict equilibrium partitioning and Setchenow coefficients of a suite of oxidation products of α-pinene ozonolysis in an aerosol that is assumed to separate into an organic-enriched phase and an electrolyte-enriched aqueous phase. The predicted coefficients are used to estimate the phase distribution of the organic compounds, water and ammonium sulfate, the resulting phase composition, and the SOA yield. Four scenarios that differ in terms of organic loading, liquid water content, and chemical aging are compared. The organic compounds partition preferentially to the organic phase rather than the aqueous phase for the studied aerosol scenarios, partially due to the salting-out effect. Extremely low volatile organic compounds are predicted to be the dominant species in the organic aerosols at low loadings and an important component at higher loadings. The highest concentration of oxidation products in the condensed phase is predicted for a scenario assuming the presence of non-phase-separated cloud droplets. Partitioning into an organic aerosol phase composed of the oxidation products is predicted to be similar to partitioning into a phase composed of a single organic surrogate molecule, suggesting that the calculation procedure can be simplified without major loss of accuracy. COSMOtherm is shown to produce results that are comparable to those obtained using group contribution methods. COSMOtherm is likely to have a much larger application domain than those group contribution methods because
NASA Technical Reports Server (NTRS)
Grams, G. W.; SHARDANAND
1972-01-01
The inherent errors of applying terrestrial atmospheric ozone distribution studies to the atmosphere of other planets are discussed. Limitations associated with some of the earlier treatments of photochemical equilibrium distributions of ozone in planetary atmospheres are described. A technique having more universal application is presented. Ozone concentration profiles for the Martian atmosphere based on the results of the Mariner 4 radio occultation experiment and the more recent results with Mariner 6 and Mariner 7 have been calculated using this approach.
Lattice Boltzmann model for generalized nonlinear wave equations
NASA Astrophysics Data System (ADS)
Lai, Huilin; Ma, Changfeng
2011-10-01
In this paper, a lattice Boltzmann model is developed to solve a class of the nonlinear wave equations. Through selecting equilibrium distribution function and an amending function properly, the governing evolution equation can be recovered correctly according to our proposed scheme, in which the Chapman-Enskog expansion is employed. We validate the algorithm on some problems where analytic solutions are available, including the second-order telegraph equation, the nonlinear Klein-Gordon equation, and the damped, driven sine-Gordon equation. It is found that the numerical results agree well with the analytic solutions, which indicates that the present algorithm is very effective and can be used to solve more general nonlinear problems.
Gauss Quadratures - the Keystone of Lattice Boltzmann Models
NASA Astrophysics Data System (ADS)
Piaud, Benjamin; Blanco, Stéphane; Fournier, Richard; Ambruş, Victor Eugen; Sofonea, Victor
2014-01-01
In this paper, we compare two families of Lattice Boltzmann (LB) models derived by means of Gauss quadratures in the momentum space. The first one is the HLB(N;Qx,Qy,Qz) family, derived by using the Cartesian coordinate system and the Gauss-Hermite quadrature. The second one is the SLB(N;K,L,M) family, derived by using the spherical coordinate system and the Gauss-Laguerre, as well as the Gauss-Legendre quadratures. These models order themselves according to the maximum order N of the moments of the equilibrium distribution function that are exactly recovered. Microfluidics effects (slip velocity, temperature jump, as well as the longitudinal heat flux that is not driven by a temperature gradient) are accurately captured during the simulation of Couette flow for Knudsen number (kn) up to 0.25.
A Boltzmann treatment for the vorton excess problem
Peter, Patrick; Ringeval, Christophe E-mail: christophe.ringeval@uclouvain.be
2013-05-01
We derive and solve a Boltzmann equation governing the cosmological evolution of the number density of current carrying cosmic string loops, whose centrifugally supported equilibrium configurations are also referred to as vortons. The phase space is three-dimensional and consists of the time variable, the loop size, and a conserved quantum number. Our approach includes gravitational wave emission, a possibly finite lifetime for the vortons and works with any initial loop distribution and for any loop production function. We then show how our results generalize previous approaches on the vorton excess problem by tracking down the time evolution of the various sub-populations of current-carrying loops in a string network.
Local non-equilibrium thermodynamics
Jinwoo, Lee; Tanaka, Hajime
2015-01-01
Local Shannon entropy lies at the heart of modern thermodynamics, with much discussion of trajectory-dependent entropy production. When taken at both boundaries of a process in phase space, it reproduces the second law of thermodynamics over a finite time interval for small scale systems. However, given that entropy is an ensemble property, it has never been clear how one can assign such a quantity locally. Given such a fundamental omission in our knowledge, we construct a new ensemble composed of trajectories reaching an individual microstate, and show that locally defined entropy, information, and free energy are properties of the ensemble, or trajectory-independent true thermodynamic potentials. We find that the Boltzmann-Gibbs distribution and Landauer's principle can be generalized naturally as properties of the ensemble, and that trajectory-free state functions of the ensemble govern the exact mechanism of non-equilibrium relaxation. PMID:25592077
Local non-equilibrium thermodynamics.
Jinwoo, Lee; Tanaka, Hajime
2015-01-01
Local Shannon entropy lies at the heart of modern thermodynamics, with much discussion of trajectory-dependent entropy production. When taken at both boundaries of a process in phase space, it reproduces the second law of thermodynamics over a finite time interval for small scale systems. However, given that entropy is an ensemble property, it has never been clear how one can assign such a quantity locally. Given such a fundamental omission in our knowledge, we construct a new ensemble composed of trajectories reaching an individual microstate, and show that locally defined entropy, information, and free energy are properties of the ensemble, or trajectory-independent true thermodynamic potentials. We find that the Boltzmann-Gibbs distribution and Landauer's principle can be generalized naturally as properties of the ensemble, and that trajectory-free state functions of the ensemble govern the exact mechanism of non-equilibrium relaxation. PMID:25592077
Lattice Boltzmann model for a steady radiative transfer equation.
Yi, Hong-Liang; Yao, Feng-Ju; Tan, He-Ping
2016-08-01
A complete lattice Boltzmann model (LBM) is proposed for the steady radiative transfer equation (RTE). The RTE can be regarded as a pure convection equation with a source term. To derive the expressions for the equilibrium distribution function and the relaxation time, an artificial isotropic diffusion term is introduced to form a convection-diffusion equation. When the dimensionless relaxation time has a value of 0.5, the lattice Boltzmann equation (LBE) is exactly applicable to the original steady RTE. We also perform a multiscale analysis based on the Chapman-Enskog expansion to recover the macroscopic RTE from the mesoscopic LBE. The D2Q9 model is used to solve the LBE, and the numerical results obtained by the LBM are comparable to the results obtained by other methods or analytical solutions, which demonstrates that the proposed model is highly accurate and stable in simulating multidimensional radiative transfer. In addition, we find that the convergence rate of the LBM depends on the transport properties of RTE: for diffusion-dominated RTE with a large optical thickness, the LBM shows a second-order convergence rate in space, while for convection-dominated RTE with a small optical thickness, a lower convergence rate is observed. PMID:27627417
Lattice Boltzmann Stokesian dynamics.
Ding, E J
2015-11-01
Lattice Boltzmann Stokesian dynamics (LBSD) is presented for simulation of particle suspension in Stokes flows. This method is developed from Stokesian dynamics (SD) with resistance and mobility matrices calculated using the time-independent lattice Boltzmann algorithm (TILBA). TILBA is distinguished from the traditional lattice Boltzmann method (LBM) in that a background matrix is generated prior to the calculation. The background matrix, once generated, can be reused for calculations for different scenarios, thus the computational cost for each such subsequent calculation is significantly reduced. The LBSD inherits the merits of the SD where both near- and far-field interactions are considered. It also inherits the merits of the LBM that the computational cost is almost independent of the particle shape. PMID:26651812
Lattice Boltzmann Stokesian dynamics
NASA Astrophysics Data System (ADS)
Ding, E. J.
2015-11-01
Lattice Boltzmann Stokesian dynamics (LBSD) is presented for simulation of particle suspension in Stokes flows. This method is developed from Stokesian dynamics (SD) with resistance and mobility matrices calculated using the time-independent lattice Boltzmann algorithm (TILBA). TILBA is distinguished from the traditional lattice Boltzmann method (LBM) in that a background matrix is generated prior to the calculation. The background matrix, once generated, can be reused for calculations for different scenarios, thus the computational cost for each such subsequent calculation is significantly reduced. The LBSD inherits the merits of the SD where both near- and far-field interactions are considered. It also inherits the merits of the LBM that the computational cost is almost independent of the particle shape.
Boltzmann H function and entropy in the plasma sheet
NASA Astrophysics Data System (ADS)
Kaufmann, Richard L.; Paterson, William R.
2009-09-01
Boltzmann's H function was evaluated using 10 years of 1-min distribution functions. These results were used to study the long-term averaged spatial distributions of four entropy parameters. The average entropy density sa(x), where a = i for ions and a = e for electrons, increased when moving Earthward or toward the flanks. The magnitudes of these entropy changes were similar for ions and electrons. The entropy per unit flux tube Sf,a(x) decreased when moving Earthward or toward midnight. The spatial changes of sa(x) and of Sf,a(x) were attributed primarily to variations of the particle density na(x) and of the particle content of each unit flux tube Nf(x), respectively. A dimensionless parameter (S/Ncv)H,i that is proportional to the average entropy per ion increased when moving Earthward or toward midnight near the neutral sheet. The dimensionless parameter proportional to the entropy per ion that would exist in a plasma with the measured pressure and density but with a Maxwellian distribution function (S/Ncv)P,i was also calculated. Differences between (S/Ncv)P,i and (S/Ncv)H,i showed that the plasma was closer to equilibrium near the neutral sheet at x = -10 RE than at x = -28 RE. These gradients of the entropy per ion and of the deviations from equilibrium suggest that nonadiabatic processes and particle scattering are significant throughout the region studied.
NASA Astrophysics Data System (ADS)
Matsuyama, Akinobu; Aiba, Nobuyuki; Yagi, Masatoshi
2015-11-01
An axisymmetric MHD equilibrium model is studied to allow the inclusion of both beam inertia and energy spectrum for runaway electron beam. Following kinetic-MHD hybrid approach, we evaluate the RE beam current from the integrals of the RE distribution function. The distribution function is here evaluated by a relativistic guiding-center trace code ETC-Rel, where we have implemented the effects of collisions, radiations, and exponential growth into the code. Because to directly treat the Dreicer mechanism in particle simulations is time consuming, the primary RE source is modeled by a Monte-Carlo weighing scheme taking into account the instantaneous generation rate. This paper applies ETC-Rel to the parametric study of the MHD equilibrium with different RE beam parameters. Kinetic effects on the MHD equilibrium appears, e.g., as enhanced Shafranov shifts due to the inertia of highly relativistic electrons. A kinetic modification to the equilibrium becomes significant if the contribution of the beam inertia - being increased with the total electron mass of multi-MeV RE populations - becomes large enough to affect the radial force balance. This work was supported in part by MEXT KAKENHI Grant No. 23561009 and 26820404.
NASA Astrophysics Data System (ADS)
Kulchytskyy, Bohdan; Andriyash, Evgeny; Amin, Mohammed; Melko, Roger
The field of machine learning has been revolutionized by the recent improvements in the training of deep networks. Their architecture is based on a set of stacked layers of simpler modules. One of the most successful building blocks, known as a restricted Boltzmann machine, is an energetic model based on the classical Ising Hamiltonian. In our work, we investigate the benefits of quantum effects on the learning capacity of Boltzmann machines by extending its underlying Hamiltonian with a transverse field. For this purpose, we employ exact and stochastic training procedures on data sets with physical origins.
NASA Astrophysics Data System (ADS)
Biagetti, Matteo; Desjacques, Vincent; Kehagias, Alex; Racco, Davide; Riotto, Antonio
2016-04-01
Dark matter halos are the building blocks of the universe as they host galaxies and clusters. The knowledge of the clustering properties of halos is therefore essential for the understanding of the galaxy statistical properties. We derive an effective halo Boltzmann equation which can be used to describe the halo clustering statistics. In particular, we show how the halo Boltzmann equation encodes a statistically biased gravitational force which generates a bias in the peculiar velocities of virialized halos with respect to the underlying dark matter, as recently observed in N-body simulations.
Theory and lattice Boltzmann simulations of rapidly oscillating flows: Application to nanofluidics
NASA Astrophysics Data System (ADS)
Colosqui, Carlos Esteban
This dissertation investigates oscillating flows of simple fluids in a wide range of frequency and density variation. Understanding these flows is of fundamental importance for state-of-the-art engineering applications such as nano/microelectromechanical systems (N/MEMS) operating in diverse fluidic environments. A fluid weakly perturbed from thermodynamic equilibrium relaxes back to equilibrium state on a time scale tau ˜lambda/ cs. In the low frequency limit tauo << 1 a predominantly viscous behavior is accurately predicted by the Navier-Stokes equations. We demonstrate that in the high-frequency limit tauo >> 1 a novel transition to viscoelastic flow occurs. The viscoelastic flow of simple fluids is mathematically described via analytical solution of the Boltzmann-BGK equation of kinetic gas theory. Theoretically predicted phenomena are numerically simulated via kinetic-based approaches (i.e. high-order lattice Boltzmann methods) and observed in physical experiments with nanomechanical resonators operating over a wide pressure (10-1 ≤ p ≤ 103torr) and frequency range (103 ≤ o ≤ 108sec-1). The principal results of this dissertation are: (1) We state that the limit tauo >> 1 of high-frequency flows is governed by a damped-wave equation (i.e. the telegraph equation); (2) A universal function for the energy dissipation rate in the entire range of frequency and pressure variation is derived and favorably compared against experimental and numerical data; (3) The derived extended hydrodynamic description predicts a novel viscoelastic transition of simple fluids that holds fundamental importance from both scientific and technological standpoints. (4) We prove that high-order lattice Boltzmann schemes employing finite Hermite expansions of the Boltzmann distribution and proper regularization procedures can accurately represent macroscopic physics of the investigated flows. (5) Qualitative and quantitative agreement between theory, simulation and experiment
Equilibrium fluctuation theorems compatible with anomalous response
NASA Astrophysics Data System (ADS)
Velazquez, L.; Curilef, S.
2010-12-01
Previously, we have derived a generalization of the canonical fluctuation relation between heat capacity and energy fluctuations C = β2langδU2rang, which is able to describe the existence of macrostates with negative heat capacities C < 0. In this work, we extend our previous results for an equilibrium situation with several control parameters to account for the existence of states with anomalous values in other response functions. Our analysis leads to the derivation of three different equilibrium fluctuation theorems: the fundamental and the complementary fluctuation theorems, which represent the generalization of two fluctuation identities already obtained in previous works, and the associated fluctuation theorem, a result that has no counterpart in the framework of Boltzmann-Gibbs distributions. These results are applied to study the anomalous susceptibility of a ferromagnetic system, in particular, the case of the 2D Ising model.
NASA Astrophysics Data System (ADS)
Xu, Dazhi; Cao, Jianshu
2016-08-01
The concept of polaron, emerged from condense matter physics, describes the dynamical interaction of moving particle with its surrounding bosonic modes. This concept has been developed into a useful method to treat open quantum systems with a complete range of system-bath coupling strength. Especially, the polaron transformation approach shows its validity in the intermediate coupling regime, in which the Redfield equation or Fermi's golden rule will fail. In the polaron frame, the equilibrium distribution carried out by perturbative expansion presents a deviation from the canonical distribution, which is beyond the usual weak coupling assumption in thermodynamics. A polaron transformed Redfield equation (PTRE) not only reproduces the dissipative quantum dynamics but also provides an accurate and efficient way to calculate the non-equilibrium steady states. Applications of the PTRE approach to problems such as exciton diffusion, heat transport and light-harvesting energy transfer are presented.
Boltzmann-Electron Model in Aleph.
Hughes, Thomas Patrick; Hooper, Russell
2014-11-01
We apply the Boltzmann-electron model in the electrostatic, particle-in-cell, finite- element code Aleph to a plasma sheath. By assuming a Boltzmann energy distribution for the electrons, the model eliminates the need to resolve the electron plasma fre- quency, and avoids the numerical "grid instability" that can cause unphysical heating of electrons. This allows much larger timesteps to be used than with kinetic electrons. Ions are treated with the standard PIC algorithm. The Boltzmann-electron model re- quires solution of a nonlinear Poisson equation, for which we use an iterative Newton solver (NOX) from the Trilinos Project. Results for the spatial variation of density and voltage in the plasma sheath agree well with an analytic model
Equivalence of equations describing trace element distribution during equilibrium partial melting
NASA Technical Reports Server (NTRS)
Consolmagno, G. J.; Drake, M. J.
1976-01-01
It is shown that four equations used for calculating the evolution of trace-element abundances during equilibrium partial melting are mathematically equivalent. The equations include those of Hertogen and Gijbels (1976), Shaw (1970), Schilling (1971), and O'Nions and Clarke (1972). The general form to which all these equations reduce is presented, and an analysis is performed to demonstrate their mathematical equivalence. It is noted that the utility of the general equation flows from the nature of equilibrium (i.e., the final state is independent of the path by which that state is attained).
Endo, Satoshi
2016-08-10
Narcosis occurs as a result of the accumulation of chemicals in the phospholipid membrane. The toxic threshold concentration in the membrane is thought to be relatively constant across different chemicals and species. Hence, estimating chemical concentrations in the membrane is expected to reduce the variability of narcotic critical body residue (CBR) data. In this study, a high quality CBR dataset for three aquatic species reported recently in the literature was evaluated with the internal equilibrium distribution concept. The raw wet-weight-based CBR values were converted to membrane-weight-based CBR values by assuming that the chemical is distributed in storage lipids, membranes, proteins, and water according to the respective equilibrium partition coefficients. Several sets of partition coefficients were compared for this analysis. The results were consistent with the notion that the use of a structural protein instead of serum albumin as a surrogate for the body protein fraction could reduce the variability of CBRs. Partition coefficients predicted by polyparameter linear free energy relationships (PP-LFERs) reduced the variability of CBRs as much as or even more than experimental partition coefficients did. It is suggested that CBR data for chemicals with larger structural diversity and biological species with more distinct compositions are needed to evaluate further the equilibrium distribution concept and the constant membrane threshold hypothesis. PMID:27136717
Full Boltzmann equations for leptogenesis including scattering
Hahn-Woernle, F.; Plümacher, M.; Wong, Y.Y.Y. E-mail: pluemi@mppmu.mpg.de
2009-08-01
We study the evolution of a cosmological baryon asymmetry produced via leptogenesis by means of the full classical Boltzmann equations, without the assumption of kinetic equilibrium and including all quantum statistical factors. Beginning with the full mode equations, we derive the usual equations of motion for the right-handed neutrino number density and integrated lepton asymmetry, and show explicitly the impact of each assumption on these quantities. For the first time, we investigate also the effects of scattering of the right-handed neutrino with the top quark to leading order in the Yukawa couplings by means of the full Boltzmann equations. We find that in our full Boltzmann treatment the final lepton asymmetry can be suppressed by as much as a factor of ∼ 1.5 in the weak wash-out regime (K ∼< 1), compared to the usual integrated approach which assumes kinetic equilibrium and neglects quantum statistics. This suppression is in contrast with the enhancement seen in some previous studies that considered only decay and inverse decay of the right-handed neutrino. However, this suppression quickly decreases as we increase K. In the strong wash-out regime (K ∼> 1), the full Boltzmann treatment and the integrated approach give nearly identical final lepton asymmetries (within 10% of each other at K > 3). Finally, we show that the opposing effects of quantum statistics on decays/inverse decays and the scattering processes tend to reduce the net importance of scattering on leptogenesis in the full treatment compared to the integrated approach.
NASA Astrophysics Data System (ADS)
Tang, M.; Zhang, S.; Liu, Y.
2015-12-01
Several important equilibrium Si isotope fractionation factors among minerals, organic molecules and the H4SiO4 solution are complemented to facilitate explanation of distributions of Si isotope in the Earth's surface environments. The results reveal that heavy Si isotopes will be significantly enriched in the secondary silicate minerals in comparison to aqueous H4SiO4. On the contrary, quadra-coordinated organosilicon complexes are enriched in light silicon isotope relative to the solution. The extent of 28Si-enrichment in hyper-coordinated organosilicon complexes is found the largest. In addition, the large kinetic isotope effect associated with the polymerization of monosilicic acid and dimer is calculated and the result supports previous statement that highly 28Si-enrichment in the formation of amorphous quartz precursor contributes to the discrepancy between theoretical calculations and field observations. With equilibrium Si isotope fractionation factors provided here, Si isotope distributions in many surface systems of the Earth can be explained. For example, the change of bulk soil δ30Si can be predicted as a concave pattern with respect to weathering degree, with the minimum value where allophane completely dissolves and the total amount of sesqui-oxides and poorly crystalline minerals reaches its maximum. When well-crystallized clays start to precipitate from pore solutions under equilibrium conditions, the bulk soil δ30Si will increase again and reach a constant value. Similarly, the precipitation of crystalline smectite and the dissolution of poorly crystalline kaolinite may explain δ30Si variations in the ground water profile. Equilibrium Si isotope fractionations among quadra-coordinated organosilicon complexes and the H4SiO4 solution may also shed the light on the Si isotope distributions in Si-accumulating plants.
Boltzmann's H theorem for systems with frictional dissipation.
Bizarro, João P S
2011-03-01
By use of Boltzmann's equation to describe an ensemble of particles under the influence of a friction force, Boltzmann's H theorem is refined to explicitly include frictional dissipation, the accompanying fluctuations being modeled via an added diffusive, Fokker-Planck term. If the friction force per particle mass is proportional to velocity, as is the case with viscous drag with a friction coefficient γ, Boltzmann's H theorem for the time rate of change of the quantity H reads dH/dt ≤ γ. The classical formulation stating that H can never increase is thus replaced by the statement that H cannot increase at a rate higher than γ, a general result but of particular relevance when fluctuations are negligible and the system is far from equilibrium. When the particles are not far from thermal equilibrium, an alternative, more suitable expression emerges which can be written in the form of a Clausius inequality. PMID:21517545
Polar-coordinate lattice Boltzmann modeling of compressible flows
NASA Astrophysics Data System (ADS)
Lin, Chuandong; Xu, Aiguo; Zhang, Guangcai; Li, Yingjun; Succi, Sauro
2014-01-01
We present a polar coordinate lattice Boltzmann kinetic model for compressible flows. A method to recover the continuum distribution function from the discrete distribution function is indicated. Within the model, a hybrid scheme being similar to, but different from, the operator splitting is proposed. The temporal evolution is calculated analytically, and the convection term is solved via a modified Warming-Beam (MWB) scheme. Within the MWB scheme a suitable switch function is introduced. The current model works not only for subsonic flows but also for supersonic flows. It is validated and verified via the following well-known benchmark tests: (i) the rotational flow, (ii) the stable shock tube problem, (iii) the Richtmyer-Meshkov (RM) instability, and (iv) the Kelvin-Helmholtz instability. As an original application, we studied the nonequilibrium characteristics of the system around three kinds of interfaces, the shock wave, the rarefaction wave, and the material interface, for two specific cases. In one of the two cases, the material interface is initially perturbed, and consequently the RM instability occurs. It is found that the macroscopic effects due to deviating from thermodynamic equilibrium around the material interface differ significantly from those around the mechanical interfaces. The initial perturbation at the material interface enhances the coupling of molecular motions in different degrees of freedom. The amplitude of deviation from thermodynamic equilibrium around the shock wave is much higher than those around the rarefaction wave and material interface. By comparing each component of the high-order moments and its value in equilibrium, we can draw qualitatively the main behavior of the actual distribution function. These results deepen our understanding of the mechanical and material interfaces from a more fundamental level, which is indicative for constructing macroscopic models and other kinds of kinetic models.
The standard map: From Boltzmann-Gibbs statistics to Tsallis statistics.
Tirnakli, Ugur; Borges, Ernesto P
2016-01-01
As well known, Boltzmann-Gibbs statistics is the correct way of thermostatistically approaching ergodic systems. On the other hand, nontrivial ergodicity breakdown and strong correlations typically drag the system into out-of-equilibrium states where Boltzmann-Gibbs statistics fails. For a wide class of such systems, it has been shown in recent years that the correct approach is to use Tsallis statistics instead. Here we show how the dynamics of the paradigmatic conservative (area-preserving) stan-dard map exhibits, in an exceptionally clear manner, the crossing from one statistics to the other. Our results unambiguously illustrate the domains of validity of both Boltzmann-Gibbs and Tsallis statistical distributions. Since various important physical systems from particle confinement in magnetic traps to autoionization of molecular Rydberg states, through particle dynamics in accelerators and comet dynamics, can be reduced to the standard map, our results are expected to enlighten and enable an improved interpretation of diverse experimental and observational results. PMID:27004989
The standard map: From Boltzmann-Gibbs statistics to Tsallis statistics
Tirnakli, Ugur; Borges, Ernesto P.
2016-01-01
As well known, Boltzmann-Gibbs statistics is the correct way of thermostatistically approaching ergodic systems. On the other hand, nontrivial ergodicity breakdown and strong correlations typically drag the system into out-of-equilibrium states where Boltzmann-Gibbs statistics fails. For a wide class of such systems, it has been shown in recent years that the correct approach is to use Tsallis statistics instead. Here we show how the dynamics of the paradigmatic conservative (area-preserving) stan-dard map exhibits, in an exceptionally clear manner, the crossing from one statistics to the other. Our results unambiguously illustrate the domains of validity of both Boltzmann-Gibbs and Tsallis statistical distributions. Since various important physical systems from particle confinement in magnetic traps to autoionization of molecular Rydberg states, through particle dynamics in accelerators and comet dynamics, can be reduced to the standard map, our results are expected to enlighten and enable an improved interpretation of diverse experimental and observational results. PMID:27004989
Lattice Boltzmann morphodynamic model
NASA Astrophysics Data System (ADS)
Zhou, Jian Guo
2014-08-01
Morphological change due to sediment transport is a common natural phenomenon in real flows. It involves complex processes of erosion and deposition such as those along beaches and in river beds, imposing a strong strain on human beings. Studying and understanding morphodynamic evolution are essential to protect living environment. Although there are conventional numerical methods like finite difference method and finite volume method for forecast of morphological change by solving flow and morphodynamic equations, the methods are too complex/inefficient to be applied to a real large scale problem. To overcome this, a lattice Boltzmann method is developed to simulate morphological evolution under flows. It provides an alternative way of studying morphodynamics at the full advantages of the lattice Boltzmann methodology. The model is verified by applications to the evolution of one and two dimensional sand dunes under shallow water flows.
NASA Astrophysics Data System (ADS)
Brun-Battistini, Dominique; Mondragon-Suarez, Jose Humberto; Sandoval-Villalbazo, Alfredo; García-Perciante, Ana Laura
2015-11-01
In 1936, Richard C. Tolman showed that in thermodynamic equilibrium a temperature gradient can be compensated by a gravitational potential gradient. In reference, in a linearized gravity approximation, Tolman's law was extended for inhomogeneous non-equilibrium systems, suggesting that the contribution of the gravitational field to heat flow can be seen as a cross effect. In this work this contribution to the heat flux for a dilute simple fluid in an isotropic Schwarzschild metric is analyzed. In this case, the effect of the field is contained in the covariant derivative, such that the molecules follow geodesics. The results show that the effect of the field on the heat flux does not vanish, in contrast with what is suggested by other authors. The authors acknowledge support from CONACyT through grant CB2011/167563.
Classical non-Markovian Boltzmann equation
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.
Fermion particle production in semiclassical Boltzmann-Vlasov transport theory
Dawson, John F.; Mihaila, Bogdan; Cooper, Fred
2009-07-01
We present numerical solutions of the semiclassical Boltzmann-Vlasov equation for fermion particle-antiparticle production by strong electric fields in boost-invariant coordinates in (1+1) and (3+1) dimensional QED. We compare the Boltzmann-Vlasov results with those of recent quantum field theory calculations and find good agreement. We conclude that extending the Boltzmann-Vlasov approach to the case of QCD should allow us to do a thorough investigation of how backreaction affects recent results on the dependence of the transverse momentum distribution of quarks and antiquarks on a second Casimir invariant of color SU(3)
Equilibrium distributions and relaxation times in gaslike economic models: an analytical derivation.
Calbet, Xavier; López, José-Luis; López-Ruiz, Ricardo
2011-03-01
A step-by-step procedure to derive analytically the exact dynamical evolution equations of the probability density functions (PDFs) of well-known kinetic wealth exchange economic models is shown. This technique gives a dynamical insight into the evolution of the PDF, for example, allowing the calculation of its relaxation times. Their equilibrium PDFs can also be calculated by finding its stationary solutions. This gives as a result an integro-differential equation, which can be solved analytically in some cases and numerically in others. This should provide some guidance into the type of PDFs that can be derived from particular economic agent exchange rules or, for that matter, any other kinetic model of gases with particular collision physics. PMID:21517559
Zhao, Yong; Abreu, Eladio; Kim, Jinyong; Stadler, Guido; Eskiocak, Ugur; Terns, Michael P.; Terns, Rebecca M.; Shay, Jerry W.; Wright, Woodring E.
2011-01-01
SUMMARY Specific information about how telomerase acts in vivo is necessary for understanding telomere dynamics in human tumor cells. Our results imply that under homeostatic telomere length-maintenance conditions only one molecule of telomerase acts at each telomere during every cell division and processively adds ~60 nt to each end. In contrast, multiple molecules of telomerase act at each telomere when telomeres are elongating (non-equilibrium conditions). Telomerase extension is less processive during the first few weeks following the reversal of long-term treatment with the telomerase inhibitor GRN163L, a time when Cajal bodies fail to deliver telomerase RNA to telomeres. This result implies that processing of telomerase by Cajal bodies may affect its processivity. Overexpressed telomerase is also less processive than the endogenously expressed telomerase. These findings reveal two major distinct extension modes adopted by telomerase in vivo. PMID:21549308
Equilibrium distributions and relaxation times in gaslike economic models: An analytical derivation
NASA Astrophysics Data System (ADS)
Calbet, Xavier; López, José-Luis; López-Ruiz, Ricardo
2011-03-01
A step-by-step procedure to derive analytically the exact dynamical evolution equations of the probability density functions (PDFs) of well-known kinetic wealth exchange economic models is shown. This technique gives a dynamical insight into the evolution of the PDF, for example, allowing the calculation of its relaxation times. Their equilibrium PDFs can also be calculated by finding its stationary solutions. This gives as a result an integro-differential equation, which can be solved analytically in some cases and numerically in others. This should provide some guidance into the type of PDFs that can be derived from particular economic agent exchange rules or, for that matter, any other kinetic model of gases with particular collision physics.
NASA Astrophysics Data System (ADS)
Wiegel, Frederik W.; Perelson, Alan S.
1982-12-01
When placed in suspension red blood cells adhere face-to-face and form long, cylindrical, and sometimes branched structures called rouleaux. We use methods developed in statistical mechanics to compute various statistical properties describing the size and shape of rouleaux in thermodynamic equilibrium. This leads to analytical expressions for (1) the average number of rouleaux consisting of n cells and having m branch points; (2) the average number of cells per rouleau; (3) the average number of branch points per rouleau; and (4) the number of rouleaux with n cells in a system containing a total of N cells. We also derive asymptotic formulas that simplify these analytic expressions, and present numerical comparisons of the exact and asymptotic results.
Accurate deterministic solutions for the classic Boltzmann shock profile
NASA Astrophysics Data System (ADS)
Yue, Yubei
The Boltzmann equation or Boltzmann transport equation is a classical kinetic equation devised by Ludwig Boltzmann in 1872. It is regarded as a fundamental law in rarefied gas dynamics. Rather than using macroscopic quantities such as density, temperature, and pressure to describe the underlying physics, the Boltzmann equation uses a distribution function in phase space to describe the physical system, and all the macroscopic quantities are weighted averages of the distribution function. The information contained in the Boltzmann equation is surprisingly rich, and the Euler and Navier-Stokes equations of fluid dynamics can be derived from it using series expansions. Moreover, the Boltzmann equation can reach regimes far from the capabilities of fluid dynamical equations, such as the realm of rarefied gases---the topic of this thesis. Although the Boltzmann equation is very powerful, it is extremely difficult to solve in most situations. Thus the only hope is to solve it numerically. But soon one finds that even a numerical simulation of the equation is extremely difficult, due to both the complex and high-dimensional integral in the collision operator, and the hyperbolic phase-space advection terms. For this reason, until few years ago most numerical simulations had to rely on Monte Carlo techniques. In this thesis I will present a new and robust numerical scheme to compute direct deterministic solutions of the Boltzmann equation, and I will use it to explore some classical gas-dynamical problems. In particular, I will study in detail one of the most famous and intrinsically nonlinear problems in rarefied gas dynamics, namely the accurate determination of the Boltzmann shock profile for a gas of hard spheres.
On the Solution of a Boltzmann System for Gas Mixtures
NASA Astrophysics Data System (ADS)
Sotirov, Alexander; Yu, Shih-Hsien
2010-02-01
We study the Boltzmann equation for a mixture of two gases in one space dimension with initial condition of one gas near vacuum and the other near a Maxwellian equilibrium state. A qualitative-quantitative mathematical analysis is developed to study this mass diffusion problem based on the Green’s function of the Boltzmann equation for the single species hard sphere collision model in Liu andYu (Commun Pure Appl Math 57:1543-1608, 2004). The cross-species resonance of the mass diffusion and the diffusion-sound wave is investigated. An exponentially sharp global solution is obtained.
Spherical Harmonic Expansion Method for Coupled Electron-Phonon Boltzmann Transport
NASA Astrophysics Data System (ADS)
Santia, Marco; Albrecht, John
2014-03-01
Thermoelectric transport modeling often relies on independent Boltzmann transport equations (BTEs) for electrons and phonons which work best near equilibrium (linearized) and steady-state. Device design relies heavily on this baseline approximation. Monte Carlo methods can allow for complex physical interactions (e.g., anharmonicity) but their stochastic nature has practical limits. Distribution functions with wide disparities in population (e.g., ratios >108 between majority and minority carriers.[1]) are a computational challenge. We present a coupled BTE solver based on a k-space spherical harmonic expansion (SHE) of the distribution functions and eigenstates of electrons and phonons. The method is deterministic and allows for detailed treatments of scattering processes, yet ameliorates the issues with population disparity within phase space. We set the formalism and examine the accuracy of the SHE for phonon band structures, calculate scattering rates determined within that representation, and compare our preliminary results for distribution statistics in control examples such as thermal conductivity and drift velocity.
U.S. stock market interaction network as learned by the Boltzmann machine
Borysov, Stanislav S.; Roudi, Yasser; Balatsky, Alexander V.
2015-12-07
Here, we study historical dynamics of joint equilibrium distribution of stock returns in the U.S. stock market using the Boltzmann distribution model being parametrized by external fields and pairwise couplings. Within Boltzmann learning framework for statistical inference, we analyze historical behavior of the parameters inferred using exact and approximate learning algorithms. Since the model and inference methods require use of binary variables, effect of this mapping of continuous returns to the discrete domain is studied. The presented results show that binarization preserves the correlation structure of the market. Properties of distributions of external fields and couplings as well as themore » market interaction network and industry sector clustering structure are studied for different historical dates and moving window sizes. We demonstrate that the observed positive heavy tail in distribution of couplings is related to the sparse clustering structure of the market. We also show that discrepancies between the model’s parameters might be used as a precursor of financial instabilities.« less
[Welding arc temperature field measurements based on Boltzmann spectrometry].
Si, Hong; Hua, Xue-Ming; Zhang, Wang; Li, Fang; Xiao, Xiao
2012-09-01
Arc plasma, as non-uniform plasma, has complicated energy and mass transport processes in its internal, so plasma temperature measurement is of great significance. Compared with absolute spectral line intensity method and standard temperature method, Boltzmann plot measuring is more accurate and convenient. Based on the Boltzmann theory, the present paper calculates the temperature distribution of the plasma and analyzes the principle of lines selection by real time scanning the space of the TIG are measurements. PMID:23240385
Boltzmann Solver with Adaptive Mesh in Velocity Space
Kolobov, Vladimir I.; Arslanbekov, Robert R.; Frolova, Anna A.
2011-05-20
We describe the implementation of direct Boltzmann solver with Adaptive Mesh in Velocity Space (AMVS) using quad/octree data structure. The benefits of the AMVS technique are demonstrated for the charged particle transport in weakly ionized plasmas where the collision integral is linear. We also describe the implementation of AMVS for the nonlinear Boltzmann collision integral. Test computations demonstrate both advantages and deficiencies of the current method for calculations of narrow-kernel distributions.
Lund, Steven M.; Friedman, Alex; Bazouin, Guillaume
2011-01-10
A one-dimensional Vlasov-Poisson model for sheet beams is reviewed and extended to provide a simple framework for analysis of space-charge effects. Centroid and rms envelope equations including image charge effects are derived and reasonable parameter equivalences with commonly employed 2D transverse models of unbunched beams are established. This sheet beam model is then applied to analyze several problems of fundamental interest. A sheet beam thermal equilibrium distribution in a continuous focusing channel is constructed and shown to have analogous properties to two- d three-dimensional thermal equilibrium models in terms of the equilibrium structure and Deybe screening properties. The simpler formulation for sheet beams is exploited to explicitly calculate the distribution of particle oscillation frequencies within a thermal equilibrium beam. It is shown that as space-charge intensity increases, the frequency distribution becomes broad, suggesting that beams with strong space-charge can have improved stability.
Bean, William T.; Stafford, Robert; Butterfield, H. Scott; Brashares, Justin S.
2014-01-01
Species distributions are known to be limited by biotic and abiotic factors at multiple temporal and spatial scales. Species distribution models, however, frequently assume a population at equilibrium in both time and space. Studies of habitat selection have repeatedly shown the difficulty of estimating resource selection if the scale or extent of analysis is incorrect. Here, we present a multi-step approach to estimate the realized and potential distribution of the endangered giant kangaroo rat. First, we estimate the potential distribution by modeling suitability at a range-wide scale using static bioclimatic variables. We then examine annual changes in extent at a population-level. We define “available” habitat based on the total suitable potential distribution at the range-wide scale. Then, within the available habitat, model changes in population extent driven by multiple measures of resource availability. By modeling distributions for a population with robust estimates of population extent through time, and ecologically relevant predictor variables, we improved the predictive ability of SDMs, as well as revealed an unanticipated relationship between population extent and precipitation at multiple scales. At a range-wide scale, the best model indicated the giant kangaroo rat was limited to areas that received little to no precipitation in the summer months. In contrast, the best model for shorter time scales showed a positive relation with resource abundance, driven by precipitation, in the current and previous year. These results suggest that the distribution of the giant kangaroo rat was limited to the wettest parts of the drier areas within the study region. This multi-step approach reinforces the differing relationship species may have with environmental variables at different scales, provides a novel method for defining “available” habitat in habitat selection studies, and suggests a way to create distribution models at spatial and temporal scales
Bean, William T; Stafford, Robert; Butterfield, H Scott; Brashares, Justin S
2014-01-01
Species distributions are known to be limited by biotic and abiotic factors at multiple temporal and spatial scales. Species distribution models, however, frequently assume a population at equilibrium in both time and space. Studies of habitat selection have repeatedly shown the difficulty of estimating resource selection if the scale or extent of analysis is incorrect. Here, we present a multi-step approach to estimate the realized and potential distribution of the endangered giant kangaroo rat. First, we estimate the potential distribution by modeling suitability at a range-wide scale using static bioclimatic variables. We then examine annual changes in extent at a population-level. We define "available" habitat based on the total suitable potential distribution at the range-wide scale. Then, within the available habitat, model changes in population extent driven by multiple measures of resource availability. By modeling distributions for a population with robust estimates of population extent through time, and ecologically relevant predictor variables, we improved the predictive ability of SDMs, as well as revealed an unanticipated relationship between population extent and precipitation at multiple scales. At a range-wide scale, the best model indicated the giant kangaroo rat was limited to areas that received little to no precipitation in the summer months. In contrast, the best model for shorter time scales showed a positive relation with resource abundance, driven by precipitation, in the current and previous year. These results suggest that the distribution of the giant kangaroo rat was limited to the wettest parts of the drier areas within the study region. This multi-step approach reinforces the differing relationship species may have with environmental variables at different scales, provides a novel method for defining "available" habitat in habitat selection studies, and suggests a way to create distribution models at spatial and temporal scales
NASA Astrophysics Data System (ADS)
Liu, Y.; He, H. T.; Zhu, C.
2014-12-01
Several important equilibrium Si isotope fractionation factors are calculated here. We use a so-called volume-variable-cluster-model (VVCM) method for solids and the "water-droplet" method for aqueous species for isotope fractionation calculation at the same quantum chemistry level. The calculation results show that several silicate minerals, such as quartz, feldspar, kaolinite, etc., all enrich heavy Si isotopes relative to aqueous H4SiO4 and can be up to 3.3‰ at 25°C, different from most field observations. Meanwhile stable organosilicon complexes can enrich even lighter Si isotopes than aqueous H4SiO4. For explaining the difference between the calculation results and field observations, we calculate the kinetic isotope effect (KIE) associated with the formation of amorphous silica, and find that amorphous silica will enrich extremely light Si isotopes. From amorphous silica to crystalline quartz, the structural adjustment & transition needs getting rid of small amount of Si to re-organize the structure. Light Si isotopes will be preferentially lost and let the final crystalline quartz with a little bit more heavy Si isotopes. However, such late-stage Si heavy isotope enrichment cannot erase the total isotopic signal, crystalline quartz still inherit much light Si isotopic composition from amorphous quartz. That is the reason for the discrepancy between the calculation results and the field observations, because the formation of amorphous quartz is under a non-equilibrium process but theoretical calculations are for equilibrium isotope fractionations. With accurate equilibrium fractionation factors provided here, Si isotope distributions in earth surface environments including soil, groundwater and plants can be further interpreted. We find that δ30Si variations in soil are mainly driven by secondary minerals precipitation and adsorption. Also, bulk soil δ30Si maybe have a parabolic distribution with soil age, with a minimum value at where allophane is
Kaita, R.; Kozub, T.; Logan, N.; Majeski, R.; Menard, J.; Zakharov, L.
2010-12-10
The lithium tokamak experiment LTX is a modest-sized spherical tokamak R0=0.4 m and a =0.26 m designed to investigate the low-recycling lithium wall operating regime for magnetically confined plasmas. LTX will reach this regime through a lithium-coated shell internal to the vacuum vessel, conformal to the plasma last-closed-flux surface, and heated to 300-400 oC. This structure is highly conductive and not axisymmetric. The three-dimensional nature of the shell causes the eddy currents and magnetic fields to be three-dimensional as well. In order to analyze the plasma equilibrium in the presence of three-dimensional eddy currents, an extensive array of unique magnetic diagnostics has been implemented. Sensors are designed to survive high temperatures and incidental contact with lithium and provide data on toroidal asymmetries as well as full coverage of the poloidal cross-section. The magnetic array has been utilized to determine the effects of nonaxisymmetric eddy currents and to model the start-up phase of LTX. Measurements from the magnetic array, coupled with two-dimensional field component modeling, have allowed a suitable field null and initial plasma current to be produced. For full magnetic reconstructions, a three-dimensional electromagnetic model of the vacuum vessel and shell is under development.
Berzak, L.; Jones, A. D.; Kaita, R.; Kozub, T.; Logan, N.; Majeski, R.; Menard, J.; Zakharov, L.
2010-10-15
The lithium tokamak experiment (LTX) is a modest-sized spherical tokamak (R{sub 0}=0.4 m and a=0.26 m) designed to investigate the low-recycling lithium wall operating regime for magnetically confined plasmas. LTX will reach this regime through a lithium-coated shell internal to the vacuum vessel, conformal to the plasma last-closed-flux surface, and heated to 300-400 deg. C. This structure is highly conductive and not axisymmetric. The three-dimensional nature of the shell causes the eddy currents and magnetic fields to be three-dimensional as well. In order to analyze the plasma equilibrium in the presence of three-dimensional eddy currents, an extensive array of unique magnetic diagnostics has been implemented. Sensors are designed to survive high temperatures and incidental contact with lithium and provide data on toroidal asymmetries as well as full coverage of the poloidal cross-section. The magnetic array has been utilized to determine the effects of nonaxisymmetric eddy currents and to model the start-up phase of LTX. Measurements from the magnetic array, coupled with two-dimensional field component modeling, have allowed a suitable field null and initial plasma current to be produced. For full magnetic reconstructions, a three-dimensional electromagnetic model of the vacuum vessel and shell is under development.
Thermal equilibrium of non-neutral plasma in dipole magnetic field
Sato, N.; Kasaoka, N.; Yoshida, Z.
2015-04-15
Self-organization of a long-lived structure is one of the remarkable characteristics of macroscopic systems governed by long-range interactions. In a homogeneous magnetic field, a non-neutral plasma creates a “thermal equilibrium,” which is a Boltzmann distribution on a rigidly rotating frame. Here, we study how a non-neutral plasma self-organizes in inhomogeneous magnetic field; as a typical system, we consider a dipole magnetic field. In this generalized setting, the plasma exhibits its fundamental mechanism that determines the relaxed state. The scale hierarchy of adiabatic invariants is the determinant; the Boltzmann distribution under the topological constraint by the robust adiabatic invariants (hence, the homogeneous distribution with respect to the fragile invariant) is the relevant relaxed state, which turns out to be a rigidly rotating clump of particles (just same as in a homogeneous magnetic field), while the density is no longer homogeneous.
The temperature and size distribution of large water clusters from a non-equilibrium model
Gimelshein, N.; Gimelshein, S.; Pradzynski, C. C.; Zeuch, T.; Buck, U.
2015-06-28
A hybrid Lagrangian-Eulerian approach is used to examine the properties of water clusters formed in neon-water vapor mixtures expanding through microscale conical nozzles. Experimental size distributions were reliably determined by the sodium doping technique in a molecular beam machine. The comparison of computed size distributions and experimental data shows satisfactory agreement, especially for (H{sub 2}O){sub n} clusters with n larger than 50. Thus validated simulations provide size selected cluster temperature profiles in and outside the nozzle. This information is used for an in-depth analysis of the crystallization and water cluster aggregation dynamics of recently reported supersonic jet expansion experiments.
The temperature and size distribution of large water clusters from a non-equilibrium model.
Gimelshein, N; Gimelshein, S; Pradzynski, C C; Zeuch, T; Buck, U
2015-06-28
A hybrid Lagrangian-Eulerian approach is used to examine the properties of water clusters formed in neon-water vapor mixtures expanding through microscale conical nozzles. Experimental size distributions were reliably determined by the sodium doping technique in a molecular beam machine. The comparison of computed size distributions and experimental data shows satisfactory agreement, especially for (H2O)n clusters with n larger than 50. Thus validated simulations provide size selected cluster temperature profiles in and outside the nozzle. This information is used for an in-depth analysis of the crystallization and water cluster aggregation dynamics of recently reported supersonic jet expansion experiments. PMID:26133426
Generalizing the Boltzmann equation in complex phase space.
Zadehgol, Abed
2016-08-01
In this work, a generalized form of the BGK-Boltzmann equation is proposed, where the velocity, position, and time can be represented by real or complex variables. The real representation leads to the conventional BGK-Boltzmann equation, which can recover the continuity and Navier-Stokes equations. We show that the complex representation yields a different set of equations, and it can also recover the conservation and Navier-Stokes equations, at low Mach numbers, provided that the imaginary component of the macroscopic mass can be neglected. We briefly review the Constant Speed Kinetic Model (CSKM), which was introduced in Zadehgol and Ashrafizaadeh [J. Comp. Phys. 274, 803 (2014)JCTPAH0021-999110.1016/j.jcp.2014.06.053] and Zadehgol [Phys. Rev. E 91, 063311 (2015)PLEEE81539-375510.1103/PhysRevE.91.063311]. The CSKM is then used as a basis to show that the complex-valued equilibrium distribution function of the present model can be identified with a simple singularity in the complex phase space. The virtual particles, in the present work, are concentrated on virtual "branes" which surround the computational nodes. Employing the Cauchy integral formula, it is shown that certain variations of the "branes," in the complex phase space, do not affect the local kinetic states. This property of the new model, which is referred to as the "apparent jumps" in the present work, is used to construct new models. The theoretical findings have been tested by simulating three benchmark flows. The results of the present simulations are in excellent agreement with the previous results reported by others. PMID:27627421
ON QUIET-TIME SOLAR WIND ELECTRON DISTRIBUTIONS IN DYNAMICAL EQUILIBRIUM WITH LANGMUIR TURBULENCE
Zaheer, S.; Yoon, P. H.
2013-10-01
A recent series of papers put forth a self-consistent theory of an asymptotically steady-state electron distribution function and Langmuir turbulence intensity. The theory was developed in terms of the κ distribution which features Maxwellian low-energy electrons and a non-Maxwellian energetic power-law tail component. The present paper discusses a generalized κ distribution that features a Davydov-Druyvesteyn type of core component and an energetic power-law tail component. The physical motivation for such a generalization is so that the model may reflect the influence of low-energy electrons interacting with low-frequency kinetic Alfvénic turbulence as well as with high-frequency Langmuir turbulence. It is shown that such a solution and the accompanying Langmuir wave spectrum rigorously satisfy the balance requirement between the spontaneous and induced emission processes in both the particle and wave kinetic equations, and approximately satisfy the similar balance requirement between the spontaneous and induced scattering processes, which are nonlinear. In spite of the low velocity modification of the electron distribution function, it is shown that the resulting asymptotic velocity power-law index α, where f{sub e} ∼ v {sup –α} is close to the average index observed during the quiet-time solar wind condition, i.e., α ∼ O(6.5) whereas α{sub average} ∼ 6.69, according to observation.
NASA Astrophysics Data System (ADS)
Garrett, T. J.; Alva, S.; Glenn, I. B.; Krueger, S. K.
2015-12-01
There are two possible approaches for parameterizing sub-grid cloud dynamics in a coarser grid model. The most common is to use a fine scale model to explicitly resolve the mechanistic details of clouds to the best extent possible, and then to parameterize these behaviors cloud state for the coarser grid. A second is to invoke physical intuition and some very general theoretical principles from equilibrium statistical mechanics. This approach avoids any requirement to resolve time-dependent processes in order to arrive at a suitable solution. The second approach is widely used elsewhere in the atmospheric sciences: for example the Planck function for blackbody radiation is derived this way, where no mention is made of the complexities of modeling a large ensemble of time-dependent radiation-dipole interactions in order to obtain the "grid-scale" spectrum of thermal emission by the blackbody as a whole. We find that this statistical approach may be equally suitable for modeling convective clouds. Specifically, we make the physical argument that the dissipation of buoyant energy in convective clouds is done through mixing across a cloud perimeter. From thermodynamic reasoning, one might then anticipate that vertically stacked isentropic surfaces are characterized by a power law dlnN/dlnP = -1, where N(P) is the number clouds of perimeter P. In a Giga-LES simulation of convective clouds within a 100 km square domain we find that such a power law does appear to characterize simulated cloud perimeters along isentropes, provided a sufficient cloudy sample. The suggestion is that it may be possible to parameterize certain important aspects of cloud state without appealing to computationally expensive dynamic simulations.
Noronha, Jorge; Denicol, Gabriel S.
2015-12-30
In this paper we obtain an analytical solution of the relativistic Boltzmann equation under the relaxation time approximation that describes the out-of-equilibrium dynamics of a radially expanding massless gas. This solution is found by mapping this expanding system in flat spacetime to a static flow in the curved spacetime AdS2 Ⓧ S2. We further derive explicit analytic expressions for the momentum dependence of the single-particle distribution function as well as for the spatial dependence of its moments. We find that this dissipative system has the ability to flow as a perfect fluid even though its entropy density does not matchmore » the equilibrium form. The nonequilibrium contribution to the entropy density is shown to be due to higher-order scalar moments (which possess no hydrodynamical interpretation) of the Boltzmann equation that can remain out of equilibrium but do not couple to the energy-momentum tensor of the system. Furthermore, in this system the slowly moving hydrodynamic degrees of freedom can exhibit true perfect fluidity while being totally decoupled from the fast moving, nonhydrodynamical microscopic degrees of freedom that lead to entropy production.« less
Schuck, Peter; Gillis, Richard B.; Besong, Tabot M.D.; Almutairi, Fahad; Adams, Gary G.; Rowe, Arthur J.; Harding, Stephen E.
2014-01-01
Sedimentation equilibrium (analytical ultracentrifugation) is one of the most inherently suitable methods for the determination of average molecular weights and molecular weight distributions of polymers, because of its absolute basis (no conformation assumptions) and inherent fractionation ability (without the need for columns or membranes and associated assumptions over inertness). With modern instrumentation it is also possible to run up to 21 samples simultaneously in a single run. Its application has been severely hampered because of difficulties in terms of baseline determination (incorporating estimation of the concentration at the air/solution meniscus) and complexity of the analysis procedures. We describe a new method for baseline determination based on a smart-smoothing principle and built into the highly popular platform SEDFIT for the analysis of the sedimentation behavior of natural and synthetic polymer materials. The SEDFIT-MSTAR procedure – which takes only a few minutes to perform - is tested with four synthetic data sets (including a significantly non-ideal system) a naturally occurring protein (human IgG1) and two naturally occurring carbohydrate polymers (pullulan and λ–carrageenan) in terms of (i) weight average molecular weight for the whole distribution of species in the sample (ii) the variation in “point” average molecular weight with local concentration in the ultracentrifuge cell and (iii) molecular weight distribution. PMID:24244936
On the Non-Equilibrium Population Distribution of E-Methanol in Dark Clouds
NASA Astrophysics Data System (ADS)
Wollman, Emma
2007-12-01
The goal of this project was to determine the typical distribution of rotational level populations in the k=0 ladder of E-methanol in dark clouds in order to provide another observational test for theoretical models of pumping. We used our own observations of several sources with the 12-m ARO telescope on Kitt Peak as well as the published observational results by Slysh et al. (1999). The relative level populations (excitation temperatures) were determined from the measured intensity ratios of a series of the J(0)-J(-1) transitions of E-methanol under the assumption of spontaneous, optically thin emission. We observed the J(0)-J(-1) lines in six sources: W75N, DR21N, DR21, and three positions at DR21OH. The J=1 to J=5 lines were observed for all sources and the J=7 line was observed for W75N, DR21N, and one position in DR21OH. We also used Slysh et al.'s results for the J=1 through 4 lines in 52 sources, for the J=5 line in 50 sources, for the J=6 line in 15 sources, and for the J=7 and 8 lines in 2 sources. We determined the excitation temperatures of the involved levels in the k=0 ladder relative to the 1(0) level for each source and averaged the results over the sources. The average excitation temperatures demonstrate strong evidence of overcooling in the k=0 ladder - the excitation temperature increases linearly with increasing energy, from 8 K to 35 K. Our observations confirm this tendency of overcooling. We will discuss the agreement of these results with the predictions of the current models of methanol pumping. The author thanks the technical staff of the 12-m ARO telescope for help with the observations. This project was supported by the NSF/REU grant AST-0354056 and the Nantucket Maria Mitchell Association.
Jacxsens, L; Devlieghere, F; Debevere, J
2002-03-01
The impact of temperature fluctuations in a simulated cold distribution chain, typical of commercial practice, was investigated on both the microbial and sensorial quality of equilibrium modified atmosphere (EMA) packaged minimally processed vegetables. The internal O2 concentration of the designed packages could be predicted for the different steps of the simulated distribution chain by applying an integrated mathematical system. The internal atmosphere in the packages remained in its aerobic range during storage in the chain due to the application of high permeable packaging films for O2 and CO2. Spoilage microorganisms were proliferating fast on minimally processed bell peppers and lettuce. Yeasts showed to be the shelf-life limiting group. Visual properties limited the sensorial shelf-life. Listeria monocytogenes was able to multiply on cucumber slices, survived on minimally processed lettuce and decreased in number on bell peppers due to the combination of low pH and refrigeration. Aeromonas caviae was multiplying on both cucumber slices and mixed lettuce, but was as well inhibited by the low pH of bell peppers. Storage temperature control was found to be of paramount importance for the microbial (spoilage and safety) and sensorial quality evaluation of EMA-packaged minimally processed vegetables. PMID:11934040
Information geometry of Boltzmann machines.
Amari, S; Kurata, K; Nagaoka, H
1992-01-01
A Boltzmann machine is a network of stochastic neurons. The set of all the Boltzmann machines with a fixed topology forms a geometric manifold of high dimension, where modifiable synaptic weights of connections play the role of a coordinate system to specify networks. A learning trajectory, for example, is a curve in this manifold. It is important to study the geometry of the neural manifold, rather than the behavior of a single network, in order to know the capabilities and limitations of neural networks of a fixed topology. Using the new theory of information geometry, a natural invariant Riemannian metric and a dual pair of affine connections on the Boltzmann neural network manifold are established. The meaning of geometrical structures is elucidated from the stochastic and the statistical point of view. This leads to a natural modification of the Boltzmann machine learning rule. PMID:18276427
Boltzmann Fluctuations in Numerical Simulations of Nonequilibrium Lattice Threshold Systems
Rundle, J.B.; Klein, W.; Gross, S.; Turcotte, D.L.
1995-08-21
Nonequilibrium threshold systems such as slider blocks are now used to model a variety of dynamical systems, including earthquake faults, driven neural networks, and sliding charge density waves. We show that for general mean field models driven at low rates fluctuations in the internal energy field are characterized by Boltzmann statistics. Numerical simulations confirm this prediction. Our results indicate that mean field models can be effectively treated as equilibrium systems.
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.
Chen, Yunjie; Roux, Benoît
2015-01-14
A family of hybrid simulation methods that combines the advantages of Monte Carlo (MC) with the strengths of classical molecular dynamics (MD) consists in carrying out short non-equilibrium MD (neMD) trajectories to generate new configurations that are subsequently accepted or rejected via an MC process. In the simplest case where a deterministic dynamic propagator is used to generate the neMD trajectories, the familiar Metropolis acceptance criterion based on the change in the total energy ΔE, min[1, exp( − βΔE)], guarantees that the hybrid algorithm will yield the equilibrium Boltzmann distribution. However, the functional form of the acceptance probability is more complex when the non-equilibrium switching process is generated via a non-deterministic stochastic dissipative propagator coupled to a heat bath. Here, we clarify the conditions under which the Metropolis criterion remains valid to rigorously yield a proper equilibrium Boltzmann distribution within hybrid neMD-MC algorithm.
Ling, G N; Niu, Z; Ochsenfeld, M
1993-01-01
We determined the equilibrium distribution of twenty-one nonmetabolized nonelectrolytes in frog muscle cells. In all cases, plots of the equilibrium intracellular concentrations of a solute in the cell water against the external concentrations of the solute yielded straight lines in agreement with the prediction of such a rectilinear plot by the polarized multilayer (PM) theory. The slopes of these straight lines yield the equilibrium distribution coefficients or q-value of that solute. It was shown that, again in agreement with the PM theory, the q-values of fourteen nonelectrolytes vary with the molecular volumes of the nonelectrolytes, obeying the "size rule", i.e., the larger the solute, the lower its q-value. The q-values of the remaining seven nonelectrolytes also decrease with their molecular volumes but on a separate curve. These q-value vs. molecular volume plots (q-v plots) show strong resemblance to similar q-v plots of solutes in dialysis sacs containing proteins and polymers assuming the fully-extended conformation (extrovert models) but no, or only weak, resemblance to q-v plots of solutions containing native globular proteins (introvert models). These findings also support the PM theory, according to which some protein(s) pervasively present in cells are in the fully-extended conformation; and that these fully extended cell protein(s) polarize(s) in multilayers all or virtually all cell water. The relationship between the q-values of the nonelectrolytes and the solutes' respective molecular volume are described by two sets of theoretical curves, calculated from an equation introduced in the preceding paper. Both curves were computed on the basis of the same exclusion intensity (Uvp = 126 cal/mole). This factor measures the extra water-to-water interaction of the polarized water which acts to keep solute out of the cell water in degree according to the size of the solute. The two curves are computed on the basis of two different values of U(s), which
Crystallographic Lattice Boltzmann Method.
Namburi, Manjusha; Krithivasan, Siddharth; Ansumali, Santosh
2016-01-01
Current approaches to Direct Numerical Simulation (DNS) are computationally quite expensive for most realistic scientific and engineering applications of Fluid Dynamics such as automobiles or atmospheric flows. The Lattice Boltzmann Method (LBM), with its simplified kinetic descriptions, has emerged as an important tool for simulating hydrodynamics. In a heterogeneous computing environment, it is often preferred due to its flexibility and better parallel scaling. However, direct simulation of realistic applications, without the use of turbulence models, remains a distant dream even with highly efficient methods such as LBM. In LBM, a fictitious lattice with suitable isotropy in the velocity space is considered to recover Navier-Stokes hydrodynamics in macroscopic limit. The same lattice is mapped onto a cartesian grid for spatial discretization of the kinetic equation. In this paper, we present an inverted argument of the LBM, by making spatial discretization as the central theme. We argue that the optimal spatial discretization for LBM is a Body Centered Cubic (BCC) arrangement of grid points. We illustrate an order-of-magnitude gain in efficiency for LBM and thus a significant progress towards feasibility of DNS for realistic flows. PMID:27251098
Crystallographic Lattice Boltzmann Method
Namburi, Manjusha; Krithivasan, Siddharth; Ansumali, Santosh
2016-01-01
Current approaches to Direct Numerical Simulation (DNS) are computationally quite expensive for most realistic scientific and engineering applications of Fluid Dynamics such as automobiles or atmospheric flows. The Lattice Boltzmann Method (LBM), with its simplified kinetic descriptions, has emerged as an important tool for simulating hydrodynamics. In a heterogeneous computing environment, it is often preferred due to its flexibility and better parallel scaling. However, direct simulation of realistic applications, without the use of turbulence models, remains a distant dream even with highly efficient methods such as LBM. In LBM, a fictitious lattice with suitable isotropy in the velocity space is considered to recover Navier-Stokes hydrodynamics in macroscopic limit. The same lattice is mapped onto a cartesian grid for spatial discretization of the kinetic equation. In this paper, we present an inverted argument of the LBM, by making spatial discretization as the central theme. We argue that the optimal spatial discretization for LBM is a Body Centered Cubic (BCC) arrangement of grid points. We illustrate an order-of-magnitude gain in efficiency for LBM and thus a significant progress towards feasibility of DNS for realistic flows. PMID:27251098
Crystallographic Lattice Boltzmann Method
NASA Astrophysics Data System (ADS)
Namburi, Manjusha; Krithivasan, Siddharth; Ansumali, Santosh
2016-06-01
Current approaches to Direct Numerical Simulation (DNS) are computationally quite expensive for most realistic scientific and engineering applications of Fluid Dynamics such as automobiles or atmospheric flows. The Lattice Boltzmann Method (LBM), with its simplified kinetic descriptions, has emerged as an important tool for simulating hydrodynamics. In a heterogeneous computing environment, it is often preferred due to its flexibility and better parallel scaling. However, direct simulation of realistic applications, without the use of turbulence models, remains a distant dream even with highly efficient methods such as LBM. In LBM, a fictitious lattice with suitable isotropy in the velocity space is considered to recover Navier-Stokes hydrodynamics in macroscopic limit. The same lattice is mapped onto a cartesian grid for spatial discretization of the kinetic equation. In this paper, we present an inverted argument of the LBM, by making spatial discretization as the central theme. We argue that the optimal spatial discretization for LBM is a Body Centered Cubic (BCC) arrangement of grid points. We illustrate an order-of-magnitude gain in efficiency for LBM and thus a significant progress towards feasibility of DNS for realistic flows.
Three-dimensional lattice Boltzmann model for magnetic reconnection.
Mendoza, M; Muñoz, J D
2008-02-01
We develop a three-dimensional (3D) lattice Boltzmann model that recovers in the continuous limit the two-fluids theory for plasmas, and consequently includes the generalized Ohm's law. The model reproduces the magnetic reconnection process just by giving the right initial equilibrium conditions in the magnetotail, without any assumption on the resistivity in the diffusive region. In this model, the plasma is handled similar to two fluids with an interaction term, each one with distribution functions associated to a cubic lattice with 19 velocities (D3Q19). The electromagnetic fields are considered as a third fluid with an external force on a cubic lattice with 13 velocities (D3Q13). The model can simulate either viscous fluids in the incompressible limit or nonviscous compressible fluids, and successfully reproduces both the Hartmann flow and the magnetic reconnection in the magnetotail. The reconnection rate in the magnetotail obtained with this model lies between R=0.062 and R=0.073, in good agreement with the observations. PMID:18352154
Treatment of moving boundaries in lattice-Boltzmann simulations.
NASA Astrophysics Data System (ADS)
Indireshkumar, K.; Pal, A.; Brasseur, J. G.
2000-11-01
We consider the treatment of moving boundaries with the lattice-Boltzmann (LB) technique, where the treatment of the boundary often does not precisely conserve mass and spurious fluctuations in density/pressure result from boundary motion through fixed grids. First, we applied the extrapolation method proposed by Chen et. al.(S. Y. Chen, D. Martinez, and R Mei, Phys. Fluids) 8, 2527 (1996) to incompressible flow induced by the movement of a piston in a 2D ``cylinder'' with mass flow out of or into the cylinder. In these simulations, the velocity of the boundary nodes is set equal to the (known) velocity of the boundary (piston) in the equilibrium distribution function (Method I). In a second set of simulations, the boundary node velocities are obtained by interpolating between interior nodes and the boundary, thus including the effect of boundary position more precisely (Method II). Comparison of LB predictions with simulations using FIDAP show pressure agreement to witnin 2 %. The total mass is conserved to within 0.1% with Method I and improves to within 0.02 % using method II. Spurious fluctuations in density/pressure due to boundary movement is about 0.9% with Method I, which improves significantly to about 0.3% with Method II. The application of these simple techniques to more complex geometries and wall (and fluid) motions in a stomach during gastric emptying will be presented.
Three-dimensional lattice Boltzmann model for magnetic reconnection
Mendoza, M.; Munoz, J. D.
2008-02-15
We develop a three-dimensional (3D) lattice Boltzmann model that recovers in the continuous limit the two-fluids theory for plasmas, and consequently includes the generalized Ohm's law. The model reproduces the magnetic reconnection process just by giving the right initial equilibrium conditions in the magnetotail, without any assumption on the resistivity in the diffusive region. In this model, the plasma is handled similar to two fluids with an interaction term, each one with distribution functions associated to a cubic lattice with 19 velocities (D3Q19). The electromagnetic fields are considered as a third fluid with an external force on a cubic lattice with 13 velocities (D3Q13). The model can simulate either viscous fluids in the incompressible limit or nonviscous compressible fluids, and successfully reproduces both the Hartmann flow and the magnetic reconnection in the magnetotail. The reconnection rate in the magnetotail obtained with this model lies between R=0.062 and R=0.073, in good agreement with the observations.
Observation of a Persistent Non-Equilibrium State in an Extremely Isotropic Harmonic Potential
NASA Astrophysics Data System (ADS)
Lobser, D. S.
Ludwig Boltzmann made tremendously important contributions to the problem of con- necting macroscopic, empirical phenomena with microscopic, atomistic dynamics. At the end of the nineteenth century, Boltzmann was confronted with various strong objections to his work. For example, Boltzmann's atomistic explanations presuppose the reality of atoms, a notion that was vigorously rejected in some circles [14, 38]. Then too, there was the critique by Loschmidt that Boltzmann's H-theorem, put forth as a microscopic explanation for the Second Law of Thermodynamics, could hardly account for irreversible physics when the individual two-atom collisions were each reversible [18, 42]. Still intriguing today is the existence of special cases of the Boltzmann equation in which time-varying distributions of atoms re- sist the imperative of equilibration, even in the presence of collisions. Boltzmann discussed such situations in a paper dedicated to responding to Loschmidt's critique [7, 4]. Perhaps Boltzmann's motivation was to enumerate special cases where his famous H value does not relax as it should, and by enumerating them, point out their nonnaturalness, their artificiality. Damping, or relaxation to equilibrium, of a time-invariant phase-space distribution, is an all-but universal result predicted by the Boltzmann equation. Such improbable systems of atoms have only very recently been realized experimentally. Kinoshita et al. [36] experimentally confirmed that atoms constrained to move in a quasi one-dimensional potential, an atomistic Newtons cradle, exhibit vastly suppressed relaxation. Chevy et al. [15] observed long-lived breathe-mode oscillations in highly elongated but still 3D geometries. Perhaps one of the more interesting cases is the vanishing damping of the monopole breathe-mode oscillation in a spherically symmetric harmonic oscillator [29], where a cloud of atoms experiences undamped temperature oscillations, causing the cloud to expand and contract as if it
NASA Astrophysics Data System (ADS)
Chau, J. F.; Or, D.; Jones, S.; Sukop, M.
2004-05-01
Liquid distribution in unsaturated porous media under different gravitational forces and resulting gaseous diffusion coefficients were investigated to enhance understanding of plant growth conditions in microgravity. Different fluid behavior in plant growth media under microgravity conditions as compared to earth presents a challenge to plant growth in long duration space exploration missions. Our primary objective was to provide qualitative description and quantitative measures of the role of reduced gravity on hydraulic and gaseous transport properties in simulated porous media. We implemented a multi-phase lattice Boltzmann code for equilibrium distribution of liquid in an idealized two-dimensional porous medium under microgravity and "normal" gravity conditions. The information was then used to provide boundary conditions for simulation of gaseous diffusion through the equilibrium domains (considering diffusion through liquid phase negligibly small). The models were tested by comparison with several analytical solutions to the diffusion equation, with excellent results. The relative diffusion coefficient for both series of simulations (with and without gravity) as functions of air-filled porosity was in good agreement with established models of Millington-Quirk. Liquid distribution under earth's gravity featured increased water content at the lower part of the medium relative to the distribution in reduced gravity, which resulted in decreased gas diffusion through a vertically oriented column of a porous medium. Simulation results for larger domains under various orientations will be presented.
Koga, S.; Shibata, T.; Terasaki, R.; Kameyama, N.; Hatayama, A.; Bacal, M.; Tsumori, K.
2012-02-15
In negative ion sources for the neutral beam injection, it is important to calculate H atom flux onto the plasma grid (PG) surface for the evaluation of H{sup -} production on the PG surface. We have developed a neutral (H{sub 2} molecules and H atoms) transport code. In the present study, the neutral transport code is applied to the analysis of the H{sub 2} and H transport in a NIFS-R and D ion source in order to calculate the flux onto the PG surface. Taking into account non-equilibrium feature of the electron energy distribution function (EEDF), i.e., the fast electron component, we have done the neutral transport simulation. The results suggest that the precise evaluation of the EEDF, especially in the energy range 15 eV < E < 30 eV is important for the dissociation rate of H{sub 2} molecules by the electron impact collision and the resultant H atom flux on the PG.
A new lattice Boltzmann model for solving the coupled viscous Burgers’ equation
NASA Astrophysics Data System (ADS)
Lai, Huilin; Ma, Changfeng
2014-02-01
In this paper, a new lattice Boltzmann model for the coupled nonlinear system of viscous Burgers’ equation is proposed by using the double evolutionary equations. Through selecting equilibrium distribution functions and amending functions properly, the governing evolution system can be recovered correctly according to our proposed scheme, in which the Chapman-Enskog expansion is employed. The effects of space and time resolutions on the accuracy and stability of the model are numerically investigated in detail. The numerical solutions for various initial and boundary conditions are calculated and validated against analytic solutions or other numerical solutions reported in previous studies. It is found that the numerical results agree well with the analytic solutions, which indicates the potential of the present algorithm for solving the coupled nonlinear system of viscous Burgers’ equation.
Lattice Boltzmann method for the age concentration equation in shallow water
NASA Astrophysics Data System (ADS)
Liu, Haifei; Ding, Yu; Wang, Hongda; Zhang, Jie
2015-10-01
Water age is a critical parameter in reflecting the extent of water exchange. It represents the time that water parcels or contaminants are transported from source to current positions. In this study, an equilibrium distribution function for water age concentration is proposed within the Eulerian framework based on the existing theory of water age, and it can recover the age concentration equation. In addition, the lattice Boltzmann method for water age in the Lagrangian procedures is developed. This method also enables the Lagrangian age to be fundamentally simulated under computationally expensive conditions. In numerical tests, cubic and circular reservoirs with narrow inflow-outflow boundaries are used to verify the applicability of the model. Finally, the proposed approaches are applied to the Baiyangdian Lake, the biggest freshwater lake in northern China. The result is compared with that acquired by the Environmental Fluid Dynamic Code (EFDC) as well.
A modified multiple-relaxation-time lattice Boltzmann model for convection-diffusion equation
NASA Astrophysics Data System (ADS)
Huang, Rongzong; Wu, Huiying
2014-10-01
A modified lattice Boltzmann model with multiple relaxation times (MRT) for the convection-diffusion equation (CDE) is proposed. By modifying the relaxation matrix, as well as choosing the corresponding equilibrium distribution function properly, the present model can recover the CDE with anisotropic diffusion coefficient with no deviation term even when the velocity vector varies generally with space or time through the Chapman-Enskog analysis. This model is firstly validated by simulating the diffusion of a Gaussian hill, which demonstrates it can handle the anisotropic diffusion problem correctly. Then it is adopted to calculate the longitudinal dispersion coefficient of the Taylo-Aris dispersion. Numerical results show that the present model can further reduce the numerical error under the condition of non-zero velocity vector, especially when the dimensionless relaxation time is relatively large.
Privacy-Preserving Restricted Boltzmann Machine
Li, Yu
2014-01-01
With the arrival of the big data era, it is predicted that distributed data mining will lead to an information technology revolution. To motivate different institutes to collaborate with each other, the crucial issue is to eliminate their concerns regarding data privacy. In this paper, we propose a privacy-preserving method for training a restricted boltzmann machine (RBM). The RBM can be got without revealing their private data to each other when using our privacy-preserving method. We provide a correctness and efficiency analysis of our algorithms. The comparative experiment shows that the accuracy is very close to the original RBM model. PMID:25101139
Privacy-preserving restricted boltzmann machine.
Li, Yu; Zhang, Yuan; Ji, Yue
2014-01-01
With the arrival of the big data era, it is predicted that distributed data mining will lead to an information technology revolution. To motivate different institutes to collaborate with each other, the crucial issue is to eliminate their concerns regarding data privacy. In this paper, we propose a privacy-preserving method for training a restricted boltzmann machine (RBM). The RBM can be got without revealing their private data to each other when using our privacy-preserving method. We provide a correctness and efficiency analysis of our algorithms. The comparative experiment shows that the accuracy is very close to the original RBM model. PMID:25101139
Determining Equilibrium Position For Acoustical Levitation
NASA Technical Reports Server (NTRS)
Barmatz, M. B.; Aveni, G.; Putterman, S.; Rudnick, J.
1989-01-01
Equilibrium position and orientation of acoustically-levitated weightless object determined by calibration technique on Earth. From calibration data, possible to calculate equilibrium position and orientation in presence of Earth gravitation. Sample not levitated acoustically during calibration. Technique relies on Boltzmann-Ehrenfest adiabatic-invariance principle. One converts resonant-frequency-shift data into data on normalized acoustical potential energy. Minimum of energy occurs at equilibrium point. From gradients of acoustical potential energy, one calculates acoustical restoring force or torque on objects as function of deviation from equilibrium position or orientation.
Bhowmick, Amiya Ranjan; Bandyopadhyay, Subhadip; Rana, Sourav; Bhattacharya, Sabyasachi
2016-01-01
The stochastic versions of the logistic and extended logistic growth models are applied successfully to explain many real-life population dynamics and share a central body of literature in stochastic modeling of ecological systems. To understand the randomness in the population dynamics of the underlying processes completely, it is important to have a clear idea about the quasi-equilibrium distribution and its moments. Bartlett et al. (1960) took a pioneering attempt for estimating the moments of the quasi-equilibrium distribution of the stochastic logistic model. Matis and Kiffe (1996) obtain a set of more accurate and elegant approximations for the mean, variance and skewness of the quasi-equilibrium distribution of the same model using cumulant truncation method. The method is extended for stochastic power law logistic family by the same and several other authors (Nasell, 2003; Singh and Hespanha, 2007). Cumulant truncation and some alternative methods e.g. saddle point approximation, derivative matching approach can be applied if the powers involved in the extended logistic set up are integers, although plenty of evidence is available for non-integer powers in many practical situations (Sibly et al., 2005). In this paper, we develop a set of new approximations for mean, variance and skewness of the quasi-equilibrium distribution under more general family of growth curves, which is applicable for both integer and non-integer powers. The deterministic counterpart of this family of models captures both monotonic and non-monotonic behavior of the per capita growth rate, of which theta-logistic is a special case. The approximations accurately estimate the first three order moments of the quasi-equilibrium distribution. The proposed method is illustrated with simulated data and real data from global population dynamics database. PMID:26561778
NASA Astrophysics Data System (ADS)
Nomura, Yasunori
2015-10-01
Understanding the observed arrow of time is equivalent, under general assumptions, to explaining why Boltzmann brains do not overwhelm ordinary observers. It is usually thought that this provides a condition on the decay rate of every cosmologically accessible de Sitter vacuum, and that this condition is determined by the production rate of Boltzmann brains calculated using semiclassical theory built on each such vacuum. We argue, based on a recently developed picture of microscopic quantum gravitational degrees of freedom, that this thinking needs to be modified. In particular, depending on the structure of the fundamental theory, the decay rate of a de Sitter vacuum may not have to satisfy any condition except for the one imposed by the Poincaré recurrence. The framework discussed here also addresses the question of whether a Minkowski vacuum may produce Boltzmann brains.
Zermelo, Boltzmann, and the recurrence paradox
NASA Astrophysics Data System (ADS)
Steckline, Vincent S.
1983-10-01
The papers exchanged by Ludwig Boltzmann and Ernst Zermelo concerning the recurrence paradox are summarized. The historical context of the paradox, Zermelo's proof of the paradox, his opinions of its consequences, Boltzmann's reply, and the ensuing discussion are described.
NASA Astrophysics Data System (ADS)
Hong, Y.; Kirschbaum, D. B.; Fukuoka, H.
2011-12-01
The key to advancing the predictability of rainfall-triggered landslides is to use physically based slope-stability models that simulate the dynamical response of the subsurface moisture to spatiotemporal variability of rainfall in complex terrains. An early warning system applying such physical models has been developed to predict rainfall-induced shallow landslides over Java Island in Indonesia and Honduras. The prototyped early warning system integrates three major components: (1) a susceptibility mapping or hotspot identification component based on a land surface geospatial database (topographical information, maps of soil properties, and local landslide inventory etc.); (2) a satellite-based precipitation monitoring system (http://trmm.gsfc.nasa.gov) and a precipitation forecasting model (i.e. Weather Research Forecast); and (3) a physically-based, rainfall-induced landslide prediction model SLIDE (SLope-Infiltration-Distributed Equilibrium). The system utilizes the modified physical model to calculate a Factor of Safety (FS) that accounts for the contribution of rainfall infiltration and partial saturation to the shear strength of the soil in topographically complex terrains. The system's prediction performance has been evaluated using a local landslide inventory. In Java Island, Indonesia, evaluation of SLIDE modeling results by local news reports shows that the system successfully predicted landslides in correspondence to the time of occurrence of the real landslide events. Further study of SLIDE is implemented in Honduras where Hurricane Mitch triggered widespread landslides in 1998. Results shows within the approximately 1,200 square kilometers study areas, the values of hit rates reached as high as 78% and 75%, while the error indices were 35% and 49%. Despite positive model performance, the SLIDE model is limited in the early warning system by several assumptions including, using general parameter calibration rather than in situ tests and neglecting
NASA Astrophysics Data System (ADS)
Izzo, Dario; Petazzi, Lorenzo
2006-08-01
We present a satellite path planning technique able to make identical spacecraft aquire a given configuration. The technique exploits a behaviour-based approach to achieve an autonomous and distributed control over the relative geometry making use of limited sensorial information. A desired velocity is defined for each satellite as a sum of different contributions coming from generic high level behaviours: forcing the final desired configuration the behaviours are further defined by an inverse dynamic calculation dubbed Equilibrium Shaping. We show how considering only three different kind of behaviours it is possible to acquire a number of interesting formations and we set down the theoretical framework to find the entire set. We find that allowing a limited amount of communication the technique may be used also to form complex lattice structures. Several control feedbacks able to track the desired velocities are introduced and discussed. Our results suggest that sliding mode control is particularly appropriate in connection with the developed technique.
The Poisson-Boltzmann model for tRNA
Gruziel, Magdalena; Grochowski, Pawel; Trylska, Joanna
2008-01-01
Using tRNA molecule as an example, we evaluate the applicability of the Poisson-Boltzmann model to highly charged systems such as nucleic acids. Particularly, we describe the effect of explicit crystallographic divalent ions and water molecules, ionic strength of the solvent, and the linear approximation to the Poisson-Boltzmann equation on the electrostatic potential and electrostatic free energy. We calculate and compare typical similarity indices and measures, such as Hodgkin index and root mean square deviation. Finally, we introduce a modification to the nonlinear Poisson-Boltzmann equation, which accounts in a simple way for the finite size of mobile ions, by applying a cutoff in the concentration formula for ionic distribution at regions of high electrostatic potentials. We test the influence of this ionic concentration cutoff on the electrostatic properties of tRNA. PMID:18432617
Navier-Stokes Dynamics by a Discrete Boltzmann Model
NASA Technical Reports Server (NTRS)
Rubinstein, Robet
2010-01-01
This work investigates the possibility of particle-based algorithms for the Navier-Stokes equations and higher order continuum approximations of the Boltzmann equation; such algorithms would generalize the well-known Pullin scheme for the Euler equations. One such method is proposed in the context of a discrete velocity model of the Boltzmann equation. Preliminary results on shock structure are consistent with the expectation that the shock should be much broader than the near discontinuity predicted by the Pullin scheme, yet narrower than the prediction of the Boltzmann equation. We discuss the extension of this essentially deterministic method to a stochastic particle method that, like DSMC, samples the distribution function rather than resolving it completely.
Lattice Boltzmann method simulations of Stokes number effects on particle motion in a channel flow
NASA Astrophysics Data System (ADS)
Zhang, Lenan; Jebakumar, Anand Samuel; Abraham, John
2016-06-01
In a recent experimental study by Lau and Nathan ["Influence of Stokes number on the velocity and concentration distributions in particle-laden jets," J. Fluid Mech. 757, 432 (2014)], it was found that particles in a turbulent pipe flow tend to migrate preferentially toward the wall or the axis depending on their Stokes number (St). Particles with a higher St (>10) are concentrated near the axis while those with lower St (<1) move toward the walls. Jebakumar et al. ["Lattice Boltzmann method simulations of Stokes number effects on particle trajectories in a wall-bounded flow," Comput. Fluids 124, 208 (2016)] have carried out simulations of a particle in a laminar channel flow to investigate this behavior. In their work, they report a similar behavior where particles with low St migrate toward the wall and oscillate about a mean position near the wall while those with high St oscillate about the channel center plane. They have explained this behavior in terms of the Saffman lift, Magnus lift, and wall repulsion forces acting on the particle. The present work extends the previous work done by Jebakumar et al. and aims to study the behavior of particles at intermediate St ranging from 10 to 20. It is in this range where the equilibrium position of the particle changes from near the wall to the axis and the particle starts oscillating about the axis. The Lattice Boltzmann method is employed to carry out this study. It is shown that the change in mean equilibrium position is related to increasing oscillations of the particle with mean position near the wall which results in the particle moving past the center plane to the opposite side. The responsible mechanisms are explained in detail.
Noronha, Jorge; Denicol, Gabriel S.
2015-12-30
In this paper we obtain an analytical solution of the relativistic Boltzmann equation under the relaxation time approximation that describes the out-of-equilibrium dynamics of a radially expanding massless gas. This solution is found by mapping this expanding system in flat spacetime to a static flow in the curved spacetime AdS_{2} Ⓧ S_{2}. We further derive explicit analytic expressions for the momentum dependence of the single-particle distribution function as well as for the spatial dependence of its moments. We find that this dissipative system has the ability to flow as a perfect fluid even though its entropy density does not match the equilibrium form. The nonequilibrium contribution to the entropy density is shown to be due to higher-order scalar moments (which possess no hydrodynamical interpretation) of the Boltzmann equation that can remain out of equilibrium but do not couple to the energy-momentum tensor of the system. Furthermore, in this system the slowly moving hydrodynamic degrees of freedom can exhibit true perfect fluidity while being totally decoupled from the fast moving, nonhydrodynamical microscopic degrees of freedom that lead to entropy production.
Perfect fluidity of a dissipative system: Analytical solution for the Boltzmann equation in AdS2⊗S2
NASA Astrophysics Data System (ADS)
Noronha, Jorge; Denicol, Gabriel S.
2015-12-01
In this paper we obtain an analytical solution of the relativistic Boltzmann equation under the relaxation time approximation that describes the out-of-equilibrium dynamics of a radially expanding massless gas. This solution is found by mapping this expanding system in flat spacetime to a static flow in the curved spacetime AdS2⊗S2 . We further derive explicit analytic expressions for the momentum dependence of the single-particle distribution function as well as for the spatial dependence of its moments. We find that this dissipative system has the ability to flow as a perfect fluid even though its entropy density does not match the equilibrium form. The nonequilibrium contribution to the entropy density is shown to be due to higher-order scalar moments (which possess no hydrodynamical interpretation) of the Boltzmann equation that can remain out of equilibrium but do not couple to the energy-momentum tensor of the system. Thus, in this system the slowly moving hydrodynamic degrees of freedom can exhibit true perfect fluidity while being totally decoupled from the fast moving, nonhydrodynamical microscopic degrees of freedom that lead to entropy production.
Interface-capturing lattice Boltzmann equation model for two-phase flows
NASA Astrophysics Data System (ADS)
Lou, Qin; Guo, Zhaoli
2015-01-01
In this work, an interface-capturing lattice Boltzmann equation (LBE) model is proposed for two-phase flows. In the model, a Lax-Wendroff propagation scheme and a properly chosen equilibrium distribution function are employed. The Lax-Wendroff scheme is used to provide an adjustable Courant-Friedrichs-Lewy (CFL) number, and the equilibrium distribution is presented to remove the dependence of the relaxation time on the CFL number. As a result, the interface can be captured accurately by decreasing the CFL number. A theoretical expression is derived for the chemical potential gradient by solving the LBE directly for a two-phase system with a flat interface. The result shows that the gradient of the chemical potential is proportional to the square of the CFL number, which explains why the proposed model is able to capture the interface naturally with a small CFL number, and why large interface error exists in the standard LBE model. Numerical tests, including a one-dimensional flat interface problem, a two-dimensional circular droplet problem, and a three-dimensional spherical droplet problem, demonstrate that the proposed LBE model performs well and can capture a sharp interface with a suitable CFL number.
Multi-Species Thermal Lattice Boltzmann Models
NASA Astrophysics Data System (ADS)
Wah, Darren; Vahala, George; Vahala, Linda; Pavlo, Pavol; Carter, Jonathan
1998-11-01
Thermal Lattice Boltzmann models (TLBM) are ideal for simulating nonlinear macroscopic conservation systems because of their inherent parallelizeability (nearly all operations are purely local). The TLBM solves a linear BGK-like kinetic equation so that the standard nonlinear convective terms in the standard fluid codes are now replaced by a simple shift operator (linear advection) at the kinetic level. Here we extend our previous TLBM to handle a two-species system, utilizing the models of Morse (1964),Greene (1973) and Kotelnikov & Montgomery (1997). Each kinetic equation now has 2 BGK-like relaxation terms : the first is due to self-collisions and the other is due to different- species collisions. The relaxation rates used are appropriate for electron-ion collisions. Certain constraints can be imposed on the relaxed distribution functions so that the cross-species momentum and energy evolutions relax at the rate determined from the full nonlinear Boltzmann integral collision operator. Ionization and recombination processes will also be examined. Both hexagonal and octagonal lattices are studied.
Reduction of the temperature jump in the immersed boundary-thermal lattice Boltzmann method
NASA Astrophysics Data System (ADS)
Seta, Takeshi; Hayashi, Kosuke; Tomiyama, Akio
2015-11-01
We analytically and numerically investigate the boundary errors computed by the immersed boundary-thermal lattice Boltzmann method (IB-TLBM) with the two-relaxation-time (TRT) collision operator. In the linear collision operator of the TRT, we decompose the distribution function into symmetric and antisymmetric components and define the relaxation parameters for each part. We derive the theoretical relation between the relaxation parameters for the symmetric and antisymmetric parts of the distribution function so as to eliminate the temperature jump. The simple TRT collision operator succeeds in reducing the temperature jump occurring at the high relaxation time in the IB-TLBM calculation. The porous plate problem numerically and analytically demonstrate that the velocity squared terms should be neglected in the equilibrium distribution function in order to eliminate the effect of the advection velocity on the temperature jump in the IB-TLBMs. The passive scalar model without the velocity squared terms more accurately calculates the incompressible temperature equation in the IB-TLBMs, compared to the double distribution model, which is based on the relation of the distribution function gk = (ek - u)2fk / 2 . We apply the passive scalar model without the velocity squared terms to the simulation of the natural convection between a hot circular cylinder and a cold square enclosure. The proposed method adequately sets the boundary values and provides reasonable average Nusselt numbers and maximum absolute values of the stream function.
Podolsky electromagnetism at finite temperature: Implications on the Stefan-Boltzmann law
Bonin, C. A.; Bufalo, R.; Pimentel, B. M.; Zambrano, G. E. R.
2010-01-15
In this work we study Podolsky electromagnetism in thermodynamic equilibrium. We show that a Podolsky mass-dependent modification to the Stefan-Boltzmann law is induced and we use experimental data to limit the possible values for this free parameter.
Boltzmann equation and hydrodynamic fluctuations.
Colangeli, Matteo; Kröger, Martin; Ottinger, Hans Christian
2009-11-01
We apply the method of invariant manifolds to derive equations of generalized hydrodynamics from the linearized Boltzmann equation and determine exact transport coefficients, obeying Green-Kubo formulas. Numerical calculations are performed in the special case of Maxwell molecules. We investigate, through the comparison with experimental data and former approaches, the spectrum of density fluctuations and address the regime of finite Knudsen numbers and finite frequencies hydrodynamics. PMID:20364972
NASA Astrophysics Data System (ADS)
Braziewicz, J.; Polasik, M.; Słabkowska, K.; Majewska, U.; Banaś, D.; Jaskóła, M.; Korman, A.; Kozioł, K.; Kretschmer, W.; Choinski, J.
2010-08-01
In the present work, an alternative approach for the evaluation of the equilibrium K-, L-, and M-shell ionizations and the mean charge state q¯ for projectiles passing through various targets has been proposed. The approach is based on measured K x-ray energy shifts and line intensity ratios and utilizes the theoretical analysis of projectile spectra using multiconfiguration Dirac-Fock calculations. It was applied for the satellite and hypersatellite K lines in the x-ray spectra emitted by sulfur projectiles passing with energies of 9.6-122 MeV through carbon, aluminium, titanium, and iron targets, recorded by a Si(Li) detector. It was found that only in the high projectile energy region there was significant dependence of mean equilibrium K-shell ionization on the target atomic number. The equilibrium L-shell ionization rises with the increase of sulfur energy until 32 MeV, but for higher energies the changes are very weak. The equilibrium M-shell ionization changes very weakly for low projectile energy while for higher energies this ionization is practically constant. For each target, the estimated value of q¯ rises with the increase of the sulfur energy value. The dependence of the sulfur charge state on the target atomic number was discussed by taking into account the cross sections for ionization, decay, and electron capture processes. The data were compared with the experimental data measured by other authors and with the predictions of Shima’s and Schiwietz and Grande’s semiempirical formulas. The presented good agreement points out that this alternative approach delivers quantitative results.
Braziewicz, J.; Majewska, U.; Banas, D.; Polasik, M.; Slabkowska, K.; Koziol, K.; Jaskola, M.; Korman, A.; Kretschmer, W.; Choinski, J.
2010-08-15
In the present work, an alternative approach for the evaluation of the equilibrium K-, L-, and M-shell ionizations and the mean charge state q for projectiles passing through various targets has been proposed. The approach is based on measured K x-ray energy shifts and line intensity ratios and utilizes the theoretical analysis of projectile spectra using multiconfiguration Dirac-Fock calculations. It was applied for the satellite and hypersatellite K lines in the x-ray spectra emitted by sulfur projectiles passing with energies of 9.6-122 MeV through carbon, aluminium, titanium, and iron targets, recorded by a Si(Li) detector. It was found that only in the high projectile energy region there was significant dependence of mean equilibrium K-shell ionization on the target atomic number. The equilibrium L-shell ionization rises with the increase of sulfur energy until 32 MeV, but for higher energies the changes are very weak. The equilibrium M-shell ionization changes very weakly for low projectile energy while for higher energies this ionization is practically constant. For each target, the estimated value of q rises with the increase of the sulfur energy value. The dependence of the sulfur charge state on the target atomic number was discussed by taking into account the cross sections for ionization, decay, and electron capture processes. The data were compared with the experimental data measured by other authors and with the predictions of Shima's and Schiwietz and Grande's semiempirical formulas. The presented good agreement points out that this alternative approach delivers quantitative results.
A note on a Discrete Boltzmann Equation with multiple collisions
NASA Astrophysics Data System (ADS)
Oliveira, Filipe; Soares, Ana Jacinta
2008-05-01
We compute a non-trivial explicit solution for the one-dimensional plane 6-velocity discrete Boltzmann model with multiple collisions introduced in [E. Longo, R. Monaco, On the discrete kinetic theory with multiple collisions: Plane six-velocity and unsteady Couette flow, in: Muntz, et al. (Eds.), The Proceedings of Rarefied Gas Dynamics, in: AIAA Publ., vol. 118, 1989, pp. 118-130] which asymptotically connects two particular equilibrium states. We prove that such a solution exists provided that a suitable condition on the differential elastic cross sections holds.
NASA Astrophysics Data System (ADS)
Zhang, Jianying; Yan, Guangwu
2016-04-01
A lattice Boltzmann model for solving the (2+1) dimensional cubic-quintic complex Ginzburg-Landau equation (CQCGLE) is proposed. Different from the classic lattice Boltzmann models, this lattice Boltzmann model is based on uniformly distributed lattice points in a two-dimensional space, and the evolution of the model is about a spatial axis rather than time. The algorithm provides advantages similar to the lattice Boltzmann method in that it is easily adapted to complex Ginzburg-Landau equations. Numerical results reproduce the phenomena of the fusion of necklace-ring pattern and the effect of non-linearity on the soliton in the CQCGLE.
NASA Astrophysics Data System (ADS)
Zhang, Jianying; Yan, Guangwu
2015-12-01
A spatiotemporal lattice Boltzmann model for solving the three-dimensional cubic-quintic complex Ginzburg-Landau equation (CQCGLE) is proposed. Different from the classic lattice Boltzmann models, this lattice Boltzmann model is based on uniformly distributed lattice points in a three-dimensional spatiotemporal space, and the evolution of the model is about a spatial axis rather than time. The algorithm possesses advantages similar to the lattice Boltzmann method in that it is easily adapted to complex Ginzburg-Landau equations. Examples show that the model reproduces the phenomena in the CQCGLE accurately.
NASA Astrophysics Data System (ADS)
Feldman, Michal; Tennenholtz, Moshe
We introduce partition equilibrium and study its existence in resource selection games (RSG). In partition equilibrium the agents are partitioned into coalitions, and only deviations by the prescribed coalitions are considered. This is in difference to the classical concept of strong equilibrium according to which any subset of the agents may deviate. In resource selection games, each agent selects a resource from a set of resources, and its payoff is an increasing (or non-decreasing) function of the number of agents selecting its resource. While it has been shown that strong equilibrium exists in resource selection games, these games do not possess super-strong equilibrium, in which a fruitful deviation benefits at least one deviator without hurting any other deviator, even in the case of two identical resources with increasing cost functions. Similarly, strong equilibrium does not exist for that restricted two identical resources setting when the game is played repeatedly. We prove that for any given partition there exists a super-strong equilibrium for resource selection games of identical resources with increasing cost functions; we also show similar existence results for a variety of other classes of resource selection games. For the case of repeated games we identify partitions that guarantee the existence of strong equilibrium. Together, our work introduces a natural concept, which turns out to lead to positive and applicable results in one of the basic domains studied in the literature.
Lattice-Boltzmann Simulation of Tablet Disintegration
NASA Astrophysics Data System (ADS)
Jiang, Jiaolong; Sun, Ning; Gersappe, Dilip
Using the lattice-Boltzmann method, we developed a 2D model to study the tablet disintegration involving the swelling and wicking mechanisms. The surface area and disintegration profile of each component were obtained by tracking the tablet structure in the simulation. Compared to pure wicking, the total surface area is larger for swelling and wicking, which indicates that the swelling force breaks the neighboring bonds. The disintegration profiles show that the tablet disintegrates faster than pure wicking, and there are more wetted active pharmaceutical ingredient particles distributed on smaller clusters. Our results indicate how the porosity would affect the disintegration process by changing the wetting area of the tablet as well as by changing the swelling force propagation.
Boltzmann babies in the proper time measure
Bousso, Raphael; Bousso, Raphael; Freivogel, Ben; Yang, I-Sheng
2007-12-20
After commenting briefly on the role of the typicality assumption in science, we advocate a phenomenological approach to the cosmological measure problem. Like any other theory, a measure should be simple, general, well defined, and consistent with observation. This allows us to proceed by elimination. As an example, we consider the proper time cutoff on a geodesic congruence. It predicts that typical observers are quantum fluctuations in the early universe, or Boltzmann babies. We sharpen this well-known youngness problem by taking into account the expansion and open spatial geometry of pocket universes. Moreover, we relate the youngness problem directly to the probability distribution for observables, such as the temperature of the cosmic background radiation. We consider a number of modifications of the proper time measure, but find none that would make it compatible with observation.
Boltzmann babies in the proper time measure
NASA Astrophysics Data System (ADS)
Bousso, Raphael; Freivogel, Ben; Yang, I.-Sheng
2008-05-01
After commenting briefly on the role of the typicality assumption in science, we advocate a phenomenological approach to the cosmological measure problem. Like any other theory, a measure should be simple, general, well defined, and consistent with observation. This allows us to proceed by elimination. As an example, we consider the proper time cutoff on a geodesic congruence. It predicts that typical observers are quantum fluctuations in the early universe, or Boltzmann babies. We sharpen this well-known youngness problem by taking into account the expansion and open spatial geometry of pocket universes. Moreover, we relate the youngness problem directly to the probability distribution for observables, such as the temperature of the cosmic background radiation. We consider a number of modifications of the proper time measure, but find none that would make it compatible with observation.
Harleston, H. Instituto de Ciencias Nucleares, Universidad Nacional Autonoma de Mexico, Distrito Federal ); Vishniac, E.T. )
1992-06-15
The Arnowitt-Deser-Misner formalism is used to write the Einstein-Boltzmann coupled system of equations. The sources of gravitational field are represented by ordinary matter described by a perfect-fluid approximation together with a particle gas described by a phase-space distribution function obeying the general-relativistic Boltzmann transport equation. Through the use of the Liouville operator in phase space, we obtain a form of the Boltzmann equation that makes it very amenable for numerical treatment. The resulting system of equations can be used for the numerical study of either massless or massive particles interacting with ordinary matter.
Equilibrium distribution of dissolved sulphur species in water at 25°C and 1 atm total pressure
Garrels, R.M.; Naeser, C.R.
1958-01-01
The Eh-pH diagrams for the equilibrium concentrations in aqueous solution at 25°C of native sulphur and all the various sulphur-containing ions and acids from which the ions are produced have been constructed for systems having a total sulphur concentration of 0.1 molar. The composite of these diagrams indicates that elemental sulphur, H2S, HS− HSO4− and SO4 are the species that predominate in the environments that might be found in nature. This indication is in agreement with the composition of all sulphur-containing minerals.
Niiyama, Tomoaki; Shimizu, Yasushi; Kobayashi, Taizo R; Okushima, Teruaki; Ikeda, Kensuke S
2009-05-01
We investigate numerically and analytically the effects of conservation of total translational and angular momentum on the distribution of kinetic energy among particles in microcanonical particle systems with small number of degrees of freedom, specifically microclusters. Molecular dynamics simulations of microclusters with constant total energy and momenta, using Lennard-Jones, Morse, and Coulomb plus Born-Mayer-type potentials, show that the distribution of kinetic energy among particles can be inhomogeneous and depend on particle mass and position even in thermal equilibrium. Statistical analysis using a microcanonical measure taking into account of the additional conserved quantities gives theoretical expressions for kinetic energy as a function of the mass and position of a particle with only O(1/N;{2}) deviation from the Maxwell-Boltzmann distribution. These expressions fit numerical results well. Finally, we propose an intuitive interpretation for the inhomogeneity of the kinetic energy distributions. PMID:19518410
Spectrum structure and behavior of the Vlasov-Maxwell-Boltzmann system without angular cutoff
NASA Astrophysics Data System (ADS)
Huang, Yongting
2016-02-01
The spectrum structure and behavior of the Vlasov-Maxwell-Boltzmann (VMB) system with physical angular non-cutoff intermolecular collisions are studied in this paper. The analysis shows the effect of the Lorentz force induced by the electro-magnetic field leads to some different spectrum structure from the non-cutoff Boltzmann equation. The spectrum structure in high frequency, quite different from the VMB system with angular cutoff assumption, also illustrates the hyperbolic structure of the Maxwell equation. Furthermore, the large time behavior and optimal convergence rates to the equilibrium of the non-cutoff VMB system are established on the spectrum analysis.
Yoshikawa, S.
1981-08-01
A straight, helical plasma equilibrium equation is solved numerically for a plasma with a helical magnetic axis. As is expected, by a suitable choice of the plasma boundary, the vacuum configuration is made line ..integral.. dl/B stable. As the plasma pressure increases, the line ..integral.. dl/B criterion will improve (again as expected). There is apparently no limit on the plasma ..beta.. from the equilibrium consideration. Thus helical-axis stellarator ..beta.. will presumably be limited by MHD stability ..beta.., and not by equilibrium ..beta...
Lattice Boltzmann simulation of turbulence-induced flocculation of cohesive sediment
NASA Astrophysics Data System (ADS)
Zhang, Jin-Feng; Zhang, Qing-He; Maa, Jerome P.-Y.; Qiao, Guang-Quan
2013-10-01
Both the floc formation and floc breakup of cohesive sediment are affected by turbulent shear which is recognized as one of the most important parameters, and thus, on the settling and transport of cohesive sediment. In this study, the development of floc characteristics at early stage and steady-state of flocculation were investigated via a three-dimensional lattice Boltzmann numerical model for turbulence-induced flocculation. Simulations for collision and aggregation of various size particles, floc growth, and breakup in isotropic and homogenous turbulent flows with different shear stresses were conducted. Model results for the temporal evolution of floc size distribution show that the normalized floc size distributions is time-independent during early stage of flocculation, and at steady-state, shear rate has no effect on the shape of normalized floc size distribution. Furthermore, the size, settling velocity, and effective density of flocs at the non-equilibrium flocculation stage do not change significantly for shear stresses in the range 0-0.4 N m-2. The relationships between floc size and settling velocity established during floc growth stages and that during steady-states are different.
Lattice Boltzmann methods in Geosciences
NASA Astrophysics Data System (ADS)
Huber, Christian; Parmigiani, Andrea; Su, Yanqing
2014-05-01
Numerical models often offer the only possible approach to study the complex non-linear dynamics of geodynamical processes that are difficult or impossible to scale for laboratory experiments. The development of improved computer resources has allowed the emergence of large-scale parallel computations in Earth Sciences. These resources have lead to an increasing complexity in models where a greater number of adjustable parameters arise. Although the increasing number of free parameters offers a greater flexibility to fit satisfyingly the set of available constraints (e.g. geochemical, structural) it also provides new challenges in terms of the size of the parameter space and non-uniqueness of model solutions. Another significant challenge associated with state-of-the-art models is that their complexity is in general associated with the addition of parameterizations of the unresolved (small) scale processes. This trend calls for the development of complementary high-performance models to constrain the physics at small-scales where mass, momentum and energy exchanges at interfaces between different phases control the dynamics in heterogeneous media. We argue that more attention should be devoted to the development of multiphase numerical modeling at the granular (pore) scale to investigate the dynamical behavior of heterogeneous media and the emergence of feedbacks that influence the response of these media at much greater scales. The lattice Boltzmann method is a paradigm that emerged almost three decades ago. It is based on kinetic theory and follows a bottom-up approach that contrast the top-down strategy of standard methods such as Finite Volumes, FEM and Finite Differences. Lattice Boltzmann is ideally suited to handle the complex dynamics of multiphase systems at small spatial scales and is very efficient for parallel programing. In this presentation, we discuss the development of different lattice Boltzmann models developed in our group over the last years
NASA Astrophysics Data System (ADS)
Louis-Martinez, Domingo
2011-04-01
A classical (non-quantum-mechanical) relativistic ideal gas in thermodynamic equilibrium in a uniformly accelerated frame of reference is studied using Gibbs's microcanonical and grand canonical formulations of statistical mechanics. Using these methods explicit expressions for the particle, energy and entropy density distributions are obtained, which are found to be in agreement with the well known results of the relativistic formulation of Boltzmann's kinetic theory. Explicit expressions for the total entropy, total energy and rest mass of the gas are obtained. The position of the center of mass of the gas in equilibrium is found. The non-relativistic and ultrarelativistic approximations are also considered. The phase space volume of the system is calculated explicitly in the ultrarelativistic approximation.