Skip to content
This repository has been archived by the owner on Nov 7, 2019. It is now read-only.

The inflow/outflow condition cause the NaN problem #1

Open
GoogleCodeExporter opened this issue Aug 24, 2015 · 2 comments
Open

The inflow/outflow condition cause the NaN problem #1

GoogleCodeExporter opened this issue Aug 24, 2015 · 2 comments

Comments

@GoogleCodeExporter
Copy link

What steps will reproduce the problem?
1. I perform the experiment with the channel flow just according to the 
official code
2. And then I tried to modify the boundary condition(0 2 0) to (2 2 0),so that 
I can add the inflow and outflow condition.I did the modify because I found 
that if nclx==2(incompact3d.f90),then the inflow/outflow is added in the 
simulation. 
3. Oohs,sucks!  I tried to track the result with times going on, but the 
divergence U move to NaN.  Later on,I check the boundary conditions of the 
computational domain in the 《high order compact  schemes ...》:"The use of 
the Dirichlet conditions on the velocity(for no-slip or open conditions)need to 
be defined acoording to the time advancement procedure".
4. So,I vim the schemes.f90,but I don't know how to modify it,can someone give 
me some ideas?
-------------------------------------------------------
What is the expected output? What do you see instead?

   My expected output is the turbulent development process,rather than the Tollmien–Schlichting wave(http://en.wikipedia.org/wiki/Tollmien%E2%80%93Schlichting_wave) with the parallel flow assumption.
--------------------------------------------------------
What version of the product are you using? On what operating system?
  Incompact3d under ubuntu13.10 

Please provide any additional information below.


Original issue reported on code.google.com by [email protected] on 18 Apr 2014 at 4:32

@GoogleCodeExporter
Copy link
Author

Dear Ye Zhaoliang,

If you want to change the BC, then you need to change numbers of mesh nodes as 
well. If you want to use inflow/outflow then you need to define them of course 
and you need to pay attention to the flow rate. If you want to do a spataially 
evolving channel flow you will have to take a huge domain in the streamwise 
direction if you want to have a chance to get a turbulent state. If you want to 
do that I will suggest a recycling procedure.

Sylvain Laizet

Original comment by [email protected] on 2 May 2014 at 7:23

@GoogleCodeExporter
Copy link
Author

Thanks for your attention.
I have three questions about your answers.
1: what does you mean about the flow rate? the range of the flow rate?

2: Take a huge domain in the streamwise direction -----> 128? 1535? Maybe you 
can give us a story about the huge domain.

3: Please forgive my ignorence  and satisfy my thirsty for knowledge. What does 
it mean the recycling procedure?

Original comment by [email protected] on 3 May 2014 at 5:41

Sign up for free to subscribe to this conversation on GitHub. Already have an account? Sign in.
Projects
None yet
Development

No branches or pull requests

3 participants