SUMMER SCHOOL - KINETIC EQUATIONS LECTURE II: CONFINED CLASSICAL TRANSPORT Shanghai, 2011.
|
|
- Harold Long
- 6 years ago
- Views:
Transcription
1 SUMMER SCHOOL - KINETIC EQUATIONS LECTURE II: CONFINED CLASSICAL TRANSPORT Shanghai, C. Ringhofer ringhofer@asu.edu, math.la.asu.edu/ chris
2 OVERVIEW Quantum and classical description of many particle systems in confined geometries and periodic media. L1 Introduction. Transport pictures (Hamilton, Schrödinger, Heisenberg, Wigner). Thermodynamics and entropy. Two different spatial scales. Induce two different time scales. Macroscopic (simpler) description on the large scale. Retain the microscopic transport information on the small scale. L2 Semiclassical transport in narrow geometries. Free energy models. Treating complex geometries. Effective mass approximations in thin tubes. L3 Quantum transport in narrow geometries. Thin plates, narrow tubes. Sub-band modeling. L4 Homogenization in unstructured media and networks.
3 CLASSICAL TRANSPORT WITH BGK OPERATOR t f + [E, f ] c = Q[f ] f = f (X, P, t): X: position P: momentum [E, f ] c : classical commutator [E, f ] c = P E X f X E P f Q[f ]: BGK collision operator, conserving an observable κ with an equilibrium M Q[f ] = 1 τ (M ρ[f ] f ) ρ: Lagrange multiplier functions κm ρ [f ] dxp = κf (X, P) dxp
4 OUTLINE Model transport complicated geometric structures (irregular pipes) on large time scales. Sub-band modeling for classical transport in ion channels. Physical principle: Strong confinement to a narrow region. Scattering processes partially conserving the energy. Large time asymptotics diffusion equations with free energy. Computation of transport coefficients from approximate confinement potentials. Results.
5 TRANSPORT IN NARROW PIPES Confinement produced by repulsive charges on pipe walls. Typical example: Ion channels (proteins in cell walls). Transport of ions in water. Scattering with background (water). Confinement due to repulsive charges on protein walls. Hamiltonian dynamics + collisions with a background confined geometries.
6 SCHEMATIC ION CHANNEL
7 PROTEIN 1LNQ
8 MODEL HIERARCHIES FOR ION CHANNELS 1 Molecular Dynamics. (compute the background). Large computation (Roux) (O(10 7 )sec). 2 Monte Carlo. (collisions with random background, Brownian motion or Boltzmann) (Ravaioli, Saraniti) O( sec). 3 Time scale for device function: O(10 3 sec) 4 Macroscopic models (Hydrodynamics, Diffusion etc.) (Jerome, Peskin) Loss of geometric information. 1-3: physical transport mechanisms. 4: function of the structure on large time scales.
9 GOAL Develop a model that is computationally feasible on large time scales while still incorporating more microscopic geometric information into the model. Compromise: One more free variable than standard hydrodynamics. Cross directional energy. A squeezed SHE model. Approach: Sub-band modeling. Treat transport in the confinement direction on a microcopic level. Average kinetic equations in the transport direction on long time scales.
10 CONFINED GEOMETRY Transport equation X Ω: narrow region Split: x Ω x, y Ω y ε: aspect ratio. Plates: Ω x R 2, Ω y R Pipes:Ω x R, Ω y R 2 t f + [E, f ] c = Q[f ] [E, f ] c = P E X f X E P f X = (x, y), P = (p, q), Ω = Ω x Ω y Ω y << Ω x, Ω y ε Ω x
11 SCALING DIMENSIONLESS EQUATIONS Ω y given by strong confinement in a narrow region y << x. Forces in confinement direction large compared to forces in transport direction y V >> x V. V 0 (x): gauge potential Scaling: V(x) = V 0 (x) + V 1 (εx, y) y εy, q q ε E = E xp + E yq, E yq Eyq ε E = p 2 + q 2 +V(x, y), 2m ε: aspect ratio. E xp = p 2 2m +V 0(x), E yq = q 2 2m +V 1(x, y),
12 SCALING THE COLLISIONS Scattering with a background (water) does not conserve energy, but only mass. Scattering exchanges a given amount ω of energy per event with background (water). Q[f ] = S f (P, P )δ( P 2 2 P 2 2 ± ω)f (X, P ) dp γ(p)f After re-scaling P = (p, q), q q ε Q[f ] = S f (P, P )δ( p q 2 p 2 2ε 2 2 q 2 ± ω)f (X, P ) dp γ(p)f 2ε 2 Q asymptotically conserves (locally in space) the cross directional kinetic energy q 2 2
13 CONSERVATION IN THE WEAK FORMULATION Q[f ] = S f (P, P )δ( P 2 2 P 2 2 ± ω)f (X, P ) dp γ(p)f γ(p) = S f (P, P)δ( P 2 2 P 2 ± ω) dp 2 uq[f ] dpq = [u(p, q) u(p, q )]S f (P, P )δ( q2 2ε 2 q 2 2ε 2 )f (p, q ) dpqp q if u depends only on q 2 : u(p, q 2 )Q[f ] dpq = 0 A collision with the background in the confinement direction y is a rare event compared to collisions with the background in the transport direction x.
14 SCALED (APPROXIMATE) MODEL EQUATIONS t f + {E xp, f } xp = 1 ε {E yq, f } yq + 1 ε Q[f ] = 1 ε C[f ] E xp (x, p) = p V 0(x), E yq (x, y, q) = q V 1(x, y) {E xp, f } xp = x (fp) p (f x V 0 ), {E yq, f } = y (fq) q (f y V 1 ) Thermodynamics (entropy maximizer) in p: Q[f ] = M(p) δ( p 2 2 p 2 2 )f (x, y, p, q ) dp q γ(p)f Q conserves q 2 2, and is local in (x, y) The commutator {E yq, } yq conserves E yq C conserves E yq Large Time Dynamics: (Hilbert or Chapman - Enskog) Diffusion equation with η = E yq as an additional independent variable for the density n(x, η, t).
15 LARGE TIME DYNAMICS- THE CHAPMAN - ENSKOG EXPANSION General setting for the linear case: t f + Lf = 1 ε Cf C has a set of conserved quantities κ: κ, Cf = 0, f slow dynamics for the macroscopic variable κ, f t κ, f + κ, Lf = 0 At the same time, f driven towards the kernel of C. Parameterize the kernel manifold of C by coordinates given κ, f. (κ and the kernel of C have to have the same dimension.)
16 PROJECTIONS Projection on the kernel M(ρ), parameterized by ρ. Pf = M(ρ), κ, Pf = κ, f CPf = CM = 0, f Split solution into κ, PCf = κ, M = κ, Cf = 0 PCf = 0, f f = f 0 + εf 1, f 0 = Pf, εf 1 = (id P)f t f 0 + PL(f 0 + εf 1 ) = 0, ε t f 1 + (id P)L(f 0 + εf 1 ) = Cf 1 Large time asymptotics ε 0 in the kinetic eaution t κ, M(ρ) + κ, L(M(ρ) + εf 1 ) = 0, (id P)LM(ρ) = Cf 1
17 SUMMARY Solve: (id P)LM(ρ) = Cf 1, κ, f 1 = 0 f 1 = C + M Compute the pseudo - inverse C + of C. Solve the macroscopic equation t κ, M(ρ) + κ, L(M(ρ) + εc + (id P)LM) = 0, for ρ(t). κ, LM : convection term κ, LC + (id P)LM : diffusion term
18 t κ, M(ρ) + κ, L(M(ρ) + εc + (id P)LM) = 0, Two cases: κ, LM 0: Diffusive a small perturbation of convection (Navier - Stokes regime). κ, LM = 0: ( Diffusion) t κ, M(ρ) + ε κ, LC + LM = 0, Diffusion equation on t ε time scale.
19 CONFINED BOLTZMANN t f + L[f ] = 1 ε C[f ] L[f ](x, y, p, q) = [E xp, f ] xp = p E xp x f x E xp p f C[f ](x, y, p, q) = [E yq, f ] yq +M(p) δ( q q 2 )f (x, y, p, q ) dp q γf
20 Conserved quantities: κ = φ(x, E yq (x, y, q)) φ(x, E yq (x, y, q))c[f ] dypq = 0, f Kernel: M[ρ] = ρ(x, E yq (x, y, q))e p 2 2 = n(x, E yq(x, y, q)) p 2 N(x, E yq (x, y, q)) e 2 N: density of states function: N(x, η) = δ(e yq (x, y, q) η) dyq ( κm = n) Projection: P[f ] = ρ(x, E yq )M(p), φ(x, E yq )P[f ] dxypq = φ(x, E yq )f dxypq
21 Diffusion equation on O( t ε 2 ) time scale t n + x F x + η F η = 0, F x = F x ( x n, η n), F η = F η ( x n, η n) η E yq : energy in the confinement direction. F x, F η : Fluxes in space and energy, given by x F x + η F η = δ(e yq η)lc + (id P)L n D dypq
22 PRACTICAL PROBLEM x F x + η F η = δ(e yq η)lc + (id P)L n D dypq The computation of the fluxes F x, F η requires the inversion of the augmented collision operator C = Q[ ] {E yq, } yq. No exact solution. Has to be done numerically. Solve a 2 dim(y) 1 dimensional problem for every grid point in (x, η) to compute the transport coefficients. Plates: dim(y) = 1, solve a 1-D problem at any gridpoint to compute transport coefficients for a 3-D diffusion equation. Pipes: dim(y) = 2, solve a 3-D problem at any gridpoint to compute transport coefficients for a 2-D diffusion equation. Reduce computational complexity by approximating V 1 (x, y) by a harmonic potential in y.
23 COMPUTATION OF THE PSEUDO INVERSE C + [E yq, f 1 ] + M(p) δ(e yq E yq)f 1 dypq γf 1 = L[e p 2 2 n D (x, E yq] equation for the averages of f 1 over equipotential surfaces E yq = q V 1(x, y) = const. Harmonic potential approximation: V 1 (x, y) quadratic in confinement variable y with x dependent coefficients. V 1 (x, y) = 1 2 (y b(x))t G(x)(y b(x)) Equipotential surfaces become ellipsoids in R 4. f 1 written in coordinates of the energy E yq and a three dimensional angle in R 4. Invert C exactly in two of the three angular dimensions. Solve a one dimensional problem in the azimuthal angle numerically.
24 Choice of b and G For every gridpoint in the transport direction x, solve an L 2 minimization problem for the forces in the confinement direction y. Ω y y V 1 G(x)(y b(x)) 2 dy min, x Flux computation (the inversion of C) can be carried out exactly in 2 of the 3 dimensions. Reduced to an effective 1-D problem. Variable transformation: α, β [ π, p], θ [0, π] (y, q) (E(y, q), θ, α, β) C diagonal in α, β. Legendre polynomials in cos(θ).
25 y y x y y x APPROXIMATION QUALITY A toy channel Generate random charges. Compute the exact Coulomb Potential V(x, y) corresponding to these charges, and the local quadratic approximation.
26 TRAJECTORIES Molecular dynamics trajectories (x = p, p = x V) for the exact and approximate Coulomb potential and random initial conditions y y x 4
27 THE LARGE TIME DIFFUSION SYSTEM t n + x F x + η F η = 0 F x = ND x x n N Nµ x(1 + η ) n N, F η = ND η (1 + η ) n N Nµ η x n N, N, D x, D η, µ x, µ η Functionals of G(x), b(x) and of V, (computed numerically). Yields a parabolic system for n.
28 ENTROPY AND PARABOLICITY Theorem:(C. Heitzinger, CR, CMS 11) ( ) Dx µ (implies) x 0 D η µ η t e V 0n(x,η) 2 N(x,η) dxη 0
29 Protein 1LNQ, Charges for actual potential
30 NUMERICAL RESULTS (DENSITY n)
31 NUMERICAL RESULTS (FLUXES F x, F η )
32 CONCLUSION Incorporate arbitrary geometries in semiclassical transport. Evolution on time scales comparable to functional time scales of the channel. Q: Verification? Q: Quantum transport?
SUMMER SCHOOL - KINETIC EQUATIONS LECTURE I: TRANSPORT PICTURES AND DISSIPATION Shanghai, 2011.
SUMMER SCHOOL - KINETIC EQUATIONS LECTURE I: TRANSPORT PICTURES AND DISSIPATION Shanghai, 2011. C. Ringhofer ringhofer@asu.edu, math.la.asu.edu/ chris OVERVIEW Quantum and classical description of many
More informationSUMMER SCHOOL - KINETIC EQUATIONS LECTURE IV: HOMOGENIZATION OF KINETIC FLOWS IN NETWORKS. Shanghai, 2011.
SUMMER SCHOOL - KINETIC EQUATIONS LECTURE IV: HOMOGENIZATION OF KINETIC FLOWS IN NETWORKS. Shanghai, 2011. C. Ringhofer ringhofer@asu.edu, math.la.asu.edu/ chris OVERVIEW Quantum and classical description
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 informationTransport equations for confined structures derived from the Boltzmann equation
Transport equations for confined structures derived from the Boltzmann equation Clemens Heitzinger School of Mathematical and Statistical Sciences, Arizona State University, and Institute for Analysis
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 informationModels of collective displacements: from microscopic to macroscopic description
Models of collective displacements: from microscopic to macroscopic description Sébastien Motsch CSCAMM, University of Maryland joint work with : P. Degond, L. Navoret (IMT, Toulouse) SIAM Analysis of
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 informationCONTINUUM MODELS FOR FLOWS ON COMPLEX NETWORKS
CONTINUUM MODELS FOR FLOWS ON COMPLEX NETWORKS C. Ringhofer ringhofer@asu.edu, math.la.asu.edu/ chris OUTLINE 1 Introduction: Network Examples: Some examples of real world networks + theoretical structure.
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 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 limit and numerical solutions in a system of self-pr
Hydrodynamic limit and numerical solutions in a system of self-propelled particles Institut de Mathématiques de Toulouse Joint work with: P.Degond, G.Dimarco, N.Wang 1 2 The microscopic model The macroscopic
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 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 informationMathematical modelling of collective behavior
Mathematical modelling of collective behavior Young-Pil Choi Fakultät für Mathematik Technische Universität München This talk is based on joint works with José A. Carrillo, Maxime Hauray, and Samir Salem
More information(Super) Fluid Dynamics. Thomas Schaefer, North Carolina State University
(Super) Fluid Dynamics Thomas Schaefer, North Carolina State University Hydrodynamics Hydrodynamics (undergraduate version): Newton s law for continuous, deformable media. Fluids: Gases, liquids, plasmas,...
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 informationEntropy and irreversibility in gas dynamics. Joint work with T. Bodineau, I. Gallagher and S. Simonella
Entropy and irreversibility in gas dynamics Joint work with T. Bodineau, I. Gallagher and S. Simonella Kinetic description for a gas of hard spheres Hard sphere dynamics The system evolves under the combined
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 informationLecture No 1 Introduction to Diffusion equations The heat equat
Lecture No 1 Introduction to Diffusion equations The heat equation Columbia University IAS summer program June, 2009 Outline of the lectures We will discuss some basic models of diffusion equations and
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 informationBuoyancy-induced Flow:
Buoyancy-induced Flow: Natural Convection in a Unconfined Space If we examine the flow induced by heat transfer from a single vertical flat plat, we observe that the flow resembles that of a boundary layer.
More informationComparison of Numerical Solutions for the Boltzmann Equation and Different Moment Models
Comparison of Numerical Solutions for the Boltzmann Equation and Different Moment Models Julian Koellermeier, Manuel Torrilhon October 12th, 2015 Chinese Academy of Sciences, Beijing Julian Koellermeier,
More informationHydrodynamic Limit with Geometric Correction in Kinetic Equations
Hydrodynamic Limit with Geometric Correction in Kinetic Equations Lei Wu and Yan Guo KI-Net Workshop, CSCAMM University of Maryland, College Park 2015-11-10 1 Simple Model - Neutron Transport Equation
More informationLattice Boltzmann Method for Fluid Simulations
Lattice Boltzmann Method for Fluid Simulations Yuanxun Bill Bao & Justin Meskas April 14, 2011 1 Introduction In the last two decades, the Lattice Boltzmann method (LBM) has emerged as a promising tool
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 informationComputational Methods of Physical Problems with Mixed Scales
Computational Methods of Physical Problems with Mixed Scales Shi Jin University of Wisconsin-Madison IPAM Tutorial, March 10-12, 2008 Scales and Physical Laws Figure from E & Engquist, AMS Notice Connections
More informationNumber of pages in the question paper : 05 Number of questions in the question paper : 48 Modeling Transport Phenomena of Micro-particles Note: Follow the notations used in the lectures. Symbols have their
More informationHomework Hint. Last Time
Homework Hint Problem 3.3 Geometric series: ωs τ ħ e s= 0 =? a n ar = For 0< r < 1 n= 0 1 r ωs τ ħ e s= 0 1 = 1 e ħω τ Last Time Boltzmann factor Partition Function Heat Capacity The magic of the partition
More informationPhysical Modeling of Multiphase flow. Boltzmann method
with lattice Boltzmann method Exa Corp., Burlington, MA, USA Feburary, 2011 Scope Scope Re-examine the non-ideal gas model in [Shan & Chen, Phys. Rev. E, (1993)] from the perspective of kinetic theory
More informationElectrical Transport in Nanoscale Systems
Electrical Transport in Nanoscale Systems Description This book provides an in-depth description of transport phenomena relevant to systems of nanoscale dimensions. The different viewpoints and theoretical
More information1 The Quantum Anharmonic Oscillator
1 The Quantum Anharmonic Oscillator Perturbation theory based on Feynman diagrams can be used to calculate observables in Quantum Electrodynamics, like the anomalous magnetic moment of the electron, and
More informationQuantum Fokker-Planck models: kinetic & operator theory approaches
TECHNISCHE UNIVERSITÄT WIEN Quantum Fokker-Planck models: kinetic & operator theory approaches Anton ARNOLD with F. Fagnola (Milan), L. Neumann (Vienna) TU Vienna Institute for Analysis and Scientific
More informationApplicability of the High Field Model: An Analytical Study via Asymptotic Parameters Defining Domain Decomposition
Applicability of the High Field Model: An Analytical Study via Asymptotic Parameters Defining Domain Decomposition Carlo Cercignani,IreneM.Gamba, Joseph W. Jerome, and Chi-Wang Shu Abstract In this paper,
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 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 information9.1 System in contact with a heat reservoir
Chapter 9 Canonical ensemble 9. System in contact with a heat reservoir We consider a small system A characterized by E, V and N in thermal interaction with a heat reservoir A 2 characterized by E 2, V
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 informationOptimal Prediction for Radiative Transfer: A New Perspective on Moment Closure
Optimal Prediction for Radiative Transfer: A New Perspective on Moment Closure Benjamin Seibold MIT Applied Mathematics Mar 02 nd, 2009 Collaborator Martin Frank (TU Kaiserslautern) Partial Support NSF
More information8.333: Statistical Mechanics I Fall 2007 Test 2 Review Problems
8.333: Statistical Mechanics I Fall 007 Test Review Problems The second in-class test will take place on Wednesday 10/4/07 from :30 to 4:00 pm. There will be a recitation with test review on Monday 10//07.
More informationScattering theory I: single channel differential forms
TALENT: theory for exploring nuclear reaction experiments Scattering theory I: single channel differential forms Filomena Nunes Michigan State University 1 equations of motion laboratory Center of mass
More informationQuantum Hydrodynamic models derived from the entropy principle
Quantum Hydrodynamic models derived from the entropy principle Pierre Degond, Florian Méhats, and Christian Ringhofer Abstract. In this work, we give an overview of recently derived quantum hydrodynamic
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 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 informationEngineering. Spring Department of Fluid Mechanics, Budapest University of Technology and Economics. Large-Eddy Simulation in Mechanical
Outline Geurts Book Department of Fluid Mechanics, Budapest University of Technology and Economics Spring 2013 Outline Outline Geurts Book 1 Geurts Book Origin This lecture is strongly based on the book:
More informationTopics for the Qualifying Examination
Topics for the Qualifying Examination Quantum Mechanics I and II 1. Quantum kinematics and dynamics 1.1 Postulates of Quantum Mechanics. 1.2 Configuration space vs. Hilbert space, wave function vs. state
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 informationOne-electron Atom. (in spherical coordinates), where Y lm. are spherical harmonics, we arrive at the following Schrödinger equation:
One-electron Atom The atomic orbitals of hydrogen-like atoms are solutions to the Schrödinger equation in a spherically symmetric potential. In this case, the potential term is the potential given by Coulomb's
More informationFluid-Particles Interaction Models Asymptotics, Theory and Numerics I
Fluid-Particles Interaction Models Asymptotics, Theory and Numerics I J. A. Carrillo collaborators: T. Goudon (Lille), P. Lafitte (Lille) and F. Vecil (UAB) (CPDE 2005), (JCP, 2008), (JSC, 2008) ICREA
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 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 informationAn Effective Potential Approach to Modelling 25nm MOSFET Devices S. Ahmed 1 C. Ringhofer 2 D. Vasileska 1
An Effective Potential Approach to Modelling 25nm MOSFET Devices S. Ahmed 1 C. Ringhofer 2 D. Vasileska 1 1 Department of Electrical Engineering, 2 Department of Mathematics, Arizona State University,
More informationA HYPERBOLIC RELAXATION MODEL FOR PRODUCT FLOW IN COMPLEX PRODUCTION NETWORKS. Ali Unver, Christian Ringhofer and Dieter Armbruster
DISCRETE AND CONTINUOUS Website: www.aimsciences.org DYNAMICAL SYSTEMS Supplement 2009 pp. 790 799 A HYPERBOLIC RELAXATION MODEL FOR PRODUCT FLOW IN COMPLEX PRODUCTION NETWORKS Ali Unver, Christian Ringhofer
More informationME 595M: Monte Carlo Simulation
ME 595M: Monte Carlo Simulation T.S. Fisher Purdue University 1 Monte Carlo Simulation BTE solutions are often difficult to obtain Closure issues force the use of BTE moments Inherent approximations Ballistic
More informationG : Statistical Mechanics
G25.2651: Statistical Mechanics Notes for Lecture 1 Defining statistical mechanics: Statistical Mechanics provies the connection between microscopic motion of individual atoms of matter and macroscopically
More informationFourier Law and Non-Isothermal Boundary in the Boltzmann Theory
in the Boltzmann Theory Joint work with Raffaele Esposito, Yan Guo, Rossana Marra DPMMS, University of Cambridge ICERM November 8, 2011 Steady Boltzmann Equation Steady Boltzmann Equation v x F = Q(F,
More informationDiscrete Projection Methods for Incompressible Fluid Flow Problems and Application to a Fluid-Structure Interaction
Discrete Projection Methods for Incompressible Fluid Flow Problems and Application to a Fluid-Structure Interaction Problem Jörg-M. Sautter Mathematisches Institut, Universität Düsseldorf, Germany, sautter@am.uni-duesseldorf.de
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 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 informationThe Klein-Gordon equation
Lecture 8 The Klein-Gordon equation WS2010/11: Introduction to Nuclear and Particle Physics The bosons in field theory Bosons with spin 0 scalar (or pseudo-scalar) meson fields canonical field quantization
More informationPhys 622 Problems Chapter 5
1 Phys 622 Problems Chapter 5 Problem 1 The correct basis set of perturbation theory Consider the relativistic correction to the electron-nucleus interaction H LS = α L S, also known as the spin-orbit
More informationA quantum heat equation 5th Spring School on Evolution Equations, TU Berlin
A quantum heat equation 5th Spring School on Evolution Equations, TU Berlin Mario Bukal A. Jüngel and D. Matthes ACROSS - Centre for Advanced Cooperative Systems Faculty of Electrical Engineering and Computing
More informationLecture 2: Weak Interactions and BEC
Lecture 2: Weak Interactions and BEC Previous lecture: Ideal gas model gives a fair intuition for occurrence of BEC but is unphysical (infinite compressibility, shape of condensate...) Order parameter
More informationDifferent types of phase transitions for a simple model of alignment of oriented particles
Different types of phase transitions for a simple model of alignment of oriented particles Amic Frouvelle CEREMADE Université Paris Dauphine Joint work with Jian-Guo Liu (Duke University, USA) and Pierre
More informationREVIEW. Hamilton s principle. based on FW-18. Variational statement of mechanics: (for conservative forces) action Equivalent to Newton s laws!
Hamilton s principle Variational statement of mechanics: (for conservative forces) action Equivalent to Newton s laws! based on FW-18 REVIEW the particle takes the path that minimizes the integrated difference
More informationMultiscale hp-adaptive quiescent reactors
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 Outline 1 2 Equation (or
More informationFree energy concept Free energy approach LBM implementation Parameters
BINARY LIQUID MODEL A. Kuzmin J. Derksen Department of Chemical and Materials Engineering University of Alberta Canada August 22,2011 / LBM Workshop OUTLINE 1 FREE ENERGY CONCEPT 2 FREE ENERGY APPROACH
More informationEffective dynamics for the (overdamped) Langevin equation
Effective dynamics for the (overdamped) Langevin equation Frédéric Legoll ENPC and INRIA joint work with T. Lelièvre (ENPC and INRIA) Enumath conference, MS Numerical methods for molecular dynamics EnuMath
More informationFundamentals of Fluid Dynamics: Elementary Viscous Flow
Fundamentals of Fluid Dynamics: Elementary Viscous Flow Introductory Course on Multiphysics Modelling TOMASZ G. ZIELIŃSKI bluebox.ippt.pan.pl/ tzielins/ Institute of Fundamental Technological Research
More 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 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 informationDevice Benchmark Comparisons via Kinetic, Hydrodynamic, and High-Field Models
Device Benchmark Comparisons via Kinetic, Hydrodynamic, and High-Field Models Carlo Cercignani, Irene M. Gamba, Joseph W. Jerome, and Chi-Wang Shu Abstract This paper describes benchmark comparisons for
More informationParallel Plate Heat Exchanger
Parallel Plate Heat Exchanger Parallel Plate Heat Exchangers are use in a number of thermal processing applications. The characteristics are that the fluids flow in the narrow gap, between two parallel
More informationThe Big Picture. Thomas Schaefer. North Carolina State University
The Big Picture Thomas Schaefer North Carolina State University 1 Big Questions What is QCD? What is a Phase of QCD? What is a Plasma? What is a (perfect) Liquid? What is a wqgp/sqgp? 2 What is QCD (Quantum
More informationand in each case give the range of values of x for which the expansion is valid.
α β γ δ ε ζ η θ ι κ λ µ ν ξ ο π ρ σ τ υ ϕ χ ψ ω Mathematics is indeed dangerous in that it absorbs students to such a degree that it dulls their senses to everything else P Kraft Further Maths A (MFPD)
More informationQuantum Physics III (8.06) Spring 2007 FINAL EXAMINATION Monday May 21, 9:00 am You have 3 hours.
Quantum Physics III (8.06) Spring 2007 FINAL EXAMINATION Monday May 21, 9:00 am You have 3 hours. There are 10 problems, totalling 180 points. Do all problems. Answer all problems in the white books provided.
More informationBEC in one dimension
BEC in one dimension Tilmann John 11. Juni 2013 Outline 1 one-dimensional BEC 2 theoretical description Tonks-Girardeau gas Interaction exact solution (Lieb and Liniger) 3 experimental realization 4 conclusion
More informationIntroduction to particle physics Lecture 7
Introduction to particle physics Lecture 7 Frank Krauss IPPP Durham U Durham, Epiphany term 2009 Outline 1 Deep-inelastic scattering and the structure of protons 2 Elastic scattering Scattering on extended
More informationTime-Dependent Statistical Mechanics 1. Introduction
Time-Dependent Statistical Mechanics 1. Introduction c Hans C. Andersen Announcements September 24, 2009 Lecture 1 9/22/09 1 Topics of concern in the course We shall be concerned with the time dependent
More informationHyperbolic Models for Large Supply Chains. Christian Ringhofer (Arizona State University) Hyperbolic Models for Large Supply Chains p.
Hyperbolic Models for Large Supply Chains Christian Ringhofer (Arizona State University) Hyperbolic Models for Large Supply Chains p. /4 Introduction Topic: Overview of conservation law (traffic - like)
More informationLangevin Methods. Burkhard Dünweg Max Planck Institute for Polymer Research Ackermannweg 10 D Mainz Germany
Langevin Methods Burkhard Dünweg Max Planck Institute for Polymer Research Ackermannweg 1 D 55128 Mainz Germany Motivation Original idea: Fast and slow degrees of freedom Example: Brownian motion Replace
More informationLattice Boltzmann Method for Fluid Simulations
1 / 16 Lattice Boltzmann Method for Fluid Simulations Yuanxun Bill Bao & Justin Meskas Simon Fraser University April 7, 2011 2 / 16 Ludwig Boltzmann and His Kinetic Theory of Gases The Boltzmann Transport
More informationG : Statistical Mechanics
G25.2651: Statistical Mechanics Notes for Lecture 15 Consider Hamilton s equations in the form I. CLASSICAL LINEAR RESPONSE THEORY q i = H p i ṗ i = H q i We noted early in the course that an ensemble
More informationIntroduction Statistical Thermodynamics. Monday, January 6, 14
Introduction Statistical Thermodynamics 1 Molecular Simulations Molecular dynamics: solve equations of motion Monte Carlo: importance sampling r 1 r 2 r n MD MC r 1 r 2 2 r n 2 3 3 4 4 Questions How can
More informationNon equilibrium thermodynamics: foundations, scope, and extension to the meso scale. Miguel Rubi
Non equilibrium thermodynamics: foundations, scope, and extension to the meso scale Miguel Rubi References S.R. de Groot and P. Mazur, Non equilibrium Thermodynamics, Dover, New York, 1984 J.M. Vilar and
More informationThe Hydrogen atom. Chapter The Schrödinger Equation. 2.2 Angular momentum
Chapter 2 The Hydrogen atom In the previous chapter we gave a quick overview of the Bohr model, which is only really valid in the semiclassical limit. cf. section 1.7.) We now begin our task in earnest
More information2 Canonical quantization
Phys540.nb 7 Canonical quantization.1. Lagrangian mechanics and canonical quantization Q: How do we quantize a general system?.1.1.lagrangian Lagrangian mechanics is a reformulation of classical mechanics.
More informationand the seeds of quantisation
noncommutative spectral geometry algebra doubling and the seeds of quantisation m.s.,.,, stabile, vitiello, PRD 84 (2011) 045026 mairi sakellariadou king s college london noncommutative spectral geometry
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 informationPhysics of fusion power. Lecture 13 : Diffusion equation / transport
Physics of fusion power Lecture 13 : Diffusion equation / transport Many body problem The plasma has some 10 22 particles. No description is possible that allows for the determination of position and velocity
More informationLesson 6: Diffusion Theory (cf. Transport), Applications
Lesson 6: Diffusion Theory (cf. Transport), Applications Transport Equation Diffusion Theory as Special Case Multi-zone Problems (Passive Media) Self-shielding Effects Diffusion Kernels Typical Values
More informationQuantum Corrections for Monte Carlo Simulation
Quantum Corrections for Monte Carlo Simulation Brian Winstead and Umberto Ravaioli Beckman Institute University of Illinois at Urbana-Champaign Outline Quantum corrections for quantization effects Effective
More informationThermoacoustic Devices
Thermoacoustic Devices Peter in t panhuis Sjoerd Rienstra Han Slot 9th May 2007 Introduction Thermoacoustics All effects in acoustics in which heat conduction and entropy variations play a role. (Rott,
More informationThe 3 dimensional Schrödinger Equation
Chapter 6 The 3 dimensional Schrödinger Equation 6.1 Angular Momentum To study how angular momentum is represented in quantum mechanics we start by reviewing the classical vector of orbital angular momentum
More informationEnergy transport: convection
Outline Introduction: Modern astronomy and the power of quantitative spectroscopy Basic assumptions for classic stellar atmospheres: geometry, hydrostatic equilibrium, conservation of momentum-mass-energy,
More informationv ( ) x (b) Using Navier-Stokes equations, derive an expression
Problem 1: Microscopic Momentum Balance In sweep flow filtration (see Figure below), a pressure drop forces fluid containing particles between two porous plates and there is a transverse (cross) flow that
More informationPhase Equilibria and Molecular Solutions Jan G. Korvink and Evgenii Rudnyi IMTEK Albert Ludwig University Freiburg, Germany
Phase Equilibria and Molecular Solutions Jan G. Korvink and Evgenii Rudnyi IMTEK Albert Ludwig University Freiburg, Germany Preliminaries Learning Goals Phase Equilibria Phase diagrams and classical thermodynamics
More informationCHAPTER 6 Quantum Mechanics II
CHAPTER 6 Quantum Mechanics II 6.1 6.2 6.3 6.4 6.5 6.6 6.7 The Schrödinger Wave Equation Expectation Values Infinite Square-Well Potential Finite Square-Well Potential Three-Dimensional Infinite-Potential
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 informationComputing High Frequency Waves By the Level Set Method
Computing High Frequency Waves By the Level Set Method Hailiang Liu Department of Mathematics Iowa State University Collaborators: Li-Tien Cheng (UCSD), Stanley Osher (UCLA) Shi Jin (UW-Madison), Richard
More informationWhy Should We Be Interested in Hydrodynamics?
Why Should We Be Interested in Hydrodynamics? Li-Shi Luo Department of Mathematics and Statistics Center for Computational Sciences Old Dominion University, Norfolk, Virginia 23529, USA Email: lluo@odu.edu
More information