Skip to content

TPSA restructure (3/3): Newton defaults for linear elasticity - #7359

Draft
svenn-t wants to merge 3 commits into
OPM:masterfrom
svenn-t:tpsa_newton_restructure
Draft

TPSA restructure (3/3): Newton defaults for linear elasticity#7359
svenn-t wants to merge 3 commits into
OPM:masterfrom
svenn-t:tpsa_newton_restructure

Conversation

@svenn-t

@svenn-t svenn-t commented Aug 25, 2026

Copy link
Copy Markdown
Contributor

Change Newton to default to 1 iteration, to speed up linear elasticity solves with TPSA. Now we only do one linearization, check the error, and proceed with one linear solve if linearization error is too high. Increased linear solver tolerance to compensate for only doing 1 Newton iteration.

Also includes some minor changes and fixes:

  • Rescaling of matrix and residual; default is scaling system with average shear modulus
  • Fix for calculating average shear modulus in parallel
  • Update info for geomechanics output at endTimeStep()

Includes both PRs #7357 and #7358, which should be merged first. Draft mode until then.

- Old 7x7 block structure is broken up to several sub-matrices, purpose-made for preconditioners like Hypre.
- New structure made to interface 7x7 structure in linearizer without changes.
- One preconditioner for whole TPSA system; three sub-block preconditioner can be set by user
- Allows for Hypre, which is default if linked; else ilu0
- newton procedure now: check linearization, if error larger than tolerance, solve linear system
- lowered linear solve tolerance
- more clear fixed-stress conv. check
- bug fix: num. dofs on owned cells only
- moved calc. of info to mech. output to separate function, which is called at endTimeStep()
- option to scale linear system for better conv. check with dune solvers (preconditioners are invariant)
@svenn-t svenn-t added the manual:enhancement This is an enhancement/improvent that needs to be documented in the manual label Aug 25, 2026
@svenn-t

svenn-t commented Aug 25, 2026

Copy link
Copy Markdown
Contributor Author

jenkins build this please

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Labels

manual:enhancement This is an enhancement/improvent that needs to be documented in the manual

Projects

None yet

Development

Successfully merging this pull request may close these issues.

1 participant