Multiscale hp-adaptive quiescent reactors
|
|
- Tobias Sullivan
- 5 years ago
- Views:
Transcription
1 Multiscale hp-adaptive quiescent reactors C.E.Michoski Institute for Computational Engineering and Sciences (ICES) University of Texas (UT) at Austin Group Meeting Nov. 7, 2011
2 Outline 1 2 Equation (or normal) ordering Conservation laws 3 Fractional multistep operator splitting Law of mass action Mass diffusivty 4 The spatial discretization The semidiscrete formulation The discrete formulations Entropy enriched hp-adaptivity and stability 5
3 What are quiescent reactors? Reactors with no mixing
4 What are quiescent reactors? Generally concerned with systems (i.e. n > 1 species) of reaction-diffusion equations, which are nonlinear in the reaction term. Such systems admit discontinous solutions, with physical meaning, such as phase-changes, action-potentials, oscillatory chemical reactions,... In spatial dimension N > 2, the Hausdorff dimension of the singularity set for quadratic reactions has a known upper bound dim H V N 2 4/N. Not all reaction-diffusion systems satisfy the requisite underlying assumptions, leading to so-called unphysical solutions in the physicist s lingo.
5 What are quiescent reactors? Generally concerned with systems (i.e. n > 1 species) of reaction-diffusion equations, which are nonlinear in the reaction term. Such systems admit discontinous solutions, with physical meaning, such as phase-changes, action-potentials, oscillatory chemical reactions,... In spatial dimension N > 2, the Hausdorff dimension of the singularity set for quadratic reactions has a known upper bound dim H V N 2 4/N. Not all reaction-diffusion systems satisfy the requisite underlying assumptions, leading to so-called unphysical solutions in the physicist s lingo.
6 What are quiescent reactors? Generally concerned with systems (i.e. n > 1 species) of reaction-diffusion equations, which are nonlinear in the reaction term. Such systems admit discontinous solutions, with physical meaning, such as phase-changes, action-potentials, oscillatory chemical reactions,... In spatial dimension N > 2, the Hausdorff dimension of the singularity set for quadratic reactions has a known upper bound dim H V N 2 4/N. Not all reaction-diffusion systems satisfy the requisite underlying assumptions, leading to so-called unphysical solutions in the physicist s lingo.
7 What are quiescent reactors? Generally concerned with systems (i.e. n > 1 species) of reaction-diffusion equations, which are nonlinear in the reaction term. Such systems admit discontinous solutions, with physical meaning, such as phase-changes, action-potentials, oscillatory chemical reactions,... In spatial dimension N > 2, the Hausdorff dimension of the singularity set for quadratic reactions has a known upper bound dim H V N 2 4/N. Not all reaction-diffusion systems satisfy the requisite underlying assumptions, leading to so-called unphysical solutions in the physicist s lingo.
8 Outline 1 Equation (or normal) ordering Conservation laws 2 Equation (or normal) ordering Conservation laws 3 Fractional multistep operator splitting Law of mass action Mass diffusivty 4 The spatial discretization The semidiscrete formulation The discrete formulations Entropy enriched hp-adaptivity and stability 5
9 Where does it come from? Equation (or normal) ordering Conservation laws The idea is that physical systems must obey a normal ordering of the physics, so that each system must be a direct consequence of a more fundamental system. If not, the system gets labelled unphysical. If so, we say it follows from first principles. e.g. QFT Standard Model Quantum Mechanics Statistical mechanics Continuum Mechanics Macroscopic continuum e.g. abstract mathematics canonical quantization KG-Dirac equations Schrödinger equation Boltzmann Equation Navier-Stokes equations Shallow Water equations
10 Where does it come from? Equation (or normal) ordering Conservation laws The idea is that physical systems must obey a normal ordering of the physics, so that each system must be a direct consequence of a more fundamental system. If not, the system gets labelled unphysical. If so, we say it follows from first principles. e.g. QFT Standard Model Quantum Mechanics Statistical mechanics Continuum Mechanics Macroscopic continuum e.g. abstract mathematics canonical quantization KG-Dirac equations Schrödinger equation Boltzmann Equation Navier-Stokes equations Shallow Water equations
11 Where does it come from? Equation (or normal) ordering Conservation laws The idea is that physical systems must obey a normal ordering of the physics, so that each system must be a direct consequence of a more fundamental system. If not, the system gets labelled unphysical. If so, we say it follows from first principles. e.g. QFT Standard Model Quantum Mechanics Statistical mechanics Continuum Mechanics Macroscopic continuum e.g. abstract mathematics canonical quantization KG-Dirac equations Schrödinger equation Boltzmann Equation Navier-Stokes equations Shallow Water equations
12 The species Boltzmann equations Equation (or normal) ordering Conservation laws Consider a system of n PDEs in Ω Ω (0, T ), where Ω R N such that: t f i + v i x f i = S i (f ) + C i (f ). in position x R N, species velocity v i R N, and time t (0, T ) over index i = 1,..., n. Distribution functions: n distribution functions f i = f (t, x, v i ) Reactive collisional integral C i = C i (f 1,..., f n ) Scattering collisional integral: S i = S i (f 1,..., f n )
13 The species Boltzmann equations Equation (or normal) ordering Conservation laws Consider a system of n PDEs in Ω Ω (0, T ), where Ω R N such that: t f i + v i x f i = S i (f ) + C i (f ). in position x R N, species velocity v i R N, and time t (0, T ) over index i = 1,..., n. Distribution functions: n distribution functions f i = f (t, x, v i ) Reactive collisional integral C i = C i (f 1,..., f n ) Scattering collisional integral: S i = S i (f 1,..., f n )
14 The species Boltzmann equations Equation (or normal) ordering Conservation laws Consider a system of n PDEs in Ω Ω (0, T ), where Ω R N such that: t f i + v i x f i = S i (f ) + C i (f ). in position x R N, species velocity v i R N, and time t (0, T ) over index i = 1,..., n. Distribution functions: n distribution functions f i = f (t, x, v i ) Reactive collisional integral C i = C i (f 1,..., f n ) Scattering collisional integral: S i = S i (f 1,..., f n )
15 The species Boltzmann equations Equation (or normal) ordering Conservation laws Consider a system of n PDEs in Ω Ω (0, T ), where Ω R N such that: t f i + v i x f i = S i (f ) + C i (f ). in position x R N, species velocity v i R N, and time t (0, T ) over index i = 1,..., n. Distribution functions: n distribution functions f i = f (t, x, v i ) Reactive collisional integral C i = C i (f 1,..., f n ) Scattering collisional integral: S i = S i (f 1,..., f n )
16 The formal Enskog Expansion Equation (or normal) ordering Conservation laws Then we consider the usual Enskog expansion of to linear order t f i +v i x f i = ɛ 1 S i (f )+ɛ b C i (f ), and f i = fi 0 ( 1 + ɛηi + O(ɛ 2 ) ), where ɛ is the formal expansion parameter. The perturbations η i are used to determine the respective forms of the corresponding transport coefficients. When b = 0 then we are in the so-called strong reaction regime When b = 1 we are in the Maxwellian reaction regime When b = 1 we are in the kinetic chemical equilibrium regime.
17 The formal Enskog Expansion Equation (or normal) ordering Conservation laws Then we consider the usual Enskog expansion of to linear order t f i +v i x f i = ɛ 1 S i (f )+ɛ b C i (f ), and f i = fi 0 ( 1 + ɛηi + O(ɛ 2 ) ), where ɛ is the formal expansion parameter. The perturbations η i are used to determine the respective forms of the corresponding transport coefficients. When b = 0 then we are in the so-called strong reaction regime When b = 1 we are in the Maxwellian reaction regime When b = 1 we are in the kinetic chemical equilibrium regime.
18 The formal Enskog Expansion Equation (or normal) ordering Conservation laws Then we consider the usual Enskog expansion of to linear order t f i +v i x f i = ɛ 1 S i (f )+ɛ b C i (f ), and f i = fi 0 ( 1 + ɛηi + O(ɛ 2 ) ), where ɛ is the formal expansion parameter. The perturbations η i are used to determine the respective forms of the corresponding transport coefficients. When b = 0 then we are in the so-called strong reaction regime When b = 1 we are in the Maxwellian reaction regime When b = 1 we are in the kinetic chemical equilibrium regime.
19 The formal Enskog Expansion Equation (or normal) ordering Conservation laws Then we consider the usual Enskog expansion of to linear order t f i +v i x f i = ɛ 1 S i (f )+ɛ b C i (f ), and f i = fi 0 ( 1 + ɛηi + O(ɛ 2 ) ), where ɛ is the formal expansion parameter. The perturbations η i are used to determine the respective forms of the corresponding transport coefficients. When b = 0 then we are in the so-called strong reaction regime When b = 1 we are in the Maxwellian reaction regime When b = 1 we are in the kinetic chemical equilibrium regime.
20 Equation (or normal) ordering Conservation laws Approx. representation by continuum mechanics After expanding with respect to the collisional invariants of the species Boltzmann equation (keeping terms of the same order ɛ), we arrive with the species Navier Stokes equations: t ρ i + x (ρ i (u V i )) δ b0 m i A i ( n) = 0, t (ρu) + x (ρu u) + x S = 0, t (ρe) + x ((ρe + S)u) + x Q = 0 Notice that this is not a system of reaction-diffusion equations.
21 The quiescent reactor constraints Equation (or normal) ordering Conservation laws In order to restrict to a true reaction diffusion system, the following constraints must hold: n (1) C i (f 0 ) f C i (f 0 ) i, (2) v i dv i 0, (3) x ρ i x ϑ i, (4) h i (n i m i pd ij ) 1 θ i. i Ω Constraint (1) is standard and observable at high energy. Constraint (2) assumes that the global flow velocity averages to zero over the entire domain (no appreciable mixing). Constraints (3) and (4) restrict to solutions that constrain the admissible bounds on the thermal gradients
22 The quiescent reactor constraints Equation (or normal) ordering Conservation laws In order to restrict to a true reaction diffusion system, the following constraints must hold: n (1) C i (f 0 ) f C i (f 0 ) i, (2) v i dv i 0, (3) x ρ i x ϑ i, (4) h i (n i m i pd ij ) 1 θ i. i Ω Constraint (1) is standard and observable at high energy. Constraint (2) assumes that the global flow velocity averages to zero over the entire domain (no appreciable mixing). Constraints (3) and (4) restrict to solutions that constrain the admissible bounds on the thermal gradients
23 The quiescent reactor constraints Equation (or normal) ordering Conservation laws In order to restrict to a true reaction diffusion system, the following constraints must hold: n (1) C i (f 0 ) f C i (f 0 ) i, (2) v i dv i 0, (3) x ρ i x ϑ i, (4) h i (n i m i pd ij ) 1 θ i. i Ω Constraint (1) is standard and observable at high energy. Constraint (2) assumes that the global flow velocity averages to zero over the entire domain (no appreciable mixing). Constraints (3) and (4) restrict to solutions that constrain the admissible bounds on the thermal gradients
24 The quiescent reactor constraints Equation (or normal) ordering Conservation laws In order to restrict to a true reaction diffusion system, the following constraints must hold: n (1) C i (f 0 ) f C i (f 0 ) i, (2) v i dv i 0, (3) x ρ i x ϑ i, (4) h i (n i m i pd ij ) 1 θ i. i Ω Constraint (1) is standard and observable at high energy. Constraint (2) assumes that the global flow velocity averages to zero over the entire domain (no appreciable mixing). Constraints (3) and (4) restrict to solutions that constrain the admissible bounds on the thermal gradients
25 Equation (or normal) ordering Conservation laws (the equations for...) The we arrive with a solution over (t, x) (0, T ) Ω for Ω R N chosen to satisfy: t ρ i x (D i x ρ i ) A i ( n) = 0, n n A i ( n) = m i (νir b νir f ) k fr k br r R with initial-boundary data on Ω, ρ i (t = 0) = ρ i,0, and a 1i ρ i,b + x ρ i,b (a 2i n + a 3i τ ) = a 4i, taking arbitrary functions a ji = a ji (t, x b ) for j {1, 2, 3, 4}, where n is the unit outward normal and τ the unit tangent vector. j=1 n νf jr j j=1 n νb jr j,
26 Variables and parameters Equation (or normal) ordering Conservation laws n i is the molar concentration of the i-th chemical constituent, which up to a scaling by Avogadro s constant N A is just the number density n i = N A n i. The species densities are given by ρ i = ρα i = m i n i, where α i is the mass fraction of the i-th species, and m i is the molar mass of the i-th species. The D i are the interspecies mass diffusivity coefficients, which have a complicated form. The A i ( n) is the chemical mass action, which depends on the reaction rates k r and the realtive constituents availability ν jr.
27 Variables and parameters Equation (or normal) ordering Conservation laws n i is the molar concentration of the i-th chemical constituent, which up to a scaling by Avogadro s constant N A is just the number density n i = N A n i. The species densities are given by ρ i = ρα i = m i n i, where α i is the mass fraction of the i-th species, and m i is the molar mass of the i-th species. The D i are the interspecies mass diffusivity coefficients, which have a complicated form. The A i ( n) is the chemical mass action, which depends on the reaction rates k r and the realtive constituents availability ν jr.
28 Variables and parameters Equation (or normal) ordering Conservation laws n i is the molar concentration of the i-th chemical constituent, which up to a scaling by Avogadro s constant N A is just the number density n i = N A n i. The species densities are given by ρ i = ρα i = m i n i, where α i is the mass fraction of the i-th species, and m i is the molar mass of the i-th species. The D i are the interspecies mass diffusivity coefficients, which have a complicated form. The A i ( n) is the chemical mass action, which depends on the reaction rates k r and the realtive constituents availability ν jr.
29 Variables and parameters Equation (or normal) ordering Conservation laws n i is the molar concentration of the i-th chemical constituent, which up to a scaling by Avogadro s constant N A is just the number density n i = N A n i. The species densities are given by ρ i = ρα i = m i n i, where α i is the mass fraction of the i-th species, and m i is the molar mass of the i-th species. The D i are the interspecies mass diffusivity coefficients, which have a complicated form. The A i ( n) is the chemical mass action, which depends on the reaction rates k r and the realtive constituents availability ν jr.
30 Outline 1 Equation (or normal) ordering Conservation laws 2 Equation (or normal) ordering Conservation laws 3 Fractional multistep operator splitting Law of mass action Mass diffusivty 4 The spatial discretization The semidiscrete formulation The discrete formulations Entropy enriched hp-adaptivity and stability 5
31 The form of chemical equations Equation (or normal) ordering Conservation laws The forward and backward stoichiometric coefficients of elementary reaction r N are given by ν f ir N and νb ir N, while k fr, k br R are the respective forward and backward reaction rates of reaction r. The indexing sets R r and P r as the reactant well and the product well R r N and P r N for reaction r. Then, j R r ν f jr M j k fr k br k P r ν b kr M k, r R. (1)
32 The form of chemical equations Equation (or normal) ordering Conservation laws The forward and backward stoichiometric coefficients of elementary reaction r N are given by ν f ir N and νb ir N, while k fr, k br R are the respective forward and backward reaction rates of reaction r. The indexing sets R r and P r as the reactant well and the product well R r N and P r N for reaction r. Then, j R r ν f jr M j k fr k br k P r ν b kr M k, r R. (1)
33 Conservation principles Equation (or normal) ordering Conservation laws Since the elementary reactions are balanced the conservation of atoms in the system is an immediate consequence of (1). Let a il be the l-th atom of the i-th species M i, where l A r is the indexing set A r = {1, 2,..., n atoms,r } of the distinct atoms present in each reaction r R. Then the total atom conservation is satisfied for every atom in every reaction a il νir f = a il νir b r R, l A r. i R r i P r Since the total number of atoms is conserved, so is the total mass, m i νir f = m i νir b r R. i R r i P r
34 Conservation principles, cont. Equation (or normal) ordering Conservation laws It then immediately follows that an integration by parts yields the following bulk conservation principle which is satisfied globally: d n ρ i dx = 0. (2) dt i=1 Ω Moreover, the system of chemical reactions satisfy spontaneity G 0 (up to a constant) with respect to the standard total Gibb s free reaction energy of the system, which further provides that the entropy of the system dissipates.
35 Conservation principles, cont. Equation (or normal) ordering Conservation laws It then immediately follows that an integration by parts yields the following bulk conservation principle which is satisfied globally: d n ρ i dx = 0. (2) dt i=1 Ω Moreover, the system of chemical reactions satisfy spontaneity G 0 (up to a constant) with respect to the standard total Gibb s free reaction energy of the system, which further provides that the entropy of the system dissipates.
36 Outline 1 Fractional multistep operator splitting Law of mass action Mass diffusivty 2 Equation (or normal) ordering Conservation laws 3 Fractional multistep operator splitting Law of mass action Mass diffusivty 4 The spatial discretization The semidiscrete formulation The discrete formulations Entropy enriched hp-adaptivity and stability 5
37 Fractional multistep operator splitting Law of mass action Mass diffusivty Now we consider the multiscale solutions that split over fast ρ f and slow ρ s modes. Assume that ρ = ρ f + ρ s, such that (ρ f ) t x (D f σ f ) A f = 0, and σ f x ρ f = 0, (ρ s ) t x (D s σ s ) A s = 0, and σ s x ρ s = 0. Then the splitting of the reaction part at time t provides { (ρf ) reaction modes t A f = 0, (ρ s ) t A s = 0. and the diffusion part splits as: { (ρf ) diffusion modes t x (D f σ f ) = 0, with σ f x ρ f = 0, (ρ s ) t x (D s σ s ) = 0, with σ s x ρ s = 0.
38 Composition of modes Fractional multistep operator splitting Law of mass action Mass diffusivty Note that the operation of composition here is given in the natural way, such that (at second order) D t/2 (ρ) R t (ρ) D t/2 (ρ) means we solve the system: (ρ ) t x (Dσ ) = 0, ρ (0) = ρ 0 on [0, t/2] (ρ ) t A = 0, ρ (0) = ρ 0 on [0, t], (ρ ) t x (Dσ ) = 0, ρ (0) = ρ 0 on [0, t/2], where the solution to the composition is then given by ρ (t/2) = ρ (t/2, x).
39 Consistency of multiscale splitting Fractional multistep operator splitting Law of mass action Mass diffusivty Reaction Regnant Diffusion Limited Transition Regimes Diffusion Dominated Shock Dominated
40 Outline 1 Fractional multistep operator splitting Law of mass action Mass diffusivty 2 Equation (or normal) ordering Conservation laws 3 Fractional multistep operator splitting Law of mass action Mass diffusivty 4 The spatial discretization The semidiscrete formulation The discrete formulations Entropy enriched hp-adaptivity and stability 5
41 nfanodes Fractional multistep operator splitting Law of mass action Mass diffusivty The law of mass action splits into an n-coupled system of first order ordinary differential equations (nfanodes), which have no spatial dependence. As such, developing a solution technique generally rests somewhere between: (1) finding a relatively straightforward approximate solution, and (2) analytically solving the difficult (though often soluble) system of nfanodes for a time-evolving solution to the split equation.
42 nfanodes Fractional multistep operator splitting Law of mass action Mass diffusivty The law of mass action splits into an n-coupled system of first order ordinary differential equations (nfanodes), which have no spatial dependence. As such, developing a solution technique generally rests somewhere between: (1) finding a relatively straightforward approximate solution, and (2) analytically solving the difficult (though often soluble) system of nfanodes for a time-evolving solution to the split equation.
43 nfanodes Fractional multistep operator splitting Law of mass action Mass diffusivty The law of mass action splits into an n-coupled system of first order ordinary differential equations (nfanodes), which have no spatial dependence. As such, developing a solution technique generally rests somewhere between: (1) finding a relatively straightforward approximate solution, and (2) analytically solving the difficult (though often soluble) system of nfanodes for a time-evolving solution to the split equation.
44 nfanodes Fractional multistep operator splitting Law of mass action Mass diffusivty The law of mass action splits into an n-coupled system of first order ordinary differential equations (nfanodes), which have no spatial dependence. As such, developing a solution technique generally rests somewhere between: (1) finding a relatively straightforward approximate solution, and (2) analytically solving the difficult (though often soluble) system of nfanodes for a time-evolving solution to the split equation.
45 Outline 1 Fractional multistep operator splitting Law of mass action Mass diffusivty 2 Equation (or normal) ordering Conservation laws 3 Fractional multistep operator splitting Law of mass action Mass diffusivty 4 The spatial discretization The semidiscrete formulation The discrete formulations Entropy enriched hp-adaptivity and stability 5
46 Complicated functional parameter Fractional multistep operator splitting Law of mass action Mass diffusivty The Enskog expansion tells us that the mass diffusivity (which obey Fick s second law of diffusion) also must satisfy: where p is the pressure, and D i = ρ i p 1 D ij ρi p, D ij = D ij (p, ϑ, v i, u, x p, x ϑ, x u, x v i ), correspon to the diffusivity matrix. Note, however that matrix diag(ρ i )D with diagonal components ρ i D ij is a symmetric positive semidefinite matrix.
47 Complicated functional parameter Fractional multistep operator splitting Law of mass action Mass diffusivty The Enskog expansion tells us that the mass diffusivity (which obey Fick s second law of diffusion) also must satisfy: where p is the pressure, and D i = ρ i p 1 D ij ρi p, D ij = D ij (p, ϑ, v i, u, x p, x ϑ, x u, x v i ), correspon to the diffusivity matrix. Note, however that matrix diag(ρ i )D with diagonal components ρ i D ij is a symmetric positive semidefinite matrix.
48 Outline 1 2 Equation (or normal) ordering Conservation laws 3 Fractional multistep operator splitting Law of mass action Mass diffusivty 4 The spatial discretization The semidiscrete formulation The discrete formulations Entropy enriched hp-adaptivity and stability 5 The spatial discretization The semidiscrete formulation The discrete formulations Entropy enriched hp-adaptivity and stability
49 The spatial discretization The semidiscrete formulation The discrete formulations Entropy enriched hp-adaptivity and stability Take test functions ϕ hp, ς hp, ϖ hp W k,q (Ω h, T h ) characterized by ϕ hp (x) = p ϕ i l Ni l (x), ς hp(x) = l=0 and ϖ hp (x) = p ς i l Ni l (x) l=0 p ϖ i l Ni l (x) x Ω e i, where ϕ i l, ςi l and ϖi l are the nodal values of the test functions in each Ω ei, and with the broken Sobolev space over the partition T h defined by W k,q (Ω h, T h ) = {w : w Ωei W k,q (Ω ei ), w w, Ω ei T h }. l=0
50 Find the variational form The spatial discretization The semidiscrete formulation The discrete formulations Entropy enriched hp-adaptivity and stability The time derivative terms yield: d ( ) ρhp, ς hp d ( ) dt Ω ρ, ςhp G d Ω G d ( ) and ρhp, ϕ dt hp d ( ) Ω ρ, ϕhp G dt and the mass action term gives, ( Ahp ( n), ϕ hp )Ω G ( A ( n), ϕ hp ) Ω G, and the mass diffusion term, ( x (I σ), ς hp )Ω G = Ω ei T h Ω G, Ω ei x (ς hp I σ)dx ( I σ, x ς hp ) Ω G.
51 One more integration by parts The spatial discretization The semidiscrete formulation The discrete formulations Entropy enriched hp-adaptivity and stability We approximate the first term on the right in the diffusion using the generalized diffusive flux G ij in the unified sense, such that G i = G i (I hp, σ hp, ρ hp, ς hp ), and we see that G i = j S(i) j S(i) Γ ij Gij (I hp, σ hp Γij, σ hp Γji, ρ hp Γij, ρ hp Γji, n ij ) ς hp Γij dξ Γ ij s=1 N (I hp σ) s (n ij ) s ς hp Γij dξ.
52 The spatial discretization The semidiscrete formulation The discrete formulations Entropy enriched hp-adaptivity and stability The interior diffusion term is approximated directly by ( ) ( ) H = H (I hp, σ h, ρ hp, ς hp ) = I hp σ hp, ς hp x I x ρ, ς hp x Ω G and the auxiliary numerical flux is taken so that: L i = L i ( L ij, σ hp, ρ hp, ϖ hp, ϖ hp x, n ij ) = ( ) σ hp, ϖ hp Ω ei ( + ρ hp, ϖ hp x )Ωei L (ρhp Γij, ρ hp Γji, ϖ hp Γij, n ij )dξ, where j S(i) Γ ij j S(i) Γ ij Lij (ρ hp Γij, ρ hp Γji, ϖ hp Γij, n ij )dξ j S(i) Γ ij s=1 N (ρ) s (n ij ) s ϖ hp Γij dξ. Ω G,
53 Outline 1 2 Equation (or normal) ordering Conservation laws 3 Fractional multistep operator splitting Law of mass action Mass diffusivty 4 The spatial discretization The semidiscrete formulation The discrete formulations Entropy enriched hp-adaptivity and stability 5 The spatial discretization The semidiscrete formulation The discrete formulations Entropy enriched hp-adaptivity and stability
54 Semidiscrete formulation The spatial discretization The semidiscrete formulation The discrete formulations Entropy enriched hp-adaptivity and stability For each t > 0, find the pair (ρ hp, σ hp ) such that The semidiscrete discontinuous Galerkin scheme a) ρ hp C 1 ([0, T ); S d h ), σ hp S d h, b) ρ hp (0) = Π hp ρ 0, d ( ) c) ρhp, ς hp = G + H, L = 0, dt Ω G d ( d) ρhp, ϕ dt hp = )Ω ( ) A hp ( n), ϕ G hp. Ω G
55 Outline 1 2 Equation (or normal) ordering Conservation laws 3 Fractional multistep operator splitting Law of mass action Mass diffusivty 4 The spatial discretization The semidiscrete formulation The discrete formulations Entropy enriched hp-adaptivity and stability 5 The spatial discretization The semidiscrete formulation The discrete formulations Entropy enriched hp-adaptivity and stability
56 RKSSP or RKC discrete form The spatial discretization The semidiscrete formulation The discrete formulations Entropy enriched hp-adaptivity and stability The discrete split explicit RKDG scheme a) ρ hp (0) = Π hp ρ 0, b) ρ n f hp = ρ hp(0), ( ( cd) ρ ns+1 hp, ς hp = )Ω ( G dc) ρ ns+1 hp, ϕ hp = )Ω G ρ (χ) e) Y T = Y T (D, R, t ns+1 ). ) hp, ς hp ( ρ (χ) hp, ϕ hp, L (χ 1) = 0, Ω ) G Ω G.
57 The spatial discretization The semidiscrete formulation The discrete formulations Entropy enriched hp-adaptivity and stability The point of the RKC scheme is to extend the real stability region of the Jacobians in the operator splitting while maintaining second order consistency (others exist, such as the ROCK and DuUMKA schemes, fourth order consistent),.
58 IMEX discrete form The spatial discretization The semidiscrete formulation The discrete formulations Entropy enriched hp-adaptivity and stability The discrete split IMEX DG scheme a) ρ hp (0) = Π hp ρ 0, b) ρ n f hp = ρ hp(0), ( ( ) cd) ρ ns+1 hp, ς hp = ρ )Ω (χ) hp, ς hp, L (χ 1) = 0, G Ω ( ( ) G dc) ρ ns+1 hp, ϕ hp = ρ )Ω ns hp, ϕ hp + t ns Z ( ) A hp ( n), ϕ hp, G Ω G e) Y T = Y T (D, R, t ns+1 ).
59 Implicit schemes for the mass action The spatial discretization The semidiscrete formulation The discrete formulations Entropy enriched hp-adaptivity and stability Here the implicit timestepping in (dc) is chosen such that we implement the usual back differentiation formulas (BDF(k)) of order k. Here we use Newton Krylov methods with low accuracy tolerances such that the explicit diffusion step stability is taken as the stability limiting step. By default the Krylov method used is GMRES, while the Newton iteration is based the standard Jacobian line search methods.
60 Implicit schemes for the mass action The spatial discretization The semidiscrete formulation The discrete formulations Entropy enriched hp-adaptivity and stability Here the implicit timestepping in (dc) is chosen such that we implement the usual back differentiation formulas (BDF(k)) of order k. Here we use Newton Krylov methods with low accuracy tolerances such that the explicit diffusion step stability is taken as the stability limiting step. By default the Krylov method used is GMRES, while the Newton iteration is based the standard Jacobian line search methods.
61 Implicit schemes for the mass action The spatial discretization The semidiscrete formulation The discrete formulations Entropy enriched hp-adaptivity and stability Here the implicit timestepping in (dc) is chosen such that we implement the usual back differentiation formulas (BDF(k)) of order k. Here we use Newton Krylov methods with low accuracy tolerances such that the explicit diffusion step stability is taken as the stability limiting step. By default the Krylov method used is GMRES, while the Newton iteration is based the standard Jacobian line search methods.
62 Outline 1 2 Equation (or normal) ordering Conservation laws 3 Fractional multistep operator splitting Law of mass action Mass diffusivty 4 The spatial discretization The semidiscrete formulation The discrete formulations Entropy enriched hp-adaptivity and stability 5 The spatial discretization The semidiscrete formulation The discrete formulations Entropy enriched hp-adaptivity and stability
63 L 1 -stability (TVB) in quiescent reactors The spatial discretization The semidiscrete formulation The discrete formulations Entropy enriched hp-adaptivity and stability Ends up, that the stoichiometric coefficients determine a weighting prefactor, given each reactions equilibrium constant K eq,r = k fr k 1 br, κ i = [ r R ] 1/n( r R (νb ir νf ir )) K eq,r so that the entropy satisfies: { n S R = n i ( ln κ i n i + κ i )dx + + n i=1 sup 0 t T t 0 Ω i=1 Ω n 1 i D i x n i 2 dxds } n i=1 P 0 r R + t 0 Ω m 1 i D( n)dxds n 2(e κ i ) 1 Ω + C(T ). i=1 (3)
64 p-enrichment via dioristic entropy The spatial discretization The semidiscrete formulation The discrete formulations Entropy enriched hp-adaptivity and stability Now, we use the discrete local-in-time change in entropy density, ( ) ˆρS k+1 R,Ω ei = ˆρ S k+1 R,Ω ei SR,Ω k ei to develop the p-enrichment/de-enrichment strategy, where E k+1 i = ( P s+1 (Ω k+1 e i ) if ˆρS k+1 R,Ω ei ϱs k+1 ) R,Ω s+1 ei µ s+1 (s + 1 p max) (τ 0 t w ), ( P s 1 (Ω k+1 e i ) if ˆρS k+1 R,Ω ei ϱs k+1 ) R,Ω ei s 1 µ s 1 (s 1 p min ) (τ 0 t w ), given the change in the average global entropy density ( ϱs k+1 R = ϱ S k+1 R ) S R k. The notation { } is simply a pair from the indefinite binary relations { } {<, >,, } where the choice determines the difference between a Type I and and Type II enrichment scheme.
65 The basic feature: a stabilizing center The spatial discretization The semidiscrete formulation The discrete formulations Entropy enriched hp-adaptivity and stability Denoting the global average change in entropy density by Avg ΩG ˆρS k+1 R,Ω ei, the stabilizing center at timestep k + 1 is the discrete subdomain c Ω hp comprised of the union of elements over which the change in entropy density satisfies the condition, } c = { 1 j ne Ω ej : ˆρS k+1 R,Ω ei Avg ΩG ˆρS k+1 R,Ω ei s µ s, i.
66 The entropic jump and hp-adaptivity The spatial discretization The semidiscrete formulation The discrete formulations Entropy enriched hp-adaptivity and stability The entropic jump (or the entropic flux ) across element boundaries, similaryl lead to the density of the variation in the entropic jump leading to h-adaptivity functional: ( ) ˆρ J k+1 R,Ω ei = ˆρ J k+1 R,Ω ei JR,Ω k ei, ( T A k+1 h (Ω k+1 e i ) if ˆρ J k+1 R,Ω i = ei ϱ J k+1 ) R h η h (s + 1 h max) (τ 0 t w ) ( T h0 (Ω k+1 e i ) if ˆρ J k+1 R,Ω ei ϱ J k+1 ) R h η h0 (s 1 h min ) (τ 0 t w ) where h max and h min correspond to the maximum and minimum refinement levels, the density of the global change in the entropic jump is given such that: ( ϱ J k+1 R = ϱ J k+1 R ) J R k.
67 Kinetic switch hp-adaptivity functional The spatial discretization The semidiscrete formulation The discrete formulations Entropy enriched hp-adaptivity and stability Finally, we couple the h-adaptivity functional A k+1 i to the p-enrichment functional E k+1 i such that h-adaptivity is always preferentially chosen over p-enrichment. We evaluate the simple kinetic switch functional K k+1 i = K k+1 i (A k+1 i, E k+1 i ) determined by evaluating: A k+1 i E k+1 i if T h (Ω k+1 e i ) P s (Ω k+1 e i ) S k+1 K k+1 i = R, R, R, R, A k+1 i if T h (Ω k+1 e i ) P s+1 (Ω k+1 e i ) S k+1 A k+1 i E k+1 i if T h0 (Ω k+1 e i ) P s (Ω k+1 e i ) S k+1 A k+1 i E k+1 i if T h0 (Ω k+1 e i ) P s+1 (Ω k+1 e i ) S k+1 0 otherwise.
68 Lotka-Volterra reactor with n = 3 species Any system satisfying: t ρ i x (D i x ρ i ) L i ( n) = 0 L i ( n) = λ i n i + n i n j=1 A ij n k=1 nb jk k } GLV system, where λ i is a real or complex parameter and B and A are real or complex valued square matrices, a generalized Lotka-Volterra (GLV) system.
69 Then we consider an exact case when n = 3 corresponds to three chemical species, j R r ν f jr M j k fr k br k P r ν b kr M k, r R, that splits over the mass action in its components to satisfy, dρ 1 dt = N 1 n 1 + U 1 n O 12n 1 n 2 + O 13 n 1 n 3, dρ 3 dt r R U i = m i dρ 2 dt = N 3 n 3 + O 31 n 3 n 1 + O 32 n 3 n 2 + U 3 n 2 3, where the constants are given by, ) N i = m i (νir b νf ir (1 ) {ν f ir >0} k fr 1 {ν b ir >0} k br, r R (ν b ir νf ir ) (1 {ν f ir >1} k fr 1 {ν b ir >1} k br = N 2 n 2 + O 21 n 2 n 1 + U 2 n O 23n 2 n 3, ), r R ) O ij = m i (νir b νf ir (1 ) {ν f ir,ν jr f >0}k fr 1 {ν b ir,ν jr b >0}k br,
70 Reactor with ten reactions (r = 10) and three species (n = 3) Then for the N i s being any excess/bulk/bath constituents, a reactor system comprised of ten coupled reactions that satisfies this equation is: M 1 + M 2 k f1 2M 2, M 1 + M 3 k f2 3M 3, 2M 1 + N 1 k f3 3M 1, M 2 + M 3 + N 2 k f4 2M 2, 2M 2 + N 3 k f5 M 2 + N 4, M 1 + M 3 k f6 3M 1, 2M 3 + N 5 k f7 3M 3, ν f 18 M 1 k f8 N 4, ν f 29 M 2 k f9 N 5, ν f 310 M 3 k f10 N 6. The integrals of motion then yield the following exact time-dependent solutions: n 1 = I 1 N 1 N 2 N 3 e N 1 m 1 1 t I 3 I 1 (I 1 U 1 N 2 N 3 e N 1 m 1 1 t + I 2 U 2 N 1 N 3 e N 2 m 1 2 t + I 1 I 2 U 3 N 1 N 2 e N 3 m 1 3 t ) I 2 N 1 N 2 N 3 e N 2 m 1 2 t n 2 =, I 3 I 1 (I 1 U 1 N 2 N 3 e N 1 m 1 1 t + I 2 U 2 N 1 N 3 e N 2 m 1 2 t + I 1 I 2 U 3 N 1 N 2 e N 3 m 1 3 t ) I 1 I 2 N 1 N 2 N 3 e N 3 m 1 3 t n 3 =, I 3 I 1 (I 1 U 1 N 2 N 3 e N 1 m 1 1 t + I 2 U 2 N 1 N 3 e N 2 m 1 2 t + I 1 I 2 U 3 N 1 N 2 e N 3 m 1 3 t ).
71 When D 0, we have an analytic solution Figure: The solution ρ 1 in N=2.
72 When D 0, we have an analytic solution Figure: We plot the p-convergence of ρ1, where in N = 2 we set h/10 = 1/32 and for N = 3 we have h/10 = 1/16.
73 When D 0, we have an analytic solution, cont. Figure: The h-convergence of ρ1 with N = 2, where the h-levels are defined as the h/10 values here and below.
74 When D 0, we have an analytic solution, cont. Figure: Here we show the h-convergence of the equilibrium solution, where p = 1 and p = 2 in N = 3. Again the h-levels are the h/10 values.
75 Adding back nonvanishing mass diffusion, D 0 With vanishing mass diffusion, only p-enrichment occurs. The h-adaptivity occurs along the profile of the change in entropy density over element interiors. The p-enrichment solution satisfies the minimal convergence requirements. With nonzero mass diffusion, h-adaptivity may occur. The h-adaptivity occurs along the profile of the entropic jump over interior faces. The solution satisfies the minimal convergence requirement, with p-enrichment off.
76 Adding back nonvanishing mass diffusion, D 0 With vanishing mass diffusion, only p-enrichment occurs. The h-adaptivity occurs along the profile of the change in entropy density over element interiors. The p-enrichment solution satisfies the minimal convergence requirements. With nonzero mass diffusion, h-adaptivity may occur. The h-adaptivity occurs along the profile of the entropic jump over interior faces. The solution satisfies the minimal convergence requirement, with p-enrichment off.
77 Adding back nonvanishing mass diffusion, D 0 With vanishing mass diffusion, only p-enrichment occurs. The h-adaptivity occurs along the profile of the change in entropy density over element interiors. The p-enrichment solution satisfies the minimal convergence requirements. With nonzero mass diffusion, h-adaptivity may occur. The h-adaptivity occurs along the profile of the entropic jump over interior faces. The solution satisfies the minimal convergence requirement, with p-enrichment off.
78 Adding back nonvanishing mass diffusion, D 0 With vanishing mass diffusion, only p-enrichment occurs. The h-adaptivity occurs along the profile of the change in entropy density over element interiors. The p-enrichment solution satisfies the minimal convergence requirements. With nonzero mass diffusion, h-adaptivity may occur. The h-adaptivity occurs along the profile of the entropic jump over interior faces. The solution satisfies the minimal convergence requirement, with p-enrichment off.
79 Adding back nonvanishing mass diffusion, D 0 With vanishing mass diffusion, only p-enrichment occurs. The h-adaptivity occurs along the profile of the change in entropy density over element interiors. The p-enrichment solution satisfies the minimal convergence requirements. With nonzero mass diffusion, h-adaptivity may occur. The h-adaptivity occurs along the profile of the entropic jump over interior faces. The solution satisfies the minimal convergence requirement, with p-enrichment off.
80 e.g. h-refinement minimal convergence requirement IMEX L 2 -error RKC, ɛ = 2/13, L 2 -error h p = 1 p = 2 p = 3 p = 1 p = 2 p = 3 1/ e e-06 1/ e e e e-06 1/ e e-05 hp ι s = ι h = 0.5, 2.7e-06 ι s = ι h = 0.5, 3.3e-06 Table: We give the L 2 -errors for the h and p ramping of the equilibrium solution with homogeneous diffusion added. Here we have set t = 0.01 seconds, and T = 0.5 s. The p-enrichment scheme is of Type II and the h-adaptive scheme is of Type I, and both span all h and p levels (i.e. p {1, 2, 3} h {1/16, 1/32, 1/64}).
81 k+1 The gradient σ h versus the entropic jump ρ JR,Ω h Figure: Here we show the derived variables evaluated at timestep five using t = 0.01 with hp-adaptation turned off, and h = 1/256.
82 hp-adapted, h-refinement profile Figure: Here we show an example of the h-refinement scheme that emerges after eleven timesteps of t = 0.01, where we have a Type I refinement scheme for h and a Type II enrichment scheme for p.
83 Asante sana
Dynamic p-enrichment schemes for multicomponent reactive flows
Dynamic p-enrichment schemes for multicomponent reactive flows C. Michoski X, C. Mirabito, C. Dawson, Institute for Computational Engineering and Sciences (ICES), Computational Hydraulics Group (CHG) University
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 informationAdaptive C1 Macroelements for Fourth Order and Divergence-Free Problems
Adaptive C1 Macroelements for Fourth Order and Divergence-Free Problems Roy Stogner Computational Fluid Dynamics Lab Institute for Computational Engineering and Sciences University of Texas at Austin March
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 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 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 informationDynamic p-enrichment schemes for multicomponent reactive flows
Dynamic p-enrichment schemes for multicomponent reactive flows C. Michoski a,, C. Mirabito a, C. Dawson a, D. Wirasaet b, E. J. Kubatko c, J. J. Westerink b a Institute for Computational Engineering and
More informationHYPERSONIC AERO-THERMO-DYNAMIC HEATING PREDICTION WITH HIGH-ORDER DISCONTINOUS GALERKIN SPECTRAL ELEMENT METHODS
1 / 36 HYPERSONIC AERO-THERMO-DYNAMIC HEATING PREDICTION WITH HIGH-ORDER DISCONTINOUS GALERKIN SPECTRAL ELEMENT METHODS Jesús Garicano Mena, E. Valero Sánchez, G. Rubio Calzado, E. Ferrer Vaccarezza Universidad
More informationIterative methods for positive definite linear systems with a complex shift
Iterative methods for positive definite linear systems with a complex shift William McLean, University of New South Wales Vidar Thomée, Chalmers University November 4, 2011 Outline 1. Numerical solution
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 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 informationMulti species Lattice Boltzmann Models and Applications to Sustainable Energy Systems
Multi species Lattice Boltzmann Models and Applications to Sustainable Energy Systems Pietro Asinari, PhD Dipartimento di Energetica (DENER), Politecnico di Torino (POLITO), Torino 10129, Italy e-mail:
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 informationBAE 820 Physical Principles of Environmental Systems
BAE 820 Physical Principles of Environmental Systems Estimation of diffusion Coefficient Dr. Zifei Liu Diffusion mass transfer Diffusion mass transfer refers to mass in transit due to a species concentration
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 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 informationFinite Element Solver for Flux-Source Equations
Finite Element Solver for Flux-Source Equations Weston B. Lowrie A thesis submitted in partial fulfillment of the requirements for the degree of Master of Science in Aeronautics Astronautics University
More informationThe heat equation in time dependent domains with Neumann boundary conditions
The heat equation in time dependent domains with Neumann boundary conditions Chris Burdzy Zhen-Qing Chen John Sylvester Abstract We study the heat equation in domains in R n with insulated fast moving
More informationA high-order discontinuous Galerkin solver for 3D aerodynamic turbulent flows
A high-order discontinuous Galerkin solver for 3D aerodynamic turbulent flows F. Bassi, A. Crivellini, D. A. Di Pietro, S. Rebay Dipartimento di Ingegneria Industriale, Università di Bergamo CERMICS-ENPC
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 informationSUMMER SCHOOL - KINETIC EQUATIONS LECTURE II: CONFINED CLASSICAL TRANSPORT Shanghai, 2011.
SUMMER SCHOOL - KINETIC EQUATIONS LECTURE II: CONFINED CLASSICAL TRANSPORT Shanghai, 2011. C. Ringhofer ringhofer@asu.edu, math.la.asu.edu/ chris OVERVIEW Quantum and classical description of many particle
More informationChapter 1. Governing Equations of GFD. 1.1 Mass continuity
Chapter 1 Governing Equations of GFD The fluid dynamical governing equations consist of an equation for mass continuity, one for the momentum budget, and one or more additional equations to account for
More informationNumerical Approximation of Phase Field Models
Numerical Approximation of Phase Field Models Lecture 2: Allen Cahn and Cahn Hilliard Equations with Smooth Potentials Robert Nürnberg Department of Mathematics Imperial College London TUM Summer School
More informationHydrodynamics, Thermodynamics, and Mathematics
Hydrodynamics, Thermodynamics, and Mathematics Hans Christian Öttinger Department of Mat..., ETH Zürich, Switzerland Thermodynamic admissibility and mathematical well-posedness 1. structure of equations
More informationChapter 2: Fluid Dynamics Review
7 Chapter 2: Fluid Dynamics Review This chapter serves as a short review of basic fluid mechanics. We derive the relevant transport equations (or conservation equations), state Newton s viscosity law leading
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 informationComputational Analysis of an Imploding Gas:
1/ 31 Direct Numerical Simulation of Navier-Stokes Equations Stephen Voelkel University of Notre Dame October 19, 2011 2/ 31 Acknowledges Christopher M. Romick, Ph.D. Student, U. Notre Dame Dr. Joseph
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 informationProblem Set Number 01, MIT (Winter-Spring 2018)
Problem Set Number 01, 18.377 MIT (Winter-Spring 2018) Rodolfo R. Rosales (MIT, Math. Dept., room 2-337, Cambridge, MA 02139) February 28, 2018 Due Thursday, March 8, 2018. Turn it in (by 3PM) at the Math.
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 informationCORBIS: Code Raréfié Bidimensionnel Implicite Stationnaire
CORBIS: Code Raréfié Bidimensionnel Implicite Stationnaire main ingredients: [LM (M3AS 00, JCP 00)] plane flow: D BGK Model conservative and entropic velocity discretization space discretization: finite
More informationBasic Aspects of Discretization
Basic Aspects of Discretization Solution Methods Singularity Methods Panel method and VLM Simple, very powerful, can be used on PC Nonlinear flow effects were excluded Direct numerical Methods (Field Methods)
More informationStability of Feedback Solutions for Infinite Horizon Noncooperative Differential Games
Stability of Feedback Solutions for Infinite Horizon Noncooperative Differential Games Alberto Bressan ) and Khai T. Nguyen ) *) Department of Mathematics, Penn State University **) Department of Mathematics,
More informationVariational formulation of entropy solutions for nonlinear conservation laws
Variational formulation of entropy solutions for nonlinear conservation laws Eitan Tadmor 1 Center for Scientific Computation and Mathematical Modeling CSCAMM) Department of Mathematics, Institute for
More informationLecture 9 Approximations of Laplace s Equation, Finite Element Method. Mathématiques appliquées (MATH0504-1) B. Dewals, C.
Lecture 9 Approximations of Laplace s Equation, Finite Element Method Mathématiques appliquées (MATH54-1) B. Dewals, C. Geuzaine V1.2 23/11/218 1 Learning objectives of this lecture Apply the finite difference
More informationFokker-Planck Equation with Detailed Balance
Appendix E Fokker-Planck Equation with Detailed Balance A stochastic process is simply a function of two variables, one is the time, the other is a stochastic variable X, defined by specifying: a: the
More informationHydrodynamics. Stefan Flörchinger (Heidelberg) Heidelberg, 3 May 2010
Hydrodynamics Stefan Flörchinger (Heidelberg) Heidelberg, 3 May 2010 What is Hydrodynamics? Describes the evolution of physical systems (classical or quantum particles, fluids or fields) close to thermal
More informationFinite Difference Solution of the Heat Equation
Finite Difference Solution of the Heat Equation Adam Powell 22.091 March 13 15, 2002 In example 4.3 (p. 10) of his lecture notes for March 11, Rodolfo Rosales gives the constant-density heat equation as:
More informationEquivalence between kinetic method for fluid-dynamic equation and macroscopic finite-difference scheme
Equivalence between kinetic method for fluid-dynamic equation and macroscopic finite-difference scheme Pietro Asinari (1), Taku Ohwada (2) (1) Department of Energetics, Politecnico di Torino, Torino 10129,
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 informationEntropy and Relative Entropy
Entropy and Relative Entropy Joshua Ballew University of Maryland October 24, 2012 Outline Hyperbolic PDEs Entropy/Entropy Flux Pairs Relative Entropy Weak-Strong Uniqueness Weak-Strong Uniqueness for
More informationarxiv:comp-gas/ v1 28 Apr 1993
Lattice Boltzmann Thermohydrodynamics arxiv:comp-gas/9304006v1 28 Apr 1993 F. J. Alexander, S. Chen and J. D. Sterling Center for Nonlinear Studies and Theoretical Division Los Alamos National Laboratory
More information2. The Schrödinger equation for one-particle problems. 5. Atoms and the periodic table of chemical elements
1 Historical introduction The Schrödinger equation for one-particle problems 3 Mathematical tools for quantum chemistry 4 The postulates of quantum mechanics 5 Atoms and the periodic table of chemical
More informationNumerical methods for the Navier- Stokes equations
Numerical methods for the Navier- Stokes equations Hans Petter Langtangen 1,2 1 Center for Biomedical Computing, Simula Research Laboratory 2 Department of Informatics, University of Oslo Dec 6, 2012 Note:
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 informationA recovery-assisted DG code for the compressible Navier-Stokes equations
A recovery-assisted DG code for the compressible Navier-Stokes equations January 6 th, 217 5 th International Workshop on High-Order CFD Methods Kissimmee, Florida Philip E. Johnson & Eric Johnsen Scientific
More informationAA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS
AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS 1 / 31 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS Linearization and Characteristic Relations 1 / 31 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS
More informationNumerical Methods for Problems with Moving Fronts Orthogonal Collocation on Finite Elements
Electronic Text Provided with the Book Numerical Methods for Problems with Moving Fronts by Bruce A. Finlayson Ravenna Park Publishing, Inc., 635 22nd Ave. N. E., Seattle, WA 985-699 26-524-3375; ravenna@halcyon.com;www.halcyon.com/ravenna
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 informationIntroduction to multiscale modeling and simulation. Explicit methods for ODEs : forward Euler. y n+1 = y n + tf(y n ) dy dt = f(y), y(0) = y 0
Introduction to multiscale modeling and simulation Lecture 5 Numerical methods for ODEs, SDEs and PDEs The need for multiscale methods Two generic frameworks for multiscale computation Explicit methods
More informationCaltech Ph106 Fall 2001
Caltech h106 Fall 2001 ath for physicists: differential forms Disclaimer: this is a first draft, so a few signs might be off. 1 Basic properties Differential forms come up in various parts of theoretical
More informationDELFT UNIVERSITY OF TECHNOLOGY
DELFT UNIVERSITY OF TECHNOLOGY REPORT 05-05 Efficient Methods for Solving Advection Diffusion Reaction Equations S. van Veldhuizen ISSN 1389-6520 Reports of the Department of Applied Mathematical Analysis
More informationClassical Statistical Mechanics: Part 1
Classical Statistical Mechanics: Part 1 January 16, 2013 Classical Mechanics 1-Dimensional system with 1 particle of mass m Newton s equations of motion for position x(t) and momentum p(t): ẋ(t) dx p =
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 informationScale invariant fluid dynamics for the dilute Fermi gas at unitarity
Scale invariant fluid dynamics for the dilute Fermi gas at unitarity Thomas Schaefer North Carolina State University Fluids: Gases, Liquids, Plasmas,... Hydrodynamics: Long-wavelength, low-frequency dynamics
More informationNumerical Algorithms as Dynamical Systems
A Study on Numerical Algorithms as Dynamical Systems Moody Chu North Carolina State University What This Study Is About? To recast many numerical algorithms as special dynamical systems, whence to derive
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 informationIs My CFD Mesh Adequate? A Quantitative Answer
Is My CFD Mesh Adequate? A Quantitative Answer Krzysztof J. Fidkowski Gas Dynamics Research Colloqium Aerospace Engineering Department University of Michigan January 26, 2011 K.J. Fidkowski (UM) GDRC 2011
More informationPiecewise Smooth Solutions to the Burgers-Hilbert Equation
Piecewise Smooth Solutions to the Burgers-Hilbert Equation Alberto Bressan and Tianyou Zhang Department of Mathematics, Penn State University, University Park, Pa 68, USA e-mails: bressan@mathpsuedu, zhang
More informationTowards Multi-field Inflation with a Random Potential
Towards Multi-field Inflation with a Random Potential Jiajun Xu LEPP, Cornell Univeristy Based on H. Tye, JX, Y. Zhang, arxiv:0812.1944 and work in progress 1 Outline Motivation from string theory A scenario
More informationTitle of communication, titles not fitting in one line will break automatically
Title of communication titles not fitting in one line will break automatically First Author Second Author 2 Department University City Country 2 Other Institute City Country Abstract If you want to add
More informationAppendix 1: List of symbols
Appendix 1: List of symbols Symbol Description MKS Units a Acceleration m/s 2 a 0 Bohr radius m A Area m 2 A* Richardson constant m/s A C Collector area m 2 A E Emitter area m 2 b Bimolecular recombination
More informationEnergy Barriers and Rates - Transition State Theory for Physicists
Energy Barriers and Rates - Transition State Theory for Physicists Daniel C. Elton October 12, 2013 Useful relations 1 cal = 4.184 J 1 kcal mole 1 = 0.0434 ev per particle 1 kj mole 1 = 0.0104 ev per particle
More informationCollective behavior, from particles to fields
978-0-51-87341-3 - Statistical Physics of Fields 1 Collective behavior, from particles to fields 1.1 Introduction One of the most successful aspects of physics in the twentieth century was revealing the
More informationRecent 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 informationAdaptive hierarchic transformations for dynamically p-enriched slope-limiting over discontinuous Galerkin systems of generalized equations
arxiv:1012.5826v3 [physics.comp-ph] 11 Aug 2011 Adaptive hierarchic transformations for dynamically p-enriched slope-limiting over discontinuous Galerkin systems of generalized equations C. Michoski X,
More informationSupplement: Statistical Physics
Supplement: Statistical Physics Fitting in a Box. Counting momentum states with momentum q and de Broglie wavelength λ = h q = 2π h q In a discrete volume L 3 there is a discrete set of states that satisfy
More informationLecture 9: RR-sector and D-branes
Lecture 9: RR-sector and D-branes José D. Edelstein University of Santiago de Compostela STRING THEORY Santiago de Compostela, March 6, 2013 José D. Edelstein (USC) Lecture 9: RR-sector and D-branes 6-mar-2013
More informationDiffusive Transport Enhanced by Thermal Velocity Fluctuations
Diffusive Transport Enhanced by Thermal Velocity Fluctuations Aleksandar Donev 1 Courant Institute, New York University & Alejandro L. Garcia, San Jose State University John B. Bell, Lawrence Berkeley
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 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 informationHidden properties of the Navier-Stokes equations. Double solutions. Origination of turbulence.
Theoretical Mathematics & Applications, vol.4, no.3, 2014, 91-108 ISSN: 1792-9687 (print), 1792-9709 (online) Scienpress Ltd, 2014 Hidden properties of the Navier-Stokes equations. Double solutions. Origination
More informationWeak-Strong Uniqueness of the Navier-Stokes-Smoluchowski System
Weak-Strong Uniqueness of the Navier-Stokes-Smoluchowski System Joshua Ballew University of Maryland College Park Applied PDE RIT March 4, 2013 Outline Description of the Model Relative Entropy Weakly
More informationRelaxation methods and finite element schemes for the equations of visco-elastodynamics. Chiara Simeoni
Relaxation methods and finite element schemes for the equations of visco-elastodynamics Chiara Simeoni Department of Information Engineering, Computer Science and Mathematics University of L Aquila (Italy)
More informationAccurate representation of velocity space using truncated Hermite expansions.
Accurate representation of velocity space using truncated Hermite expansions. Joseph Parker Oxford Centre for Collaborative Applied Mathematics Mathematical Institute, University of Oxford Wolfgang Pauli
More informationChapter 2. General concepts. 2.1 The Navier-Stokes equations
Chapter 2 General concepts 2.1 The Navier-Stokes equations The Navier-Stokes equations model the fluid mechanics. This set of differential equations describes the motion of a fluid. In the present work
More informationWaves in a Shock Tube
Waves in a Shock Tube Ivan Christov c February 5, 005 Abstract. This paper discusses linear-wave solutions and simple-wave solutions to the Navier Stokes equations for an inviscid and compressible fluid
More informationAnisotropic fluid dynamics. Thomas Schaefer, North Carolina State University
Anisotropic fluid dynamics Thomas Schaefer, North Carolina State University Outline We wish to extract the properties of nearly perfect (low viscosity) fluids from experiments with trapped gases, colliding
More informationLattice Boltzmann Methods for Fluid Dynamics
Lattice Boltzmann Methods for Fluid Dynamics Steven Orszag Department of Mathematics Yale University In collaboration with Hudong Chen, Isaac Goldhirsch, and Rick Shock Transient flow around a car LBGK
More informationFluid Approximations from the Boltzmann Equation for Domains with Boundary
Fluid Approximations from the Boltzmann Equation for Domains with Boundary C. David Levermore Department of Mathematics and Institute for Physical Science and Technology University of Maryland, College
More informationSpace-time Discontinuous Galerkin Methods for Compressible Flows
Space-time Discontinuous Galerkin Methods for Compressible Flows Jaap van der Vegt Numerical Analysis and Computational Mechanics Group Department of Applied Mathematics University of Twente Joint Work
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 informationTOPOLOGICAL EQUIVALENCE OF LINEAR ORDINARY DIFFERENTIAL EQUATIONS
TOPOLOGICAL EQUIVALENCE OF LINEAR ORDINARY DIFFERENTIAL EQUATIONS ALEX HUMMELS Abstract. This paper proves a theorem that gives conditions for the topological equivalence of linear ordinary differential
More informationDyson series for the PDEs arising in Mathematical Finance I
for the PDEs arising in Mathematical Finance I 1 1 Penn State University Mathematical Finance and Probability Seminar, Rutgers, April 12, 2011 www.math.psu.edu/nistor/ This work was supported in part by
More informationONSAGER S VARIATIONAL PRINCIPLE AND ITS APPLICATIONS. Abstract
ONSAGER S VARIAIONAL PRINCIPLE AND IS APPLICAIONS iezheng Qian Department of Mathematics, Hong Kong University of Science and echnology, Clear Water Bay, Kowloon, Hong Kong (Dated: April 30, 2016 Abstract
More information0.3.4 Burgers Equation and Nonlinear Wave
16 CONTENTS Solution to step (discontinuity) initial condition u(x, 0) = ul if X < 0 u r if X > 0, (80) u(x, t) = u L + (u L u R ) ( 1 1 π X 4νt e Y 2 dy ) (81) 0.3.4 Burgers Equation and Nonlinear Wave
More informationThe Lattice Boltzmann method for hyperbolic systems. Benjamin Graille. October 19, 2016
The Lattice Boltzmann method for hyperbolic systems Benjamin Graille October 19, 2016 Framework The Lattice Boltzmann method 1 Description of the lattice Boltzmann method Link with the kinetic theory Classical
More informationHyperbolic Systems of Conservation Laws. in One Space Dimension. II - Solutions to the Cauchy problem. Alberto Bressan
Hyperbolic Systems of Conservation Laws in One Space Dimension II - Solutions to the Cauchy problem Alberto Bressan Department of Mathematics, Penn State University http://www.math.psu.edu/bressan/ 1 Global
More informationLevel Set Tumor Growth Model
Level Set Tumor Growth Model Andrew Nordquist and Rakesh Ranjan, PhD University of Texas, San Antonio July 29, 2013 Andrew Nordquist and Rakesh Ranjan, PhD (University Level Set of Texas, TumorSan Growth
More informationWell-balanced central finite volume methods for the Ripa system
Well-balanced central finite volume methods for the Ripa system R. Touma a C. Klingenberg b a Lebanese American University, Beirut, Lebanon b Würzburg University, Würzburg, Germany This paper is dedicated
More informationFrom a Mesoscopic to a Macroscopic Description of Fluid-Particle Interaction
From a Mesoscopic to a Macroscopic Description of Fluid-Particle Interaction Carnegie Mellon University Center for Nonlinear Analysis Working Group, October 2016 Outline 1 Physical Framework 2 3 Free Energy
More informationLong-Term Atomistic Simulation of Hydrogen Diffusion in Metals
Long-Term Atomistic Simulation of Hydrogen Diffusion in Metals K.G. Wang 1, M. and M. Ortiz 2 Universidad de Sevilla 1 Virginia Polytechnic Institute and State University 2 California Institute of Technology
More information3.012 PS Issued: Fall 2003 Graded problems due:
3.012 PS 7 3.012 Issued: 11.18.03 Fall 2003 Graded problems due: 11.26.03 Graded problems As usual, try as much as possible to explain in a narrative form your reasoning (as opposed to writing down bare
More informationRadiative transfer equation in spherically symmetric NLTE model stellar atmospheres
Radiative transfer equation in spherically symmetric NLTE model stellar atmospheres Jiří Kubát Astronomický ústav AV ČR Ondřejov Zářivě (magneto)hydrodynamický seminář Ondřejov 20.03.2008 p. Outline 1.
More informationStopping, blooming, and straggling of directed energetic electrons in hydrogenic and arbitrary-z plasmas
Stopping, blooming, and straggling of directed energetic electrons in hydrogenic and arbitrary-z plasmas This model Monte Carlo 1 MeV e 1 MeV e C. K. Li and R. D. Petrasso MIT 47th Annual Meeting of the
More informationA High-Order Discontinuous Galerkin Method for the Unsteady Incompressible Navier-Stokes Equations
A High-Order Discontinuous Galerkin Method for the Unsteady Incompressible Navier-Stokes Equations Khosro Shahbazi 1, Paul F. Fischer 2 and C. Ross Ethier 1 1 University of Toronto and 2 Argonne National
More informationLecture Models for heavy-ion collisions (Part III): transport models. SS2016: Dynamical models for relativistic heavy-ion collisions
Lecture Models for heavy-ion collisions (Part III: transport models SS06: Dynamical models for relativistic heavy-ion collisions Quantum mechanical description of the many-body system Dynamics of heavy-ion
More informationGetting started: CFD notation
PDE of p-th order Getting started: CFD notation f ( u,x, t, u x 1,..., u x n, u, 2 u x 1 x 2,..., p u p ) = 0 scalar unknowns u = u(x, t), x R n, t R, n = 1,2,3 vector unknowns v = v(x, t), v R m, m =
More informationThe mathematical description of the motion of Atoms, Molecules & Other Particles. University of Rome La Sapienza - SAER - Mauro Valorani (2007)
The mathematical description of the motion of Atoms, Molecules Other Particles Particle Dynamics Mixture of gases are made of different entities: atoms, molecules, ions, electrons. In principle, the knowledge
More informationOn limiting for higher order discontinuous Galerkin method for 2D Euler equations
On limiting for higher order discontinuous Galerkin method for 2D Euler equations Juan Pablo Gallego-Valencia, Christian Klingenberg, Praveen Chandrashekar October 6, 205 Abstract We present an implementation
More information