Page 1 of 1

The prolem with convective outlet bounary condition

Posted: Thu Nov 12, 2015 10:42 am
by MartianDuan
Dear developer and users.

I am using the Code_Saturne 4.0.2 to run the FSI simulation to study the effect of the movement of the valve head on the flow field. The movement of the valve is defined in the usalcl.f90. To avoid the reflection from the outlet boundary, I created a convective outlet boundary conditions, which is converted from a tutorial case of Code_Saturne 3.0.*.

However, the simualtion came across the error immediately after the convective functions being turned on. Can you please have a look at the case and give me some suggestions?

The source files of the simualtiona and the 'listing' are attached. Please have a look.

Kindest regards,

Yu.

Re: The prolem with convective outlet bounary condition

Posted: Fri Nov 13, 2015 10:55 am
by Yvan Fournier
Hello,

The crash happens in your subroutine from cs_user_extra_operations.f90.

It might be due to a division by zero or an uninitialized variable.

Running a debug build would tell you which line is incorrect. Using Valgrind could also be very helpful, as the case is small enough.

Otherwise, I recommend checking the values of all terms in your formulas separately.

Regards,

Yvan