[Nek5000-users] RB convection with high Rayleigh numbers
nek5000-users at lists.mcs.anl.gov
nek5000-users at lists.mcs.anl.gov
Sun May 16 09:25:04 CDT 2010
Hi Mani,
I find that I need to use something like:
argt = (t8-time)/0.05
if (time.gt.t8) ra = 1.e7 + 0.9e8*(1-exp(argt))
to get the run to transition smoothly to the higher Ra
values, but this seems to work.
Paul
On Sun, 16 May 2010, nek5000-users at lists.mcs.anl.gov wrote:
>
> Hi Mani,
>
> Note that your other email did not go to the users list because
> of the size of the attachments. I got it though and tried it.
>
> It looks like the following steps should work ---
>
> .reduce COURANT to 2.0 instead of 5.0
> .set TORDER=2
> .set initial DT (param 12) to .001
> .reduce lx1 to 5 and lxd to 8 in your SIZE file
>
>
> I tried this and it seems to work ok, save that I ran my
> case w/ Pr=1 (param 2).
>
> Your case w/ p2=6.8 is a bit harder, so you may need to also
> use the reduced initial Ra values per your .usr file.
>
> Once you are at the target Ra, you can increase lx1, lxd, and
> Torder, then restart.
>
>
> The issue with the walls is that they trap the fluid in the
> corner. Then when the heating is turned on you get a stagnant hot spot that
> builds up until it finally goes unstable --- this happens on a timescale
> faster than the
> early-time convective timescale (which is relatively long
> until the flow really starts moving), so the explicit
> treatment of the buoyancy leads to an instability.
>
> Reducing Courant to 2.0 instead of 5.0 forces nek5000 to
> reduce dt more severely. Note that we're able to exceed
> CFL of unity for convection because of the OIFS convective
> time integrator. The explicit Boussinesq forcing, however,
> has it's own timestep restriction and tends to dominate the
> stability constraints at early times until the flow is moving.
>
> Paul
>
>
>
>
>
> On Sun, 16 May 2010, nek5000-users at lists.mcs.anl.gov wrote:
>
>> Hi Paul,
>>
>> Thank you for your time on this problem. The geometry is a 2D box with
>> no-slip velocity boundary conditions on all sides and the temperature
>> boundary conditions are insulating on the side walls and the usual fixed
>> temperatures on the top and bottom walls. I will try out your code
>> immediately.
>>
>> Regards,
>> Mani chandra
>> On 05/16/2010 09:14 AM, nek5000-users at lists.mcs.anl.gov wrote:
>>>
>>> Hi Mani,
>>>
>>> I've added an example of a case that should run out to Ra=10^8
>>> in the /benard directory of the repo (high_ray.blah is the case).
>>>
>>> I've run this with lx1=5, lxd=9 and it will get to Ra=10^8.
>>> (In reality, I ran tests with 4 separate runs, 10^k, k=5,6,7,8. I then a
>>> created "high_ray.usr" file that would march through this sequence in a
>>> single run...)
>>>
>>> Once you establish the flow at the desired Ra, you can then
>>> bump lx1 to say 7 or 8 --- that should be sufficient for Ra=10^8
>>> I think. Even lx1=5 is looking pretty good with the given mesh.
>>>
>>>
>>> Also, you can reduce the CFL in the .rea file and increase Torder
>>> to 3, so that you will increase your temporal accuracy.
>>>
>>> I initially set the parameters so that it would march quickly
>>> to Ra=10^8.
>>>
>>> The high-Ra test case uses periodic bcs and a uniform mesh in X.
>>> If you want vertical walls then you should likely use a nonuniform
>>> mesh in X. The "high_ray.box" file shows how the nonuniform mesh in Y is
>>> created. A similar set up could be used for nonuniformity
>>> in both X and Y.
>>>
>>>
>>> Paul
>>>
>>>
>>>
>>> On Sun, 16 May 2010, nek5000-users at lists.mcs.anl.gov wrote:
>>>
>>>> Hi Paul,
>>>>
>>>> Thanks for the reply. I have a doubt whether one should increase the
>>>> polynomial order or the number of elements with a fixed order for this
>>>> case. As I understand from your illustration of the rotating cone case,
>>>> increasing the polynomial order leads to a much better solution than with
>>>> the same number of points generated using more elements with a lower
>>>> order polynomial.
>>>>
>>>> Also how do I use the fast tensor-product solver? Is there some flag that
>>>> I need to turn on?
>>>>
>>>> Regards,
>>>> Mani chandra
>>>> On 05/15/2010 09:05 PM, nek5000-users at lists.mcs.anl.gov wrote:
>>>>>
>>>>> Hi Mani,
>>>>>
>>>>> I've not tried to go this high yet, but there are several std
>>>>> techniques to pushing the envelope.
>>>>>
>>>>> First, you should try to get an estimate of your expected resolution
>>>>> requirements at the target Ra. Then, guess that you might want to
>>>>> run that simulation at say, lx1=10 or 12; and build the corresponding
>>>>> mesh (number of elements in x and y); but start the compuation at
>>>>> much lower Ra, with lx1=4, say.
>>>>>
>>>>> Then work your way up in Ra, say by a factor of 10 in each successive
>>>>> run, while also increasing lx1 as needed.
>>>>>
>>>>> Note that for the flow in a box, you might find it faster to use
>>>>> the fast tensor-product solver for the pressure, if using the Pn-Pn-2
>>>>> method --- this is not always the case, but for some very difficult
>>>>> problems (e.g., highly refined meshes), it does work better.
>>>>>
>>>>> Paul
>>>>>
>>>>>
>>>>> On Sat, 15 May 2010, nek5000-users at lists.mcs.anl.gov wrote:
>>>>>
>>>>>> Hi,
>>>>>>
>>>>>> I'm trying to simulate RB convection in a 2D box of aspect ratio
>>>>>> 2:1. I'd like to simulate cases with Rayleigh numbers of the order of
>>>>>> 10^7 - 10^8. But I'm barely able to do 10^5 and that too with a time
>>>>>> step of 10^-6. For the mesh I took 8 spectral elements in the x
>>>>>> direction and 4 spectral elements in the y direction with 48 points
>>>>>> within each of them. Is there any way one can reach higher Ra numbers?
>>>>>>
>>>>>> Regards,
>>>>>> Mani chandra
>>>>>> _______________________________________________
>>>>>> Nek5000-users mailing list
>>>>>> Nek5000-users at lists.mcs.anl.gov
>>>>>> https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users
>>>>>>
>>>>> _______________________________________________
>>>>> Nek5000-users mailing list
>>>>> Nek5000-users at lists.mcs.anl.gov
>>>>> https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users
>>>>>
>>>>
>>>> _______________________________________________
>>>> Nek5000-users mailing list
>>>> Nek5000-users at lists.mcs.anl.gov
>>>> https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users
>>>>
>>> _______________________________________________
>>> Nek5000-users mailing list
>>> Nek5000-users at lists.mcs.anl.gov
>>> https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users
>>>
>>
>> _______________________________________________
>> Nek5000-users mailing list
>> Nek5000-users at lists.mcs.anl.gov
>> https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users
>>
> _______________________________________________
> Nek5000-users mailing list
> Nek5000-users at lists.mcs.anl.gov
> https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users
>
More information about the Nek5000-users
mailing list