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

Yên Bình yenbinhpy308 at gmail.com
Tue Jul 25 00:29:59 EDT 2023 

Dear Adam Reyes,

I did just as you suggested, I changed those parameters:
sim_nsubzones = 16
lrefine_min_init = 4
The initial frame does have the grid refinement but the following frames
seem to be a uniform one with this weird shape. I also try to increase the
"lrefine_min" but this keeps happening.
[image: image.png]



This is my "flash.par" set up:






*smlrho = 1.E-27smallp = 1.E-15sim_nsubzones = 16lrefine_min_init =
4lrefine_min = 4dtmin = 0.1E-20*




*sim_pAmbient = 1.E-12sim_rhoAmbient = 1.E-24sim_expEnergy = 1.E51*



*sim_rInit = 1.E4sim_xctr = 0.5sim_yctr = 0.5*


*gamma           = 1.4*





*xmin = 0.xmax = 1.E16ymin = 0.ymax = 1.E16*





*xl_boundary_type = "axisymmetric"xr_boundary_type =
"outflow"yl_boundary_type = "axisymmetric"yr_boundary_type = "outflow"*


*cfl = 0.5*

*checkpointFileIntervalTime = 0.01*


*nend            = 10000tmax            = 1.E8*

*eintSwitch     = 1.e-4*




*lrefine_max     = 8refine_var_1 = "dens"refine_var_2 = "pres"*







*# Uniform Grid specific parameters #iGridSize = 8   #global number of
gridpoints along x, excluding gcells#jGridSize = 8   #global number of
gridpoints along y, excluding gcells#kGridSize = 1iProcs = 8 #num procs in
i directionjProcs = 8 #num procs in j direction*

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://www.facebook.com/wbinhw>

<https://www.instagram.com/yenbinne/>

<https://github.com/pea9>








Vào Th 2, 24 thg 7, 2023 vào lúc 03:58 Adam Reyes <
adam.reyes at rochester.edu> đã viết:

> Hi Binh,
>
> I suspect the initial radius for the energy is too small for the grid
> refinement that you're using. The users guide suggests that this be at
> least 3.5x your smallest grid spacing, and I've also found this to be a
> good rule of thumb. It could also be helpful to increase the parameter
> "sim_nsubzones" from its default of 7, or increase "lrefine_min_init" to
> something closer to "lrefine_max" so that paramesh will refine on the
> initial energy perturbation.
> *************************************************************************
> 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 Sun, Jul 23, 2023 at 1:07 PM Yên Bình <yenbinhpy308 at gmail.com> wrote:
>
>> Hi Adam Reyes,
>>
>> Thank you for the advice.
>> I reduced the smlrho and smallp, the simulation ran fine and produced
>> checkpoint files as expected.
>>
>> However when I try to plot, nothing happens (I checked the log file).
>> That's very strange.
>>
>> Best regards,
>> 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!dtB5HBalUMYskKnEInTyvY4WO6yQeRMZEwd-LjXh3Z7eRHUmnvxgoLH-nVwLW4YX7HTnfF0Y39aXTfMudrczUcFnB94$>
>>
>>
>> <https://urldefense.com/v3/__https://www.instagram.com/yenbinne/__;!!CGUSO5OYRnA7CQ!dtB5HBalUMYskKnEInTyvY4WO6yQeRMZEwd-LjXh3Z7eRHUmnvxgoLH-nVwLW4YX7HTnfF0Y39aXTfMudrcz_cOmyx0$>
>>
>>
>> <https://urldefense.com/v3/__https://github.com/pea9__;!!CGUSO5OYRnA7CQ!dtB5HBalUMYskKnEInTyvY4WO6yQeRMZEwd-LjXh3Z7eRHUmnvxgoLH-nVwLW4YX7HTnfF0Y39aXTfMudrczvr4vIr4$>
>>
>>
>>
>>
>>
>>
>>
>>
>> Vào CN, 23 thg 7, 2023 vào lúc 16:58 Adam Reyes <adam.reyes at rochester.edu>
>> đã viết:
>>
>>> 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/20230725/29304bc9/attachment-0001.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/20230725/29304bc9/attachment-0004.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/20230725/29304bc9/attachment-0005.png>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: image.png
Type: image/png
Size: 54325 bytes
Desc: not available
URL: <http://flash.rochester.edu/pipermail/flash-users/attachments/20230725/29304bc9/attachment-0006.png>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: image.png
Type: image/png
Size: 53523 bytes
Desc: not available
URL: <http://flash.rochester.edu/pipermail/flash-users/attachments/20230725/29304bc9/attachment-0007.png>

More information about the flash-users mailing list