Skip to content

feat: Add single phase viscosity dependency on temperature#4029

Open
dkachuma wants to merge 3 commits intodevelopfrom
dkachuma/viscosity-thermal-expansion
Open

feat: Add single phase viscosity dependency on temperature#4029
dkachuma wants to merge 3 commits intodevelopfrom
dkachuma/viscosity-thermal-expansion

Conversation

@dkachuma
Copy link
Copy Markdown
Contributor

Adds a viscosity thermal expansion coefficient for single phase fluid. This introduces a simple dependency of the viscosity on temperature.

The documentation has been updated to reflect this addition.

Replaces #3951

@dkachuma dkachuma self-assigned this Apr 14, 2026
@dkachuma dkachuma added type: feature New feature or request ci: run CUDA builds Allows to triggers (costly) CUDA jobs flag: ready for review flag: requires rebaseline Requires rebaseline branch in integratedTests ci: run integrated tests Allows to run the integrated tests in GEOS CI ci: run code coverage enables running of the code coverage CI jobs labels Apr 14, 2026
Comment on lines +12 to +13
For thermal simulations, the model density and viscosity have a temperature dependency
expressed as a constant thermal expansion coefficient.
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
For thermal simulations, the model density and viscosity have a temperature dependency
expressed as a constant thermal expansion coefficient.
For thermal simulations, fluid density and viscosity are temperature-dependent, governed by a constant thermal expansion coefficient.

for isothermal cases and,

.. math::
\rho(p,T) = \rho_0 e^{c_\rho(p - p_0)} e^{-\beta_\rho (T - T_0)}
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
\rho(p,T) = \rho_0 e^{c_\rho(p - p_0)} e^{-\beta_\rho (T - T_0)}
\rho(p,T) = \rho_0 e^{c_\rho(p - p_0)} e^{ - \beta_\rho (T - T_0)}

TestDeck(
name="thermalViscosity_1d",
description='Thermal problem with temperature dependent viscosity',
partitions=((1, 1, 1)),
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
partitions=((1, 1, 1)),
partitions=((1, 1, 1),),

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 CUDA builds Allows to triggers (costly) CUDA jobs ci: run integrated tests Allows to run the integrated tests in GEOS CI flag: ready for review flag: requires rebaseline Requires rebaseline branch in integratedTests type: feature New feature or request

Projects

None yet

Development

Successfully merging this pull request may close these issues.

3 participants