[FLASH-USERS] Two sedov explosions
Carlo Graziani
carlo at oddjob.uchicago.edu
Mon Jun 11 17:29:03 EDT 2018
Might it not be better to write a local source term that adds the extra
local energy to create the second explosion near x=0.3 exactly at the
time when you want it to go off?
Carlo
On 06/11/2018 04:16 PM, Ian Padilla wrote:
> Dear all,
>
> I’m working on two identical Sedov explosions centered at x=0.3 and x=0.8, respectively, in a [1,1,1] domain box. I’d like one of them to start later than the other and I’m looking for advise on how to delay the evolution of one of them.
>
> Here what I’ve done so far:
>
> I have set Simulation_initBlock to initialize both explosions. Naively one can “pause” the evolution of BOTH explosions by writing an if-stament in Driver_evolveFlash:
>
> if( (dr_simTime .ge. 0.0 .and. dr_simTime < 0.01) .or. (dr_simTime .ge. 0.05) ) then
> call Timers_start("Hydro”)
> call Hydro(blockCount, blockList,
> & dr_simTime, dr_dt, dr_dtOld, sweepDummy)
> call Timers_stop("Hydro”)
> call Driver_driftUnk(__FILE__,__LINE__,driftUnk_flags)
> end if
>
> Since I’d like to stop the evolution of only ONE explosion (please correct me if I’m mistaken) one solution is to identify (inside Driver_evolveFlash) the blocks with center < 0.5 (using Grid_getBlkCenterCoords(blockId, blockCenter)) and remove them from the blockList array that then is passed to Hydro sequence, so that only blocks belonging to the explosion at x = 0.8 get updated in the hydro unit.
>
> I’m not sure this is the best way of doing it. I'd be very grateful if you had any tips on this.
>
> Thanks in advance,
>
> Ian Padilla
> Lund University
>
--
Carlo Graziani (630) 252-1543 (Voice)
Argonne National Laboratory -------------------------------------
carlo at mcs.anl.gov | The less a statesman amounts to, the
| more he loves the flag. -- Kin Hubbard
More information about the flash-users
mailing list