Skip to content
Snippets Groups Projects
Commit fd1779e6 authored by Vinayak Gholap's avatar Vinayak Gholap
Browse files

Added analytically calculated EM thermal source for debugging purpose,...

Added analytically calculated EM thermal source for debugging purpose, compairing strong vs weak results in SMP
parent 21ed346a
No related branches found
No related tags found
1 merge request!386Major update: Computing average IPVar over time period, rewriting last time...
...@@ -260,6 +260,8 @@ void mlawLinearElecMagTherMech::constitutive( ...@@ -260,6 +260,8 @@ void mlawLinearElecMagTherMech::constitutive(
// here magnetic, electric field and temp dependent // here magnetic, electric field and temp dependent
// if true then we solve EM problem and thus reevaluate Magnetic potential as well
// else do not consider Magnetic contributions -> TM problem
if (this->_evaluateCurlField) if (this->_evaluateCurlField)
{ {
STensorOperation::zero(H); STensorOperation::zero(H);
...@@ -400,6 +402,21 @@ void mlawLinearElecMagTherMech::constitutive( ...@@ -400,6 +402,21 @@ void mlawLinearElecMagTherMech::constitutive(
// w_AV = j_e \cdot -dadt + h \cdot dbdt // w_AV = j_e \cdot -dadt + h \cdot dbdt
ThermalEMFieldSource += (sourceVectorField(i) * ((An(i) - A0(i)) * invdh) /*+ H(i) * (Bn(i) - B0(i)) * invdh*/); ThermalEMFieldSource += (sourceVectorField(i) * ((An(i) - A0(i)) * invdh) /*+ H(i) * (Bn(i) - B0(i)) * invdh*/);
// to debug weak vs strong coupling results by
// using analytical EM thermal source in SMP only with
// A²/2 being mean value that is calculated
// in weak coupling for considered parameter values
/*
if(this->_num == 1) // if SMP
// A² cos²(2 pi f t)
// A²/2 = Mean value over one period = 3559017.23076498
ThermalEMFieldSource = 2.0 * 2227580. * std::cos(2.0 * M_PI * _freqEM * ctime) *
std::cos(2.0 * M_PI * _freqEM * ctime);
else
ThermalEMFieldSource = 0.0;
*/
  • Author Developer

    Note: I will need to update the value used depending on f, I used.

    Use when comparing the resulting thermal field in SMP with EM thermal source computed analytically in the strong coupling test example (only SMP domain with no resolution of EM fields) with EM value (A²/2) imported from weak coupling test.

  • Please register or sign in to reply
// w_AV += je \cdot -gradV // w_AV += je \cdot -gradV
if (_useFluxT) if (_useFluxT)
ThermalEMFieldSource += (sourceVectorField(i) * gradV(i)); ThermalEMFieldSource += (sourceVectorField(i) * gradV(i));
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment