Search

In this work, we establish that discontinuous Galerkin methods are capable of producing reliable approximations for a broad class of nonlinear variational problems. In particular, we demonstrate that these schemes provide essential flexibility by removing inter-element continuity while also guaranteeing convergent approximations in the quasiconvex case. Notably, quasiconvexity is the weakest form of convexity pertinent to elasticity. Furthermore, we show that in the non-convex case discrete minimisers converge to minimisers of the relaxed problem. In this case, the minimisation problem corresponds to the energy defined by the quasiconvex envelope of the original energy. Our approach covers all discontinuous Galerkin formulations known to converge for convex energies. This work addresses an open challenge in the vectorial calculus of variations: developing and rigorously justifying numerical schemes capable of reliably approximating nonlinear energy minimization problems with potentially singular solutions, which are frequently encountered in materials science.

While much attention of neural network methods is devoted to high-dimensional PDE problems, in this work we consider methods designed to work for elliptic problems on domains $\Omega \subset \mathbb{R} ^d, $ $d=1,2,3$ in association with more standard finite elements. We suggest to connect finite elements and neural network approximations through training, i.e., using finite element spaces to compute the integrals appearing in the loss functionals. This approach, retains the simplicity of classical neural network methods for PDEs, uses well established finite element tools (and software) to compute the integrals involved and it gains in efficiency and accuracy. We demonstrate that the proposed methods are stable and furthermore, we establish that the resulting approximations converge to the solutions of the PDE. Numerical results indicating the efficiency and robustness of the proposed algorithms are presented.

In the context of Discontinuous Galerkin methods, we study approximations of nonlinear variational problems associated with convex energies. We propose element-wise nonconforming finite element methods to discretize the continuous minimisation problem. Using $\Gamma$-convergence arguments we show that the discrete minimisers converge to the unique minimiser of the continuous problem as the mesh parameter tends to zero, under the additional contribution of appropriately defined penalty terms at the level of the discrete energies. We finally substantiate the feasibility of our methods by numerical examples.

Potassium sodium niobate is considered a prominent material system as a substitute for toxic lead-containing ferroelectric materials. It exhibits first-order phase transformations and ferroelectricity with potential applications ranging from energy conversion to innovative cooling technologies, hereby addressing urgent societal challenges. However, a major obstacle in the application of potassium sodium niobate is its multi-scale heterogeneity and the lack of understanding of its phase transition pathway and microstructure. This can be seen from the findings of Pop-Ghe 2021 et al. [1] which also reveal the occurrence of intermediate twinning during the phase transition. Here we show that intermediate twinning is a consequence of energy minimization. We employ a geometrically nonlinear electroelastic energy function for potassium sodium niobate, including the cubic-tetragonal-orthorhombic transformations and ferroelectricity. The construction of the minimizers is based on compatibility conditions which ensure continuous deformations and pole-free interfaces. These minimizers agree with the experimental observations, including laminates between the tetragonal variants under the cubic to tetragonal transformation, crossing twins under the tetragonal to orthorhombic transformation, intermediate twinning and spontaneous polarization. This shows how the full nonlinear electroelastic model provides a powerful tool in understanding, exploring and tailoring the electromechanical properties of complex ferroelectric ceramics.

Through modelling and simulations we show that material instabilities play a dominant role in the mechanical behavior of the fibrous collagen Extracellular Matrix (ECM), as observed when the ECM is deformed by contractile biological cells. We compare two families of fiber network models, Family 1 with stable and Family 2 with unstable force-stretch response of individual fibers in compression. The latter is a characteristic of post-buckling of beams with hierarchical structure. Our simulations reveal different compression instabilities at play in each family, namely, fiber collapse (buckling) in Family 2 and fiber element collapse (snap-through) in Family 1. These result in highly localized densification zones consisting of strongly aligned fibers emanating from individual contractile cells or joining neighbouring cells, as observed in experiments. Despite substantial differences in the response of the two families, our work underscores the importance of buckling/compression instabilities in the behavior of fibrous biological tissues, with implications in cancer invasion and metastasis.

We consider a model of energy minimization arising in the study of the mechanical behavior caused by cell contraction within a fibrous biological medium. The macroscopic model is based on the theory of non rank-one convex nonlinear elasticity for phase transitions. We study appropriate numerical approximations based on the discontinuous Galerkin treatment of higher gradients and used succesfully in numerical simulations of experiments. We show that the discrete minimizers converge in the limit to minimizers of the continuous problem. This is achieved by employing the theory of $\Gamma$-convergence of the approximate energy functionals to the continuous model when the discretization parameter tends to zero. The analysis is involved due to the structure of numerical approximations which are defined in spaces with lower regularity than the space where the minimizers of the continuous variational problem are sought. This fact leads to the development of a new approach to $\Gamma$-convergence, appropriate for discontinuous finite element discretizations, which can be applied to quite general energy minimization problems. Furthermore, the adoption of exponential terms penalising the interpenetration of matter requires a new framework based on Orlicz spaces for discontinuous Galerkin methods which is developed in this paper as well.

Through mechanical forces, biological cells remodel the surrounding collagen network, generating striking deformation patterns. Tethers-tracts of high densification and fiber alignment-form between cells, thinner bands emanate from cell clusters. While tethers facilitate cell migration and communication, how they form is unclear. Combining modeling, simulation and experiment, we show that tether formation is a densification phase transition of the extracellular matrix, caused by buckling instability of network fibers under cell-induced compression, featuring unexpected similarities with martensitic microstructures. Multiscale averaging yields a two-phase, bistable continuum energy landscape for fibrous collagen, with a densified/aligned second phase. Simulations predict strain discontinuities between the undensified and densified phase, which localizes within tethers as experimentally observed. In our experiments, active particles induce similar localized patterns as cells. This shows how cells exploit an instability to mechanically remodel the extracellular matrix simply by contracting, thereby facilitating mechanosensing, invasion and metastasis.