I want to export "joule heating" from Elmer to OpenFOAM via EOF-Library. As discussed in the Elmer tutorial, for joule heating https://github.com/ElmerCSC/elmerfem/bl ... M/case.sif, I would like to employ joule heating (qjoule) as the source term in TEqn.H file (https://github.com/anilkunwar/OpenFOAM_ ... oam/TEqn.H). So my basic query is:
1. What is the unit of joule heating as target variable in Elmer? h_q is in W/kg and S_qjoule = \rho* h is in W/m^3.
2. As I have written it as
Code: Select all
...
fvScalarMatrix TEqn
(
fvm::ddt(T)
+ fvm::div(phi, T)
- fvm::laplacian(alphaEff, T)
== qjoule*alphaEff/kappa
);
...
I suppose the unit of qjoule of openfoam and "joule heating" from Elmer have the same dimensions if the right hand side (RHS) of equation is written as in the code. In RHS of the equation, kappa is thermal conductivity (W/m K), alphaEff is thermal diffusivity (m^2/S) and qjoule is joule heating (W/m^3).
Can anyone verify whether this is correct or not?
Note: it is assumed that the temperature dependence of sigma (electrical conductivity) is addressed during Elmer simulation and thus qjoule is solved as source term not dependent on T variable.
References:
1. https://github.com/ElmerCSC/elmerfem/tr ... ToOpenFOAM
2. https://eof-library.com/
Yours Sincerely,
Anil Kunwar