I have dealt with the atmospheric for only few weeks, so I warn you, I'm still a beginner.
I am trying to run a simulation where wind blows over a mountain. The case is pretty easy to understand :
- I want a transcient simulation,
- Wind is blowing from the east,
- North and south are symmetries
- The floor is divided in 2 : ground and mountain
- Top is set to inlet, meteorological profile from data (such as advocated in the user guide)
- The domain is : 6 km of length (3 before the ountain and 3 after). - I know it is a little bit to small but I will resize it soon. -
- For now, it is dry atmosphere but I want eventually humid atmosphere
However in 3D, my calculation diverges regardless of the efforts made to stabilize the model. I checked several times my mesh which is ok (salome mesh quality report linked below).
** With automatic meteo file generated, I tried to change the inlet / outlet :
The simulation crashes when :
1) east: inlet, meteorological profile from data, Automatic inlet/outlet nature from data
west: inlet, meteorological profile from data, Automatic inlet/outlet nature from data
2) east: inlet, meteorological profile from data,
west: outlet, meteorological profile from data
However, it runs pretty well when
3) east: inlet, meteorological profile from data
west: inlet, meteorological profile from data
I do not understand why... The user guide advocates to set the (1) but it is the one wich crash the earlier....
** With User meteo file provided (the parameters are set really close to the automatic profile)
the simulation crash after 15 iterations or less, whatever my efforts are. It seems to be really unstable.
I have attached two screen : the first one is the velocity with automatic meteo file, and boundaries case 3). With the exact same parameters, the other screen is the User defined meteo file.
So I have few questions :
a) Am I doing all wrong with the inlet boundaries ? and if so, what is advocated, why is it crashing faster withe userguide recomandations ?
b) What can I do to stabilize a simulation without structured mesh (I have already tried inflations, forcing upwind scheme and blending factor).
c) Is my meteo file wrong ?
d) Why does the custom meteo file seems to make the simulation unstable regardless of the setup ?
e) Is it normal that my simulation is much much slower with meteo userfile ?
Thank you very much for the help, it will be appreciated in any case
Please find attached the meteo file, the two screens, the mesh quality report from salome, and the two similar xml.
PS : Don't mind the time setting, I have tried everything, and the constant is currently the more robust, but I tried with time variative. I tried also with Simplec, but inner iteration seems to work well too...