Skip to content

fix: SinglePhaseWell Thermal - #4126

Open
tjb-ltk wants to merge 8 commits into
developfrom
fix/byer3/singlephasewell
Open

fix: SinglePhaseWell Thermal#4126
tjb-ltk wants to merge 8 commits into
developfrom
fix/byer3/singlephasewell

Conversation

@tjb-ltk

@tjb-ltk tjb-ltk commented Aug 19, 2026

Copy link
Copy Markdown
Contributor
  1. fix numdof for shutwell matrix assembly , set to 2 for and should be 3
  2. add missing perforation rate contribution to well system only solve

Compositional wells were correct. The code was refactored a bit to avoid duplication in matrix stamping and Dof definition

…on contribution to well system when solving well only system . These are singlephase only issues
@tjb-ltk tjb-ltk self-assigned this Aug 19, 2026
@tjb-ltk tjb-ltk added ci: run CUDA builds ci: run integrated tests Allows to run the integrated tests in GEOS CI flag: no rebaseline Does not require rebaseline ci: run code coverage enables running of the code coverage CI jobs labels Aug 19, 2026
Comment thread src/coreComponents/physicsSolvers/fluidFlow/wells/SinglePhaseWell.cpp Outdated
@herve-gross
herve-gross requested a review from victorapm August 27, 2026 16:45
@herve-gross

Copy link
Copy Markdown
Contributor

@dkachuma can you take a look at this PR? Thanks.

@victorapm victorapm added ci: run device builds Allows to run the device (CUDA/HIP) builds and removed ci: run CUDA builds labels Aug 27, 2026

@dkachuma dkachuma left a comment

Copy link
Copy Markdown
Contributor

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Looks good

@@ -233,11 +233,6 @@ void CompositionalMultiphaseWell::registerWellDataOnMesh( WellElementSubRegion &
m_numComponents = fluid0.numFluidComponents();
}
// 1 pressure + NC compositions + 1 connectionRate + temp if thermal

Copy link
Copy Markdown
Contributor

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Suggested change
// 1 pressure + NC compositions + 1 connectionRate + temp if thermal

}
};

template< integer IS_THERMAL >

Copy link
Copy Markdown
Contributor

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Do we need this template argument? Won't it always be 1?

getReference< real64 >(
CompositionalMultiphaseWell::viewKeyStruct::currentTotalVolRateString() )=0.0;
getReference< real64 >( CompositionalMultiphaseWell::viewKeyStruct::currentMassRateString() )=0.0;
shutEntireWell( subRegion, dofManager, localMatrix, localRhs );

Copy link
Copy Markdown
Contributor

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Should we call resetShutInControlState after this? Looks like we were doing that before.

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

Labels

ci: run code coverage enables running of the code coverage CI jobs ci: run device builds Allows to run the device (CUDA/HIP) builds ci: run integrated tests Allows to run the integrated tests in GEOS CI flag: no rebaseline Does not require rebaseline flag: ready for review

Projects

None yet

Development

Successfully merging this pull request may close these issues.

5 participants