Getting nothing but 0's for temperature grad data

Numerical methods and mathematical models of Elmer
Post Reply
DrHu
Posts: 6
Joined: 29 May 2024, 19:13
Antispam: Yes

Getting nothing but 0's for temperature grad data

Post by DrHu »

Hoping someone might see why I am getting nothing but zeros for the Temperature gradient data in the vtk file. I'm trying to follow the example in the Solver documentation for computing heat flux and the gradients (although I just want the gradients). The solver sections in the .sif file are listed below. The temperature data in the vtk output using Paraview looks fine (same as w/o the fluxsolver). And it appears that gradient data is being ouyput, but it's all 0, and there is obvious temperature variation.I've tried a few variations, but without any luck.
--
Solver 2
Equation = Heat Equation
Variable = Temperature
Procedure = "HeatSolve" "HeatSolver"
Exec Solver = Always
Stabilize = True
Optimize Bandwidth = True
Steady State Convergence Tolerance = 1.0e-5
Nonlinear System Convergence Tolerance = 1.0e-7
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 = Iterative
Linear System Iterative Method = BiCGStab
Linear System Max Iterations = 500
Linear System Convergence Tolerance = 1.0e-10
BiCGstabl polynomial degree = 2
Linear System Preconditioning = ILU0
Linear System ILUT Tolerance = 1.0e-3
Linear System Abort Not Converged = False
Linear System Residual Output = 10
Linear System Precondition Recompute = 1
End

Solver 1
Exec Solver = After all
Equation = "flux compute"
Output Format = Vtu
Procedure = "FluxSolver" "FluxSolver"
Calculate grad =True
Target Variable = Temperature
Steady State Convergence Tolerance = 1.0e-5
Linear System Solver = Iterative
Linear System Solver = "Iterative"
Linear System Iterative Method = "cg"
Linear System Preconditioning = ILU0
Linear System Residual Output = 10
Linear System Max Iterations = Integer 500
Linear System Convergence Tolerance = 1.0e-10
End

Equation 1
Name = "Equation 1"
Active Solvers(2) = 2 1
End
---
DH
Rich_B
Posts: 432
Joined: 24 Aug 2009, 20:18

Re: Getting nothing but 0's for temperature grad data

Post by Rich_B »

Without posting a minimum working example, with geometry and meshing instructions, it's hard to say what is going on with the simulation.

Searching through the consistency tests for a sif file containing the phrase 'calculate grad', one test in particular looked appropriate. Running the test 'fluxsolver3', shows the following results for heat gradient:
grad.png
grad.png (11.68 KiB) Viewed 201 times
kevinarden
Posts: 2428
Joined: 25 Jan 2019, 01:28
Antispam: Yes

Re: Getting nothing but 0's for temperature grad data

Post by kevinarden »

I would say since the flux solver is solver 1 and the heat solver is solver 2, the flux solver is going first and the gradient is zero because no heat or temperature is calculate yet. The heat solver goes second and calculates the heat and temperature, but the flux solver was already done, when everything was zero. Either redo the order, or repeat the flux solver after the heat solver goes.
DrHu
Posts: 6
Joined: 29 May 2024, 19:13
Antispam: Yes

Re: Getting nothing but 0's for temperature grad data

Post by DrHu »

NEVERMIND...
Putting

Exec Solver = before saving !(instead of = After all)

in Solver 1, populated the gradients. Not sure why, but so much for that issue.
Post Reply