T8b: Heat transport:
historic cooling of amphitheaters

Before the tutorials:

Think about:

• What does a cooling cycle (in your fridge or in your air conditioner) generally look like?
• How can heat transport combined with T7b be used to mimic an ancient air conditioning system where the entrances to theaters had conical shapes?

During the tutorials:

To run the experiment, we first need to

• Setup the experiment
  • Use the 400×20×20 setup of T7b,
  • but make the tube much shorter, to have some of the original volume on the left.
  • Split the moving wall in two pieces, leaving the right part at rest with an hole to allow a flow through it.
  • Also, reduce the overall velocity.
  • Hint: Here a sketch of the geometry:
• You can now run the simulation with and without conduction switched on.
• Is there a difference in temperature in the left chamber?
• Have a closer look at the differences. What happens in the run without conduction?

Programming goals for T8:

Improve your ability to set up special geometries in the initial conditions.

Solutions

• Again, extending the setup script from T7b for the first test problem.
  • setup_nozzle.pro in IDL (see setup_test2)
  you should get something like shown below
• Again, you can basically use the script from T7b for looking at the results.
  • show_nozzle.pro in IDL
  you should get something like shown below
• Example using Fortran and gnuplot:
  • ifx -g -traceback -check all -fpe0 -o nozzlesetup nozzlesetup.f90
    (this uses glass.txt from T4)
  • ./nozzlesetup
  • gnuplot grid2.plt (to check that the particle positions are set up correctly)
  • compile Gadget with PERIODIC and NOGRAVITY and SPH_BND_PARTICLES and with LONG_X=30, LONG_Y=1, and LONG_Z=1.
  • run Gadget using nozzle.ic as the initial conditions file and BoxSize = 1
    in the parameter file (and TimeMax = 100, TimeBetSnapshot = 0.1, and ConductionEfficiency = 100)
  • ifx -g -traceback -check all -fpe0 -o readsnap readsnap.f90
  • for file in snap_???; do ./readsnap $file >$file.txt; done
  • gnuplot nozzle.plt
    (this requires ffmpeg, which is installed on the mpusm06 machine)
  • xine -l nozzle.mp4
    • with ConductionEfficiency = 100:
    • with ConductionEfficiency = 0: