How can I modify the equation for electric potential?
Posted: Tue Mar 31, 2020 3:05 am
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 Code_Saturne. Now in my two-temperature electric arc model I need to introduce the effective field instead of the common electric field in the equation of electric current continuity. Right now Code_Saturne solves the following equation in terms of electric potential:
div(σ∇φ)=0
Basically this equation is the diffusive term being equal zero. It is given in the section "cs elec model routine" of Code_Saturne theory guide. The further improvement of the model that I am trying to introduce consists in correcting the previous equation to the following one:
div(A(x,y,z)+σ∇φ)=0.
Does anyone have any idea how can I implement this correction? How can I add an extra term inside the divergence of the diffusive term?
Best regards,
Rodion
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 Code_Saturne. Now in my two-temperature electric arc model I need to introduce the effective field instead of the common electric field in the equation of electric current continuity. Right now Code_Saturne solves the following equation in terms of electric potential:
div(σ∇φ)=0
Basically this equation is the diffusive term being equal zero. It is given in the section "cs elec model routine" of Code_Saturne theory guide. The further improvement of the model that I am trying to introduce consists in correcting the previous equation to the following one:
div(A(x,y,z)+σ∇φ)=0.
Does anyone have any idea how can I implement this correction? How can I add an extra term inside the divergence of the diffusive term?
Best regards,
Rodion