[Nek5000-users] Initial Condition & Recycling

[nek5000-users at lists.mcs.anl.gov]

Dear Armand,


The attached .usr file seems to work ok... just a few

small modifications to the .usr file you provided.


I think the main issue was the declaration of some of the

variables in useric.


Best,


Paul


________________________________
From: Nek5000-users <nek5000-users-bounces at lists.mcs.anl.gov> on behalf of nek5000-users at lists.mcs.anl.gov <nek5000-users at lists.mcs.anl.gov>
Sent: Thursday, June 7, 2018 3:51:09 AM
To: nek5000-users at lists.mcs.anl.gov
Subject: [Nek5000-users] Initial Condition & Recycling

Dear All,


In an attempt to simulate an established turbulent flow over a flat
plate with a recycling method for the inlet, it has to my attention that
a (1/7)th Power Law would be suitable as an initial condition. However
the solver divergences and I think I know why.

Here the layout of the recycling procedure :

1. call set_inflow_fpt in userchk

2. If it is the 1st call then call setp_inflow_fpt_setup

3. call rescale_fpt (which rescales the copied field to a user specified
mean velocity ubar)

4.In the rescale_fpt routine we call get_flux_and_area which gives the
flux and area of a face with cbc == 'v   ' so a user specified velocity
boundary condition.

4. a) As the initial condition is wrong once that routine sums the flux
on the boundary it returns zero in this case and the variable scale
which is divided by the latter returns infinity !


However I set the initial conditions as in many other examples, such as
the turbInflow. Would have an insight that could lead me to a solution ?


Please find attached the related files.


Sincerely,

Armand, ONERA - The French Aerospace Lab.


-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20180607/935a7da1/attachment-0001.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: z.usr
Type: application/octet-stream
Size: 31771 bytes
Desc: z.usr
URL: <http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20180607/935a7da1/attachment-0001.obj>