Tests of a convergent numerical scheme for nonlinear Boltzmann-like models with chemical reactions: two, three and four species
|
|
- Patrick Manning
- 6 years ago
- Views:
Transcription
1 Tests of a convergent numerical scheme for nonlinear Boltzmann-like models with chemical reactions: two, three and four species Dorin Marinescu and Cecil P. Grünfeld Abstract We review recent tests of a convergent numerical scheme for space-homogeneous, nonlinear Boltzmann-like equations for fluids with two, three and four species of gas particles, respectively, with elastic and reactive binary collisions. The tests refer both to simpler exactly solvable models with elastic collisions, and to more complicated examples, involving chemical reactions, which present only exact expressions for certain low-order moments of the one-particle distribution functions. The results of the numerical tests appear to be in very good agreement with theoretical predictions. Mathematics Subject Classification: 45K5, 76M5, 76P5, 65R, 8C4. Key words: numerical scheme, numerical test, theoretical prediction, non-linear Boltzmann-like equations, energy, operators, Dirac measures. 1 Introduction An accurate, computational method has been recently introduced [7], in order to solve numerically nonlinear Boltzmann-like equations [4], [5], [6] for reacting gas mixtures. The method represents a generalization of a rigorous numerical scheme proposed in [1], [], which combines analytical and stochastic techniques to obtain numerical solutions of the classical Boltzmann equation. As compared to the method of [1], [], the numerical techniques developed in [7] have to face increased mathematical difficulties, essentially due to the presence of several species of gas particles and occurrence of reactions. The usefulness of models with exact solutions for numerical tests is well known. In this respect, the simplest example is represented by the Krook-Wu model [8], corresponding to a system of Boltzmann equations describing nonreactive gas mixtures. The numerical tests on the Krook-Wu equations have the merit to provide some insight in the effectiveness of the numerical method in dealing with models with Proceedings of The 3-rd International Colloquium Mathematics in Engineering and Numerical Physics October 7-9, 4, Bucharest, Romania, pp c Balkan Society of Geometers, Geometry Balkan Press 5.
2 Tests of a convergent numerical scheme 161 nonlinearities, and in handling several unknown distribution functions, corresponding to several species. The scheme of [7] was tested in [1] on the exact Krook-Wu solutions of the multicomponent Boltzmann equation [8], in order to study the effectiveness of the method in applications to gas mixtures. However, the Krook-Wu equations refer only to non-reacting gases. Consequently, the results of [1] are of rather limited interest. In fact, they do not provide information on the applicability of the scheme of [7] to models with chemical reactions. More realistic evaluations, based on the confrontation between numerical and analytical results, should require the knowledge of non-trivial, exact solutions of the non-linear Boltzmann-like equations for reacting flows. So far, analytical solutions to Boltzmann models with chemical reactions, have been obtained only for linear equations [1], or equations with collision terms defined by trivial (reaction) collision laws [3]. Unfortunately, such exact, simple solutions do not serve the purposes of the aforementioned evaluations. On the other hand, one could perform useful, indirect tests on certain non-trivial Boltzmann models for reacting flows. Although the models cannot be solved analytically, they provide exactly solvable equations for the evolution of certain macroscopic variables (e.g. concentration, energy). In that case, the one-particle distribution functions of the gas species can be obtained, solving numerically the Boltzmann model, by means of the techniques developed in [7]. Then the aforementioned macroscopic quantities (at various moments) can be calculated as averages with respect to the numerical one-particle distribution functions. The resulting values can be compared with those provided by the analytical solutions of the equations for macroscopic quantities. Moreover, the computed values of the macroscopic variables can serve to study the accuracy of the numerical scheme concerning the verification of the bulk conservation properties [4], [5], [6] of the model. Results of several tests on Boltzmann-like models with chemical reactions have been recently published in [11]. In this paper we review comparatively some of the results of the above numerical tests, and also present new results concerning the situation when exothermal reactions involving three species are included in the models investigated in [11]. Our aim is to give a first answer to the important question concerning the sensitivity of the method with respect to the change of the nature of the binary interactions and the number of components of the model, the main conclusion being that the accuracy of the numerical seems not to be affected by such changes (in the regime of reactions involving few particles). Analytical Considerations In this paper, we are interested in the kinetics of space-homogeneous gas mixtures consisting of N 4 chemical species, of gas particles with binary elastic collisions and binary reactions. The gas reactions conserve the total mass, momentum and energy of the reaction partners in the reaction channels [4]. We assume that each species X i is composed of point particles with mass m i,and one-internal state characterized by internal energy E i, 1 i 4. We consider the situation when the following kinds of gas interactions are possible:
3 16 D. Marinescu and C. P. Grünfeld (a) nonreactive (elastic) collisions (.1) X i + X j X i + X j, for 1 i, j 4 and {i, j} = {3, 4}; (b) reactions of the form (.) X 1 + X X 3 + X 4 and the reverse reaction (.3) X 3 + X 4 X 1 + X Notice that reactions of the form X 1 + X X 3 can be considered as particularization of (.), (.3) to the case when X 3 = X 4. To be precise, without loss of generality, we suppose reaction (.) to be endothermal (E 1 + E E 3 + E 4 ). This means that reaction (.3) is exothermal. In this case, the weak form of the Boltzmann-like model introduced in [5], reads (.4) ( φ, f ) i = (φ, P i (f)) (φ, S i (f)), ( ) 1 i 4; φ C b (R 3 ). t for the unknown one-particle distribution functions f i = f i (t, v) (depending on t and velocity v ) and associated to the species X i, 1 i 4. Here f := (f 1,..., f 4 ) (.5) (φ, P i (f)) = j,k,l=1 D kl;ij S φ(v ij;kl ) r ij;kl (v, w, n) f k (v)f l (w) dvdwdn, and (.6) (φ, S i (f)) = φ(v) r kl;ij (v, w, n) f i (v)f j (w) dvdwdn, j,k,l=1 D ij;kl S where C b (R 3 ) denotes the space of continuous, bounded, real functions on R 3, S := {n R 3 n = 1}, (.7) with (.8) (.9) D ij;kl := {(v, w) R 3 R 3 Wij;kl (v, w) }, W ij;kl (v, w) := m i m j (m i + m j ) v w + E i + E j E k E l, v ij;kl = v ij;kl (v, w, n) := m iv + m j w m i + m j + + [ ] 1/ m l m k (m i + m j ) W ij;kl (v, w) n,
4 Tests of a convergent numerical scheme 163 Finally the collision law r ij;kl (v, w, n) characterizes the reaction X i + X j X k + X l. According to the hypotheses of the model, in formulae (.4), (.5) and (.6), only the collision laws associated to reactions (.1), (.), (.3) do not vanish identically. It is known [5], [7] that the Cauchy problem associated to (.4) has, in some sense, unique positive global solutions (for details the reader is referred to [5], [6], [7]). Following [1], [11], we consider collision laws similar to those introduced by Krook and Wu for nonreactive gaseous mixtures [8]. Assumption: (.1) const when W1;34 (v, w, n) <, r 1;1 (v, w, n) = if W1;34 (v, w, n) ; the functions r 1;34, r 34;1 are nonnegative constants (on their domains); similarly r ij;ij is nonnegative and constant if {i, j} = {1, } or {i, j} = {3, 4}. Under hypothesis (.1), isotropic initial conditions lead to isotropic solutions, i.e., f i (v) = f i (v) for 1 i 4 (where v := v is the modulus of the velocity v). In that case, one can provide a more explicit form for Eq.(.4). Indeed, define (.11) F i (v) := 4πv f i (v), for 1 i 4. Denote the concentration of the species i by (.1) I i := F i (v)dv, for 1 i 4. Let λ kl;ij := 4π r ij;kl (1 i, j, k, l 4). Using suitable changes of variables in (.5), (.6), one obtains (.13) (φ, P i (f)) = λ kl;ij I k I l φ(ṽ) F k (v)f l (w) dv dw dζ dη, D kl;ij j,k,l=1 (.14) (φ, S i (f)) = λ ij;kl I i I j φ(v) F i (v)f j (w) dv dw dζ dη, D ij;kl j,k,l=1 where (.15) D kl;ij := {(v, w, ζ, η) R + [, 1) ṽ kl;ij R} The the post collision velocity of the species i, ṽ kl;ij is given by (.16) ṽ kl;ij = ṽ kl;ij (v, w, ζ, η) := [ V kl + gm k m l ρ kl;ij + ( mk m l ) 1/ ρ kl;ijv kl(η 1) ] 1/,
5 164 D. Marinescu and C. P. Grünfeld where (.17) V kl = V kl (v, w, ζ) := (m kv + m l w + m k m l vw(ζ 1)) 1/ /(m k + m l ), (.18) ρ kl;ij = ρ kl;ij (v, w, ζ) := [m k m l (v + w + vw(ζ 1))/(m k + m l ) + +(E k + E l E i E j )/(m k + m l )] 1/. The simple particularization of the above model concerns the non-reactive gas with two species [1]: N =, E i = and λ kl;ij =, for k i and l j. Let λ ij;ij = λ ij, m ik = m i /(m i + m k ) (.19) p 1 = λ λ 1 µ(3 µ), p = λ 11 λ 1 µ(3 µ), with µ = 4m 1 m /(m 1 + m ). The constants λ ki satisfy the condition (.) λ 1 /I 1 = λ 1 /I, The Krook-Wu like solutions [8] are determined if either condition on parameters p 1 = p or condition µ (λ 1 /p 1 λ 1 /p ) = 1 is fulfilled. Then the general solution has the form: (.1) f i (t, v) = Φ(v; α i (t))[1 3p i R(t) + Q i (t)v ], i = 1,, where (.) Φ(v; α) = (πα) 3/ exp( v /α), Q i (t) = m i p i R(t)/ξ(t), α i (t) = ξ(t)/m i, R(t) = A/[A exp(a(t t )) B], ξ(t) = (I 1 + I )/[I 1 + I + 8π(I 1 p 1 + I p ) R(t)], with A and B constants defined by: (.3) A = 4π[λ 11 + λ 1 µ(3 µp /p 1 )]/6, B = 16π [λ 11 p 1 + λ 1 µ(3 µ)p ]/3. On the other hand, if reactions are present (under condition (.1)), no nontrivial analytical solution of the model is known. However, as mentioned in the Introduction, in some particular cases, indirect tests can be done, because starting from (.4) one can solve exactly the equations for chemical concentrations and energies. Indeed, the internal energy can be expressed exactly in terms of concentrations. (.4) E int (t) = I i (t)e i. i=1
6 Tests of a convergent numerical scheme 165 By the global conservation of energy, the total (kinetic + internal) energy of the gas (.5) E(t) = i=1 m i v i (t) + I i (t)e i = E(), From (.5) the total kinetic energy of the particles can be written (.6) (.7) E cin (t) = i=1 m i v i (t) = i=1 i=1 m i v i () + i=1[i i I i (t)]e i. Moreover, if species X 1 and X interact only by elastic collisions, then I i = λ I 3 I 4, i = 1, I i = λ I 3 I 4, i = 3, 4 with λ := λ 34;1. These equations can be solved immediately. Let Ii := I i(), 1 i 4. If I3 I4 then, the solutions of (.7) are I i (t) = Ii + I 3 I4 [exp (I4 λt) exp (I3 λt) ] I4 exp (I 4 λt) I 3 exp, i = 1, (I 3 λt) (.8) Ii I i (t) = (I 4 I3 ) exp (Ii λt) I4 exp (I 4 λt) I 3 exp, i = 3, 4, (I 3 λt) and when I3 = I4, the solutions of (.7) are I i (t) = Ii + (I 3 ) λt (1 + I, i = 1, 3 λt) (.9) Ii I i (t) = (1 + Ii, i = 3, 4. λt) 3 Numerical Method and Results One considers some interval [, T ] and some given step < t < T. The -discretized version of (.4) has the following form (3.3) (φ, F p+1 i ) = (φ, Q p ij ), where {i, j} = {1, } or {i, j} = {3, 4}, p {, 1,..., [[T/ t]] 1} ([[x]] denoting the integer part of the real positive number x) and φ C b (R 3 ). The operators Q p ij in (3.3) are defined through the following measures. For {i, j} = {1, } and {k, l} = {3, 4} let (3.31) dh ij (v, w, ζ, η) := 1 J ij F i (v)f j (w) dv dw dζ dη on D ij;kl, dh ij (v, w, ζ, η) := 1 J ij F i (v)f j (w) dv dw dζ dη on R + [, 1] \ D ij;kl.
7 166 D. Marinescu and C. P. Grünfeld where (3.3) J ij := J ij := F i (v)f j (w) dv dw dζ dη, D ij;kl F i (v)f j (w) dv dw dζ dη. R + [,1] \D ij;kl For {i, j} = {1, }, let (3.33) dh ij (v, w, ζ, η) := 1 I i I j F i (v)f j (w) dv dw dζ dη on R + [, 1). The operators Q p ij have the following form. (3.34) [ ( ) ] (φ, Q p ij ) := mk M p λ ik;ik t I p k λ ii;ii t I p i I p i (φ, F p i )+ k=1 tot [( ) mj + M p λ ij;ij t J p ij (φ, Hp ij )+ tot ( ) ] mj + λ ij;34 t J p ij(φ, H p ij) I p j Ip i + M p tot [ ] + λ ii;ii ( φ ii;ii, H p ii )Ip i + λ ik;ik ( φ ik;ik, H p ik )Ip k t I p i + k=1 [ + λ ij;ij J p ij ( φ ij;ij, H p ij )Ip j Ip i + λ 34;ij( φ ] 34;ij, H p 34 )Ip 3 Ip 4 t; when {i, j} = {1, } and (φ, Q p ij ) := [ k=1 ( mk M p tot ) ] λ ik;ik t I p k λ ii;ii t I p i I p i (φ, F p i )+ (3.35) [ ] + λ ii;ii ( φ ii;ii, H p ii ) Ip i + λ ik;ik ( φ ik;ik, H p ik ) Ip k t I p i + k=1 +λ 1;ij J p 1( φ 1;ij, H p 1) t I p 1 Ip. if {i, j} = {3, 4}. Note that, in (3.34) and (3.35) the quantities I p i, Hp ij, Hp ij, J p ij, J p ij, H p ij are defined as in (.1), (3.31), (3.3), (3.33), simply replacing F i and F j by F p i and F p j
8 Tests of a convergent numerical scheme 167 respectively. Moreover, (3.36) M p tot := m i I p i i=1 and φ ij;kl := φ ṽ ij;kl, with ṽ ij;kl given by (.16). We approximate the initial data Fi (v)dv by sums of Dirac measures (3.37) dµ i,n (v) = 1 N N δ(v vi,n)dv n=1 concentrated on N points vi,n, (such that dµ i,n (v) converges to F i (v)dv as N in the weak sense of measures) and we approximate the Lebesgue measure on the unit square dζdη by sums of Dirac measures concentrated on N points 1 N N n=1 δ(ζ ζ n )δ(η η n )dζdη. As one knows [7], the product measures in the r.h.s. of (3.3), yields a power-like increasing computational effort. The decrease of the computational effort (preserving the convergence of the scheme) is accomplished by random selections, applying Theorems 7, 8 of [7]. In all our tests we used the mixed congruential method to generate sequences {ω n } n N of pseudo-random numbers. The elements ω n = z n /b, where z n are given recursively by z n = λz n 1 +r (modb). In this relation b > 1, λ and r are fixed natural numbers and λ is relative prime with b. The initialization is made with some integer z < b. Here b = 3 1 3, r = , λ = Each test starts with an arbitrary positive z < b. The approximation of the initial data by sums of Dirac measures was made by means of the Hamersley-Van der Corput sequences (see [9] for details). The numerical results correspond to several tests as follows. In the case of the model with elastic collisions which presents exact solutions, a typical result can be found in [1]: m 1 = 5, m = 1, n 1 =, t = ln(3p + B/A)/A, T = 6. (This choice of t ensures the positivity of f i given by (.1) for each t.) In a first experiment we have obtained for each i = 1, ten numerical solutions for N = 15 approximation points, at T = 6 and J = 16 iteration steps. Denote G i (t, v) := 4π v u f i (t, u)du, for i = 1,. The numerical results are represented in Figure 1 and Figure. In the case of the model with chemical reactions, our experiments consisted of three series of tests performed on the model introduced in Section, corresponding to the caricatures of three kinds of chemical reactions (the below quantities are expressed in conventional, dimensionless units): Case A: three distinct chemical species without endothermal reactions (only processes (.1) and (.3) with X 3 = X 4 occur) e.g. gas mixtures consisting of N, O and NO, at sufficiently low temperature, undergoing binary elastic collisions, as well as reactive collision inducing the reverse reaction in the process N + O NO. Masses : m 1 = , m = , m 3 = m 4 = Initial concentrations: I 1 =, I =.1, I 3 = I 4 =.45.
9 168 D. Marinescu and C. P. Grünfeld 1, G 1 1, G,8,8,6 exact solution,6 exact solution,4,4,,,,,5,5,75 1, 1,5 1,5 1,75, v,,,5 1, 1,5,,5 3, 3,5 4, v Figure 1: Exact solution of the first specie and ten numerical solutions for N = 15, at T = 6 and J = 16. Figure : Exact solution of the second specie and ten numerical solutions for N = 15, at T = 6 and J = 16. Internal energies: E 1 = , E = , E 3 = E 4 =.97. Values of λ: λ 11;11 = 1.56, λ 1;1 = λ 33;33 = λ 34;1 = λ 44;44 = 1.334, λ 1;34 =, λ 13;13 = λ 14;14 = 1.69, λ ; = 1.459, λ 3;3 = λ 4;4 = Simulation T = 16. Initialization function in the -discretized equations: (3.38) F i (v) = 4πv I i ( mi ) 3/ ( exp mi v ), 1 i 4 π Case B: four distinct species, without endothermic reaction (only processes (.1), (.3) occur and X 1,...,X 4 are different) e.g. gas mixtures consisting of HF, HCl, F Cl and F, at sufficiently low temperature, undergoing binary elastic collisions, as well as reactive collision inducing the direct reaction in the process HF + F Cl F + HCl. Masses: m 1 = , m = , m 3 = m 4 = Initial concentrations: I 1 =., I =.8, I 3 =.4, I 4 =.5. Internal energies: E 1 = 4.475, E = 9.364, E 3 = 4.764, E 4 =.751. Values of λ: λ 11;11 = 1.58, λ 1;1 = λ 34;1 =.173, λ 1;34 =, λ 13;13 = 1.691, λ 14;14 = , λ ; =.9564, λ 3;3 =.458, λ 4;4 =.655, λ 33;33 = 1.96, λ 44;44 =.547. Simulation T = 8.5. Initialization function in the -discretized equations is the same as in Case A, (3.38) Case C: four distinct species (.1), (.), (.3) occur and X 1,...,X 4 are different) e.g. gas mixtures consisting of BrF, F Cl, BrCl, F with elastic collisions and all reactions of the form in the process BrF + F Cl F + BrCl. Masses: m 1 = , m = , m 3 = , m 4 = Initial concentrations: I 1 =.35, I =.6, I 3 =.5, I 4 =. Internal energies: E 1 = , E = , E 3 =.38, E 4 = Values of λ: λ 11;11 = 3.6, λ 1;1 = λ 34;1 = 3.79, λ 1;34 = 3.79, λ 13;13 =.781, λ 14;14 = 4.1, λ ; =.9564, λ 3;3 =.458, λ 4;4 = , λ 33;33 = 1.96, λ 44;44 = Simulation T = 1.75.
10 Tests of a convergent numerical scheme 169 (3.39) Initialization function in the -discretized equations: F i (v) = 4πv I i ( mi π ) 3/ exp ( m iv ), 1 i total concentration.8 1 kinetic energy O 5 concentration.6.4 N energy total energy -5 internal energy. NO Figure 3: The evolution of the concentrations in Case A for 46 points and 85 iteration steps Figure 4: The evolution of the energies in Case A for 46 points and 85 iteration steps total concentration.8 1 kinetic energy FCl 5 concentration.6.4 HF energy total energy -5 HCl internal energy. F Figure 5: The evolution of the concentrations in Case B for 46 points and 85 iteration steps Figure 6: The evolution of the energies in Case B for 46 points and 85 iteration steps. The support of the initial datum is included in [, ] when Fi is given by (3.38). It is included in [, 6] if Fi is given by (3.39). For details concerning the approximation of the initial data (3.38), (3.39), the interested reader is referred to [9], [1]. The above results present the evolution of concentrations and energies for several tests. The numerical results are indicated by dots. Each dot corresponds to a number given as the arithmetic mean of the values obtained as results of six simulations (corresponding to identical physical conditions). In addition to the numerical results, the values of the known exact solutions of the equations for the above macroscopic quantities are represented as lines, in the aforementioned figures.
11 17 D. Marinescu and C. P. Grünfeld 1 15 total concentration.8 1 kinetic energy 5 concentration.6.4 FCl BrF energy total energy -5 internal energy. BrCl Figure 7: The evolution of the concentrations in Case C for 46 points and 55 iteration steps. F Figure 8: The evolution of the energies in Case C for 46 points and 55 iteration steps total concentration.8 1 kinetic energy O 5 concentration.6.4 N energy total energy -5 internal energy. NO Figure 9: The evolution of the concentrations in Case A for 46 points and 9 iteration steps with non-constant step Figure 1: The evolution of the energies in Case A for 46 points and 9 iteration steps with non-constant step total concentration.8 1 kinetic energy FCl 5 concentration.6.4 HF energy total energy -5 HCl internal energy. F Figure 11: The evolution of the concentrations in Case B for 46 points and 41 iteration steps Figure 1: The evolution of the energies in Case B for 46 points and 41 iteration steps.
12 Tests of a convergent numerical scheme 171 Moreover, Figures 4, 6, 8, 1 and 1 detail the values of the kinetic energy, calculated by two different methods. Specifically, the first method evaluates the kinetic energy (kinetic energy 1, represented in our figures by circles) at each iteration step, as difference between the total energy (at t = ) and the numerical value of internal energy (expressed in terms of concentrations as in (.6)). The second method yields the kinetic energy (kinetic energy, represented in our figures by points) as an average with respect to the one-particle distribution functions. Figures 3-6 present the evolution of concentrations and energies (simulated quantities and exact solutions) in Cases A and B for 85 iteration steps and 46 concentration points. Figures 7 8 correspond to Case C for 46 concentration points and 55 iteration steps. Figures 9-1 correspond to Case A where we used 46 concentration points and a non-constant step ( t = 1/11 for the first 48 step, t = /11 for the next 4 steps and t = 4/11, for the last steps). Figures 11 1 correspond to Case B for 46 concentration points and 41 iteration steps. A few remarks and critical considerations are in order: The analysis of the results shows a good agreement between numerical and exact values. The results described in Figures 3 8 are a little bit finer than those obtained in [11]. Figures 9-1 have no correspondent in literature. The discretized method used in [7] is a first order Euler method. Our tests with non-constant -step, Figures 9-1, suggest that, a high-order method seems to improve the precision, without increasing the computational effort of the numerical scheme. This could be very useful for computations in space non-homogeneous cases. It appears from the analysis of Figures 11-1 that one can still obtain satisfactory results using less concentration points. However, we recall that, according to Theorem 1 in [7], the convergence of the numerical scheme, presented in this paper, requires the number of concentration points (of the measures) to be dependent on the length of -step of the iteration. Indeed, taking the same number of concentration points, a very small iteration -step may provide more inaccurate values of some estimated quantity (at a given moment) than a greater one. This numerical phenomenon was confirmed by tests in [1]. The main source of errors in the computation of the one-particle distribution function comes from the use of repeated random selections. Here it should be remarked that when the kinetic energy is calculated by means of the numeric values of concentrations ( kinetic energy 1 ), the probabilistic selection appears only in the evaluation of the integrals (3.3). Then the estimation of kinetic energy 1 is expected to be more accurate than for kinetic energy. Indeed, Figure 1 shows clearly that kinetic energy 1 is closer to the exact solution for energy than kinetic energy. Finally, we recall that the numerical method of [7] can be applied to more general models, with non constant collision kernels. However, in such situations, it is difficult to find exact solutions. Then the tests should rather be limited to comparing the data provided by the numerical scheme with those obtained from some relevant experiments.
13 17 D. Marinescu and C. P. Grünfeld References [1] H. Babovsky, On a simulation scheme for the Boltzmann equation, European Journal of Mechanics B/Fluids, 8, 41 (1989). [] H. Babovsky and R. Illner, A convergence proof for nambu s simulation method for the full Boltzmann equation, SIAM J. Num. Anal., 6, 45 (1989). [3] V. C. Boffi and G. Spiga, Extended kinetic theory for gas mixtures in the presence of removal and regeneration effects, Z. Angew. Math. Phys., 37, 7 (1986). [4] C. P. Grünfeld, On a class of kinetic equations for reacting gas mixtures with multiple collisions, C.R. Acad. Sci. Paris Serie I, 316, 953 (1993). [5] C. P. Grünfeld, Nonlinear kinetic models with chemical reactions, in N. Bellomo and M. Pulvirenti, editors, Modeling in Applied Sciences, a Kinetic Theory Approach (Birkhäuser, Boston, 1999). [6] C. P. Grünfeld and E. Georgescu, On a class of kinetic equations for reacting gas mixtures, Matematicheskaia Fizika, Analiz, Geometria,, 48 (1995). [7] C. P. Grünfeld and D. Marinescu, On the numerical simulation of a class of reactive Boltzmann type equations, Transp. Theory Stat. Phys., 6, 87 (1997). [8] M. Krook and T. T. Wu, Exact solutions of the Boltzmann equations for multicomponents systems, Phys. Rev. Lett., 38, 991 (1977). [9] L. Kuipers and H. Niederreiter, Uniform Distribution of Sequences (John Wiley & Sons, 1974). [1] D. Marinescu and A. Espesset, Tests of a simulation method for a system of Boltzmann equations, Computers & Mathematics with Applications, 4(6-7), 85 (). [11] D. Marinescu, A. Espesset and C.P. Grünfeld, Tests of a Simulation Method for Boltzmann-like Models with Chemical Reactions, Journal of Computational Physics, 175, 5,. [1] G. Spiga, Rigorous solution to the extended kinetic equations for homogeneous gas mixtures, in G. Toscani, V. C. Boffi, and S. Rionero, editors, Mathematical Aspects of Fluid and Plasma Dynamics, volume 146 of Lect. Notes in Phys. (Springer-Verlag, Berlin, 1991). Dorin Marinescu Romanian Academy, Institute of Applied Mathematics, Calea 13 Septembrie, No. 13, Sector 5, Bucharest P.O. Box 1-4, 77 dorinm@ima.ro Cecil P. Grünfeld Institute of Space Sciences, INFLPR, P.O. Box MG-3, Bucharest-Măgurele grunfeld@venus.nipne.ro
Recent advances in kinetic theory for mixtures of polyatomic gases
Recent advances in kinetic theory for mixtures of polyatomic gases Marzia Bisi Parma University, Italy Conference Problems on Kinetic Theory and PDE s Novi Sad (Serbia), September 25 27, 2014 M. Bisi,
More informationModelling and numerical methods for the diffusion of impurities in a gas
INERNAIONAL JOURNAL FOR NUMERICAL MEHODS IN FLUIDS Int. J. Numer. Meth. Fluids 6; : 6 [Version: /9/8 v.] Modelling and numerical methods for the diffusion of impurities in a gas E. Ferrari, L. Pareschi
More informationStochastic Particle Methods for Rarefied Gases
CCES Seminar WS 2/3 Stochastic Particle Methods for Rarefied Gases Julian Köllermeier RWTH Aachen University Supervisor: Prof. Dr. Manuel Torrilhon Center for Computational Engineering Science Mathematics
More informationGlobal Weak Solution to the Boltzmann-Enskog equation
Global Weak Solution to the Boltzmann-Enskog equation Seung-Yeal Ha 1 and Se Eun Noh 2 1) Department of Mathematical Science, Seoul National University, Seoul 151-742, KOREA 2) Department of Mathematical
More informationGeneralized Korn s inequality and conformal Killing vectors
Generalized Korn s inequality and conformal Killing vectors arxiv:gr-qc/0505022v1 4 May 2005 Sergio Dain Max-Planck-Institut für Gravitationsphysik Am Mühlenberg 1 14476 Golm Germany February 4, 2008 Abstract
More informationKinetic relaxation models for reacting gas mixtures
Kinetic relaxation models for reacting gas mixtures M. Groppi Department of Mathematics and Computer Science University of Parma - ITALY Main collaborators: Giampiero Spiga, Giuseppe Stracquadanio, Univ.
More informationKinetic models of Maxwell type. A brief history Part I
. A brief history Part I Department of Mathematics University of Pavia, Italy Porto Ercole, June 8-10 2008 Summer School METHODS AND MODELS OF KINETIC THEORY Outline 1 Introduction Wild result The central
More informationAn introduction to Birkhoff normal form
An introduction to Birkhoff normal form Dario Bambusi Dipartimento di Matematica, Universitá di Milano via Saldini 50, 0133 Milano (Italy) 19.11.14 1 Introduction The aim of this note is to present an
More informationAnalysis of the Reaction Rate of Binary Gaseous Mixtures near from the Chemical Equilibrium
Analysis of the Reaction Rate of Binary Gaseous Mixtures near from the Chemical Equilibrium Adriano W. da Silva May 30, 2013 Abstract A binary gaseous mixture with reversible reaction of type A + A = B
More informationarxiv:cond-mat/ v1 [cond-mat.stat-mech] 11 Dec 2002
arxiv:cond-mat/0212257v1 [cond-mat.stat-mech] 11 Dec 2002 International Journal of Modern Physics B c World Scientific Publishing Company VISCOUS FINGERING IN MISCIBLE, IMMISCIBLE AND REACTIVE FLUIDS PATRICK
More informationMonte Carlo methods for kinetic equations
Monte Carlo methods for kinetic equations Lecture 2: Monte Carlo simulation methods Lorenzo Pareschi Department of Mathematics & CMCS University of Ferrara Italy http://utenti.unife.it/lorenzo.pareschi/
More informationCurvature-homogeneous spaces of type (1,3)
Curvature-homogeneous spaces of type (1,3) Oldřich Kowalski (Charles University, Prague), joint work with Alena Vanžurová (Palacky University, Olomouc) Zlatibor, September 3-8, 2012 Curvature homogeneity
More informationDecomposition of Waveguides Propagating in Piezoelectric Crystals subject to Initial Fields
Proceedings of the Fifth Workshop on Mathematical Modelling of Environmental and Life Sciences Problems Constanţa, Romania, September, 006, pp. 191 00 Decomposition of Waveguides Propagating in Piezoelectric
More informationLow Variance Particle Simulations of the Boltzmann Transport Equation for the Variable Hard Sphere Collision Model
Low Variance Particle Simulations of the Boltzmann Transport Equation for the Variable Hard Sphere Collision Model G. A. Radtke, N. G. Hadjiconstantinou and W. Wagner Massachusetts Institute of Technology,
More informationDSMC Simulation of Binary Rarefied Gas Flows between Parallel Plates and Comparison to Other Methods
Simulation of Binary Rarefied Gas Flows between Parallel Plates and Comparison to Other Methods L. Szalmas Department of Mechanical Engineering, University of Thessaly, Pedion Areos, Volos 38334, Greece
More informationCH.9. CONSTITUTIVE EQUATIONS IN FLUIDS. Multimedia Course on Continuum Mechanics
CH.9. CONSTITUTIVE EQUATIONS IN FLUIDS Multimedia Course on Continuum Mechanics Overview Introduction Fluid Mechanics What is a Fluid? Pressure and Pascal s Law Constitutive Equations in Fluids Fluid Models
More informationIsometric elastic deformations
Isometric elastic deformations Fares Al-Azemi and Ovidiu Calin Abstract. This paper deals with the problem of finding a class of isometric deformations of simple and closed curves, which decrease the total
More informationDyadic diaphony of digital sequences
Dyadic diaphony of digital sequences Friedrich Pillichshammer Abstract The dyadic diaphony is a quantitative measure for the irregularity of distribution of a sequence in the unit cube In this paper we
More informationPhysical models for plasmas II
Physical models for plasmas II Dr. L. Conde Dr. José M. Donoso Departamento de Física Aplicada. E.T.S. Ingenieros Aeronáuticos Universidad Politécnica de Madrid Physical models,... Plasma Kinetic Theory
More informationLattice Bhatnagar Gross Krook model for the Lorenz attractor
Physica D 154 (2001) 43 50 Lattice Bhatnagar Gross Krook model for the Lorenz attractor Guangwu Yan a,b,,liyuan a a LSEC, Institute of Computational Mathematics, Academy of Mathematics and System Sciences,
More informationHydrodynamic Limits for the Boltzmann Equation
Hydrodynamic Limits for the Boltzmann Equation François Golse Université Paris 7 & Laboratoire J.-L. Lions golse@math.jussieu.fr Academia Sinica, Taipei, December 2004 LECTURE 2 FORMAL INCOMPRESSIBLE HYDRODYNAMIC
More informationELASTIC SCATTERING OF X-RAYS AND Γ-RAYS BY 2s ELECTRONS IN IONS AND NEUTRAL ATOMS *
Romanian Reports in Physics, Vol. 64, No. 4, P. 986 996, 0 ELASTIC SCATTERING OF X-RAYS AND Γ-RAYS BY s ELECTRONS IN IONS AND NEUTRAL ATOMS * K. KARIM, M. L. MUNTEANU, S. SPÂNULESCU,, C. STOICA University
More informationFluid Dynamics from Kinetic Equations
Fluid Dynamics from Kinetic Equations François Golse Université Paris 7 & IUF, Laboratoire J.-L. Lions golse@math.jussieu.fr & C. David Levermore University of Maryland, Dept. of Mathematics & IPST lvrmr@math.umd.edu
More informationINTRODUCTION TO PADÉ APPROXIMANTS. E. B. Saff Center for Constructive Approximation
INTRODUCTION TO PADÉ APPROXIMANTS E. B. Saff Center for Constructive Approximation H. Padé (1863-1953) Student of Hermite Histhesiswon French Academy of Sciences Prize C. Hermite (1822-1901) Used Padé
More informationCOMBINED EFFECTS FOR A STATIONARY PROBLEM WITH INDEFINITE NONLINEARITIES AND LACK OF COMPACTNESS
Dynamic Systems and Applications 22 (203) 37-384 COMBINED EFFECTS FOR A STATIONARY PROBLEM WITH INDEFINITE NONLINEARITIES AND LACK OF COMPACTNESS VICENŢIU D. RĂDULESCU Simion Stoilow Mathematics Institute
More informationOn Weak Solutions to the Linear Boltzmann Equation with Inelastic Coulomb Collisions
On Weak Solutions to the Linear Boltzmann Equation with Inelastic Coulomb Collisions Rolf Pettersson epartment of Mathematics, Chalmers University of Technology, SE-412 96 Göteborg, Sweden Abstract. This
More informationFrom Boltzmann Equations to Gas Dynamics: From DiPerna-Lions to Leray
From Boltzmann Equations to Gas Dynamics: From DiPerna-Lions to Leray C. David Levermore Department of Mathematics and Institute for Physical Science and Technology University of Maryland, College Park
More informationOverview of Accelerated Simulation Methods for Plasma Kinetics
Overview of Accelerated Simulation Methods for Plasma Kinetics R.E. Caflisch 1 In collaboration with: J.L. Cambier 2, B.I. Cohen 3, A.M. Dimits 3, L.F. Ricketson 1,4, M.S. Rosin 1,5, B. Yann 1 1 UCLA Math
More informationHypocoercivity for kinetic equations with linear relaxation terms
Hypocoercivity for kinetic equations with linear relaxation terms Jean Dolbeault dolbeaul@ceremade.dauphine.fr CEREMADE CNRS & Université Paris-Dauphine http://www.ceremade.dauphine.fr/ dolbeaul (A JOINT
More informationKinetic/Fluid micro-macro numerical scheme for Vlasov-Poisson-BGK equation using particles
Kinetic/Fluid micro-macro numerical scheme for Vlasov-Poisson-BGK equation using particles Anaïs Crestetto 1, Nicolas Crouseilles 2 and Mohammed Lemou 3. The 8th International Conference on Computational
More informationIncomplete exponential sums over finite fields and their applications to new inversive pseudorandom number generators
ACTA ARITHMETICA XCIII.4 (2000 Incomplete exponential sums over finite fields and their applications to new inversive pseudorandom number generators by Harald Niederreiter and Arne Winterhof (Wien 1. Introduction.
More informationComputational Solutions for the Korteweg devries Equation in Warm Plasma
COMPUTATIONAL METHODS IN SCIENCE AND TECHNOLOGY 16(1, 13-18 (1 Computational Solutions for the Korteweg devries Equation in Warm Plasma E.K. El-Shewy*, H.G. Abdelwahed, H.M. Abd-El-Hamid. Theoretical Physics
More informationLecture 5: Kinetic theory of fluids
Lecture 5: Kinetic theory of fluids September 21, 2015 1 Goal 2 From atoms to probabilities Fluid dynamics descrines fluids as continnum media (fields); however under conditions of strong inhomogeneities
More informationLong time behavior of non autonomous Fokker Planck equations and the cooling of granular gases
Long time behavior of non autonomous Fokker Planck equations and the cooling of granular gases Bertrand Lods Politecnico di Torino, Dipartimento di Matematica, Corso Duca degli Abruzzi, 1129 Torino, Italy.
More informationOn Some Mean Value Results for the Zeta-Function and a Divisor Problem
Filomat 3:8 (26), 235 2327 DOI.2298/FIL6835I Published by Faculty of Sciences and Mathematics, University of Niš, Serbia Available at: http://www.pmf.ni.ac.rs/filomat On Some Mean Value Results for the
More informationHilbert Sixth Problem
Academia Sinica, Taiwan Stanford University Newton Institute, September 28, 2010 : Mathematical Treatment of the Axioms of Physics. The investigations on the foundations of geometry suggest the problem:
More informationNon-linear Wave Propagation and Non-Equilibrium Thermodynamics - Part 3
Non-linear Wave Propagation and Non-Equilibrium Thermodynamics - Part 3 Tommaso Ruggeri Department of Mathematics and Research Center of Applied Mathematics University of Bologna January 21, 2017 ommaso
More informationc 2012 Society for Industrial and Applied Mathematics
SIAM J. NUMER. ANAL. Vol. 50, No. 4, pp. 849 860 c 0 Society for Industrial and Applied Mathematics TWO RESULTS CONCERNING THE STABILITY OF STAGGERED MULTISTEP METHODS MICHELLE GHRIST AND BENGT FORNBERG
More information1.3 Molecular Level Presentation
1.3.1 Introduction A molecule is the smallest chemical unit of a substance that is capable of stable, independent existence. Not all substances are composed of molecules. Some substances are composed of
More information2 P. L'Ecuyer and R. Simard otherwise perform well in the spectral test, fail this independence test in a decisive way. LCGs with multipliers that hav
Beware of Linear Congruential Generators with Multipliers of the form a = 2 q 2 r Pierre L'Ecuyer and Richard Simard Linear congruential random number generators with Mersenne prime modulus and multipliers
More informationFundamentals of Fluid Dynamics: Elementary Viscous Flow
Fundamentals of Fluid Dynamics: Elementary Viscous Flow Introductory Course on Multiphysics Modelling TOMASZ G. ZIELIŃSKI bluebox.ippt.pan.pl/ tzielins/ Institute of Fundamental Technological Research
More informationSupport Vector Machine Classification via Parameterless Robust Linear Programming
Support Vector Machine Classification via Parameterless Robust Linear Programming O. L. Mangasarian Abstract We show that the problem of minimizing the sum of arbitrary-norm real distances to misclassified
More informationUne approche hypocoercive L 2 pour l équation de Vlasov-Fokker-Planck
Une approche hypocoercive L 2 pour l équation de Vlasov-Fokker-Planck Jean Dolbeault dolbeaul@ceremade.dauphine.fr CEREMADE CNRS & Université Paris-Dauphine http://www.ceremade.dauphine.fr/ dolbeaul (EN
More informationChapter 3 Stress, Strain, Virtual Power and Conservation Principles
Chapter 3 Stress, Strain, irtual Power and Conservation Principles 1 Introduction Stress and strain are key concepts in the analytical characterization of the mechanical state of a solid body. While stress
More informationDRIVING FORCE IN SIMULATION OF PHASE TRANSITION FRONT PROPAGATION
Chapter 1 DRIVING FORCE IN SIMULATION OF PHASE TRANSITION FRONT PROPAGATION A. Berezovski Institute of Cybernetics at Tallinn Technical University, Centre for Nonlinear Studies, Akadeemia tee 21, 12618
More informationTHE CENTRAL INTERTWINING LIFTING AND STRICT CONTRACTIONS
PROCEEDINGS OF THE AMERICAN MATHEMATICAL SOCIETY Volume 124, Number 12, December 1996, Pages 3813 3817 S 0002-9939(96)03540-X THE CENTRAL INTERTWINING LIFTING AND STRICT CONTRACTIONS RADU GADIDOV (Communicated
More informationSimilarities and differences:
How does the system reach equilibrium? I./9 Chemical equilibrium I./ Equilibrium electrochemistry III./ Molecules in motion physical processes, non-reactive systems III./5-7 Reaction rate, mechanism, molecular
More informationExtremal Solutions of Differential Inclusions via Baire Category: a Dual Approach
Extremal Solutions of Differential Inclusions via Baire Category: a Dual Approach Alberto Bressan Department of Mathematics, Penn State University University Park, Pa 1682, USA e-mail: bressan@mathpsuedu
More informationChapter 2 A Brief Introduction to the Mathematical Kinetic Theory of Classical Particles
Chapter 2 A Brief Introduction to the Mathematical Kinetic Theory of Classical Particles 2.1 Plan of the Chapter The contents of this chapter are motivated by the second key question posed in Section 1.3,
More informationModelling of low-temperature plasmas: kinetic and transport mechanisms. L.L. Alves
Modelling of low-temperature plasmas: kinetic and transport mechanisms L.L. Alves llalves@tecnico.ulisboa.pt Instituto de Plasmas e Fusão Nuclear Instituto Superior Técnico, Universidade de Lisboa Lisboa,
More informationarxiv:quant-ph/ v1 22 Aug 2005
Conditions for separability in generalized Laplacian matrices and nonnegative matrices as density matrices arxiv:quant-ph/58163v1 22 Aug 25 Abstract Chai Wah Wu IBM Research Division, Thomas J. Watson
More informationRandom Averaging. Eli Ben-Naim Los Alamos National Laboratory. Paul Krapivsky (Boston University) John Machta (University of Massachusetts)
Random Averaging Eli Ben-Naim Los Alamos National Laboratory Paul Krapivsky (Boston University) John Machta (University of Massachusetts) Talk, papers available from: http://cnls.lanl.gov/~ebn Plan I.
More informationYAN GUO, JUHI JANG, AND NING JIANG
LOCAL HILBERT EXPANSION FOR THE BOLTZMANN EQUATION YAN GUO, JUHI JANG, AND NING JIANG Abstract. We revisit the classical ork of Caflisch [C] for compressible Euler limit of the Boltzmann equation. By using
More informationOn the Boltzmann equation: global solutions in one spatial dimension
On the Boltzmann equation: global solutions in one spatial dimension Department of Mathematics & Statistics Colloque de mathématiques de Montréal Centre de Recherches Mathématiques November 11, 2005 Collaborators
More informationf(x n ) [0,1[ s f(x) dx.
ACTA ARITHMETICA LXXX.2 (1997) Dyadic diaphony by Peter Hellekalek and Hannes Leeb (Salzburg) 1. Introduction. Diaphony (see Zinterhof [13] and Kuipers and Niederreiter [6, Exercise 5.27, p. 162]) is a
More informationFundamentals of Linear Elasticity
Fundamentals of Linear Elasticity Introductory Course on Multiphysics Modelling TOMASZ G. ZIELIŃSKI bluebox.ippt.pan.pl/ tzielins/ Institute of Fundamental Technological Research of the Polish Academy
More informationFourier Law and Non-Isothermal Boundary in the Boltzmann Theory
in the Boltzmann Theory Joint work with Raffaele Esposito, Yan Guo, Rossana Marra DPMMS, University of Cambridge ICERM November 8, 2011 Steady Boltzmann Equation Steady Boltzmann Equation v x F = Q(F,
More informationThe Boltzmann Equation and Its Applications
Carlo Cercignani The Boltzmann Equation and Its Applications With 42 Illustrations Springer-Verlag New York Berlin Heidelberg London Paris Tokyo CONTENTS PREFACE vii I. BASIC PRINCIPLES OF THE KINETIC
More informationVariational Theory of Solitons for a Higher Order Generalized Camassa-Holm Equation
International Journal of Mathematical Analysis Vol. 11, 2017, no. 21, 1007-1018 HIKAI Ltd, www.m-hikari.com https://doi.org/10.12988/ijma.2017.710141 Variational Theory of Solitons for a Higher Order Generalized
More informationRionero s critical perturbations method via weighted energy for stability of convective motions in anisotropic porous media
1 Rionero s critical perturbations method via weighted energy for stability of convective motions in anisotropic porous media F. Capone,a, M. Gentile,a, A. A. Hill,b a Department of Mathematics and Applications
More informationData on the Velocity Slip and Temperature Jump Coefficients
Data on the Velocity Slip and Temperature Jump Coefficients Felix Sharipov Departamento de Física, Universidade Federal do Paraná, Caixa Postal 19044, 81531-990 Curitiba, Brazil email: sharipov@fisica.ufpr.br;
More informationHamiltonian and Non-Hamiltonian Reductions of Charged Particle Dynamics: Diffusion and Self-Organization
NNP2017 11 th July 2017 Lawrence University Hamiltonian and Non-Hamiltonian Reductions of Charged Particle Dynamics: Diffusion and Self-Organization N. Sato and Z. Yoshida Graduate School of Frontier Sciences
More informationLecture II: The Lorentz gas
Lecture II: The Lorentz gas Applications of measure rigidity: from number theory to statistical mechanics Simons Lectures, Stony Brook October 203 Jens Marklof University of Bristol http://www.maths.bristol.ac.uk
More informationAnalytical Mechanics: Elastic Deformation
Analytical Mechanics: Elastic Deformation Shinichi Hirai Dept. Robotics, Ritsumeikan Univ. Shinichi Hirai (Dept. Robotics, Ritsumeikan Univ.) Analytical Mechanics: Elastic Deformation 1 / 60 Agenda Agenda
More informationOn the Midpoint Method for Solving Generalized Equations
Punjab University Journal of Mathematics (ISSN 1016-56) Vol. 40 (008) pp. 63-70 On the Midpoint Method for Solving Generalized Equations Ioannis K. Argyros Cameron University Department of Mathematics
More informationThe dynamics of small particles whose size is roughly 1 µmt or. smaller, in a fluid at room temperature, is extremely erratic, and is
1 I. BROWNIAN MOTION The dynamics of small particles whose size is roughly 1 µmt or smaller, in a fluid at room temperature, is extremely erratic, and is called Brownian motion. The velocity of such particles
More informationKinetic Theory of Gases
Kinetic Theory of Gases Modern Physics September 7 and 12, 2016 1 Intro In this section, we will relate macroscopic properties of gases (like Pressure, Temperature) to the behavior of the microscopic components
More informationASYMPTOTIC THEORY FOR WEAKLY NON-LINEAR WAVE EQUATIONS IN SEMI-INFINITE DOMAINS
Electronic Journal of Differential Equations, Vol. 004(004), No. 07, pp. 8. ISSN: 07-669. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu (login: ftp) ASYMPTOTIC
More informationTWO CENTER SHELL MODEL WITH WOODS-SAXON POTENTIALS
Romanian Reports in Physics, Vol. 59, No. 2, P. 523 531, 2007 Dedicated to Prof. Dorin N. Poenaru s 70th Anniversary TWO CENTER SHELL MODEL WITH WOODS-SAXON POTENTIALS M. MIREA Horia Hulubei National Institute
More informationSimplified Hyperbolic Moment Equations
Simplified Hyperbolic Moment Equations Julian Koellermeier and Manuel Torrilhon Abstract Hyperbolicity is a necessary property of model equations for the solution of the BGK equation to achieve stable
More informationA simple kinetic equation of swarm formation: blow up and global existence
A simple kinetic equation of swarm formation: blow up and global existence Miroslaw Lachowicz, Henryk Leszczyński, Martin Parisot To cite this version: Miroslaw Lachowicz, Henryk Leszczyński, Martin Parisot.
More informationDerivatives of Harmonic Bergman and Bloch Functions on the Ball
Journal of Mathematical Analysis and Applications 26, 1 123 (21) doi:1.16/jmaa.2.7438, available online at http://www.idealibrary.com on Derivatives of Harmonic ergman and loch Functions on the all oo
More informationThe existence of Burnett coefficients in the periodic Lorentz gas
The existence of Burnett coefficients in the periodic Lorentz gas N. I. Chernov and C. P. Dettmann September 14, 2006 Abstract The linear super-burnett coefficient gives corrections to the diffusion equation
More informationNEW LIFTS OF SASAKI TYPE OF THE RIEMANNIAN METRICS
NEW LIFTS OF SASAKI TYPE OF THE RIEMANNIAN METRICS Radu Miron Abstract One defines new elliptic and hyperbolic lifts to tangent bundle T M of a Riemann metric g given on the base manifold M. They are homogeneous
More informationEntropic structure of the Landau equation. Coulomb interaction
with Coulomb interaction Laurent Desvillettes IMJ-PRG, Université Paris Diderot May 15, 2017 Use of the entropy principle for specific equations Spatially Homogeneous Kinetic equations: 1 Fokker-Planck:
More informationLARGE TIME BEHAVIOR OF THE RELATIVISTIC VLASOV MAXWELL SYSTEM IN LOW SPACE DIMENSION
Differential Integral Equations Volume 3, Numbers 1- (1), 61 77 LARGE TIME BEHAVIOR OF THE RELATIVISTIC VLASOV MAXWELL SYSTEM IN LOW SPACE DIMENSION Robert Glassey Department of Mathematics, Indiana University
More informationEntropy-dissipation methods I: Fokker-Planck equations
1 Entropy-dissipation methods I: Fokker-Planck equations Ansgar Jüngel Vienna University of Technology, Austria www.jungel.at.vu Introduction Boltzmann equation Fokker-Planck equations Degenerate parabolic
More informationFlux Limiter Lattice Boltzmann for Compressible Flows
Commun. Theor. Phys. 56 (2011 333 338 Vol. 56, No. 2, August 15, 2011 Flux Limiter Lattice Boltzmann for Compressible Flows CHEN Feng (íô, 1 XU Ai-Guo (ÆÇÁ, 2, ZHANG Guang-Cai (¾, 2 and LI Ying-Jun (Ó
More informationREGULARIZATION AND BOUNDARY CONDITIONS FOR THE 13 MOMENT EQUATIONS
1 REGULARIZATION AND BOUNDARY CONDITIONS FOR THE 13 MOMENT EQUATIONS HENNING STRUCHTRUP ETH Zürich, Department of Materials, Polymer Physics, CH-8093 Zürich, Switzerland (on leave from University of Victoria,
More informationLattice Boltzmann Method
3 Lattice Boltzmann Method 3.1 Introduction The lattice Boltzmann method is a discrete computational method based upon the lattice gas automata - a simplified, fictitious molecular model. It consists of
More information12. MHD Approximation.
Phys780: Plasma Physics Lecture 12. MHD approximation. 1 12. MHD Approximation. ([3], p. 169-183) The kinetic equation for the distribution function f( v, r, t) provides the most complete and universal
More informationThe First Principle Calculation of Green Kubo Formula with the Two-Time Ensemble Technique
Commun. Theor. Phys. (Beijing, China 35 (2 pp. 42 46 c International Academic Publishers Vol. 35, No. 4, April 5, 2 The First Principle Calculation of Green Kubo Formula with the Two-Time Ensemble Technique
More informationMathematical modelling of collective behavior
Mathematical modelling of collective behavior Young-Pil Choi Fakultät für Mathematik Technische Universität München This talk is based on joint works with José A. Carrillo, Maxime Hauray, and Samir Salem
More informationON THE BOUNDEDNESS BEHAVIOR OF THE SPECTRAL FACTORIZATION IN THE WIENER ALGEBRA FOR FIR DATA
ON THE BOUNDEDNESS BEHAVIOR OF THE SPECTRAL FACTORIZATION IN THE WIENER ALGEBRA FOR FIR DATA Holger Boche and Volker Pohl Technische Universität Berlin, Heinrich Hertz Chair for Mobile Communications Werner-von-Siemens
More information14. Energy transport.
Phys780: Plasma Physics Lecture 14. Energy transport. 1 14. Energy transport. Chapman-Enskog theory. ([8], p.51-75) We derive macroscopic properties of plasma by calculating moments of the kinetic equation
More informationA fast spectral method for the inelastic Boltzmann collision operator and application to heated granular gases
A fast spectral method for the inelastic Boltzmann collision operator and application to heated granular gases Jingwei Hu a, Zheng Ma a a Department of Mathematics, Purdue University 150 N. University
More informationIBM Research Report. Conditions for Separability in Generalized Laplacian Matrices and Diagonally Dominant Matrices as Density Matrices
RC23758 (W58-118) October 18, 25 Physics IBM Research Report Conditions for Separability in Generalized Laplacian Matrices and Diagonally Dominant Matrices as Density Matrices Chai Wah Wu IBM Research
More informationChemnitz Scientific Computing Preprints
Arnd Meyer The Koiter shell equation in a coordinate free description CSC/1-0 Chemnitz Scientific Computing Preprints ISSN 1864-0087 Chemnitz Scientific Computing Preprints Impressum: Chemnitz Scientific
More informationLaplace Transform of Spherical Bessel Functions
Laplace Transform of Spherical Bessel Functions arxiv:math-ph/01000v 18 Jan 00 A. Ludu Department of Chemistry and Physics, Northwestern State University, Natchitoches, LA 71497 R. F. O Connell Department
More informationNon-equilibrium mixtures of gases: modeling and computation
Non-equilibrium mixtures of gases: modeling and computation Lecture 1: and kinetic theory of gases Srboljub Simić University of Novi Sad, Serbia Aim and outline of the course Aim of the course To present
More informationFluid dynamics for a vapor-gas mixture derived from kinetic theory
IPAM Workshop II The Boltzmann Equation: DiPerna-Lions Plus 20 Years (IPAM-UCLA, April 15-17, 2009) Fluid dynamics for a vapor-gas mixture derived from kinetic theory Kazuo Aoki Department of Mechanical
More informationVibrational degrees of freedom in the Total Collision Energy DSMC chemistry model
Vibrational degrees of freedom in the Total Collision Energy DSMC chemistry model Mark Goldsworthy, Michael Macrossan Centre for Hypersonics, School of Engineering, University of Queensland, Brisbane,
More informationAverage Orders of Certain Arithmetical Functions
Average Orders of Certain Arithmetical Functions Kaneenika Sinha July 26, 2006 Department of Mathematics and Statistics, Queen s University, Kingston, Ontario, Canada K7L 3N6, email: skaneen@mast.queensu.ca
More informationClassical inequalities for the Boltzmann collision operator.
Classical inequalities for the Boltzmann collision operator. Ricardo Alonso The University of Texas at Austin IPAM, April 2009 In collaboration with E. Carneiro and I. Gamba Outline Boltzmann equation:
More informationONSAGER S RECIPROCAL RELATIONS AND SOME BASIC LAWS
Journal of Computational and Applied Mechanics, Vol. 5., No. 1., (2004), pp. 157 163 ONSAGER S RECIPROCAL RELATIONS AND SOME BASIC LAWS József Verhás Department of Chemical Physics, Budapest University
More informationKinetic Transport Theory
Lecture notes on Kinetic Transport Theory Christian Schmeiser 1 Preface Kinetic transport equations are mathematical descriptions of the dynamics of large particle ensembles in terms of a phase space (i.e.
More informationRELATIONS BETWEEN THE ARRHENIUS ACTIVATION ENERGY AND THRESHOLD ENERGY FOR SIMPLE MODELS OF THE REACTIVE CROSS SECTIONS IN A DILUTE GAS
Vol. 37 (2006) CT PHYSIC POLONIC B No 6 RELTIONS BETWEEN THE RRHENIUS CTIVTION ENERGY ND THRESHOLD ENERGY FOR SIMPLE MODELS OF THE RECTIVE CROSS SECTIONS IN DILUTE GS.S. Cukrowski Institute of Physical
More informationTHE LAPLACE TRANSFORM OF THE FOURTH MOMENT OF THE ZETA-FUNCTION. Aleksandar Ivić
THE LAPLACE TRANSFORM OF THE FOURTH MOMENT OF THE ZETA-FUNCTION Aleksandar Ivić Univ. Beograd. Publ. Elektrotehn. Fak. Ser. Mat. (2), 4 48. Abstract. The Laplace transform of ζ( 2 +ix) 4 is investigated,
More informationInternal Stabilizability of Some Diffusive Models
Journal of Mathematical Analysis and Applications 265, 91 12 (22) doi:1.16/jmaa.21.7694, available online at http://www.idealibrary.com on Internal Stabilizability of Some Diffusive Models Bedr Eddine
More informationROTATIONAL STABILITY OF A CHARGED DIELEC- TRIC RIGID BODY IN A UNIFORM MAGNETIC FIELD
Progress In Electromagnetics Research Letters, Vol. 11, 103 11, 009 ROTATIONAL STABILITY OF A CHARGED DIELEC- TRIC RIGID BODY IN A UNIFORM MAGNETIC FIELD G.-Q. Zhou Department of Physics Wuhan University
More information