Skip to content

contact-miniapp #348

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.

Sign up for GitHub

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

Merged
tzanio merged 8 commits into from
Dec 8, 2025
Merged

contact-miniapp #348

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
8 commits
Select commit Hold shift + click to select a range
0cafb54
contact-miniapp
psocratis Oct 20, 2025
f1373e4
amgf paper entry
psocratis Oct 20, 2025
8211ef6
Update src/examples.md
psocratis Dec 4, 2025
7b9157b
Update src/examples.md
psocratis Dec 4, 2025
be034d6
Merge branch 'mfem-4.9' into contact-miniapp
psocratis Dec 4, 2025
279e2d2
Fix verb agreement in miniapp documentation
psocratis Dec 4, 2025
0206b2a
Update src/examples.md
psocratis Dec 4, 2025
9542058
reviewer's comments
psocratis Dec 8, 2025
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
17 changes: 17 additions & 0 deletions src/examples.md
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -2346,7 +2346,23 @@ _This miniapp has only a parallel ([contact-patch-test.cpp](https://github.com/m
<div style="clear:both;"/></div>
<br></div>

<div id="contact" markdown="1">
## Contact miniapp

<a href="../img/gallery/workshop25/contact.mp4" target="_blank">
<img class="floatright" width="250" src="../img/gallery/workshop25/contact.png">
</a>

This [miniapp](https://github.com/mfem/mfem/blob/master/miniapps/contact/contact.cpp) demonstrates how to use MFEM in combination with [Tribol](https://github.com/LLNL/Tribol)
to simulate frictionless contact between deformable solids. The contact gap function and its Jacobian are evaluated using Tribol’s mortar segment-to-segment method (see [Puso and Laursen, 2004](https://doi.org/10.1016/j.cma.2003.10.010)). An **Interior-Point (IP)** optimization solver is used to solve an inequality-constrained minimization problem enforcing a non-penetration condition. Linear systems arising in the IP iterations are solved using **Conjugate Gradient (CG)** preconditioned with the [AMG with Filtering (AMGF)](https://arxiv.org/abs/2505.18576) solver.

Several benchmark problems are provided, including the two-block compression, ironing and beam-sphere tests. These examples illustrate how MFEM and Tribol can be combined to build large-scale scalable simulations of contact mechanics problems.
Copy link
Member

@acfisher acfisher Dec 5, 2025

Choose a reason for hiding this comment

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

Are there any limitations to variety of Tribol contact methods that someone could use this miniapp to start with? It would be good to call that out.

Copy link
Member Author

@psocratis psocratis Dec 8, 2025

Choose a reason for hiding this comment

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

The main limitation is that the contact is frictionless, as already mentioned.


_This miniapp has only a parallel ([contact.cpp](https://github.com/mfem/mfem/blob/master/miniapps/contact/contact.cpp)) implementation. For more details, please see the documentation in [miniapps/contact/README.md](https://github.com/mfem/mfem/blob/master/miniapps/contact/README.md).
**We recommend viewing Example 2 and the Tribol miniapp before viewing this miniapp.**_

<div style="clear:both;"/></div>
<br></div>

<!-- ------------------------------------------------------------------------- -->

Expand Down Expand Up @@ -2492,6 +2508,7 @@ function update()
+ showElement("spde", (diffusion || nonlocal || stochastic) && h1 && galerkin && amg)
+ showElement("dpgminiapp", (diffusion || convectiondiffusion || maxwell || acoustics || wave) && (group2) && (dpg || amr) && (amg || ams || ads || pcg))
+ showElement("tribol", elasticity && h1 && galerkin && (superlu || minres || jacobi))
+ showElement("contact", elasticity && h1 && galerkin && (pcg || amg))

; // ...end of expression

Expand Down
1 change: 1 addition & 0 deletions src/publications.md
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -11,6 +11,7 @@
1. T. Dzanic, T. Kolev, K. Mittal, [A method for bounding high-order finite element functions: Applications to mesh validity and bounds-preserving limiters](https://www.sciencedirect.com/science/article/pii/S0021999125007260), *Journal of Computational Physics*, 545, 114444, **2026**. Also available as [arXiv:2504.11688](https://arxiv.org/abs/2504.11688).

#### 2025
1. S. Petrides, T. Hartland, T. Kolev, C. S. Lee, M. Puso, J. Solberg, E. B. Chin, J. Wang, C. Petra, [AMG with Filtering: An Efficient Preconditioner for Interior Point Methods in Large-Scale Contact Mechanics Optimization](https://arxiv.org/abs/2505.18576), *in review*, **2025**. Also available as [arXiv:2505.18576](https://arxiv.org/abs/2505.18576).
1. K. Mittal, J.-L. Barrera, T. Kolev, M. Schmidt, V. Tomov, [High-Order Mesh hr-Adaptivity for Surface Fitting to Implicit Geometries](https://link.springer.com/chapter/10.1007/978-3-031-76988-7_18). *Lecture Notes in Computational Science and Engineering*, **2025**.
1. T. Kolev, B. Lazarov, K. Mittal, M. Schmidt, V. Tomov, [PDE-Constrained High-Order Mesh Optimization](https://arxiv.org/pdf/2507.01917), *to appear*, **2025**.
1. J. Andrej, T. Kolev, B. S. Lazarov, [Scalable Analysis and Design Using Automatic Differentiation](https://arxiv.org/abs/2506.00746), *to appear*, **2025**. Also available as [arXiv:2506.00746](https://arxiv.org/abs/2506.00746).
Expand Down