Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Add plane to stasis constraint and remove pressure constraints #1193

Merged
merged 13 commits into from
Jul 15, 2024

Conversation

AmishgaAlphonius
Copy link
Collaborator

@AmishgaAlphonius AmishgaAlphonius commented Jul 11, 2024

Description

Following the implementation of temperature-dependent solid domain constraint for both single-phase and two-phase flows, a need for an additional condition was identified. Indeed, in some simulations, splatters may happen, and with the previous implementation, splatter could freeze "mid-air" which is non-physical and therefore undesired. This PR adds a plane to divide the problem domain in two parts and only the one in the opposite direction to the normal vector is considered for stasis constraints to be applied. The plane is defined through a point and an outward pointing normal vector.

During this PR, a bug regarding the pressure constraints was detected and fixed. Previously, the pressure DOFs were also constrained, but this resulted in an ill-posed problem when sources terms in the momentum equation were pressure-dependent. Therefore, as of this PR, the constraints on pressure DOFs were removed.

Testing

Previous tests with solid domain constraints were updated due to the pressure DOFs that are no more constrained.
Additionally, two new application tests were added for the new plane restriction feature (1 for single fluid flows and one for VOF).

Documentation

The "Constrain Stasis" page in the parameter documentation was updated with the new features and disabled ones.

Miscellaneous (will be removed when merged)

Checklist (will be removed when merged)

See this page for more information about the pull request process.

Code related list:

  • All in-code documentation related to this PR is up to date (Doxygen format)
  • Lethe documentation is up to date
  • New feature has unit test(s) (preferred) or application test(s), and restart files are in the generator folder
  • The branch is rebased onto master
  • Changelog (CHANGELOG.md) is up to date
  • Code is indented with indent-all and .prm files (examples and tests) with prm-indent

Pull request related list:

  • Labels are applied
  • There are at least 2 reviewers (or 1 if small feature) excluding the responsible for the merge
  • If this PR closes an issue or is related to a project, it is linked in the "Projects" or "Development" section
  • If the fix is temporary, an issue is opened
  • The PR description is cleaned and ready for merge

@AmishgaAlphonius AmishgaAlphonius added Bug Something isn't working Enhancement New feature or request Quick fix labels Jul 11, 2024
@AmishgaAlphonius AmishgaAlphonius self-assigned this Jul 11, 2024
@AmishgaAlphonius AmishgaAlphonius added Ready for review Needs more reviewers This pull request needs more review before a merge is possible labels Jul 11, 2024
@AmishgaAlphonius AmishgaAlphonius marked this pull request as ready for review July 11, 2024 20:09
@AmishgaAlphonius AmishgaAlphonius changed the title Add plane to stasis constraint Add plane to stasis constraint and remove pressure constraints Jul 12, 2024
Copy link
Collaborator

@oguevremont oguevremont left a comment

Choose a reason for hiding this comment

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

Good work!

include/core/parameters.h Outdated Show resolved Hide resolved
include/solvers/navier_stokes_base.h Outdated Show resolved Hide resolved
include/solvers/navier_stokes_base.h Show resolved Hide resolved
include/solvers/navier_stokes_base.h Outdated Show resolved Hide resolved
source/solvers/navier_stokes_base.cc Outdated Show resolved Hide resolved
source/solvers/navier_stokes_base.cc Outdated Show resolved Hide resolved
doc/source/parameters/cfd/constrain_stasis.rst Outdated Show resolved Hide resolved
doc/source/parameters/cfd/constrain_stasis.rst Outdated Show resolved Hide resolved
CHANGELOG.md Outdated Show resolved Hide resolved
Copy link
Contributor

@blaisb blaisb left a comment

Choose a reason for hiding this comment

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

Other than Olivier's comment, I am fine with this Pr.

include/solvers/navier_stokes_base.h Outdated Show resolved Hide resolved
include/solvers/navier_stokes_base.h Show resolved Hide resolved
@blaisb
Copy link
Contributor

blaisb commented Jul 14, 2024

can you rebase? Then I would merge afterwards

@blaisb blaisb added Reviewed and ready to merge and removed Ready for review Needs more reviewers This pull request needs more review before a merge is possible labels Jul 14, 2024
@blaisb blaisb requested a review from oguevremont July 14, 2024 18:08
@AmishgaAlphonius AmishgaAlphonius force-pushed the add-plane-to-stasis-constraint branch from 753557e to 6989b50 Compare July 15, 2024 02:20
@AmishgaAlphonius AmishgaAlphonius requested a review from hepap July 15, 2024 14:09
@blaisb
Copy link
Contributor

blaisb commented Jul 15, 2024

I'll wait until @oguevremont and @hepap accept the PR before merging

include/solvers/navier_stokes_base.h Outdated Show resolved Hide resolved
@AmishgaAlphonius AmishgaAlphonius force-pushed the add-plane-to-stasis-constraint branch from 52ef1b4 to 653027c Compare July 15, 2024 15:21
Copy link
Collaborator

@hepap hepap left a comment

Choose a reason for hiding this comment

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

Looks very good to me :) Can't wait to try it!

@blaisb blaisb merged commit c41f675 into master Jul 15, 2024
8 checks passed
@blaisb blaisb deleted the add-plane-to-stasis-constraint branch July 15, 2024 19:13
M-Badri pushed a commit to M-Badri/lethe that referenced this pull request Sep 29, 2024
…-polymtl#1193)

Description
Following the implementation of temperature-dependent solid domain constraint for both single-phase and two-phase flows, a need for an additional condition was identified. Indeed, in some simulations, splatters may happen, and with the previous implementation, splatter could freeze "mid-air" which is non-physical and therefore undesired. This PR adds a plane to divide the problem domain in two parts and only the one in the opposite direction to the normal vector is considered for stasis constraints to be applied. The plane is defined through a point and an outward pointing normal vector.

During this PR, a bug regarding the pressure constraints was detected and fixed. Previously, the pressure DOFs were also constrained, but this resulted in an ill-posed problem when sources terms in the momentum equation were pressure-dependent. Therefore, as of this PR, the constraints on pressure DOFs were removed.

Testing
Previous tests with solid domain constraints were updated due to the pressure DOFs that are no more constrained.
Additionally, two new application tests were added for the new plane restriction feature (1 for single fluid flows and one for VOF).

Documentation
The "Constrain Stasis" page in the parameter documentation was updated with the new features and disabled ones.



Former-commit-id: c41f675
blaisb pushed a commit that referenced this pull request Sep 30, 2024
Description
Following the implementation of temperature-dependent solid domain constraint for both single-phase and two-phase flows, a need for an additional condition was identified. Indeed, in some simulations, splatters may happen, and with the previous implementation, splatter could freeze "mid-air" which is non-physical and therefore undesired. This PR adds a plane to divide the problem domain in two parts and only the one in the opposite direction to the normal vector is considered for stasis constraints to be applied. The plane is defined through a point and an outward pointing normal vector.

During this PR, a bug regarding the pressure constraints was detected and fixed. Previously, the pressure DOFs were also constrained, but this resulted in an ill-posed problem when sources terms in the momentum equation were pressure-dependent. Therefore, as of this PR, the constraints on pressure DOFs were removed.

Testing
Previous tests with solid domain constraints were updated due to the pressure DOFs that are no more constrained.
Additionally, two new application tests were added for the new plane restriction feature (1 for single fluid flows and one for VOF).

Documentation
The "Constrain Stasis" page in the parameter documentation was updated with the new features and disabled ones.



Former-commit-id: c41f675
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
Projects
None yet
Development

Successfully merging this pull request may close these issues.

4 participants