Heating Power by Body ID

Post processing utility for Elmer
Post Reply
Edmund
Posts: 9
Joined: 10 May 2013, 11:28
Antispam: Yes

Heating Power by Body ID

Post by Edmund »

Hi,

running the StatCurrentSolver the Total heating power is calculated by default.
I want to put out this power separated by body-ID.

How to do this?

Kind Regards.

Edmund
Edmund
Posts: 9
Joined: 10 May 2013, 11:28
Antispam: Yes

Re: Heating Power by Body ID

Post by Edmund »

I tried savescalars (Operator 1 = int Variable 1 = joule heating) and Target Variable 1 = String Powercalc by executing the savescalar solver in an extra equation.
Without success. The output is still the overall Power calculated in a different way for all Bodies in an extra variable.

How can I tell the savescalars solver to compute it only in a body of choice?

Code: Select all

Body 1
  Target Bodies(1) = 3
  Name = "mat2"
  Equation = 1
  Material = 1
End

Body 2
  Target Bodies(1) = 4
  Name = "mat2"
  Equation = 2
  Material = 2
End

Solver 1
  Equation = Static Current Conduction
  Calculate Volume Current = True
  Variable = -dofs 1 Potential
  Procedure = "StatCurrentSolve" "StatCurrentSolver"
  Exec Solver = Always
  Stabilize = True
  Bubbles = False
  Calculate Joule Heating = True
  Lumped Mass Matrix = False
  Optimize Bandwidth = True
  Steady State Convergence Tolerance = 1.0e-6
  Nonlinear System Convergence Tolerance = 1.0e-9
  Nonlinear System Max Iterations = 20
  Nonlinear System Newton After Iterations = 3
  Nonlinear System Newton After Tolerance = 1.0e-3
  Nonlinear System Relaxation Factor = 1
  Linear System Solver = Direct
!  Linear System Iterative Method = Mumps
  Linear System Iterative Method = Umfpack
!  Linear System Solver = Iterative
!  Linear System Iterative Method = BiCGStab
  Linear System Max Iterations = 1000
  Linear System Convergence Tolerance = 1.0e-10
  Linear System Preconditioning = ILU0
  Linear System ILUT Tolerance = 1.0e-3
  Linear System Abort Not Converged = True
  Linear System Residual Output = 1
  Linear System Precondition Recompute = 1
End

Solver 2
Exec Solver = String "after timestep"
exec interval =  1
Equation = String "ResultOutput"
Procedure = File "ResultOutputSolve" "ResultOutputSolver"
Output File Name = String "paraver." !or any other output file name of your choice
Save Geometry Ids = Logical True
Output Format = String "vtu"
Vtu Format = Logical True
Binary Output = Logical True !binaryformatisthedefault
End

Solver 3
  Equation = SaveScalars
  Operator 1 = int
  Variable 1 = joule heating
  Procedure = "SaveData" "SaveScalars"
  Target Variable 1 = String Powercalc
  Exec Solver = Always
  Filename = File "scalar_data.dat"
!  variable 2 = Powercalc
End

Equation 1
  Name = "Equation 1"
  Active Solvers(2) = 1 2
End

Equation 2
  Name = "Equation 2"
  Active Solvers(3) = 1 2 3
End


raback
Site Admin
Posts: 4812
Joined: 22 Aug 2009, 11:57
Antispam: Yes
Location: Espoo, Finland
Contact:

Re: Heating Power by Body ID

Post by raback »

Hi

Perhaps you could introduce in Material section

Code: Select all

MyCoeff = Real 1.0  ! or 0.0 in some other material
and the give for SaveScalars

Code: Select all

Coefficient 1 = String MyCoeff
Not that elegant but could work.

Even better, if you compute the "Joule Field" then your coefficient could be the real conductivity where you want it and you would not get the errors introduced by the non-continuous.

Now this is just out of memory, little bit busy with some deadlines...

-Peter
Edmund
Posts: 9
Joined: 10 May 2013, 11:28
Antispam: Yes

Re: Heating Power by Body ID

Post by Edmund »

Yes it works.
Thanks again.
Post Reply