scholarly journals On a nonlocal Cahn–Hilliard model permitting sharp interfaces

Author(s):  
Olena Burkovska ◽  
Max Gunzburger

A nonlocal Cahn–Hilliard model with a non-smooth potential of double-well obstacle type that promotes sharp interfaces in the solution is presented. To capture long-range interactions between particles, a nonlocal Ginzburg–Landau energy functional is defined which recovers the classical (local) model as the extent of nonlocal interactions vanish. In contrast to the local Cahn–Hilliard problem that always leads to diffuse interfaces, the proposed nonlocal model can lead to a strict separation into pure phases of the substance. Here, the lack of smoothness of the potential is essential to guarantee the aforementioned sharp-interface property. Mathematically, this introduces additional inequality constraints that, in a weak formulation, lead to a coupled system of variational inequalities which at each time instance can be restated as a constrained optimization problem. We prove the well-posedness and regularity of the semi-discrete and continuous in time weak solutions, and derive the conditions under which pure phases are admitted. Moreover, we develop discretizations of the problem based on finite element methods and implicit–explicit time-stepping methods that can be realized efficiently. Finally, we illustrate our theoretical findings through several numerical experiments in one and two spatial dimensions that highlight the differences in features of local and nonlocal solutions and also the sharp interface properties of the nonlocal model.

Author(s):  
Ryuta Imai ◽  
Naoki Kasui ◽  
Masayuki Yamada ◽  
Koji Hada ◽  
Hiroyuki Fujiwara

In this paper, we propose a smoothing scheme for seismic wave propagation simulation. The proposed scheme is based on a diffusionized wave equation with the fourth-order spatial derivative term. So, the solution requires higher regularity in the usual weak formulation. Reducing the diffusionized wave equation to a coupled system of diffusion equations yields a mixed FEM to ease the regularity. We mathematically explain how our scheme works for smoothing. We construct a semi-implicit time integration scheme and apply it to the wave equation. This numerical experiment reveals that our scheme is effective for filtering short wavelength components in seismic wave propagation simulation.


1997 ◽  
Vol 8 (6) ◽  
pp. 623-637 ◽  
Author(s):  
LUCA BILLI

The motivation for this work arises from the study of the processes involved in the manufacturing of a class of composite materials, in particular, those that are obtained by injecting a resin through a porous preform. A one-dimensional model that describes the non-isothermal filtration of an incompressible fluid is presented, and it also includes the possibility of curing, i.e. the polymerization of the penetrating resin. It comes out a fully coupled system consisting of the heat diffusion equation, Darcy's law and an equation related to the kinetics of the chemical reaction. The system is regarded as a free boundary problem for the heat equation with non-constant discontinuous coefficients. Its weak formulation is studied and the local existence of solutions is proved.


2014 ◽  
Vol 26 (07) ◽  
pp. 1450014 ◽  
Author(s):  
Antoine Levitt

The multiconfiguration Dirac–Fock (MCDF) model uses a linear combination of Slater determinants to approximate the electronic N-body wave function of a relativistic molecular system, resulting in a coupled system of nonlinear eigenvalue equations, the MCDF equations. In this paper, we prove the existence of solutions of these equations in the weakly relativistic regime. First, using a new variational principle as well as the results of Lewin on the multiconfiguration non-relativistic model, and Esteban and Séré on the single-configuration relativistic model, we prove the existence of critical points for the associated energy functional, under the constraint that the occupation numbers are not too small. Then, this constraint can be removed in the weakly relativistic regime, and we obtain non-constrained critical points, i.e. solutions of the multiconfiguration Dirac–Fock equations.


Author(s):  
Robert Nourgaliev ◽  
Meng-Sing Liou ◽  
Theo Theofanous

We assess the state of the art in numerical prediction of interfacial instabilities due to shear in layered flows of viscous fluids. Basic ingredients of this assessment include linear stability analysis results for both miscible (diffuse) and immiscible (sharp) interfaces, the physics and resolution requirements of the critical layer, and convergence properties of the relevant numerical schemes. Behaviors of physically and numerically diffuse interfaces are contrasted and the case for a sharp interface treatment for reliable predictions of this class of flows is made.


2018 ◽  
Vol 24 (5) ◽  
pp. 1530-1555 ◽  
Author(s):  
CJ van Duijn ◽  
Andro Mikelić ◽  
Thomas Wick

In this paper, we present a full phase-field model for a fluid-driven fracture in a nonlinear poroelastic medium. The nonlinearity arises in the Biot equations when the permeability depends on porosity. This extends previous work (see Mikelić et al. Phase-field modeling of a fluid-driven fracture in a poroelastic medium. Comput Geosci 2015; 19: 1171–1195), where a fully coupled system is considered for the pressure, displacement, and phase field. For the extended system, we follow a similar approach: we introduce, for a given pressure, an energy functional, from which we derive the equations for the displacement and phase field. We establish the existence of a solution of the incremental problem through convergence of a finite-dimensional Galerkin approximation. Furthermore, we construct the corresponding Lyapunov functional, which is related to the free energy. Computational results are provided that demonstrate the effectiveness of this approach in treating fluid-driven fracture propagation. Specifically, our numerical findings confirm differences with test cases using the linear Biot equations.


