For problems in the calculus of cariations that exhibit the Lavrentiev phenomenon, it is known that the \textit{repulsion property} holds, that is, if one approximates the global minimizer in these problems by smooth functions, then the approximate energies will blow up. Thus standard numerical schemes, like the finite element method, may fail when applied directly to these type of problems. In this paper we prove that the repulsion property holds for variational problems in three dimensional elasticity that exhibit cavitation. In addition we propose a numerical scheme that circumvents the repulsion property, which is an adaptation of the Modica and Mortola functional for phase transitions in liquids, in which the phase function is coupled to the mechanical part of the stored energy functional, via the determinant of the deformation gradient. We show that the corresponding approximations by this method satisfy the lower bound $\Gamma$--convergence property in the multi-dimensional non--radial case. The convergence to the actual cavitating minimizer is established for a spherical body, in the case of radial deformations, and for the case of an elastic fluid without assuming radial symmetry.
翻译:对于在变分计算中出现拉夫连提耶夫现象的问题,已知存在排斥特性,即如果用平滑函数来逼近这些问题的全局极小值,则近似能量将会突增。因此,标准的数值方案(如有限元方法)直接应用于这些类型的问题可能会失败。在本文中,我们证明对于展现出空洞现象的三维弹性波动的变分问题,排斥特性成立。此外,我们提出了一种绕过排斥特性的数值方案,该方案是液体相变中Modica和Mortola函数(函数的相转移)的一种适应,其中相函数通过形变梯度行列式与存储能量函数的机械部分耦合。我们展示了通过该方法得到的相应近似满足非径向多维情况下的$\Gamma$-收敛性。对于球体的情况下,在径向形变的情况下推导了到实际空洞化极小值的收敛性,并且对于不假设径向对称的弹性流体的情况,也推导了其收敛性。