Skip to content
This repository was archived by the owner on Sep 26, 2025. It is now read-only.

WIP: Add mixed form 2D Darcy problem example#82

Closed
alexanderianblair wants to merge 18 commits intomainfrom
alexanderianblair/darcy-example
Closed

WIP: Add mixed form 2D Darcy problem example#82
alexanderianblair wants to merge 18 commits intomainfrom
alexanderianblair/darcy-example

Conversation

@alexanderianblair
Copy link
Contributor

@alexanderianblair alexanderianblair commented Feb 20, 2025

Adds a test solving a 2D Darcy problem, based on MFEM Example 5.
Depends on #78

@alexanderianblair alexanderianblair marked this pull request as draft February 20, 2025 19:20
@codecov
Copy link

codecov bot commented Feb 20, 2025

Codecov Report

All modified and coverable lines are covered by tests ✅

Files with missing lines Coverage Δ
...bcs/MFEMVectorFEFunctionBoundaryFluxIntegratedBC.C 100.00% <100.00%> (ø)
src/kernels/MFEMDomainLFKernel.C 100.00% <100.00%> (ø)
src/kernels/MFEMMixedBilinearFormKernel.C 100.00% <100.00%> (ø)
src/kernels/MFEMMixedScalarCurlKernel.C 100.00% <100.00%> (ø)
src/kernels/MFEMMixedVectorGradientKernel.C 100.00% <100.00%> (ø)
src/kernels/MFEMVectorFEDivergenceKernel.C 100.00% <100.00%> (ø)

... and 1 file with indirect coverage changes

🚀 New features to boost your workflow:
  • ❄️ Test Analytics: Detect flaky tests, report on failures, and find test suite problems.

@alexanderianblair alexanderianblair force-pushed the alexanderianblair/darcy-example branch from 3a37987 to 4f73150 Compare February 26, 2025 18:24
@nmnobre nmnobre force-pushed the alexanderianblair/darcy-example branch 3 times, most recently from 94deb43 to ea7110c Compare April 4, 2025 16:46
@nmnobre
Copy link
Collaborator

nmnobre commented Apr 4, 2025

While helping @bryceshirley to rebase this onto main, we realised that the tests were failing due to a bug introduced in #84. I've applied the fix here to confirm all tests are still passing with @bryceshirley's changes, but I've also created a separate PR, #88, to get the fix in asap (at which point we should rebase this again). EDIT: Done.

@nmnobre nmnobre force-pushed the alexanderianblair/darcy-example branch from ea7110c to 57a78cd Compare April 4, 2025 21:01
Comment on lines +30 to +43
\begin{equation}
(\mathbf{K} \mathbf{u}, \mathbf{v})_\Omega + (\nabla p, \mathbf{v})_\Omega - (\nabla \cdot \mathbf{u}, q)_\Omega
= (\mathbf{f}, \mathbf{v})_\Omega
\quad \forall \mathbf{v} \in V, \forall q \in W
\end{equation}

where

\begin{equation}
\begin{split}
\mathbf{u} \in H(\mathrm{div})(\Omega) &: \mathbf{u} \cdot \hat n = \mathbf{g} \quad \text{on} \quad \partial \Omega \\
p \in L^2(\Omega) &: p = p_0 \quad \text{on} \quad \Gamma_\mathrm{D}
\end{split}
\end{equation}
Copy link
Contributor Author

Choose a reason for hiding this comment

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

This weak form here isn't consistent with the actual formulation used in the example (the choice of spaces and which gradient terms are present don't agree)

@nmnobre nmnobre force-pushed the alexanderianblair/darcy-example branch from 57a78cd to 5f549d1 Compare May 19, 2025 14:56
@nmnobre
Copy link
Collaborator

nmnobre commented Jul 22, 2025

Closing as completed via idaholab/moose#30637.

@nmnobre nmnobre closed this Jul 22, 2025
Sign up for free to subscribe to this conversation on GitHub. Already have an account? Sign in.

Labels

None yet

Projects

None yet

Development

Successfully merging this pull request may close these issues.

3 participants