[FLASH-USERS] [EXT] Re: Sedov running time problem

[Adam Reyes]

Hi Binh,

I think the density inside the rarefied region is going below the value of
the runtime parameter "smlrho", which defaults to 1e-10. You can try
lowering that in your flash.par.
*************************************************************************
Adam Reyes

Code Group Leader, Flash Center for Computational Science
Research Scientist, Dept. of Physics and Astronomy
University of Rochester
River Campus: Bausch and Lomb Hall, 369
500 Wilson Blvd. PO Box 270171, Rochester, NY 14627
Email adam.reyes at rochester.edu
Web https://flash.rochester.edu
 (he / him / his)

*************************************************************************


On Fri, Jul 21, 2023 at 3:53 PM Yên Bình <yenbinhpy308 at gmail.com> wrote:

> Dear Adam Reyes,
>
> I can change the explosion energy to 10^51 erg but when I try to minimize
> the initial radius to below 10^14, the simulation cracks and sometimes
> produces weird results like in the photo.
>
> It does not seem to be related to other runtime parameters, it thinks it
> is because of the grid ? I try to set nxb and nyb to higher value like 200
> but it output this error:
>
> *Abort(1) on node 0 (rank 0 in comm 0): application called
> MPI_Abort(MPI_COMM_WORLD, 1) - process 0 *
>
> How can I set the initial radius to 10^4 like a star core before an
> explosion?
>
> Thank you for your kind support,
>
> Binh
>
> [image: image.png]
>
>
>
> Nguyen Thi Yen Binh
>
>
> Summer Student, Academia Sinica Institute of Astronomy and Astrophysics
>
> Space Science and Satellite Technology, University of Science and
> Technology of Ha Noi
>
> binhnty.bi12-060 at st.usth.edu.vn
>
> +84367533986
>
> Tuy Hoa, Phu Yen, Viet Nam
>
>
> <https://urldefense.com/v3/__https://www.facebook.com/wbinhw__;!!CGUSO5OYRnA7CQ!eHzY5sAvt9JVcQA5cJKZc6ry3jGnKBHVycUBNVDfU0KfAwSYw4QTVhFLbkxC7sqZjy9eKhHdhKa2njTzIDK0EINBdi0$>
>
>
> <https://urldefense.com/v3/__https://www.instagram.com/yenbinne/__;!!CGUSO5OYRnA7CQ!eHzY5sAvt9JVcQA5cJKZc6ry3jGnKBHVycUBNVDfU0KfAwSYw4QTVhFLbkxC7sqZjy9eKhHdhKa2njTzIDK0HytT4xk$>
>
>
> <https://urldefense.com/v3/__https://github.com/pea9__;!!CGUSO5OYRnA7CQ!eHzY5sAvt9JVcQA5cJKZc6ry3jGnKBHVycUBNVDfU0KfAwSYw4QTVhFLbkxC7sqZjy9eKhHdhKa2njTzIDK0Q3WW6Mc$>
>
>
>
>
>
>
>
>
> Vào Th 4, 19 thg 7, 2023 vào lúc 16:43 Adam Reyes <
> adam.reyes at rochester.edu> đã viết:
>
>> Hi Binh,
>>
>> You probably want to play with the runtime parameters "sim_rInit" and
>> "sim_expEnergy"  in your flash.par. These control the energy of the
>> explosion and the initial radius in which it is deposited for the initial
>> conditions. You'll want to scale these to your much larger domain.
>>
>> The FLASH setup of the problem including the many runtime parameters as
>> well as the dimensionless analytic solution are described here
>> <https://flash.rochester.edu/site/flashcode/user_support/flash4_ug_4p7/node187.html#SECTION010114000000000000000> in
>> the user guide.
>>
>> *************************************************************************
>> Adam Reyes
>>
>> Code Group Leader, Flash Center for Computational Science
>> Research Scientist, Dept. of Physics and Astronomy
>> University of Rochester
>> River Campus: Bausch and Lomb Hall, 369
>> 500 Wilson Blvd. PO Box 270171, Rochester, NY 14627
>> Email adam.reyes at rochester.edu
>> Web https://flash.rochester.edu
>>  (he / him / his)
>>
>> *************************************************************************
>>
>> On Wed, Jul 19, 2023 at 7:44 AM Yên Bình <yenbinhpy308 at gmail.com> wrote:
>>
>>> Dear FLASH Users,
>>>
>>> I have a problem with setting up the code for Sedov. I try to set up the
>>> simulation box to 10^16 x 10^16:
>>>
>>> *sim_xctr = 5.0e15*
>>> *sim_yctr = 5.0e15*
>>>
>>>
>>>
>>> *xmin = 0.xmax = 1.0e16ymin = 0.ymax = 1.0e16*
>>>
>>> But when I try to run the simulation, the plot file does not have any
>>> sign of an explosion (just a blank). I think it's because the box is too
>>> big so I need to set the right value of the running time:
>>>
>>>
>>>
>>>
>>>
>>>
>>>
>>>
>>>
>>>
>>>
>>>
>>> *# checkpoint file output parameterscheckpointFileIntervalTime =
>>> 0.01checkpointFileIntervalStep = 0checkpointFileNumber = 0# plotfile output
>>> parametersplotfileIntervalTime = 1plotfileIntervalStep = 0plotfileNumber =
>>> 0nend            = 10000tmax            = 1.0e8*
>>>
>>> But even if I set the time very high (10^8), nothing happens.
>>>
>>> I really appreciate any help to solve this problem. Thank you in advance.
>>>
>>> Binh.
>>>
>>>
>>> [image: image.png]
>>>
>>>
>>> Nguyen Thi Yen Binh
>>>
>>> College Student
>>>
>>> Space Science and Satellite Technology
>>>
>>> University of Science and Technology of Ha Noi
>>>
>>> binhnty.bi12-060 at st.usth.edu.vn
>>>
>>> +84367533986
>>>
>>> Tuy Hoa, Phu Yen, Viet Nam
>>>
>>>
>>> <https://urldefense.com/v3/__https://www.facebook.com/wbinhw__;!!CGUSO5OYRnA7CQ!eHzY5sAvt9JVcQA5cJKZc6ry3jGnKBHVycUBNVDfU0KfAwSYw4QTVhFLbkxC7sqZjy9eKhHdhKa2njTzIDK0EINBdi0$>
>>>
>>>
>>> <https://urldefense.com/v3/__https://www.instagram.com/yenbinne/__;!!CGUSO5OYRnA7CQ!eHzY5sAvt9JVcQA5cJKZc6ry3jGnKBHVycUBNVDfU0KfAwSYw4QTVhFLbkxC7sqZjy9eKhHdhKa2njTzIDK0HytT4xk$>
>>>
>>>
>>> <https://urldefense.com/v3/__https://github.com/pea9__;!!CGUSO5OYRnA7CQ!eHzY5sAvt9JVcQA5cJKZc6ry3jGnKBHVycUBNVDfU0KfAwSYw4QTVhFLbkxC7sqZjy9eKhHdhKa2njTzIDK0Q3WW6Mc$>
>>>
>>>
>>>
>>>
>>>
>>>
>>> _______________________________________________
>>> flash-users mailing list
>>> flash-users at flash.rochester.edu
>>>
>>> For list info, including unsubscribe:
>>> https://flash.rochester.edu/mailman/listinfo/flash-users
>>>
>>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://flash.rochester.edu/pipermail/flash-users/attachments/20230723/ff4a1663/attachment.htm>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: image.png
Type: image/png
Size: 57168 bytes
Desc: not available
URL: <http://flash.rochester.edu/pipermail/flash-users/attachments/20230723/ff4a1663/attachment.png>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: image.png
Type: image/png
Size: 50864 bytes
Desc: not available
URL: <http://flash.rochester.edu/pipermail/flash-users/attachments/20230723/ff4a1663/attachment-0001.png>