Repository navigation
Use the element-average temperature for CTETRA4 thermal stresses - #345
Merged
victorkemp merged 1 commit intoOct 19, 2026
Merged
Conversation
The 4-node tetrahedron has a constant strain, so the only thermal strain consistent with it is the one from the element-average temperature. Since the solid thermal recovery change, CTETRA4 used the temperature at each recovery point instead, which put large spurious thermal stresses at the corners of any element with a temperature gradient across it. CTETRA4 recovery points now use the element-average temperature; CTETRA10 and the other solids keep the temperature at the point, which is consistent with their strain fields.
This file contains hidden or bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
Sign up for free
to join this conversation on GitHub.
Already have an account?
Sign in to comment
2 participants
Add this suggestion to a batch that can be applied as a single commit.This suggestion is invalid because no changes were made to the code.Suggestions cannot be applied while the pull request is closed.Suggestions cannot be applied while viewing a subset of changes.Only one suggestion per line can be applied in a batch.Add this suggestion to a batch that can be applied as a single commit.Applying suggestions on deleted lines is not supported.You must change the existing code in this line in order to create a valid suggestion.Outdated suggestions cannot be applied.This suggestion has been applied or marked resolved.Suggestions cannot be applied from pending reviews.Suggestions cannot be applied on multi-line comments.Suggestions cannot be applied while the pull request is queued to merge.Suggestion cannot be applied right now. Please check back later.
This undoes the previous change for only the CTETRA4, these are constant strain so can never represent any of the thermal strains induced by using the temp at the recovery point leading to erroneous corner stresses.
All other linear elements show some improvement, including the CPENTA6, so we keep those on the new path.