<div dir="ltr">Hi Clio and Adam,<div><br></div><div> Thanks for the information. </div><div><br></div><div> I think in order to compute the base flow I need the ###.box file as well as a different ###.usr file as I think the boundary and initial conditions for the ###.usr file in the short_tests/dfh_cav directory is applicable for the perturbations only. Please send me the ###.box and the corresponding ###.usr file applicable for the base flow for dfh_cav case. </div><div><br></div><div>Swarandeep<br><div><br></div><div> </div></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Thu, Aug 17, 2017 at 7:03 PM, <span dir="ltr"><<a href="mailto:nek5000-users@lists.mcs.anl.gov" target="_blank">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">Hi Swarandeep,<br>
<br>
The baseflow has been computed with a Newton solver (which is not available in the nek examples). The additional code had been implemented in the version 1040 of nek5000, and there is no .par file for that. However, in general, if your case is laminar, to compute a baseflow, you just need to evaluate the steady state of your nonlinear case. For that, you can use the same .par file, and just change to the nonlinear solver (perturbationmodes = 0, perturbations = no, solveBaseflow = yes).<br>
<br>
The non dimensionalization is based on Hellums and Churchill (1962).<br>
<br>
The equations are:<br>
<br>
Pr (du/dt + (U*nabla)u + (u*nabla)U) = 1/sqrt(Ra)laplacian(u) -grad p + e_g theta<br>
d theta/dt + (u*nabla)Theta + (U*nabla)theta = 1/sqrt(Ra)laplacian(theta)<br>
div(u)=0<br>
<br>
where u, theta are respectively the velocity and temperature of the perturbation, U, Theta the velocity and temperature of the baseflow.<br>
<br>
You can find more information in the paper "Saglietti, Clio, et al. "Adjoint optimization of natural convection problems: differentially heated cavity." Theoretical and Computational Fluid Dynamics (2016): 1-17." on which the example is based.<br>
<br>
Regards,<br>
<br>
clio<br>
<br>
<br>
Quoting <a href="mailto:nek5000-users@lists.mcs.anl.gov" target="_blank">nek5000-users@lists.mcs.anl.go<wbr>v</a>:<br>
<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div><div class="h5">
Hi Adam,<br>
<br>
Thanks for the clarification. Can you please send me a sample .par file for<br>
the baseflow for the case of "dfh_cav". Also is it possible to get the<br>
equations solved for this problem(dfh_cav). I would like to know the way<br>
that this has been made dimensionless.<br>
<br>
For your kind information I briefly describe my case which I have been used<br>
to solve with the regular non-linear framework using older versions (.rea<br>
files). My setup is as follows:<br>
Geometry is a 3D cylindrical annullus.<br>
Temperature are fixed at inner and outer cylinders.<br>
Gravity is acting in -z direction<br>
There is a back ground rotation which is constant.<br>
<br>
Please let me know if any other info is required. Thanks for the help.<br>
<br>
Swarandeep<br>
<br>
<br>
<br>
<br>
On Thu, Aug 17, 2017 at 1:48 PM, <<a href="mailto:nek5000-users@lists.mcs.anl.gov" target="_blank">nek5000-users@lists.mcs.anl.g<wbr>ov</a>> wrote:<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">
Hi Swarandeep,<br>
<br>
Yes. You generate ###.rea in the standard way. In the old version it was<br>
txt and contained both mesh information and runtime parameter. One could<br>
split mesh information (###.re2, binary file) and runtime parameters<br>
(###.rea with negative number of elements or currently ###.par) using<br>
reatore2 tool. The map file is generated with genmap (can be txt or<br>
binary). So the procedure is more or less the same. If you have ###.rea<br>
with mesh you split it and generate ###.re2. Next you generate ###.map with<br>
genmap and finally ###.par by hand. Be careful to use tolls from the same<br>
release, as they've been changed significantly (e.g. binary ###.map file).<br>
Regards<br>
Adam<br>
<br>
<br>
On 16/08/17 18:52, <a href="mailto:nek5000-users@lists.mcs.anl.gov" target="_blank">nek5000-users@lists.mcs.anl.go<wbr>v</a> wrote:<br>
<br>
Hi Adam,<br>
<br>
Thank you very much for the help.<br>
<br>
Does the ###.map, ###.box file and from them the ###.rea file<br>
(in the dfh_cav case: the ###.re2 files) are created as done for the<br>
earlier versions i.e., by using genbox and genmap. Please clarify.<br>
<br>
Swarandeep<br>
<br>
<br>
<br>
On Wed, Aug 16, 2017 at 5:35 PM, <<a href="mailto:nek5000-users@lists.mcs.anl.gov" target="_blank">nek5000-users@lists.mcs.anl.g<wbr>ov</a>> wrote:<br>
<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Hi,<br>
<br>
The simplest way would be to look into existing example under<br>
Nek5000/short_test/dfh_cav/<br>
It is a differentially heated cavity case and covers both direct and<br>
adjoint linear solver. It is compatible with the current git version and<br>
uses ###.par file instead of ###.rea.<br>
Regards<br>
Adam<br>
<br>
<br>
<br>
<br>
<br>
On 11/08/17 11:05, <a href="mailto:nek5000-users@lists.mcs.anl.gov" target="_blank">nek5000-users@lists.mcs.anl.go<wbr>v</a> wrote:<br>
<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
This information might be outdated (old version) but previously the base<br>
flow was read in using the restart option in the rea file. The disturbance<br>
you can read in yourself in userchk.<br>
<br>
Philipp<br>
<br>
On 2017-08-11 10:49, <a href="mailto:nek5000-users@lists.mcs.anl.gov" target="_blank">nek5000-users@lists.mcs.anl.go<wbr>v</a> wrote:<br>
<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Hi Philipp,<br>
<br>
That was very much helpful. I have a steady base flow and<br>
temperature state. Please let me know where to provide the base flow.<br>
<br>
<br>
Thank you<br>
Swarandeep<br>
______________________________<wbr>_________________<br>
Nek5000-users mailing list<br>
<a href="mailto:Nek5000-users@lists.mcs.anl.gov" target="_blank">Nek5000-users@lists.mcs.anl.go<wbr>v</a><br>
<a href="https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users" rel="noreferrer" target="_blank">https://lists.mcs.anl.gov/mail<wbr>man/listinfo/nek5000-users</a><br>
<br>
______________________________<wbr>_________________<br>
</blockquote>
Nek5000-users mailing list<br>
<a href="mailto:Nek5000-users@lists.mcs.anl.gov" target="_blank">Nek5000-users@lists.mcs.anl.go<wbr>v</a><br>
<a href="https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users" rel="noreferrer" target="_blank">https://lists.mcs.anl.gov/mail<wbr>man/listinfo/nek5000-users</a><br>
<br>
</blockquote>
<br>
<br>
______________________________<wbr>_________________<br>
Nek5000-users mailing list<br>
<a href="mailto:Nek5000-users@lists.mcs.anl.gov" target="_blank">Nek5000-users@lists.mcs.anl.go<wbr>v</a><br>
<a href="https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users" rel="noreferrer" target="_blank">https://lists.mcs.anl.gov/mail<wbr>man/listinfo/nek5000-users</a><br>
<br>
</blockquote>
<br>
<br>
<br>
______________________________<wbr>_________________<br></div></div>
Nek5000-users mailing <a href="mailto:listNek5000-users@lists.mcs.an">listNek5000-users@lists.mcs.an</a><wbr>l.govhttps://<a href="http://lists.mcs.anl.gov/mailman/listinfo/nek5000-users" rel="noreferrer" target="_blank">lists.mcs.anl.gov<wbr>/mailman/listinfo/nek5000-user<wbr>s</a><span class=""><br>
<br>
<br>
<br>
______________________________<wbr>_________________<br>
Nek5000-users mailing list<br>
<a href="mailto:Nek5000-users@lists.mcs.anl.gov" target="_blank">Nek5000-users@lists.mcs.anl.go<wbr>v</a><br>
<a href="https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users" rel="noreferrer" target="_blank">https://lists.mcs.anl.gov/mail<wbr>man/listinfo/nek5000-users</a><br>
<br>
<br>
</span></blockquote></blockquote><div class="HOEnZb"><div class="h5">
<br>
<br>
______________________________<wbr>_________________<br>
Nek5000-users mailing list<br>
<a href="mailto:Nek5000-users@lists.mcs.anl.gov" target="_blank">Nek5000-users@lists.mcs.anl.go<wbr>v</a><br>
<a href="https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users" rel="noreferrer" target="_blank">https://lists.mcs.anl.gov/mail<wbr>man/listinfo/nek5000-users</a><br>
</div></div></blockquote></div><br></div>