Numerical analysis of the Biot equations coupled to frictional contact mechanics

Technische Samenvatting: Numerieke Analyse van de Biot-vergelijkingen gekoppeld aan wrijvingscontactmechanica

Probleemformulering
Dit artikel behandelt de wiskundige modellering en numerieke benadering van een porovisco-elastisch medium verzadigd met een vloeistof, onderhevig aan wrijvingscontact met een rigide obstakel. Het fysieke systeem wordt beheerst door de quasi-statische Biot-vergelijkingen, die lineaire elasticiteit (specifiek het Kelvin-Voigt visco-elastische model) koppelen aan Darcy-stroming. De contactmechanica worden gemodelleerd via twee specifieke randvoorwaarden op een deel van de rand ($\Gamma_c$):

1. Normale compliantie: Een constitutieve relatie waarbij de normale contactkracht proportioneel is aan de penetratiediepte, wat de niet-gladde Signorini-conditie vermijdt.
2. Coulomb-wrijving: Een wet die de tangentiële beweging regelt, waarbij het lichaam blijft staan (stick) als de tangentiële kracht onder een wrijvingsgrens ligt (proportioneel aan de normale kracht) en glijdt (slip) indien dit niet het geval is.

Het resulterende wiskundige model is een systeem dat een lineaire partiële differentiaalvergelijking (de massabalansvergelijking) koppelt met een niet-lineaire variationele ongelijkheid (de impulsbalansvergelijking met contactvoorwaarden). De primaire uitdaging ligt in de koppeling van stroming en mechanica, de tijdsafhankelijkheid van het probleem, en de niet-lineariteit geïntroduceerd door de wrijvingscontactvoorwaarden.

Methodologie
De auteurs stellen een volledig discrete numerieke methode voor en analyseren deze voor dit probleem:

• Ruimtelijke discretisatie: Conforme eindelementen worden gebruikt. De verplaatsingsruimte wordt benaderd met Lagrangiaanse elementen van orde $k_1$, en de drukruimte met elementen van orde $k_2$.
• Temporele discretie: De impliciete Euler-methode wordt toegepast voor de tijdstapverwerking.
• Strategie voor het bewijs van existentie: Om de existentie en uniciteit van de discrete oplossing vast te stellen, hanteren de auteurs een innovatieve theoretische aanpak. Eerst beschouwen zij een vereenvoudigd probleem waarbij de contacttermen "vertraagd" zijn (wat het wrijvingsmodel effectief reduceert tot het Tresca-model). Voor dit vereenvoudigde probleem maken zij gebruik van het fixed-stress splitting schema — doorgaans een praktisch algoritme voor het ontkoppelen van stroming en mechanica — als een theoretisch instrument om de existentie van een oplossing te bewijzen via een contractie-argument. Vervolgens wordt een fixed-point argument uit de contactmechanica-literatuur toegepast om het volledige Coulomb-wrijvingsmodel te herstellen.
• Foutanalyse: De afleiding van a priori foutschattingen combineert technieken uit de bestaande literatuur over de Biot-vergelijkingen (Girault et al., Phillips en Wheeler) en tijd-afhankelijke variationele ongelijkheden (Sofonea en Han). Een cruciale analytische stap is het aantonen dat de koppelingsterm tussen stroming en mechanica wegvalt in de stabiliteits- en foutanalyse, waardoor deze de schattingen niet verslechtert.

Belangrijkste Resultaten

1. Welbepaaldheid: Het artikel stelt de existentie en uniciteit van de discrete oplossing vast, mits de tijdstapgrootte voldoet aan een specifieke voorwaarde gerelateerd aan de coërciviteitsconstanten van de bilineaire vormen en de Lipschitz-constante van de contactfunctionaal.
2. Stabiliteit: Er wordt een stabiliteitsschatting afgeleid, die aantoont dat de discrete oplossing onafhankelijk van de roostergrootte ($h$) en tijdstap ($\Delta t$) begrensd blijft, ervan uitgaande dat de tijdstap voldoende klein is om de niet-lineariteit van de contacttermen te behandelen.
3. A Priori Foutschattingen: De auteurs leiden foutschattingen af in de $H^1$-norm voor verplaatsing en de $L^2$-norm voor druk. De theoretische convergentiesnelheden zijn:
   • Tijd: Eerste-orde convergentie ($O(\Delta t)$).
   • Ruimte: De druk convergeert met optimale orde ($O(h^s)$). De convergentieorde van de verplaatsing wordt met een factor $\sqrt{h}$ verminderd ten opzichte van standaard elasticiteitsproblemen, wat resulteert in een orde van ongeveer $O(h^{r-0.5})$ voor lineaire elementen. Deze reductie wordt toegeschreven aan de contactvoorwaarden die de regulariteit van de oplossing verlagen.
4. Numerieke Experimenten: Een numeriek experiment is uitgevoerd met de FEniCSx-software. De resultaten bevestigden de theoretische voorspellingen:
   • Eerste-orde convergentie in de tijd werd waargenomen voor zowel druk als verplaatsing.
   • Ruimtelijke convergentie voor druk was optimaal.
   • Ruimtelijke convergentie voor verplaatsing was licht verminderd (geobserveerde snelheden rond de 0,8 voor lineaire elementen), wat consistent is met de theoretische verwachting van verminderde regulariteit door contact.

Betekenis en Claims
Het artikel beweert een specifiek onderzoeksgat te overbruggen in de numerieke analyse van porositeit. Terwijl er uitgebreide literatuur bestaat over elastische lichamen met contact en over porositeit zonder contact, was de combinatie van tijd-afhankelijke Biot-vergelijkingen met wrijvingscontactvoorwaarden (specifiek normale compliantie en Coulomb-wrijving) voorheen niet geanalyseerd met a priori foutschattingen.

De auteurs benadrukken dat hun werk de behapbare analyse van visco-elastische contactproblemen (voorheen vastgesteld voor niet-poreuze media) uitbreidt naar de gekoppelde porovisco-elastische setting. De betekenis van het werk ligt in:

• Het bieden van een rigoureuze wiskundige basis voor de numerieke simulatie van gekoppelde stroming en wrijvingscontact in poreuze media.
• Het aantonen dat de koppelingsterm de afleiding van foutschattingen niet belemmert wanneer technieken uit de porositeitsleer en de contactmechanica worden gecombineerd.
• Het aanbieden van een volledig discrete methode met bewezen stabiliteit en convergentiesnelheden, gevalideerd door numerieke experimenten.

Het artikel blijft bescheiden wat betreft toepassingen en richt zich strikt op de wiskundige analyse en verificatie van de voorgestelde numerieke methode, in plaats van het voorstellen van nieuwe fysieke toepassingen of toekomstige uitbreidingen buiten de reikwijdte van de afgeleide schattingen.