Hello,
Are you following the tutorial? The expected results are shown at the end of the pdf.
Could you post the listing file.
Regards,
Luciano
Search found 284 matches
- Wed Nov 09, 2016 4:45 am
- Forum: code_saturne usage
- Topic: Heated square cavity flow Tutorial 3 in Code Saturne 4.0
- Replies: 11
- Views: 11019
- Mon Nov 07, 2016 3:06 pm
- Forum: code_saturne usage
- Topic: free surface not moving
- Replies: 21
- Views: 15502
Re: free surface not moving
Hello Olivier, I was debuging the alelav.f90 subroutine and I found that may be there is a bug in the following line, because to get the type of boundary condition at the boundary face ifac is used itypfb(ifac). if (ialtyb(ifac) .eq. ifresf) then so if you change to if (itypfb(ifac).eq.ifresf) then ...
- Thu Nov 03, 2016 5:55 pm
- Forum: code_saturne usage
- Topic: Heated square cavity flow Tutorial 3 in Code Saturne 4.0
- Replies: 11
- Views: 11019
Re: Heated square cavity flow Tutorial 3 in Code Saturne 4.0
Hello,
CFL -> Courant–Friedrichs–Lewy (CFL) condition.
https://en.wikipedia.org/wiki/Courant%E ... _condition
Regards,
Luciano
CFL -> Courant–Friedrichs–Lewy (CFL) condition.
https://en.wikipedia.org/wiki/Courant%E ... _condition
Regards,
Luciano
- Thu Nov 03, 2016 4:01 pm
- Forum: code_saturne usage
- Topic: Heated square cavity flow Tutorial 3 in Code Saturne 4.0
- Replies: 11
- Views: 11019
Re: Heated square cavity flow Tutorial 3 in Code Saturne 4.0
Hello, You need a temperature reference state for the Boussinesq approximation. The tutorial mention that "The reference state is taken at the cold wall temperature", so the Tref sould be 293.15K. For convergence checking as mentioned in the tutorial you can " Verify that the residual...
- Thu Nov 03, 2016 1:31 pm
- Forum: code_saturne usage
- Topic: free surface not moving
- Replies: 21
- Views: 15502
Re: free surface not moving
Hello Olivier, Please, next time post the updated subroutine Whe you mentioned that you "digging" the alelav.f90, did you debug or printed the values of pimpv(i) for any face??. Because I debugging (printing) the data in alelav.f90 and I don't get any face that satisfied the condition if (...
- Thu Nov 03, 2016 1:20 pm
- Forum: code_saturne usage
- Topic: Heated square cavity flow Tutorial 3 in Code Saturne 4.0
- Replies: 11
- Views: 11019
Re: Heated square cavity flow Tutorial 3 in Code Saturne 4.0
Hello, I was checking you case and the problem it is with the value of the air expasion coeff that you compute. With the data that you defined you get a \beta=3.3224e-12 , so something in your formula is wrong. If you set \beta=0.0033 the flow start moving. Also, if you look at the listing file all ...
- Tue Nov 01, 2016 10:34 pm
- Forum: code_saturne usage
- Topic: free surface not moving
- Replies: 21
- Views: 15502
Re: free surface not moving
Hello olivier, I was looking at your problem. First of all, I couldn't run but I found some issues in the user subroutine. 1) You have a if (.false.) then before the boundary selection. 2) call getfbr('5', nlelt, lstelt) any of your baundaries has this ('5') indentifier. 3) The axis of your mesh are...
- Tue Nov 01, 2016 10:15 pm
- Forum: code_saturne usage
- Topic: Heated square cavity flow Tutorial 3 in Code Saturne 4.0
- Replies: 11
- Views: 11019
Re: Heated square cavity flow Tutorial 3 in Code Saturne 4.0
Hello In CS 4.0.5 you can use the Mathematical Expression Editor for a variable density \rho(T) and write the Boussinesq Model with the correspondig thermal expansion coefficient of the fluid that you want to simulate. The following is an example: #air expansion coefficient (1/K) b=0.00343; #Ref tem...
- Tue Nov 01, 2016 4:10 pm
- Forum: Discussion
- Topic: Hexa dominant mesh
- Replies: 3
- Views: 6895
Re: Hexa dominant mesh
Hello Yvan, I have compile the converter and it works. I can convert snappy meshes (with and without boundary layers) to CS, but no all the genereted meshes are valid for CS. Some quality parameters in SHM has to be tunned in order to obtain a valid mesh, but once you get the appropiated parameters ...
- Fri Oct 28, 2016 7:04 pm
- Forum: Discussion
- Topic: Hexa dominant mesh
- Replies: 3
- Views: 6895
Re: Hexa dominant mesh
Hello Yvan, I had seen this converter but the last version has 2 years and I thought it wouldn't work, so I write a python script to convert VTKtoMED (polyhedron meshes) but it is not as robust as I want, for some meshes appears negatives volume. I will try to compile the converter and comment the r...