<div dir="ltr"><div dir="ltr">On Fri, Sep 4, 2020 at 11:09 AM Thibault Bridel-Bertomeu <<a href="mailto:thibault.bridelbertomeu@gmail.com">thibault.bridelbertomeu@gmail.com</a>> wrote:<br></div><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">Dear all,<div><br></div><div dir="auto">I am trying to figure out a working set of options for TS tuto ex11 involving the 3D Euler Linear Wave test. I cannot seem to go anywhere as it keeps crashing saying the physics did not compute the wave speed ... anyone could advise ? </div></blockquote><div><br></div><div>Does the shock test run for you?</div><div><br></div><div>  test:<br>    suffix: shock_0<br>    requires: p4est !single !complex<br>    args: -ufv_vtk_interval 0 -monitor density,energy -f -grid_size 2,1 -grid_bounds -1,1.,0.,1 -bc_wall 1,2,3,4 -dm_type p4est -dm_forest_partition_overlap 1 -dm_forest_maximum_refinement 6 -dm_forest_minimum_refinement 2 -dm_forest_initial_refinement 2 -ufv_use_amr -refine_vec_tagger_box 0.5,inf -coarsen_vec_tagger_box 0,1.e-2 -refine_tag_view -coarsen_tag_view -physics euler -eu_type iv_shock -ufv_cfl 10 -eu_alpha 60. -grid_skew_60 -eu_gamma 1.4 -eu_amach 2.02 -eu_rho2 3. -petscfv_type leastsquares -petsclimiter_type minmod -petscfv_compute_gradients 0 -ts_max_time 0.5 -ts_ssp_type rks2 -ts_ssp_nstages 10 -ufv_vtk_basename ${wPETSC_DIR}/ex11<br> <br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="auto">It is not however completely about ex11. I am trying to figure out how to set up periodic boundary conditions for a box mesh. I tried to get some infos here and there following examples : I got that you have to set DM_BOUNDARY_PERIODIC for the dimensions that you want in the array that you pass to DM Create Box Mesh, but i am not clear whether I have to do anything else after. Do I have to call DM Set Periodicity ? I also saw sometimes calls to DMLocalizeCoordinates ? ... ?</div></blockquote><div><br></div><div>DMSetPeriodicity() will be called by CreateBoxMesh(). You will need to call DMLocalizeCoordinates().</div><div> </div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="auto">And are there any constraints regarding the number of cells in the periodic dimensions ? I would like to just have 1 cell in Y dimension for instance to have a pseudo 1D computation, is that ok ?</div></blockquote><div><br></div><div>Yes. 1 cell will not work right now. We now have everything in place to make it work, but have not put the code into CreateBoxMesh() to do it. For now, you can just</div><div>put 3 cells in this direction. If you make an Issue, I will fix it as soon as I can.</div><div><br></div><div>  Thanks,</div><div><br></div><div>     Matt</div><div> </div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="auto">Thank you for your help,</div><div dir="auto"><br></div><div dir="auto">Thibault</div>-- <br><div dir="ltr"><div dir="ltr"><div><div dir="ltr"><div dir="ltr"><div><div><div><div><div>Thibault Bridel-Bertomeu<br>—<br></div></div></div></div>Eng, MSc, PhD</div><div>Research Engineer</div><div>CEA/CESTA</div><div>33114 LE BARP</div><div>Tel.: (+33)557046924</div><div>Mob.: (+33)611025322<br></div><div>Mail: <a href="mailto:thibault.bridelbertomeu@gmail.com" target="_blank">thibault.bridelbertomeu@gmail.com</a><br></div></div></div></div></div></div>
</blockquote></div><br clear="all"><div><br></div>-- <br><div dir="ltr" class="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div>What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.<br>-- Norbert Wiener</div><div><br></div><div><a href="http://www.cse.buffalo.edu/~knepley/" target="_blank">https://www.cse.buffalo.edu/~knepley/</a><br></div></div></div></div></div></div></div></div>