dolfin team mailing list archive

[Victor <flyaway1212@xxxxxxxxx>]

Dear developers and users of DOLFIN.
I am considering using DOLFIN to solve the equations of
magnetohydrodynamics (MHD).
The equation looks like this:

dU/dt + grad(F)=0

where U is a vector of physical quantities, dU/dt is a partial
derivative and F is flux, which is a non-linear function of U. It is a
non-linear hyperbolic system of equations similar to Euler equations.

The method I am trying to implement is an explicit Runge-Kutta
Discontinuous Galerkin (RKDG) method by Cockburn et al. So the
discretized equation should look like this:

M dU/dt = R(U)

where U is the solution vector, M is a mass matrix, dU/dt is a full time
derivative and R(U) is a function of U.

So my question is as follows: Is it possible to use the DOLFIN ODE
solver to solve this kind of system or am I saying complete nonsense?
Has anybody done any work with similar equations or methods using DOLFIN?


Also, most of the ODE demos in DOLFIN 0.8.1 won't compile on my machine.
Is there interest in bug reports?

Any advise or pointers would be greatly appreciated.

Sincerely,
Victor Prosolin.
Graduate student,
University of Calgary.