Hi Paul,<div><br></div><div>Sorry for the late reply.</div><div><br></div><div>That would be great. I will get back once I start with 3D runs.</div><div><br></div><div>Regards,</div><div>Pradeep<br><br><div class="gmail_quote">
On Tue, Feb 21, 2012 at 7:13 AM, <span dir="ltr"><<a href="mailto:nek5000-users@lists.mcs.anl.gov">nek5000-users@lists.mcs.anl.gov</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<br>
Hi Pradeep,<br>
<br>
There was one thing that I'd intended to mention regarding mesh motion.<br>
<br>
The file that I set up uses the most general ALE formulation,<br>
where the mesh velocity is computed through an elasticity solver.<br>
The performance for this solver has not been optimized, but there<br>
are several tricks we can use to make it faster. At the moment,<br>
I would guess that well over half the time is spent in this solver,<br>
which could be an issue if you plan to look at 3D problems. As<br>
you get into production runs we should examine the performance and<br>
optimize if necessary.<br>
<br>
Regards,<br>
<br>
Paul<div class="im"><br>
<br>
<br>
On Thu, 16 Feb 2012, <a href="mailto:nek5000-users@lists.mcs.anl.gov" target="_blank">nek5000-users@lists.mcs.anl.<u></u>gov</a> wrote:<br>
<br>
</div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div class="im">
Hi Paul,<br>
<br>
Thank you so much for your help!<br>
<br>
Regards,<br>
Pradeep<br>
<br>
On Thu, Feb 16, 2012 at 7:28 AM, <<a href="mailto:nek5000-users@lists.mcs.anl.gov" target="_blank">nek5000-users@lists.mcs.anl.<u></u>gov</a>> wrote:<br>
<br>
</div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div class="im">
<br>
PS -- an example of the results can be found here:<br>
<br></div>
<a href="http://www.mcs.anl.gov/~fischer/ocyl.**gif" target="_blank">www.mcs.anl.gov/~fischer/ocyl.<u></u>**gif</a><<a href="http://www.mcs.anl.gov/%7Efischer/ocyl.gif" target="_blank">http://www.mcs.anl.gov/%<u></u>7Efischer/ocyl.gif</a>><br>
<br>
<br>
<br>
<br>
On Thu, 16 Feb 2012, nek5000-users@lists.mcs.anl.**<u></u>gov<<a href="mailto:nek5000-users@lists.mcs.anl.gov" target="_blank">nek5000-users@lists.mcs.<u></u>anl.gov</a>>wrote:<br>
<br>
<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div class="im">
Hi Pradeep,<br>
<br>
Attached is a case that appears to work. Hope this helps.<br>
<br>
Paul<br>
<br>
<br></div>
On Wed, 15 Feb 2012, nek5000-users@lists.mcs.anl.**<u></u>gov<<a href="mailto:nek5000-users@lists.mcs.anl.gov" target="_blank">nek5000-users@lists.mcs.<u></u>anl.gov</a>>wrote:<br>
<br>
Thanks Paul and Josh,<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div class="im">
<br>
I shall start by making the modification in the rea file and trying visit<br>
2.3.2<br>
<br>
Regards,<br>
Pradeep<br>
<br></div>
On Wed, Feb 15, 2012 at 9:39 AM, <nek5000-users@lists.mcs.anl.*<u></u>*gov<<a href="mailto:nek5000-users@lists.mcs.anl.gov" target="_blank">nek5000-users@lists.mcs.<u></u>anl.gov</a>>><br>
wrote:<br>
<br>
Pradeep,<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div class="im">
<br>
Which version of Visit are you using? We have experienced similar<br>
issues with version 2.4 when viewing 2D meshes. If this is the case,<br>
try using 2.3.2 and see if you have the same problem.<br>
<br>
Josh<br>
<br></div>
On Tue, Feb 14, 2012 at 8:08 AM, <nek5000-users@lists.mcs.anl.*<u></u>*gov<<a href="mailto:nek5000-users@lists.mcs.anl.gov" target="_blank">nek5000-users@lists.mcs.<u></u>anl.gov</a>>><div><div class="h5"><br>
wrote:<br>
<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Hi Neks,<br>
<br>
I'm trying to simulate 2D flow over an oscillating cylinder using<br>
moving<br>
mesh, and tried to set up a test case by looking at the peris example.<br>
I<br>
however am getting vanishing Jacobians after the first few timesteps.<br>
<br>
I am modifying the the boundary velocity of the cylinder in userbc<br>
<br>
</blockquote>
instead<br>
<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
of userchk (as in the peris example). Would it be possible for you to<br>
<br>
</blockquote>
look<br>
<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
at my .usr file and tell me what I am doing wrong?<br>
<br>
I have attached the .rea, .usr, and SIZE with this mail.<br>
<br>
Further, when I try to view the mesh of the output (or the xyz fld<br>
outputfile) in visit, visit crashes with the following error:<br>
VisIt's viewer exited abnormally! Aborting the Graphical User<br>
Interface.<br>
<br>
Lastly, would it be possible to simulate moving walls with large<br>
deformations (i.e. where multiple elements get deformed), where the<br>
deformation is known ahead of time? (eg. a flapping aerofoil)<br>
<br>
Thank you for your help!<br>
<br>
Regards,<br>
Pradeep<br>
<br>
Extract from logfile:<br>
<br>
60 U-PRES gmres: 100 7.5345E+42 1.0000E-02 4.0062E+57<br>
6.7556E-03 1.5093E-02<br>
ERROR: alphad .le. 0 in econjp -1.31175831444004495E+179 16<br>
ERROR: alphad .le. 0 in econjp -1.31175831444004495E+179 16<br>
ERROR: alphad .le. 0 in econjp -1.31175831444004495E+179 16<br>
ERROR: alphad .le. 0 in econjp -1.31175831444004495E+179 16<br>
60 DNORM, DIVEX 1.18120867781751263E+043<br>
7.53448962619089522E+042<br>
60 6.0000E-02 1.6697E-02 Fluid done<br>
DT/DTCFL/DTFS/DTINIT 0.635E-45 0.635E-45 0.000E+00<br>
0.100E-02<br>
Step 61, t= 6.0800000E-02, DT= 8.0000000E-04, C=******* 2.2667E-01<br>
1.8287E-02<br>
Solving for fluid<br>
<br>
<br>
3 ERROR: Vanishing Jacobian near 2th node of element<br>
<br>
</blockquote>
3.<br>
<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<br>
<br>
2 ERROR: Vanishing Jacobian near 2th node of element<br>
<br>
</blockquote>
1.<br>
<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
2.27348348709230109E+081 -5.12478988977734591E+080<br>
-2.60483379627720531E+079 3.21465902039155941E+039<br>
3 xyz: -5.30056E+39 -5.84642E+39<br>
3 xyz: -5.24850E+39 -6.93109E+39<br>
<br>
<br>
3 ERROR: Vanishing Jacobian near 2th node of element<br>
<br>
</blockquote>
4.<br>
<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
2 xyz: 1.25061E+38 4.64957E+39<br>
9.79493533554062390E+080 -3.01548938509156252E+080<br>
2 xyz: -2.46589E-08 -5.64092E-01<br>
<br>
<br>
3 xyz: -5.30056E+39 -5.84642E+39<br>
3 xyz: -5.24850E+39 -6.93109E+39<br>
2 ERROR: Vanishing Jacobian near 3th node of element<br>
<br>
</blockquote>
2.<br>
<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<br>
<br>
3 ERROR: Vanishing Jacobian near 2th node of element<br>
<br>
</blockquote>
5.<br>
<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
1.43125817410005932E+039 -6.03607926692375119E+038<br>
5.54918579647991880E+039 -1.67302572815141788E+039<br>
3 xyz: -5.30056E+39 -5.84642E+39<br>
2 xyz: -2.46589E-08 -5.64092E-01<br>
3 xyz: -5.24850E+39 -6.93109E+39<br>
<br>
<br>
3 ERROR: Vanishing Jacobian near 2th node of element<br>
<br>
</blockquote>
6.<br>
<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
2 xyz: -3.07794E-08 -7.04112E-01<br>
<br>
<br>
2 ERROR: Vanishing Jacobian near 7th node of element<br>
<br>
</blockquote>
7.<br>
<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
8.35923896949543198E+038 -3.55455289325888722E+038<br>
3 xyz: -5.30056E+39 -5.84642E+39<br>
1.21747337192549314E+039 -1.79454406120891328E+038<br>
3 xyz: -5.24850E+39 -6.93109E+39<br>
2 xyz: -5.66434E-08 -1.29581E+00<br>
2 xyz: -6.27639E-08 -1.43583E+00<br>
<br>
<br>
2 ERROR: Vanishing Jacobian near 9th node of element<br>
<br>
</blockquote>
8.<br>
<br>
</div></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div><div class="h5">
1.17802538305091326E+039 -4.21489851051369387E+038<br>
2 xyz: -6.55671E-08 -1.49996E+00<br>
2 xyz: -3.25016E+38 4.93026E+39<br>
call outfld: ifpsco: F<br>
<br>
61 6.0800E-02 Write checkpoint:<br>
<br>
61 6.0800E-02 OPEN:<br>
xyzcyl2.fld01<br>
Jac error 1, setting p66=4, ifxyo=t<br>
<br>
call exitt: dying ...<br>
<br>
<br></div></div>
______________________________<u></u>**_________________<br>
Nek5000-users mailing list<br>
Nek5000-users@lists.mcs.anl.**<u></u>gov <<a href="mailto:Nek5000-users@lists.mcs.anl.gov" target="_blank">Nek5000-users@lists.mcs.anl.<u></u>gov</a>><br>
<a href="https://lists.mcs.anl.gov/**mailman/listinfo/nek5000-users" target="_blank">https://lists.mcs.anl.gov/**<u></u>mailman/listinfo/nek5000-users</a><u></u><<a href="https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users" target="_blank">https://lists.mcs.anl.gov/<u></u>mailman/listinfo/nek5000-users</a><u></u>><br>
<br>
<br>
</blockquote><div class="im">
<br>
<br>
--<br>
Josh Camp<br>
<br>
"All that is necessary for the triumph of evil is that good men do<br>
nothing" -- Edmund Burke<br></div>
______________________________<u></u>**_________________<br>
Nek5000-users mailing list<br>
Nek5000-users@lists.mcs.anl.**<u></u>gov <<a href="mailto:Nek5000-users@lists.mcs.anl.gov" target="_blank">Nek5000-users@lists.mcs.anl.<u></u>gov</a>><br>
<a href="https://lists.mcs.anl.gov/**mailman/listinfo/nek5000-users" target="_blank">https://lists.mcs.anl.gov/**<u></u>mailman/listinfo/nek5000-users</a><u></u><<a href="https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users" target="_blank">https://lists.mcs.anl.gov/<u></u>mailman/listinfo/nek5000-users</a><u></u>><br>
<br>
<br>
</blockquote></blockquote>
______________________________<u></u>**_________________<br>
</blockquote>
Nek5000-users mailing list<br>
Nek5000-users@lists.mcs.anl.**<u></u>gov <<a href="mailto:Nek5000-users@lists.mcs.anl.gov" target="_blank">Nek5000-users@lists.mcs.anl.<u></u>gov</a>><br>
<a href="https://lists.mcs.anl.gov/**mailman/listinfo/nek5000-users" target="_blank">https://lists.mcs.anl.gov/**<u></u>mailman/listinfo/nek5000-users</a><u></u><<a href="https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users" target="_blank">https://lists.mcs.anl.gov/<u></u>mailman/listinfo/nek5000-users</a><u></u>><br>
<br>
</blockquote>
<br>
</blockquote><div class="HOEnZb"><div class="h5">
______________________________<u></u>_________________<br>
Nek5000-users mailing list<br>
<a href="mailto:Nek5000-users@lists.mcs.anl.gov" target="_blank">Nek5000-users@lists.mcs.anl.<u></u>gov</a><br>
<a href="https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users" target="_blank">https://lists.mcs.anl.gov/<u></u>mailman/listinfo/nek5000-users</a><br>
</div></div></blockquote></div><br></div>