|
[Sponsors] |
October 17, 2008, 12:56 |
1D Lagrangian code
|
#1 |
Guest
Posts: n/a
|
Hi,
I am trying to write a 1D Lagrangian hydrocode, my lagrangian co-ordinate is the mass of the cell. My equations are: dv/dt-du/dh=0 du/dt+dp/dh=0 d/dt(e+1/2u^2)+d/dh(pu)=0 Where v is specific volume, p is pressure, u is velocity and e is internal energy/unit mass and I use an EoS to calculate pressure from e and v. There are a few things which I need a but of clarification on. Firstly, do I use the initial mass distribution throughout the code? Will a predictor-corrector method work as a numerical method? When looking at the Sod shock tube problem what artificial viscosity method should I be using? I get large spikes near the shock front and oscillations near the expansion wave. Any advice? |
|
October 17, 2008, 15:08 |
Re: 1D Lagrangian code
|
#2 |
Guest
Posts: n/a
|
Do a google search on "Neumann artificial viscosity". He was one of the first people to create a working Lagrangian hydrocode (if I remember correctly).
|
|
October 17, 2008, 15:15 |
Re: 1D Lagrangian code
|
#3 |
Guest
Posts: n/a
|
I know many flux limiter methods which I can use, I have tried a few of these and they don't seem to work for some reason. Is there a non-shock problem I can try my code on to see if the actual code is working properly?
Mat |
|
October 18, 2008, 09:32 |
Re: 1D Lagrangian code
|
#4 |
Guest
Posts: n/a
|
Just out of curiosity, how are you applying a flux limiter in a Lagrangian scheme where there is no advective flux? Or are you referring to a slope limiter? Either way, take a look at the Neumann artificial viscosity - as I recall it basically uses a pressure switch to control a 2nd-order pseudo-viscous term, and is fairly straight-forward to implement.
|
|
October 18, 2008, 17:24 |
Re: 1D Lagrangian code
|
#5 |
Guest
Posts: n/a
|
I have been reading the book, "difference methods for initial value problems." where the switch is velocity. I am using a flux limiter from a paper on artificial viscosity in more than 1 dimension, it was in the journal of computational physics. I can't remember the authors but it combines a quadratic and linear viscosity.
I can send you my code if you wish, it can be run in matlab/octave with ease. Mat |
|
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
The FOAM Documentation Project - SHUT-DOWN | holger_marschall | OpenFOAM | 242 | March 7, 2013 12:30 |
How to make code run in parallel? | cwang5 | OpenFOAM Programming & Development | 1 | May 30, 2011 04:47 |
Open Source Vs Commercial Software | MechE | OpenFOAM | 28 | May 16, 2011 11:02 |
Design Integration with CFD? | John C. Chien | Main CFD Forum | 19 | May 17, 2001 15:56 |
public CFD Code development | Heinz Wilkening | Main CFD Forum | 38 | March 5, 1999 11:44 |