[FLASH-USERS] freezing laseslab target until temperature threshold

[Ong Jian Fuh]

Hi Tomer,

I also try to freeze the grid and get the same problem as you. The 
target freeze all the way even if the temperature has raise above 1000K. 
It has been very blur how this can be done as suggested in previous threads.

Regards,
Ong





Ong Jian Fuh
Extreme Light Infrastructure - Nuclear Physics
"Horia Hulubei" National Institute for Physics and Nuclear Engineering
30 Reactorului Street, 077125 Magurele (Ilfov), Romania
www.eli-np.ro

On 07/11/2022 10:51 pm, Tomer Catabi wrote:
> Hello Flash users and developers,
>
> I am currently imulating laser induced plasma from a solid target, and 
> am using the laserslab simulation as a template. I noticed that the 
> plasma expands quicker than I would expect, especially in simulation 
> where I compute the power to be above 1.0e6 W. As a result I am trying 
> to freeze my target until it gets to some temperature threshold.
>
> I have defined the BDRY variable and am trying to implement a 
> Simulation_adjustEvolution file into my folder, but the target seems 
> to stay frozen for the entire simulation even though the temperature 
> has surpassed the threshold.
> The weird thing is that if I compute a value that is smaller than the 
> initial temperature the target does not freeze (which makes me think 
> the code itself is ok and I am missing something else).
>
> My Simulation_adjustEvolution loop is:
>
> #ifdef BDRY_VAR
>
>   do lb = 1, blkcnt
>      call Grid_getBlkIndexLimits(blklst(lb),blkLimits,blkLimitsGC)
>      if (NDIM > 1) call Grid_getBlkBC(blklst(lb),bcs)
>      call Grid_getBlkPtr(blklst(lb), blkPtr)
>
>      do k = blkLimits(LOW,KAXIS), blkLimits(HIGH,KAXIS)
>         do j = blkLimits(LOW,JAXIS), blkLimits(HIGH,JAXIS)
>            do i = blkLimits(LOW,IAXIS), blkLimits(HIGH,IAXIS)
>
>               !  hardwired temperature threshold of 1000 K
>                if (blkPtr(TEMP_VAR,i,j,k) > 1.e3) then
>                   blkPtr(BDRY_VAR,i,j,k) = -1.0
>                end if
>
>            enddo
>         end do
>      end do
>
>      call Grid_releaseBlkPtr(blklst(lb), blkPtr)
>
>   end do
>
> #endif
>
> My pulse length is 43 ps (I am using 4 sections to describe it) and I 
> am running a 1d simulation.
>
> I would really appreciate any input on the matter.
>
> Kind regards,
> Tomer
>
> _______________________________________________
> flash-users mailing list
> flash-users at flash.rochester.edu
>
> For list info, including unsubscribe:
> https://flash.rochester.edu/mailman/listinfo/flash-users