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

1 comment:

Unknown said...

Hi Rhys,

I noticed your comment on my blog -- do you want to go for a beer? I am at LLNL too this summer. You can find my contact info at my webpage:


Just give me a call/email.


Subscribe Subscribe to The Return of Agent Zlerich