c 2012 International Press FAST COMMUNICATION
|
|
- Mark Page
- 5 years ago
- Views:
Transcription
1 COMMU. MATH. SCI. Vol. 1, o. 3, pp c 1 International Press FAST COMMUICATIO A FAST SPECTRAL ALGORITHM FOR THE QUATUM BOLTZMA COLLISIO OPERATOR JIGWEI HU AD LEXIG YIG Abstract. This paper introduces a fast spectral algorithm for the quantum Boltzmann collision operator. In the usual spectral framework, one of the terms in the operator cannot be evaluated efficiently. The new approach is based on the fundamental property of the exponential function which allows one to construct a new decomposition of the collision kernel to speed up the computation. umerical results in -D and 3-D for both the Bose gas and the Fermi gas are presented to illustrate the accuracy and efficiency of the method. Key words. Quantum Boltzmann equation, fast spectral method. AMS subject classifications. 35Q, 65M7. 1. Introduction The quantum Boltzmann equation describes the time evolution of a dilute Bose gas or Fermi gas. It was first formulated by ordheim [13] and Uehling and Uhlenbeck [17] from the classical Boltzmann equation by heuristic arguments. If f(t,x,v) is the phase space distribution function of time t, position x, and particle velocity v, then the equation reads f t +v xf=q q (f), x Ω R dx, v R dv. (1.1) The quantum collision operator Q q is given by Q q (f)(v)= B(v v,ω)[f f (1±θ f)(1±θ f ) R dv S ff (1±θ f )(1±θ f )]dωdv, (1.) where θ =h dv and h is the rescaled Planck constant. In this paper, the upper sign will always correspond to the Bose gas while the lower sign to the Fermi gas. f, f, f, and f are the shorthand notations for f(t,x,v), f(t,x,v ), f(t,x,v ), and f(t,x,v ) respectively. (v,v ) and (v,v ) are the velocities before and after collision. They are related by the following parametrization: v = v+v v = v+v + v v ω, v v ω, (1.3) Received: September 7, 11; accepted (in revised version): October 14, 11. Communicated by Shi Jin. J.H. would like to thank KAUST for their generous support. L.Y. was supported in part by the SF under CAREER grant DMS Institute for Computational Engineering and Sciences (ICES), The University of Texas at Austin, 1 University Station C, Austin, TX 7871, USA (hu@ices.utexas.edu). Department of Mathematics and ICES, The University of Texas at Austin, 1 University Station C1, Austin, TX 7871, USA (lexing@math.utexas.edu). 989
2 99 FAST ALGORITHM FOR QUATUM BOLTZMA COLLISIO OPERATOR where ω is the unit vector along v v. The collision kernel B is a nonnegative function that only depends on v v and cosθ=w (v v )/ v v. From now on we will write it as B( v v,cosθ). The quantum collision operator Q q satisfies the conservation of mass, momentum, and energy, Q q (f)dv= Q q (f)vdv= Q q (f) v dv=, (1.4) R dv R dv R dv and Boltzmann s H-theorem R dv ln ( ) f Q q (f)dv. (1.5) 1±θ f Moreover, the equality in (1.5) holds if and only if f attains the local equilibrium (the quantum Maxwellian): M q = 1 1, (1.6) θ z 1 e (v u) T 1 where z is the fugacity, T is the temperature, and u is the macroscopic velocity. (1.6) is the well-known Bose-Einstein ( - ) and Fermi-Dirac ( + ) distributions. For further details about the derivation of M q and the properties of the quantum Boltzmann equation we refer to [9, 5]. The numerical challenge of solving the Boltzmann equation mainly comes from the multidimensional structure of the collision integral. For the classical Boltzmann operator, many different approaches are available: for example, the direct simulation Monte Carlo (DSMC) methods [1, 1], the discrete velocity models (DVM) [16,, 4], and more recently the spectral methods [3, 14, 15, 11, 6, 7]. On the contrary, very few studies have been conducted in the quantum case. See [8] and references therein for the Monte Carlo simulations and [1] for a fast deterministic scheme for the energy space boson Boltzmann equation. In [5] the authors constructed a spectral method for the quantum Boltzmann collision operator (1.) following the framework of [11, 6]. However, due to its cubic structure, one of the terms in Q q cannot be evaluated by the fast Fourier transform (FFT), which results in a computational cost of O(M dv log) (M is the number of points in angular direction; is the number of points in each v direction). In this work, we improve the above method and develop a fast algorithm for the quantum Boltzmann operator based on a simple observation of the exponential function. The numerical cost is reduced to O(M dv+1 log). The rest of the paper is organized as follows: In the next section we first describe the spectral method in [5] and then present the new fast algorithm. umerical results in both -D and 3-D are shown in Section 3 to illustrate the accuracy and efficiency of the method. Finally the concluding remarks are given in Section 4.. The spectral methods for the quantum Boltzmann collision operator Webeginbyabriefdescriptionofthespectralmethodin[5]. Itservesasthestarting point of our fast algorithm and as a comparison in the numerical tests. For simplicity, we drop the spatial dependence and consider the homogeneous quantum Boltzmann equation: f t =Q q(f), v R dv. (.1)
3 where J.-W. HU AD L.-X. YIG First attempt. The quantum collision operator Q q can be recast as Q q =Q c ±θ (Q 1 +Q Q 3 Q 4 ), (.) Q c (f)(v)= B( v v,cosθ)[f f ff ]dωdv (.3) R dv S is the classical collision operator. The cubic terms Q 1 Q 4 are Q 1 (f)(v)= B( v v,cosθ)f f f dωdv, R dv S Q (f)(v)= B( v v,cosθ)f f fdωdv, R dv S (.4) Q 3 (f)(v)= B( v v,cosθ)ff f dωdv, R dv S Q 4 (f)(v)= B( v v,cosθ)ff f dωdv. R dv S In order to perform the Fourier transform, we periodize the function f on the domain D L =[ L,L] dv (L is chosen such that L 3+ S, where B S is an approximation of the support of f(v) [15]). Using a Carleman like representation [11], one can transform the above operators to Q c (f)(v)= B(x,y)δ(x y)[f(v+x)f(v+y) f(v+x+y)f(v)]dxdy (.5) B R B R and Q 1 (f)(v)= B(x,y)δ(x y)f(v+x)f(v+y)f(v+x+y)dxdy, B R B R Q (f)(v)= B(x,y)δ(x y)f(v+x)f(v+y)f(v)dxdy, B R B R (.6) Q 3 (f)(v)= B(x,y)δ(x y)f(v+x)f(v+x+y)f(v)dxdy, B R B R Q 4 (f)(v)= B(x,y)δ(x y)f(v+y)f(v+x+y)f(v)dxdy, B R B R where B R is the truncation of the collision integral which satisfies R S; B(x,y) ( ) x = B + y,1 x ( x + y ) dv x + y = B( x, y ). (.7) We approximate f by a truncated Fourier series (k is the d v -dimensional index): f(v) k= ˆf k e i π L k v, ˆfk = 1 f(v)e i π L k v dv. (.8) (L) dv D L
4 99 FAST ALGORITHM FOR QUATUM BOLTZMA COLLISIO OPERATOR Plugging it into (.1), one can get the following equation for ˆf k by the orthogonality of the basis e i π L k v, k Z dv : dˆf k dt = ˆQ q,k = ˆQ c,k ±θ ( ˆQ 1,k + ˆQ,k ˆQ 3,k ˆQ 4,k ), (.9) where ˆQ c,k = l,m= l+m=k [β(l,m) β(m,m)] ˆf l ˆfm, (.1) and ˆQ 1,k = β(l+n,m+n)ˆf l ˆfm ˆfn, ˆQ,k = β(l,m)ˆf l ˆfm ˆfn, ˆQ 3,k = l,m,n= l,m,n= β(l+m,m)ˆf l ˆfm ˆfn, ˆQ4,k = l,m,n= l,m,n= β(m,l+m)ˆf l ˆfm ˆfn, (.11) with the kernel mode β defined by π i β(l,m) = B( x, y )δ(x y)e L l x e i π L m y dxdy. (.1) B R B R The classical part was treated in [11] by introducing a decomposition for β: β(l,m) α p (l)α p(m), (.13) where α p, α p are some functions to be given below. ow with the same decomposition the remaining four cubic quantum terms can be approximated as ˆQ 1,k n= l,m= l+m=k n α p (l+n)α p(m+n)ˆf l ˆfm ˆf n = n= ĝ k n (n)ˆf n, (.14)
5 ˆQ,k ˆQ 3,k ˆQ 4,k J.-W. HU AD L.-X. YIG 993 n= n= n= l,m= l+m=k n α p (k n) α p (l)α p(m)ˆf l ˆfm ˆf n, (.15) α p(k n) l,m= l+m=k n l,m= l+m=k n α p(m)ˆf l ˆfm ˆf n, (.16) α p (m)ˆf l ˆfm ˆf n. (.17) Thus, with the above arrangements, Q, Q 3, and Q 4 (all in double-convolution forms) can be evaluated by the FFT in O(M dv log) operations. For instance, to compute Q 3, we first treat ˆf l and α p(m)ˆf m as two functions, then terms inside the bracket of (.16) form a convolution and can be evaluated as multiplication in the original space via the FFT. The resulting function (with subindex k n) multiplied by α p (k n) together with ˆf n forms a convolution again and can be computed by the FFT as well. However, the fast algorithm cannot be applied to Q 1, since ĝ k n (n) itself depends on n (although the terms inside the bracket are a convolution). We compute Q 1 directly at the moment and the total computational cost is dominated by this part, which is O(M dv log) The decomposition of β. We now discuss the decomposition (.13) of β. Following [11], (.1) becomes β(l,m)= 1 δ(e e )dede 4 e,e S R R B( ρ, ρ )e i π L ρ(l e) e i π L ρ (m e ) ρ dv ρ dv dρdρ (.18) in polar coordinates. We make the decoupling assumption that B( x, y ) = a( x )b( y ), (.19) then β(l,m)= 1 dede 4 e,e S,e e ( )( R ) R a( ρ ) ρ dv e i π L ρ(l e) dρ b( ρ ) ρ dv e i π L ρ (m e ) dρ = 1 φ(l e)φ (m e )dede, (.) 4 e,e S,e e where φ(s) =. R φ (s). = R a( ρ ) ρ dv e i π L ρs dρ, (.1) b( ρ ) ρ dv e i π L ρs dρ. (.)
6 994 FAST ALGORITHM FOR QUATUM BOLTZMA COLLISIO OPERATOR In the -D case, β(l,m) = π so we can approximate it by φ(l e θ )φ (m e θ+ π )dθ, e θ=(cosθ,sinθ), (.3) β(l,m) π M φ(l e θp )φ (m e θp+ π ), (.4) where M is the number of equally spaced points in [,π]. In the 3-D case, β(l,m) = φ(l e)ψ( Π e (m) )de, ψ(s) =. π φ (scosη)sinηdη, (.5) e S + where Π e is the orthogonal projection onto the plane e. Parametrizing e as e =(sinϕcosθ,sinϕsinθ,cosϕ), θ,ϕ [,π], we can approximate β by ( π β(l,m) M) M M is the same as in -D. p,q=1 φ(l e θp,ϕ q )ψ( Π e θp,ϕq (m) )sinϕ q, (.6).. A fast algorithm. Based on previous discussions we see that the decomposition (.13) is good for the fast evaluation of Q c (see [11]), Q, Q 3, and Q 4. To improve the computation of Q 1, we need to seek other decompositions. Suppose we can write β as follows: β(l,m) J exp(iα j (l))α j(m), (.7) j=1 where α j is assumed to be linear in l. Then one can compute Q 1 efficiently by using the fundamental property of the exponential function: e a+b =e a e b. In fact, (.14) becomes ˆQ 1,k = β(l+n,m+n)ˆf l ˆfm ˆfn = l,m,n= J j=1l,m,n= J j=1l= exp(iα j (l))exp(iα j (n))α j(m+n)ˆf l ˆfn ˆfm exp(iα j (l))α j(k l) m,n= m+n=k l exp(iα j (n))ˆf n ˆfm ˆf l, (.8) so it is again in the double-convolution form and can be evaluated by the FFT.
7 J.-W. HU AD L.-X. YIG 995 Fortunately, the decomposition (.7) can be readily obtained thanks to the special form of φ(s) (.1). Specifically, we still adopt the discretizations (.4) and (.6), but instead of precomputing φ(s), we approximate it by a numerical quadrature rule: J [ ( φ(s) a( ρ j ) ρ j dv exp i π )] L ρ js w j, (.9) j=1 where ρ j,w j, j=1,...j are the integration points and weights in [,R]. Therefore, in the -D case, β(l,m) π M in the 3-D case, J j=1 ( π M β(l,m) M) [ ( a( ρ j )exp i π L ρ ) j( ) l eθp φ (m e θp+ ] w π) j ; (.3) J p,q=1j=1 [ ( a( ρ j ) ρ j exp i π L ρ ( ) ) j l eθp,ϕ q ψ( Π e θp,ϕq (m) )sinϕ q ] w j. (.31) In our implementation, we use the Gauss quadrature or composite Gauss quadrature for the integration depending on the properties of a( ρ ) ρ dv. For reasonably nice a( ρ ) ρ dv, the integrand of φ(l e) oscillates on the scale of O() for l [, 1], so J should be taken as O(). Thus the total number of terms in the decomposition of β is equal to M J=O(M ). We are ready to summarize the final algorithm for computing Q q (f). Given f(v) at a fixed time t, 1. Compute ˆf k by the FFT.. Evaluate Q c, Q, Q 3, and Q 4 by the method presented in Section.1, with β(l,m) given by the formulas (.4) and (.6). 3. Evaluate Q 1 by the method presented in this subsection, with β(l,m) given by the formulas (.3) and (.31). 4. Construct Q q through (.). The main work load still comes from step 3, which is now reduced to O(M dv+1 log). 3. umerical tests In this section we provide some numerical results of our new algorithm. We first test its performance on the steady state, i.e. compute Q q (M q ) and check the accuracy. Fixing the macroscopic quantities: density ρ = 1, temperature T = 1, macrovelocity u=, one can adjust θ to get M q that lies in different physical regimes [9]. We consider the -D Maxwellian molecules and 3-D hard sphere models which both satisfythedecouplingassumption(.19), as B isaconstantineithercase(withoutloss of generality, assume a=b=1). The computational domain for v is taken as [ 8,8] dv. The Gauss-Legendre quadrature rule with J = O() is employed to integrate φ(s). In the 3-D case, due to the kink ρ in the integrand we divide [,R] to [,] and [,R], and apply the Gauss quadrature on each interval. The implementation is in MATLAB and all the numerical results are obtained on a desktop computer with.8ghz CPU. Further acceleration can be achieved by careful implementation in C or Fortran.
8 996 FAST ALGORITHM FOR QUATUM BOLTZMA COLLISIO OPERATOR D Maxwellian molecules. In this case ( π ) φ(s)=φ (s)=rsinc L Rs, (3.1) where Sinc(x) = sinx/x. Let θ =.1 (Planck constant h=.1), then z Bose =1.593e 3, z Fermi = 1.598e 3. Here z 1 implies that the quantum effect is very small. Either the Bose gas or Fermi gas should behave like a classical gas. When we increase θ, say θ =3 (h=3), then z Bose =.7613, z Fermi = ow the difference between the quantum gases and classical gas is evident (Figure 3.1) Planck constant h=3 Bose Einstein distribution Fermi Dirac distribution classical Maxwellian Fig The cross section of the Maxwellians when θ =3. Tables3.1and3.reportthemaxnormsofQ q (M q )generatedondifferentmeshes. For comparison we also list the values computed by the spectral method without the speedup strategy for Q 1 (we call it direct method). Clearly our new algorithm gives accuracy comparable to the direct method. Table 3.3 shows the average running time of both methods for computing Q 1. The speedup factor is more than 1 in =18 case. fast, Bose (J=) direct, Bose fast, Fermi (J=) direct, Fermi e e e e e e e e e e e e e e e e-1 Table 3.1. Q q(m q) L in the classical regime (θ =.1 ) computed by the fast algorithm and the direct method. M =8 in all the simulations D hard sphere models. In this case φ(s)=r [ ( π ) Sinc L Rs Sinc ( π )] ( π ) L Rs, ψ(s)=r Sinc L Rs. (3.) We test the fast algorithm in (1) the classical regime: θ =.3 3 (Planck constant h=.3), which corresponds to z Bose =1.7133e 3, z Fermi =1.7154e 3; and () the
9 J.-W. HU AD L.-X. YIG 997 fast, Bose (J=) direct, Bose fast, Fermi (J=) direct, Fermi e e-6.874e e e e e e e-9.414e e-11 Table 3.. Q q(m q) L in the quantum regime (θ =3 ) computed by the fast algorithm and the direct method. M =8 in all the simulations. fast algorithm (J = ) direct method 16.4s.36s 3.9s.3s 64.55s 8.91s 18.43s s Table 3.3. The average running time of the fast algorithm and the direct method for computing Q 1. M =8 in all the simulations. fast, Bose (J=) direct, Bose fast, Fermi (J=) direct, Fermi e e e e e e e e e e e e-1 Table 3.4. Q q(m q) L in the classical regime (θ =.3 3 ) computed by the fast algorithm and the direct method. M =8 in all the simulations. fast, Bose (J=) direct, Bose fast, Fermi (J=) direct, Fermi e e e e e e e e-11 Table 3.5. Q q(m q) L in the quantum regime (θ =.6 3 for Bose gas; θ =3 3 for Fermi gas) computed by the fast algorithm and the direct method. M =8 in all the simulations. quantum regime: θ =.6 3 (h=.6) for the Bose gas, z Bose =.7474; θ =3 3 (h=3) for the Fermi gas, z Fermi = The numerical results of both methods are given in Tables In the 3-D case, the fast algorithm provides a huge saving of the computational time, while the direct method is too slow Relaxation to equilibrium. In this subsection, we would like to test the numerical performance of the algorithm on a general distribution f, instead of M q. Consider the homogeneous Equation (.1) for a -D Bose gas of Maxwellian
10 998 FAST ALGORITHM FOR QUATUM BOLTZMA COLLISIO OPERATOR fast algorithm (J = ) direct method 8.8s 5.4s 16.99s 36.3s s 1s s 18 19s Table 3.6. The average running time of the fast algorithm and the direct method for computing Q 1. M =8 in all the simulations f M q L 1 relative entropy H time time Fig. 3.. The time evolution of f M q L 1 and the relative entropy H for a -D Bose gas. 1.1 density time.1 macro velocity 1.1 internal energy time time Fig The time evolution of density, macro-velocity, and specific internal energy for a -D Bose gas.
11 molecules, with the nonequilibrium initial data f (v)= ρ ( exp ( v u ) +exp 4πT T J.-W. HU AD L.-X. YIG 999 ( v+u T )), (3.3) where ρ =1, T =3/8, and u =(1,1/). The fast algorithm for Q q and the fourthorder Runge-Kutta method for the time derivative are applied to solve the equation. Here the Planck constant h is chosen as 3 such that z= The computational domain is [ 8,8]. Here =64=J and M=8. Figure 3. shows the time evolution of f M q L 1 and the relative entropy H, from which we clearly observe the relaxation to equilibrium and the entropy decay. Furthermore, the mass, momentum, and energy are also preserved very well (Figure 3.3). 4. Conclusion We improved the spectral method in [5] and developed a fast algorithm for the quantum Boltzmann collision operator. The special form of the collision kernel in the spatial domain enables us to accelerate the computation. With the fast Fourier transform, the algorithm complexity is brought down to O(M dv+1 log) as opposed to O(M dv log) of the previous method. umerical results showed that the new algorithm preserves the spectral accuracy as well as offers a significant speedup over the original approach. REFERECES [1] G.A. Bird, Molecular Gas Dynamics and the Direct Simulation of Gas Flows, Clarendon Press, Oxford, [] A.V. Bobylev, A. Palczewski, and J. Schneider, On approximation of the Boltzmann equation by discrete velocity models, C.R. Acad. Sci. Paris Sér. I Math., 3, , [3] A.V. Bobylev and S. Rjasanow, Fast deterministic method of solving the Boltzmann equation for hard spheres, European J. Mech. B Fluids, 18, , [4] C. Buet, A discrete velocity scheme for the Boltzmann operator of rarefied gas dynamics, Transp. Theory Stat. Phys., 5, 33 6, [5] F. Filbet, J.W. Hu, and S. Jin, A numerical scheme for the quantum Boltzmann equation with stiff collision terms, ESAIM-Math. Model. umer. Anal., 46, , 1. [6] F. Filbet, C. Mouhot, and L. Pareschi, Solving the Boltzmann equation in log, SIAM J. Sci. Comput., 8, , 6. [7] I.M. Gamba and S.H. Tharkabhushanam, Spectral-Lagrangian methods for collisional models of non-equilibrium statistical states, J. Comput. Phys., 8, 1 36, 9. [8] A.L. Garcia and W. Wagner, Direct simulation Monte Carlo method for the Uehling-Uhlenbeck- Boltzmann equation, Phys. Rev. E, 68, 5673, 3. [9] J.W. Hu and S. Jin, On kinetic flux vector splitting schemes for quantum Euler equations, Kinet. Relat. Models, 4, , 11. [1] P. Markowich and L. Pareschi, Fast, conservative and entropic numerical methods for the Bosonic Boltzmann equation, umerische Math., 99, 59 53, 5. [11] C. Mouhot and L. Pareschi, Fast algorithms for computing the Boltzmann collision operator, Math. Comput., 75, , 6. [1] K. anbu, Direct simulation scheme derived from the Boltzmann equation. I. Monocomponent gases, J. Phys. Soc. Japan, 5, 4 49, [13] L.W. ordheim, On the kinetic method in the new statistics and its application in the electron theory of conductivity, Proc. R. Soc. London, Ser. A, 119, , 198. [14] L. Pareschi and B. Perthame, A Fourier spectral method for homogeneous Boltzmann equations, Transp. Theory Stat. Phys., 5, , [15] L. Pareschi and G. Russo, umerical solution of the Boltzmann equation I. Spectrally accurate approximation of the collision operator, SIAM J. umer. Anal., 37, ,. [16] F. Rogier and J. Schneider, A direct method for solving the Boltzmann equation, Transp. Theory Stat. Phys., 3, , [17] E.A. Uehling and G.E. Uhlenbeck, Transport phenomena in Einstein-Bose and Fermi-Dirac gases I, Phys. Rev., 43, , 1933.
A 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 informationhal , version 1-18 Sep 2013
Author manuscript, published in "Multiscale Modeling & Simulation 0 (202 792-87" ON DETERMINISTIC APPROXIMATION OF THE BOLTZMANN EQUATION IN A BOUNDED DOMAIN FRANCIS FILBET Abstract. In this paper we present
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 informationA conservative spectral method for the Boltzmann equation with anisotropic scattering and the grazing collisions limit
Work supported by NSF grants DMS-636586 and DMS-1217254 Jeff Haack (UT Austin) Conservative Spectral Boltzmann 7 Jun 213 1 / 21 A conservative spectral method for the Boltzmann equation with anisotropic
More informationA Conservative Discontinuous Galerkin Scheme with O(N 2 ) Operations in Computing Boltzmann Collision Weight Matrix
A Conservative Discontinuous Galerkin Scheme with O(N 2 ) Operations in Computing Boltzmann Collision Weight Matrix Irene M. Gamba, and Chenglong Zhang ICES, The University of Texas at Austin, 201 E. 24th
More informationMonte Carlo methods for kinetic equations
Monte Carlo methods for kinetic equations Lecture 4: Hybrid methods and variance reduction Lorenzo Pareschi Department of Mathematics & CMCS University of Ferrara Italy http://utenti.unife.it/lorenzo.pareschi/
More informationHybrid and Moment Guided Monte Carlo Methods for Kinetic Equations
Hybrid and Moment Guided Monte Carlo Methods for Kinetic Equations Giacomo Dimarco Institut des Mathématiques de Toulouse Université de Toulouse France http://perso.math.univ-toulouse.fr/dimarco giacomo.dimarco@math.univ-toulouse.fr
More informationFast methods for the Boltzmann collision integral
Fast methods for the Boltzmann collision integral Clément Mouhot Lorenzo Pareschi Abstract In this note we present methods for the development of fast numerical schemes for the Boltzmann collision integral.
More informationExponential methods for kinetic equations
Exponential methods for kinetic equations Lorenzo Pareschi Department of Mathematics & CMCS University of Ferrara, Italy http://utenti.unife.it/lorenzo.pareschi/ lorenzo.pareschi@unife.it Joint research
More informationAccepted Manuscript. Deterministic numerical solutions of the Boltzmann equation using the fast spectral method
Accepted Manuscript Deterministic numerical solutions of the Boltzmann equation using the fast spectral method Lei Wu, Craig White, Thomas J. Scanlon, Jason M. Reese, Yonghao Zhang PII: S1-9991(13)37-6
More informationReview of last lecture I
Review of last lecture I t f + v x f = Q γ (f,f), (1) where f (x,v,t) is a nonegative function that represents the density of particles in position x at time t with velocity v. In (1) Q γ is the so-called
More informationA Numerical Scheme for the Quantum Fokker-Planck-Landau Equation Efficient in the Fluid Regime
Commun. Comput. Phys. doi:.48/cicp.4.9a Vol., No. 5, pp. 54-56 November A Numerical Scheme for the Quantum Fokker-Planck-Landau Equation Efficient in the Fluid Regime Jingwei Hu, Shi Jin,3, and Bokai Yan
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 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 informationHypocoercivity and Sensitivity Analysis in Kinetic Equations and Uncertainty Quantification October 2 nd 5 th
Hypocoercivity and Sensitivity Analysis in Kinetic Equations and Uncertainty Quantification October 2 nd 5 th Department of Mathematics, University of Wisconsin Madison Venue: van Vleck Hall 911 Monday,
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 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 informationNumerical methods for kinetic equations
Numerical methods for kinetic equations Lecture 6: fluid-kinetic coupling and hybrid methods Lorenzo Pareschi Department of Mathematics and Computer Science University of Ferrara, Italy http://www.lorenzopareschi.com
More informationNumerical methods for plasma physics in collisional regimes
J. Plasma Physics (15), vol. 81, 358116 c Cambridge University Press 1 doi:1.117/s3778176 1 Numerical methods for plasma physics in collisional regimes G. Dimarco 1,Q.Li,L.Pareschi 1 and B. Yan 3 1 Department
More informationExponential moments for the homogeneous Kac equation
Exponential moments for the homogeneous Kac equation Maja Tasković University of Pennsylvania Young Researchers Workshop: Stochastic and deterministic methods in kinetic theory, Duke University, November
More informationDevelopment of a 3D Spectral Boltzmann Solver for Plasma Modeling
Development of a 3D Spectral Boltzmann Solver for Plasma Modeling IEPC-9-38 Presented at the 3 st International Electric Propulsion Conference, University of Michigan, Ann Arbor, Michigan, USA Jacques
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 informationExponential tail behavior for solutions to the homogeneous Boltzmann equation
Exponential tail behavior for solutions to the homogeneous Boltzmann equation Maja Tasković The University of Texas at Austin Young Researchers Workshop: Kinetic theory with applications in physical sciences
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 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 informationarxiv: v1 [math.na] 25 Oct 2018
Multi-scale control variate methods for uncertainty quantification in kinetic equations arxiv:80.0844v [math.na] 25 Oct 208 Giacomo Dimarco and Lorenzo Pareschi October 26, 208 Abstract Kinetic equations
More informationFluid Equations for Rarefied Gases
1 Fluid Equations for Rarefied Gases Jean-Luc Thiffeault Department of Applied Physics and Applied Mathematics Columbia University http://plasma.ap.columbia.edu/~jeanluc 21 May 2001 with E. A. Spiegel
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 informationExponential Runge-Kutta for inhomogeneous Boltzmann equations with high order of accuracy
Exponential Runge-Kutta for inhomogeneous Boltzmann equations with high order of accuracy Qin Li, Lorenzo Pareschi Abstract We consider the development of exponential methods for the robust time discretization
More informationKinetic theory of gases
Kinetic theory of gases Toan T. Nguyen Penn State University http://toannguyen.org http://blog.toannguyen.org Graduate Student seminar, PSU Jan 19th, 2017 Fall 2017, I teach a graduate topics course: same
More informationFluid Equations for Rarefied Gases
1 Fluid Equations for Rarefied Gases Jean-Luc Thiffeault Department of Applied Physics and Applied Mathematics Columbia University http://plasma.ap.columbia.edu/~jeanluc 23 March 2001 with E. A. Spiegel
More information19 th INTERNATIONAL CONGRESS ON ACOUSTICS MADRID, 2-7 SEPTEMBER 2007
19 th INTERNATIONAL CONGRESS ON ACOUSTICS MADRID, -7 SEPTEMBER 007 INVESTIGATION OF AMPLITUDE DEPENDENCE ON NONLINEAR ACOUSTICS USING THE DIRECT SIMULATION MONTE CARLO METHOD PACS: 43.5.Ed Hanford, Amanda
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 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 information(2014) : ISSN
Wu, Lei and Reese, Jason and Zhang, Yonghao (214) Solving the Boltzmann equation deterministically by the fast spectral method : Application to gas microflows. Journal of Fluid Mechanics, 746. pp. 53-84.
More informationSOLUTIONS OF THE LINEAR BOLTZMANN EQUATION AND SOME DIRICHLET SERIES
SOLUTIONS OF THE LINEAR BOLTZMANN EQUATION AND SOME DIRICHLET SERIES A.V. BOBYLEV( ) AND I.M. GAMBA( ) Abstract. It is shown that a broad class of generalized Dirichlet series (including the polylogarithm,
More informationON THE MINIMIZATION PROBLEM OF SUB-LINEAR CONVEX FUNCTIONALS. Naoufel Ben Abdallah. Irene M. Gamba. Giuseppe Toscani. (Communicated by )
Kinetic and Related Models c American Institute of Mathematical Sciences Volume, Number, doi: pp. X XX ON THE MINIMIZATION PROBLEM OF SUB-LINEAR CONVEX FUNCTIONALS Naoufel Ben Abdallah Laboratoire MIP,
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 informationSpectral - Lagrangian methods for Collisional Models of Non - Equilibrium Statistical States
Spectral - Lagrangian methods for Collisional Models of Non - Equilibrium Statistical States Irene M. Gamba Dept. of Mathematics & Institute of Computational Engineering and Sciences, University of Texas
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 informationA comparative study of no-time-counter and majorant collision frequency numerical schemes in DSMC
Purdue University Purdue e-pubs School of Aeronautics and Astronautics Faculty Publications School of Aeronautics and Astronautics 2012 A comparative study of no-time-counter and majorant collision frequency
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 informationAn Asymptotic-Preserving Monte Carlo Method for the Boltzmann Equation
An Asymptotic-Preserving Monte Carlo Method for the Boltzmann Equation Wei Ren a, Hong Liu a,, Shi Jin b,c a J C Wu Center for Aerodynamics, School of Aeronautics and Aerospace, Shanghai Jiao Tong University,
More informationAnalisis para la ecuacion de Boltzmann Soluciones y Approximaciones
Analisis para la ecuacion de Boltzmann Soluciones y Approximaciones Irene M. Gamba Department of Mathematics and ICES The University of Texas at Austin Buenos Aires, June 2012 Collaborators: R. Alonso,
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 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 informationKinetic Models and Gas-Kinetic Schemes with Rotational Degrees of Freedom for Hybrid Continuum/Kinetic Boltzmann Methods
Kinetic Models and Gas-Kinetic Schemes with Rotational Degrees of Freedom for Hybrid Continuum/Kinetic Boltzmann Methods Simone Colonia, René Steijl and George N. Barakos CFD Laboratory - School of Engineering
More informationPhysics 607 Exam 2. ( ) = 1, Γ( z +1) = zγ( z) x n e x2 dx = 1. e x2
Physics 607 Exam Please be well-organized, and show all significant steps clearly in all problems. You are graded on your work, so please do not just write down answers with no explanation! Do all your
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 informationMicro-macro methods for Boltzmann-BGK-like equations in the diffusion scaling
Micro-macro methods for Boltzmann-BGK-like equations in the diffusion scaling Anaïs Crestetto 1, Nicolas Crouseilles 2, Giacomo Dimarco 3 et Mohammed Lemou 4 Saint-Malo, 14 décembre 2017 1 Université de
More informationHolder regularity for hypoelliptic kinetic equations
Holder regularity for hypoelliptic kinetic equations Alexis F. Vasseur Joint work with François Golse, Cyril Imbert, and Clément Mouhot The University of Texas at Austin Kinetic Equations: Modeling, Analysis
More informationAnomalous transport of particles in Plasma physics
Anomalous transport of particles in Plasma physics L. Cesbron a, A. Mellet b,1, K. Trivisa b, a École Normale Supérieure de Cachan Campus de Ker Lann 35170 Bruz rance. b Department of Mathematics, University
More informationA Unified Gas-kinetic Scheme for Continuum and Rarefied Flows
A Unified Gas-kinetic Scheme for Continuum and Rarefied Flows K. Xu and J.C. Huang Mathematics Department, Hong Kong University of Science and Technology, Hong Kong Department of Merchant Marine, National
More informationRevisit to Grad s Closure and Development of Physically Motivated Closure for Phenomenological High-Order Moment Model
Revisit to Grad s Closure and Development of Physically Motivated Closure for Phenomenological High-Order Moment Model R. S. Myong a and S. P. Nagdewe a a Dept. of Mechanical and Aerospace Engineering
More informationQUANTUM MODELS FOR SEMICONDUCTORS AND NONLINEAR DIFFUSION EQUATIONS OF FOURTH ORDER
QUANTUM MODELS FOR SEMICONDUCTORS AND NONLINEAR DIFFUSION EQUATIONS OF FOURTH ORDER MARIA PIA GUALDANI The modern computer and telecommunication industry relies heavily on the use of semiconductor devices.
More information1 Phase Spaces and the Liouville Equation
Phase Spaces and the Liouville Equation emphasize the change of language from deterministic to probablistic description. Under the dynamics: ½ m vi = F i ẋ i = v i with initial data given. What is the
More informationMath Questions for the 2011 PhD Qualifier Exam 1. Evaluate the following definite integral 3" 4 where! ( x) is the Dirac! - function. # " 4 [ ( )] dx x 2! cos x 2. Consider the differential equation dx
More informationAccurate Deterministic Solutions for the Classic Boltzmann Shock Profile
City University of New York (CUNY) CUNY Academic Works Dissertations, Theses, and Capstone Projects Graduate Center 9-2015 Accurate Deterministic Solutions for the Classic Boltzmann Shock Profile Yubei
More informationMACROSCOPIC VARIABLES, THERMAL EQUILIBRIUM. Contents AND BOLTZMANN ENTROPY. 1 Macroscopic Variables 3. 2 Local quantities and Hydrodynamics fields 4
MACROSCOPIC VARIABLES, THERMAL EQUILIBRIUM AND BOLTZMANN ENTROPY Contents 1 Macroscopic Variables 3 2 Local quantities and Hydrodynamics fields 4 3 Coarse-graining 6 4 Thermal equilibrium 9 5 Two systems
More informationQuantum Hydrodynamics models derived from the entropy principle
1 Quantum Hydrodynamics models derived from the entropy principle P. Degond (1), Ch. Ringhofer (2) (1) MIP, CNRS and Université Paul Sabatier, 118 route de Narbonne, 31062 Toulouse cedex, France degond@mip.ups-tlse.fr
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 informationUne décomposition micro-macro particulaire pour des équations de type Boltzmann-BGK en régime de diffusion
Une décomposition micro-macro particulaire pour des équations de type Boltzmann-BGK en régime de diffusion Anaïs Crestetto 1, Nicolas Crouseilles 2 et Mohammed Lemou 3 La Tremblade, Congrès SMAI 2017 5
More informationAsymptotic-Preserving Exponential Methods for the Quantum Boltzmann Equation with High-Order Accuracy
J Sci Comput () 6: 74 DOI.7/s9-4-9869- Asymptotic-Preserving Eponential Methods for the Quantum Boltzmann Equation with High-Order Accuracy Jingwei Hu Qin Li Lorenzo Pareschi Received: 3 October 3 / Revised:
More informationA discontinuous Galerkin fast spectral method for the multi-species Boltzmann equation
A discontinuous Galerkin fast spectral method for the multi-species Boltzmann equation Shashank Jaiswal a, Alina A. Alexeenko a, Jingwei Hu b, a School of Aeronautics and Astronautics, Purdue University,
More informationFluctuating Hydrodynamics and Direct Simulation Monte Carlo
Fluctuating Hydrodynamics and Direct Simulation Monte Carlo Kaushi Balarishnan Lawrence Bereley Lab John B. Bell Lawrence Bereley Lab Alesandar Donev New Yor University Alejandro L. Garcia San Jose State
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 informationOn a New Diagram Notation for the Derivation of Hyperbolic Moment Models
On a New Diagram Notation for the Derivation of Hyperbolic Moment Models Julian Koellermeier, Manuel Torrilhon, Yuwei Fan March 17th, 2017 Stanford University J. Koellermeier 1 / 57 of Hyperbolic Moment
More informationDerivation of quantum hydrodynamic equations with Fermi-Dirac and Bose-Einstein statistics
Derivation of quantum hydrodynamic equations with Fermi-Dirac and Bose-Einstein statistics Luigi Barletti (Università di Firenze) Carlo Cintolesi (Università di Trieste) 6th MMKT Porto Ercole, june 9th
More informationSOME APPLICATIONS OF THE METHOD OF MOMENTS FOR THE HOMOGENEOUS BOLTZMANN AND KAC EQUATIONS
SOME APPLICATIONS OF THE METHOD OF MOMENTS FOR THE HOMOGENEOUS BOLTZMANN AND KAC EQUATIONS Laurent Desvillettes ECOLE NORMALE SUPERIEURE 45, Rue d Ulm 7530 Paris Cédex 05 March 6, 013 Abstract Using the
More informationComparison of the DSMC Method with an Exact Solution of the Boltzmann Equation
Comparison of the DSMC Method with an Exact Solution of the Boltzmann Equation J M Montanero and A Santos Departamento de Fisica, Universidad de Extremadura, Spain Abstract Results obtained from the Direct
More informationOn the Interior Boundary-Value Problem for the Stationary Povzner Equation with Hard and Soft Interactions
On the Interior Boundary-Value Problem for the Stationary Povzner Equation with Hard and Soft Interactions Vladislav A. Panferov Department of Mathematics, Chalmers University of Technology and Göteborg
More informationA Solution of the Two-dimensional Boltzmann Transport Equation
Applied Mathematical Sciences, Vol. 9, 2015, no. 147, 7347-7355 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.12988/ams.2015.510665 A Solution of the Two-dimensional Boltzmann Transport Equation F.
More informationNumerical simulation of the dynamics of the rotating dipolar Bose-Einstein condensates
Numerical simulation of the dynamics of the rotating dipolar Bose-Einstein condensates Qinglin TANG Inria, University of Lorraine Joint work with: Prof. Weizhu BAO, Daniel MARAHRENS, Yanzhi ZHANG and Yong
More informationHyperbolic Approximation of Kinetic Equations Using Quadrature-Based Projection Methods
Master Thesis Hyperbolic Approximation of Kinetic Equations Using Quadrature-Based Projection Methods Julian Köllermeier RWTH Aachen University Supervisor: Prof. Dr. Manuel Torrilhon Center for Computational
More informationNumerical Simulation of Rarefied Gases using Hyperbolic Moment Equations in Partially-Conservative Form
Numerical Simulation of Rarefied Gases using Hyperbolic Moment Equations in Partially-Conservative Form Julian Koellermeier, Manuel Torrilhon May 18th, 2017 FU Berlin J. Koellermeier 1 / 52 Partially-Conservative
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 informationSolution of time-dependent Boltzmann equation for electrons in non-thermal plasma
Solution of time-dependent Boltzmann equation for electrons in non-thermal plasma Z. Bonaventura, D. Trunec Department of Physical Electronics Faculty of Science Masaryk University Kotlářská 2, Brno, CZ-61137,
More informationStochastic equations for thermodynamics
J. Chem. Soc., Faraday Trans. 93 (1997) 1751-1753 [arxiv 1503.09171] Stochastic equations for thermodynamics Roumen Tsekov Department of Physical Chemistry, University of Sofia, 1164 Sofia, ulgaria The
More informationA hybrid method for hydrodynamic-kinetic flow - Part II - Coupling of hydrodynamic and kinetic models
A hybrid method for hydrodynamic-kinetic flow - Part II - Coupling of hydrodynamic and kinetic models Alessandro Alaia, Gabriella Puppo May 31, 2011 Abstract In this work we present a non stationary domain
More informationGlobal Maxwellians over All Space and Their Relation to Conserved Quantites of Classical Kinetic Equations
Global Maxwellians over All Space and Their Relation to Conserved Quantites of Classical Kinetic Equations C. David Levermore Department of Mathematics and Institute for Physical Science and Technology
More informationRapporto di ricerca Research report
Dipartimento di Informatica Università degli Studi di Verona Rapporto di ricerca Research report September 2009 76/2009 Spectral methods for dissipative nonlinear Schrödinger equations Laura M. Morato
More informationAn asymptotic-preserving micro-macro scheme for Vlasov-BGK-like equations in the diffusion scaling
An asymptotic-preserving micro-macro scheme for Vlasov-BGK-like equations in the diffusion scaling Anaïs Crestetto 1, Nicolas Crouseilles 2 and Mohammed Lemou 3 Saint-Malo 13 December 2016 1 Université
More information16. GAUGE THEORY AND THE CREATION OF PHOTONS
6. GAUGE THEORY AD THE CREATIO OF PHOTOS In the previous chapter the existence of a gauge theory allowed the electromagnetic field to be described in an invariant manner. Although the existence of this
More informationFokker-Planck collision operator
DRAFT 1 Fokker-Planck collision operator Felix I. Parra Rudolf Peierls Centre for Theoretical Physics, University of Oxford, Oxford OX1 3NP, UK (This version is of 16 April 18) 1. Introduction In these
More informationSEMICLASSICAL LATTICE BOLTZMANN EQUATION HYDRODYNAMICS
The Seventh Nobeyama Workshop on CFD: To the Memory of Prof. Kuwahara Sanjo Kaikan, the University of Tokyo, Japan, September. 23-24, 2009 SEMICLASSICAL LATTICE BOLTZMANN EQUATION HYDRODYNAMICS Jaw-Yen
More informationPlasmas as fluids. S.M.Lea. January 2007
Plasmas as fluids S.M.Lea January 2007 So far we have considered a plasma as a set of non intereacting particles, each following its own path in the electric and magnetic fields. Now we want to consider
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 informationEXPONENTIAL AND ALGEBRAIC RELAXATION IN KINETIC MODELS FOR WEALTH DISTRIBUTION
1 EXPONENTIAL AND ALGEBAIC ELAXATION IN KINETIC MODELS FO WEALTH DISTIBUTION B. DÜING Institut für Analysis und Scientific Computing, Technische Universität Wien, 1040 Wien, Austria. D. MATTHES and G.
More informationUncertainty Quantification and hypocoercivity based sensitivity analysis for multiscale kinetic equations with random inputs.
Uncertainty Quantification and hypocoercivity based sensitivity analysis for multiscale kinetic equations with random inputs Shi Jin University of Wisconsin-Madison, USA Shanghai Jiao Tong University,
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 informationAnalisis para la ecuacion de Boltzmann Soluciones y Approximaciones
Analisis para la ecuacion de Boltzmann Soluciones y Approximaciones Irene M. Gamba Department of Mathematics and ICES The University of Texas at Austin Mar del PLata, Julio 2012 Collaborators: R. Alonso,
More informationDirect Simulation Monte Carlo Calculation: Strategies for Using Complex Initial Conditions
Mat. Res. Soc. Symp. Proc. Vol. 731 2002 Materials Research Society Direct Simulation Monte Carlo Calculation: Strategies for Using Complex Initial Conditions Michael I. Zeifman 1, Barbara J. Garrison
More informationOSGOOD TYPE REGULARITY CRITERION FOR THE 3D NEWTON-BOUSSINESQ EQUATION
Electronic Journal of Differential Equations, Vol. 013 (013), No. 3, pp. 1 6. ISSN: 107-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu OSGOOD TYPE REGULARITY
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 informationAn improved unified gas-kinetic scheme and the study of shock structures
IMA Journal of Applied Mathematics (2011) 76, 698 711 doi:10.1093/imamat/hxr002 Advance Access publication on March 16, 2011 An improved unified gas-kinetic scheme and the study of shock structures KUN
More informationA class of asymptotic-preserving schemes for the Fokker-Planck-Landau equation
A class of asymptotic-preserving schemes for the Fokker-Planck-Landau equation Shi Jin Bokai Yan October 2, 200 Abstract We present a class of asymptotic-preserving AP) schemes for the nonhomogeneous Fokker-Planck-Landau
More informationThe inviscid limit to a contact discontinuity for the compressible Navier-Stokes-Fourier system using the relative entropy method
The inviscid limit to a contact discontinuity for the compressible Navier-Stokes-Fourier system using the relative entropy method Alexis Vasseur, and Yi Wang Department of Mathematics, University of Texas
More informationPOSITIVITY-PRESERVING DISCONTINUOUS GALERKIN SCHEMES FOR LINEAR VLASOV-BOLTZMANN TRANSPORT EQUATIONS
POSITIVITY-PRESERVING DISCONTINUOUS GALERKIN SCHEMES FOR LINEAR VLASOV-BOLTZMANN TRANSPORT EQUATIONS YINGDA CHENG, IRENE M. GAMBA, AND JENNIFER PROFT Abstract. We develop a high-order positivity-preserving
More informationUncertainty Quantification for multiscale kinetic equations with random inputs. Shi Jin. University of Wisconsin-Madison, USA
Uncertainty Quantification for multiscale kinetic equations with random inputs Shi Jin University of Wisconsin-Madison, USA Where do kinetic equations sit in physics Kinetic equations with applications
More information4. The Green Kubo Relations
4. The Green Kubo Relations 4.1 The Langevin Equation In 1828 the botanist Robert Brown observed the motion of pollen grains suspended in a fluid. Although the system was allowed to come to equilibrium,
More information