Search found 235 matches

by Luciano Garelli
Fri May 31, 2019 1:15 pm
Forum: Code_Saturne usage
Topic: VOF and ALE compatible?
Replies: 3
Views: 577

Re: VOF and ALE compatible?

Hello,

VOF and ALE can be used simultaneously. I did some test with a wave generator.

Be careful with the version that you are using becouse were a bug fixed in the ALE Bc (commit 4b879fa).

Regards,

Luciano
by Luciano Garelli
Wed May 22, 2019 3:53 pm
Forum: Code_Saturne usage
Topic: FSI internal coupling and sub-iterations
Replies: 2
Views: 479

Re: FSI internal coupling and sub-iterations

Hello, The coefficient cfopre is used to predict the forces at time (n+1), not to relax them. In modini.f90 is set cfopre = 2.0d0, then in strini.f90 is modified in case of nalimx>1 else if (nbstru.gt.0) then if (nalimx.eq.1) then write(nfecra,2020) aexxst, bexxst, cfopre else cfopre = 1.d0 write(nf...
by Luciano Garelli
Mon May 20, 2019 5:32 pm
Forum: Code_Saturne usage
Topic: Newmark-HHT algorithm in Code_Saturne
Replies: 9
Views: 795

Re: Newmark-HHT algorithm in Code_Saturne

Hello, I don't known what kind of FSI problem are you solving. The time integration is implemented like in the paper "Improved Numerical Dissipation for Time Integration Algorithms in Structural Dynamics, Hilber, H.M, Hughes,T.J.R and Talor, R.L." . The change in the parameters will modify the behav...
by Luciano Garelli
Mon May 20, 2019 4:27 pm
Forum: Code_Saturne usage
Topic: Newmark-HHT algorithm in Code_Saturne
Replies: 9
Views: 795

Re: Newmark-HHT algorithm in Code_Saturne

Hello,

To what specific set of values do you refer? alpnmk, betnmk, gamnmk?

Regards,

Luciano
by Luciano Garelli
Mon May 20, 2019 3:09 pm
Forum: Code_Saturne usage
Topic: Newmark-HHT algorithm in Code_Saturne
Replies: 9
Views: 795

Re: Newmark-HHT algorithm in Code_Saturne

The parameter alpnmk is initialized in iniini.f90 and the paremeters betnmk and gamnmk a defined in modini.f90. I think that you can copy the newmrk.f90 or strdep.f90 to the SRC directory and change the values. I don't known if the values of the alpnmk, betnmk, gamnmk can be accessed from another us...
by Luciano Garelli
Mon May 20, 2019 12:35 pm
Forum: Code_Saturne usage
Topic: Newmark-HHT algorithm in Code_Saturne
Replies: 9
Views: 795

Re: Newmark-HHT algorithm in Code_Saturne

Hello,

I don't understand very well your question, but did you check the file newmrk.f90. In this subroutine are computed the accelerations, velocities and displacements of the structure.

Regards,

Luciano
by Luciano Garelli
Thu May 16, 2019 12:47 pm
Forum: Code_Saturne usage
Topic: Boundary condition definition error
Replies: 2
Views: 455

Re: Boundary condition definition error

Hello Roberto, This error appears when some boundary face doesn't have a boundary condition imposed. You can try to select all the boundary faces that weren't selected. To do this, add a new boundary condition in the "definition of boundary region" and in the selection criteria will be appear a sele...
by Luciano Garelli
Tue May 14, 2019 2:09 pm
Forum: Code_Saturne usage
Topic: subroutine problem
Replies: 4
Views: 609

Re: subroutine problem

Hello, Did you check the documentation for coding of non-standard boundary conditions in the user guide? In your setup you start from the position #16 of the face list. Why ? do ilelt = 16, nlelt and when you set the temperature at the face you set the #16 also... rcodcl(ifac,isca(iscalt), 16 ) = 20...
by Luciano Garelli
Thu Apr 25, 2019 1:41 pm
Forum: Code_Saturne usage
Topic: Meshing warning
Replies: 4
Views: 515

Re: Meshing warning

Hello, I attach an example of the flow past an sphere. I use the geometry that you shared and I have created a new mesh. This case with a coarse mesh (aprox 115.000 cellls) took arround 13 min in a laptop to solve 500 time steps with Courant max around 7. I added the listing file in the zip. I hope ...
by Luciano Garelli
Mon Apr 22, 2019 1:34 pm
Forum: Code_Saturne usage
Topic: Meshing warning
Replies: 4
Views: 515

Re: Meshing warning

Hello, I use you mesh and the case run without diverging. The mesh is very coarse and the check quality criteria report some issues. I think that you have to reduce the time step in order to get a Courant ~ 5 or lower. You can read some guidelines at https://www.code-saturne.org/cms/documentation/BP...