31 May 2010

Summer internship

Following a 2,500 mile motorcycle trip from Austin to the Bay Area with my friend Paul, I'll be starting a summer internship at Lawrence Livermore National Labs tomorrow morning.

Instead of my usual my thesis work, I'll be looking at RANS turbulence models/numerics for Rayleigh-Taylor problems. As preparation, over the past couple of months, I've been spending some time staring at conservation law numerics classics like LeVeque's book as well as lecture notes from Caltech's AE 232 course. On the odd chance it is useful to anyone, I've put my toy code up at http://github.com/RhysU/codefour.

The Fortran 90 code solves inviscid or viscid 1D scalar conservation laws of the form u_t + f(u)_x = 0 or u_t + f(u)_x = \nu u_{xx} using WENO reconstruction and a Lax-Friedrichs flux.

It is based heavily on one provided through Caltech's AE 232 course. Compared to the course's code, the following changes have been made:

  1. Adding 3rd-order WENO reconstruction
  2. Adding 2nd and 4th order viscous capabilities
  3. Adding HDF5-based result output
  4. Adding Doxygen-based documentation
  5. Adding Mathematica-based error analysis against Hopf-Cole-based analytic solutions of the viscid Burger's equation

06 May 2010

Made GSL Committer

Subscribe Subscribe to The Return of Agent Zlerich