When DFLUX is applied to one C3D8 element in transient thermal calculations, strange things happen…
I provide the example file with just one element:
*ELEMENT, TYPE=C3D8, ELSET=EALL
*SOLID SECTION, ELSET=EALL, MATERIAL=MID1
*INITIAL CONDITIONS, TYPE=TEMPERATURE
*HEAT TRANSFER, SOLVER=ITERATIVE SCALING, DIRECT
Here the bottom face of the element is heated and the top face is heat insulated.
The results are following:
temperatures for set NALL and time 0.1000000E-04
1 3.999640E-03 2 3.999640E-03 3 3.999640E-03 4 3.999640E-03 5 -1.999640E-03 6 -1.999640E-03 7 -1.999640E-03 8 -1.999640E-03
heat generation for set NALL and time 0.1000000E-04
1 1.499820E-03 2 1.499820E-03 3 1.499820E-03 4 1.499820E-03 5 -1.499820E-03 6 -1.499820E-03 7 -1.499820E-03 8 -1.499820E-03
I didn’t expect negative temperatures at the top face. Also, heat generation values values at the insulated top face are negative. The magnitudes themselves don’t agree with the flux specified in DFLUX card.
For reference, I’ve run verification test oneel20df and the results are correct. I’ve played with different options, types and shapes of the element and apparently the cause is transient calculation.
Am I doing something wrong? Or maybe I misunderstood semantics of the DFLUX card.