<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=us-ascii">
<style type="text/css" style="display:none;"><!-- P {margin-top:0;margin-bottom:0;} --></style>
</head>
<body dir="ltr">
<div id="divtagdefaultwrapper" style="font-size:12pt;color:#000000;font-family:Calibri,Helvetica,sans-serif;" dir="ltr">
<p style="margin-top:0;margin-bottom:0"><br>
</p>
<p style="margin-top:0;margin-bottom:0">Dear Armand,</p>
<p style="margin-top:0;margin-bottom:0"><br>
</p>
<p style="margin-top:0;margin-bottom:0">The attached .usr file seems to work ok... just a few</p>
<p style="margin-top:0;margin-bottom:0">small modifications to the .usr file you provided.</p>
<p style="margin-top:0;margin-bottom:0"><br>
</p>
<p style="margin-top:0;margin-bottom:0">I think the main issue was the declaration of some of the</p>
<p style="margin-top:0;margin-bottom:0">variables in useric.</p>
<p style="margin-top:0;margin-bottom:0"><br>
</p>
<p style="margin-top:0;margin-bottom:0">Best,</p>
<p style="margin-top:0;margin-bottom:0"><br>
</p>
<p style="margin-top:0;margin-bottom:0">Paul</p>
<p style="margin-top:0;margin-bottom:0"><br>
</p>
</div>
<hr style="display:inline-block;width:98%" tabindex="-1">
<div id="divRplyFwdMsg" dir="ltr"><font face="Calibri, sans-serif" style="font-size:11pt" color="#000000"><b>From:</b> Nek5000-users <nek5000-users-bounces@lists.mcs.anl.gov> on behalf of nek5000-users@lists.mcs.anl.gov <nek5000-users@lists.mcs.anl.gov><br>
<b>Sent:</b> Thursday, June 7, 2018 3:51:09 AM<br>
<b>To:</b> nek5000-users@lists.mcs.anl.gov<br>
<b>Subject:</b> [Nek5000-users] Initial Condition & Recycling</font>
<div> </div>
</div>
<div class="BodyFragment"><font size="2"><span style="font-size:11pt;">
<div class="PlainText">Dear All,<br>
<br>
<br>
In an attempt to simulate an established turbulent flow over a flat <br>
plate with a recycling method for the inlet, it has to my attention that <br>
a (1/7)th Power Law would be suitable as an initial condition. However <br>
the solver divergences and I think I know why.<br>
<br>
Here the layout of the recycling procedure :<br>
<br>
1. call set_inflow_fpt in userchk<br>
<br>
2. If it is the 1st call then call setp_inflow_fpt_setup<br>
<br>
3. call rescale_fpt (which rescales the copied field to a user specified <br>
mean velocity ubar)<br>
<br>
4.In the rescale_fpt routine we call get_flux_and_area which gives the <br>
flux and area of a face with cbc == 'v ' so a user specified velocity <br>
boundary condition.<br>
<br>
4. a) As the initial condition is wrong once that routine sums the flux <br>
on the boundary it returns zero in this case and the variable scale <br>
which is divided by the latter returns infinity !<br>
<br>
<br>
However I set the initial conditions as in many other examples, such as <br>
the turbInflow. Would have an insight that could lead me to a solution ?<br>
<br>
<br>
Please find attached the related files.<br>
<br>
<br>
Sincerely,<br>
<br>
Armand, ONERA - The French Aerospace Lab.<br>
<br>
<br>
</div>
</span></font></div>
</body>
</html>