## Computational Fluid Dynamics and Numerical Analysis

### Lagrangian and Semi-Lagrangian for pure advection

A couple of students have been trying to grasp the beauty of Lagrangian and Semi-Lagrangian scheme for pure advection problems which arises in image matching. We have created an online course on the numerical solution of the 1D linear advection equation.

### EPDiff: Euler-Poincaré Equation in the full group of Diffeomorphisms

Shape analysis is predicated on comparison of differences relative to template coordinates in the Grenander deformable template model. Specifically shape is encoded via momentum vactors in template coordinates. These momenta are analogous to solitons, wave fronts and higher dimensional forms found in classical mechanics. Therefore one can apply methods established in computational mechanics.

It is important to develop suitable algorithms for evolution of shapes in Computational Anatomy via EPDiff. Hence the need to develop a numerical method for solving a Hamiltonian that preserves energy as well as linear and angular momentum? But this may not be possible unless adaptive timestep is used according to Ge and Marsden. So can we sacrifice energy for linear and angular momentum or vice-versa.

### Turbulence modelling

Asymptotic analysis of turbulence models can be used as an effective tool in verifying both numerical codes and turbulence models as exemplified by the following:

We have created an online course on a simple model of turbulent flow in a plane channel.

### Forward-backward parabolic equations

Not only do these equations arise in the numerical solution of separating boundary layers but also in stochastic processes in financial modelling, particle transport problems and modelling of counter-current separator. Examples in thermal boundary layer separation can be found here.

### Convection-diffusion equations with source and sink terms

The k-&epsilon turbulence model presents interesting problems in numerical analysis. First is the problem of low values of k and &epsilon in the freestream just outside of a turbulent boundary layer. Asymptotic analysis (Mohammadi & Pirroneau; Deriat) show that a lower bound for &epsilon in these regimes is necessary to guarantee non-negative solution. The uniqueness of numerical solution of k-&epsilon equations using Newton-Raphson solvers was left unanswered in my thesis. This may be due in part to the presence of nonlinear source and sink terms in the coupled convection-diffusion equations.

### Numerical Solution of Infinite Integrals of Products of Bessel Functions

A problem in auditory physiology requires computing infinite integrals involving products of two Bessel functions of the first or second kind. Such integrals occur in fluid dynamics, elasticity, electrodynamics and biophysics to name but a few applications. We have developed and tested a MATLAB toolbox called IIPBF. It was used to compute kernels in a system of Fredholm Integral equation of the second kind that arose in models of a viscous fluid jet impinging on an infinite plane (Davis et al. 2012; Davis et al. 2013).

### Recommended books

This selection of recent and old books is evidently biased:

Associate Research Professor
Center for Imaging Science and
Institute for Computational Medicine,
Department of Biomedical Engineering,
The Johns Hopkins University