[FLASH-USERS] Wiggles appearing in density profile

[Esha Kundu]

Hi Klaus,

Thanks for your suggestions. Here are few clarifications.  
          
  >>It is not clear to me what kind of configuration you have created, the
  >>information give appears contradictory: If you are using the "radflaHD"
  >>code, you should be using the hydro code under
  >>source/physics/Hydro/HydroMain/unsplit_rad/ ,
  >>and not the code under
  >>source/physics/Hydro/HydroMain/split/PPM/ ,
  >>as you message implies.  So more information would be useful, including
  >>your setup command, list of units and list of variables that are included.

I am using the PPM solver because initially when I tried with the 
physics/Hydro/HydroMain/unsplit_rad/Hydro_Unsplit 
I got the following error:

 Driver_abort called. See log file for details.
 Error message is [hy_uhd_unsplitUpdateCastroLike] Negative 3T internal energy, CHECK LOG
 Calling MPI_Abort() for shutdown in   2 seconds!
 
--------------------------------------------------------------------------
MPI_ABORT was invoked on rank 0 in communicator MPI_COMM_WORLD 
with errorcode 1.

NOTE: invoking MPI_ABORT causes Open MPI to kill all MPI processes.
You may or may not see output from other processes, depending on
exactly when Open MPI kills them.
-----------------------------------------------------------------------------

In this case (with Hydro_Unsplit ) the code ran for a while and then this error appeared which aborted the simulation. With the PPM solver the code is running properly. 

I  have modified the source/Simulation/SimulationMain/radflaHD/RadBlastWave codes to meet my requirements.  Therefore the set up command I am giving is
./setup SimulationEsha/radflaHD/mySim1T -1d -auto  +spherical -maxblocks=4000 -nxb=100 species=h1 ManualSpeciesDirectives=True  +parallelio  +splitHydro  +1t

    >>Assuming that your are indeed using
    >>source/physics/Hydro/HydroMain/unsplit_rad/ -
    >>1) Your code should also include the directory (see output from setup
    >> command)
    >>  source/physics/Hydro/HydroMain/unsplit_rad/multiTemp/
    >> and one subdirectory of
    >> source/physics/Eos/EosMain/multiTemp/MatRad3/  ;
    >>make sure that is the case.

I am simulating with only one type of species (hydrogen), therefore using
  physics/Eos/EosMain/Gamma
and not 
source/physics/Hydro/HydroMain/unsplit_rad/multiTemp/ 
&
source/physics/Eos/EosMain/multiTemp/MatRad3/.

In the manual it is mentioned that (Section 30.9) RadFlaH differ from generic 3T simulations, in that there is no physical distinction between the temperatures of electrons and ions; that is, electrons and ions are understood as forming a “matter” component. The “matter” component can have a different temperature from radiation. In my simulation the effect of radiation is not signification. So I am using the 1T and physics/Eos/EosMain/Gamma units.

    >>With that variant of Eos implementation, you will probably
    >>need something like
    >>eosModeInit = "dens_temp_equi" # or "dens_temp_gather" ,
    >>depending how you initialize your temperature(s) (note plural)
     >>in Simulation_initBlock.

 As I am simulating with a single species, I defined 
eosModeInit                    = "dens_temp"
eosMode                        = "dens_ie".

Regarding those wiggles:  I found (according to the suggestions of one my colleagues) that the wiggles are coming because of coarse spatial resolution. With high enough resolution those are gone now.

Regards,
Esha

Department of Astronomy
The Oskar Klein Centre
Stockholm University
10691 Stockholm 
Sweden