code_saturne User's Forum

Hello Luciano,

Another option to stop any momentum in the solid is to impose a high negative implicit source term for the momentum conservation equation as it is described in [Patankar, S., 1980. Numerical Heat Transfer and Fluid Flow. Hemisphere Publishing Corporation, Washington, D.C.]
The part ...

Hello Yvan,

Thanks a lot for your explanations. It is clear now. I will try to implement all these new things and get back to this thread to tell what I got.

Best regards,
Rodion

Update:
Okay, I implemented the equation for effective field in the model. In fact it didn't work as intended. It ...

Dear Yvan,

Thanks a lot for the answer. It is surprising how much I did not know about the gradient computation in Code_Saturne, even though I looked into it for a while.

In the source code I can see that cgdcel is called by gradient_s in case without potential-type values and without weighting ...

Hello Yvan,

Sorry, I meant cs_gradient_scalar or _gradient_scalar . I guess I have to use cs_gradient_scalar , right?

Also if I assume the homogeneous Neumann boundary condition, does it mean that I have to assign bc_coeff_a=0 and bc_coeff_b=1?

Furthermore, is there any way to call cs_gradient ...

Hello Yvan,

Thanks a lot for your answer.
The low-level gradient function that you are talking about, is it cs_field_gradient_scalar or cs_gradient_scalar or cs_gradient_scalar ?

I think in my case the homogeneous Neumann boundary condition should be alright. It is the zero flux boundary condition ...

Hello Yvan,

Thanks a lot for your comment.
Unfortunately I am still using Code_Saturne 5.0.9 on Windows. The transition to the version 6.0 or 6.1 and to Linux will take some more time, since I have to wait for the end of the lock-down/confinement in France to create a WSL machine on the ...

Hello Erwan,

Thank you very much for your response.

The thing is that the factor (σ*k/e)*Q(he) is known for each cell and it depends on the enthalpy (or temeprature, does not matter). But this factor does not depend on the electric potential φ, since it is a parameter of plasma and not of the ...

Hello Erwan,

Thank you very much for your response.

Basically I am trying to add the effect of electron pressure on the electric field. The current form of the equation for electric potential φ is:
div(σ∇φ)=0 (1)
Or if we write it in terms of electric field div(σ E )=0

I am trying to replace the ...

Hello everyone,

I work with simulation of electric arc plasma torches and currently I am trying to improve the theory of the electric arc model of Code_Saturne,
For that I developed a two-temperature model (one temperature fore electrons and one temperature for heavy species) of electric arc in ...

Hello evereone,

I am developing a two temperature model of electric arc in Code_Saturne. In this model I have two thermal scalars: electron enthalpy and heavy species enthalpy. These two scalars have diffferent diffusivity and different source terms. The source term for electrons is Joule power and ...