Search found 46 matches

by C0st4s
Mon Jul 06, 2020 12:19 pm
Forum: code_saturne usage
Topic: Cell_Gradient_Upwind
Replies: 3
Views: 628

Re: Cell_Gradient_Upwind

Dear Yvan, Many thanks for the reply. I saw that function, but I thought you have to have ischcp == 2 (UPWIND) for all the convective gradients in the domain. Two questions then: 1) Correct me if I am wrong, but the src/alge/convection_diffusion.c upwind_gradient function is it local for the .c file...
by C0st4s
Sun Jul 05, 2020 5:17 pm
Forum: code_saturne usage
Topic: Cell_Gradient_Upwind
Replies: 3
Views: 628

Cell_Gradient_Upwind

Hi, I know that if I want to get the gradient of a vector or scalar field in a cell volume, there is the predefined 'field_gradient_scalar' function which returns based on the least square or iterative process etc.. However, say that I want to compute the dU/dx across a cell but instead to use upwin...
by C0st4s
Thu Mar 26, 2020 7:33 pm
Forum: code_saturne usage
Topic: EpsilonBoundaryCondition
Replies: 2
Views: 838

Re: EpsilonBoundaryCondition

Hi Yvan,

Thanks for the reply,
Costas
by C0st4s
Mon Mar 23, 2020 11:35 am
Forum: code_saturne usage
Topic: EpsilonBoundaryCondition
Replies: 2
Views: 838

EpsilonBoundaryCondition

Hi, I am trying to impose an epsilon boundary condition to the near wall node I that is equal to the one proposed by Chien and Launder (1980) paper. The approach is to tweak the discretised equation of the near wall node by using a big source term and an imposed value equivalent to the equilibrium l...
by C0st4s
Sat Sep 07, 2019 10:23 pm
Forum: code_saturne usage
Topic: Backtrace_SIGSEGV signal
Replies: 2
Views: 648

Re: Backtrace_SIGSEGV signal

Dear Yvan,

Thanks a lot for the quick reply. I had a division by zero


Best regards,
Costas
by C0st4s
Fri Sep 06, 2019 2:07 pm
Forum: code_saturne usage
Topic: Backtrace_SIGSEGV signal
Replies: 2
Views: 648

Backtrace_SIGSEGV signal

Hi,

Hope you are well.
What is the easiest way to trace the array that is causing SIGSEGV signal (forbidden memory area access) intercepted! in cs? I am writing some module routines and I am trying to find the root of the error.

Best regards,
Costas
by C0st4s
Mon Aug 19, 2019 5:34 pm
Forum: code_saturne usage
Topic: Identify the type and get the value of a boundary
Replies: 6
Views: 1329

Re: Identify the type and get the value of a boundary

Hi, Thanks Yvan again! I use the approach of defining the boundary conditions in cs_user_boundary_conditions.f90 and calling that routine in my module. It works fine. Where can I find an example of the coefa and coefb approach that you suggested just to have a look? Is it similar to the one in the e...
by C0st4s
Wed Aug 14, 2019 11:05 am
Forum: code_saturne usage
Topic: Identify the type and get the value of a boundary
Replies: 6
Views: 1329

Re: Identify the type and get the value of a boundary

Dear Yvan, I call the routines of the cs_user_modules in cs_user_source_terms.f90 under cs_user_turbulence_source_terms since I modifying sources in the turbulent transport equations. The import thing is to get the face index of each boundary wall face which I am not sure if can do it with rcodcl(if...
by C0st4s
Sun Aug 11, 2019 5:55 pm
Forum: code_saturne usage
Topic: Identify the type and get the value of a boundary
Replies: 6
Views: 1329

Re: Identify the type and get the value of a boundary

Dear Yvan,

Really appreciate your reply. Though, the main issue that I face is that the boundary conditions are prescribed in the GUI (smooth wall with heat flux) and I want to call the value and the type of that in the cs_user_modules.f90.

Best regards,
Costas
by C0st4s
Fri Aug 09, 2019 12:43 pm
Forum: code_saturne usage
Topic: Identify the type and get the value of a boundary
Replies: 6
Views: 1329

Identify the type and get the value of a boundary

Hi, I have a simple question. Is it possible to get the type and the value of an assigned boundary condition inside the cs_user_modules.f90? I tried: call user_boundary_conditions(nvar, itypfb, icodcl, rcodcl) write(*,*) itypfb(ifac), icodcl(ifac, ivart), rcodcl(ifac, ivart,3) !,cvar_scalt(ifac) And...