Re: Shock Tube

From: Sander Valcke <>
Date: Fri, 28 May 2010 13:38:07 +0200

Hi Pablo,

I would say that you should either go with the boundary particles (you
have to set those up yourself in your IC file by giving those particles
an ID of 0) or with periodic boundary conditions. If the latter then you
actually have 2 shock-tube experiments (one in the center of the tube
and one where the low-density gas is wrapped into the high-density gas)
so you have to make the tube sufficiently long so those shocks have not
met each-other at the latest time you want to have information about the
shock tube. I would vote for the periodic boundary conditions as the
easiest solution, as you only have to accurately specify the
-DLONG_X/Y/Z values as required and make sure that the tube is long enough.

Note that Gadget does not do 1-D simulations out of the box so without
any extra coding you will have to do the shock tube in 2-D (with the
-DTWODIMS flag) or 3-D. Boundary particles will by required here in the
x- and y direction or the x, y and z directions respectively (unless the
shock tube is sufficiently large so the "evaporating" along the edges
has not reached the region of interest by the end of the simulation).



Pablo Araya Melo wrote:
> Hi,
> I'm trying to run a shock tube, but I am not sure of the necessary
> Makefile conditions. So far, I have been using -DNogravity,
> -DLong_x_y_z and -Dperiodic, together with the recommended ones.
> Should I also set -DSPH_BND_PARTICLES, so that the particles at the
> x-borders of the tube act as walls?
> Cheers,
> Pablo
> -----------------------------------------------------------
> If you wish to unsubscribe from this mailing, send mail to
> with a subject of: unsubscribe
> gadget-list
> A web-archive of this mailing list is available here:

| Sander Valcke                                  |
| Sterrenkundig Observatorium, Universiteit Gent |
| Krijgslaan 281 S9, B-9000 Gent, Belgium        |
| phone +32-9-264-4796  |  fax +32-9-264-4989    |
|  |
The moon is made of green cheese.
		-- John Heywood
By sending an email to ANY of my addresses you are agreeing that:
   1. I am by definition, "the intended recipient"
   2. All information in the email is mine to do with as I see fit and
      make such financial profit, political mileage, or good joke as
      it lends itself to.
   3. I may take the contents as representing the views of your company.
   4. This overrides any disclaimer or statement of confidentiality
      that may be included in your message.
Received on 2010-05-28 13:51:21

This archive was generated by hypermail 2.3.0 : 2023-01-10 10:01:31 CET