Other physics applications

Experimental html version of downloadable textbook, see http://www.tacc.utexas.edu/~eijkhout/istc/istc.html
\[ \newcommand\inv{^{-1}}\newcommand\invt{^{-t}} \newcommand\bbP{\mathbb{P}} \newcommand\bbR{\mathbb{R}} \newcommand\defined{ \mathrel{\lower 5pt \hbox{${\equiv\atop\mathrm{\scriptstyle D}}$}}} \newcommand\macro[1]{$\langle$#1$\rangle$} \] 19.1 : Lattice Boltzmann methods
19.2 : Hartree-Fock / Density Functional Theory
Back to Table of Contents

19 Other physics applications

19.1 Lattice Boltzmann methods

crumb trail: > lbm > Lattice Boltzmann methods LBMs

\footnote{This chapter owes much to the presentation in  [ChenDoolen:LBM] .} offer a different way of computing a discretized solution to a PDE , based on dividing a domain in cells. Its basic idea is then to wonder `if there is a particle in this cell, and given the force on this cell, where does he particle move to?'. Rather than reasoning about actual particles, the LBM considers the probability that there is a particle in a cell, and how that probability is updated.

Distribution function $f$ (or particle velocity distribution, just under equation 2?). For each cell $i$ \[ f_i(\mathbf{x}+\mathbf{e}_i\Delta x,t+\Delta t) = f_i(\mathbf{x},t)+\Omega_i(f(\mathbf{x},t)). \] where $\Omega_i$ is the collision operator giving the rate of change of the distribution function $f_i$.

Density: $\rho=\sum_i f_i$, momentum density: $\rho\mathbf{u}=\sum_if_i\mathbf{e}_i$.

Conservation of mass $\sum_i\Omega_i=0$ conservation of momentum $\sum_i\Omega_i\mathbf{e}_i=0$

Basic equation, second order in $\epsilon$: \[ \frac{\partial f_i}{\partial t} + \mathbf{e}_i\cdot \nabla f_i + \epsilon\left( \frac12 \mathbf{e}_e\mathbf{e}_i\colon \nabla\nabla f_i + \mathbf{e}_i\cdot\frac{\partial f_i}{\partial t} + \frac12 \frac{\partial^2f_i}{\partial t^2} \right) = \frac{\Omega_i}{\epsilon} \]

19.2 Hartree-Fock / Density Functional Theory

crumb trail: > lbm > Hartree-Fock / Density Functional Theory

Replace particle-particle interactions with computing a field: repeat until field is self-consistent.

Field is matrix-valued

Repeated diagonalization of the Kohn-sham matrix.

Back to Table of Contents