[FLASH-USERS] freezing laseslab target until temperature threshold

[Tomer Catabi]

 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
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://flash.rochester.edu/pipermail/flash-users/attachments/20221107/6038068a/attachment.htm>