2019 ◽  
Vol 40 (2) ◽  
pp. 1122-1153
Author(s):  
Rodolfo Araya ◽  
Rodolfo Rodríguez ◽  
Pablo Venegas

Abstract This paper deals with the numerical analysis of a system of second order in time partial differential equations modeling the vibrations of a coupled system that consists of an elastic solid in contact with an inviscid compressible fluid. We analyze a weak formulation with the unknowns in both media being the respective displacement fields. For its numerical approximation, we propose first a semidiscrete in space discretization based on standard Lagrangian elements in the solid and Raviart–Thomas elements in the fluid. We establish its well-posedness and derive error estimates in appropriate norms for the proposed scheme. In particular, we obtain an $\mathrm L^{\infty }(\mathrm L^2)$ optimal rate of convergence under minimal regularity assumptions of the solution, which are proved to hold for appropriate data of the problem. Then we consider a fully discrete approximation based on a family of implicit finite difference schemes in time, from which we obtain optimal error estimates for sufficiently smooth solutions. Finally, we report some numerical results, which allow us to assess the performance of the method. These results also show that the numerical solution is not polluted by spurious modes as is the case with other alternative approaches.


2014 ◽  
Vol 24 (07) ◽  
pp. 1265-1341 ◽  
Author(s):  
Elisabetta Rocca ◽  
Riccarda Rossi

In this paper, we analyze a PDE system arising in the modeling of phase transition and damage phenomena in thermoviscoelastic materials. The resulting evolution equations in the unknowns ϑ (absolute temperature), u (displacement), and χ (phase/damage parameter) are strongly nonlinearly coupled. Moreover, the momentum equation for u contains χ-dependent elliptic operators, which degenerate at the pure phases (corresponding to the values χ = 0 and χ = 1), making the whole system degenerate. That is why, we have to resort to a suitable weak solvability notion for the analysis of the problem: it consists of the weak formulations of the heat and momentum equation, and, for the phase/damage parameter χ, of a generalization of the principle of virtual powers, partially mutuated from the theory of rate-independent damage processes. To prove an existence result for this weak formulation, an approximating problem is introduced, where the elliptic degeneracy of the displacement equation is ruled out: in the framework of damage models, this corresponds to allowing for partial damage only. For such an approximate system, global-in-time existence and well-posedness results are established in various cases. Then, the passage to the limit to the degenerate system is performed via suitable variational techniques.


Author(s):  
Theo Theofanous ◽  
Svetlana Sushchikh ◽  
Robert Nourgaliev

We show that diffuse interfaces in shear flow are prone to instabilities that are of the same origin and character as the Yih-instability (for sharp interfaces between immiscible fluids). We also show that the development of these “Yih-like” instabilities towards “Yih” is monotonic, and that the agreement with Yih becomes quantitative as the limits of zero thickness and no-diffusion are approached.


2021 ◽  
Vol 87 (2) ◽  
Author(s):  
S. A. Henneberg ◽  
S. R. Hudson ◽  
D. Pfefferlé ◽  
P. Helander

Combined plasma–coil optimization approaches for designing stellarators are discussed and a new method for calculating free-boundary equilibria for multiregion relaxed magnetohydrodynmics (MRxMHD) is proposed. Four distinct categories of stellarator optimization, two of which are novel approaches, are the fixed-boundary optimization, the generalized fixed-boundary optimization, the quasi-free-boundary optimization, and the free-boundary (coil) optimization. These are described using the MRxMHD energy functional, the Biot–Savart integral, the coil-penalty functional and the virtual casing integral and their derivatives. The proposed free-boundary equilibrium calculation differs from existing methods in how the boundary-value problem is posed, and for the new approach it seems that there is not an associated energy minimization principle because a non-symmetric functional arises. We propose to solve the weak formulation of this problem using a spectral-Galerkin method, and this will reduce the free-boundary equilibrium calculation to something comparable to a fixed-boundary calculation. In our discussion of combined plasma–coil optimization algorithms, we emphasize the importance of the stability matrix.


Analysis ◽  
2017 ◽  
Vol 37 (4) ◽  
Author(s):  
Patrick Henkemeyer

AbstractWe discuss certain quantitative geometric properties of energy stationary currents describing minimal surfaces under gravitational forces. Enclosure theorems give statements about the confinement of the support of currents to certain enclosing sets on the basis that one knows something about the position of their boundaries. These results are closely related to non-existence theorems for currents with connected support. Finally, we define a weak formulation in the theory of varifolds for the curvature flow associated to this energy functional. We extend the enclosure results to the flow and discuss several comparison principles.


Sign in / Sign up

Export Citation Format

Share Document