Search

In this paper, we study a hyperbolic relaxation of the viscous Cahn--Hilliard system with zero Neumann boundary conditions. In fact, we consider a relaxation term involving the second time derivative of the chemical potential in the first equation of the system. We develop a well-posedness, continuous dependence and regularity theory for the initial-boundary value problem. Moreover, we investigate the asymptotic behavior of the system as the relaxation parameter tends to 0 and prove the convergence to the viscous Cahn--Hilliard system., Comment: Key words: Cahn--Hilliard system, hyperbolic relaxation, partial differential equations, initial-boundary value problem, well-posedness, continuous dependence, regularity, asymptotic convergence

This paper investigates an optimal control problem associated with a two-dimensional multi-species Cahn-Hilliard-Keller-Segel tumor growth model, which incorporates complex biological processes such as species diffusion, chemotaxis, angiogenesis, and nutrient consumption, resulting in a highly nonlinear system of nonlinear partial differential equations. The modeling derivation and corresponding analysis have been addressed in a previous contribution. Building on this foundation, the scope of this study involves investigating a distributed control problem with the goal of optimizing a tracking-type cost functional. This latter aims to minimize the deviation of tumor cell location from desired target configurations while penalizing the costs associated with implementing control measures, akin to introducing a suitable medication. Under appropriate mathematical assumptions, we demonstrate that sufficiently regular solutions exhibit continuous dependence on the control variable. Furthermore, we establish the existence of optimal controls and characterize the first-order necessary optimality conditions through a suitable variational inequality.

In this paper we study the optimal control of an initial-boundary value problem for the classical nonviscous Cahn-Hilliard system with zero Neumann boundary conditions. Phase field systems of this type govern the evolution of diffusive phase transition processes with conserved order parameter. For such systems, optimal control problems have been studied in the past. We focus here on the situation when the cost functional of the optimal control problem contains a sparsity-enhancing nondifferentiable term like the L1-norm. For such cases, we establish first-order necessary and second-order sufficient optimality conditions for locally optimal controls, where in the approach to second-order sufficient conditions we employ a technique introduced by E. Casas, C. Ryll and F. Tr\"oltzsch in the paper [SIAM J. Control Optim. 53 (2015), 2168-2202]. The main novelty of this paper is that this method, which has recently been successfully applied to systems of viscous Cahn-Hilliard type, can be adapted also to the classical nonviscous case. Since in the case without viscosity the solutions to the state and adjoint systems turn out to be considerably less regular than in the viscous case, numerous additional technical difficulties have to be overcome, and additional conditions have to be imposed. In particular, we have to restrict ourselves to the case when the nonlinearity driving the phase separation is regular, while in the presence of a viscosity term also nonlinearities of logarithmic type turn could be admitted. In addition, the implicit function theorem, which was employed to establish the needed differentiability properties of the control-to-state operator in the viscous case, does not apply in our situation and has to be substituted by other arguments., Comment: 40 pages. The key words are: Cahn-Hilliard equation, optimal control, sparsity, first- and second- order optimality conditions. The results contained in this paper can be compared with the corresponding ones in the paper arXiv:2402.18506. arXiv admin note: substantial text overlap with arXiv:2402.18506

Given an undirected connected graph $G = (V(G), E(G))$ on $n$ vertices, the minimum Monitoring Edge-Geodetic Set (MEG-set) problem asks to find a subset $M \subseteq V(G)$ of minimum cardinality such that, for every edge $e \in E(G)$, there exist $x,y \in M$ for which all shortest paths between $x$ and $y$ in $G$ traverse $e$. We show that, for any constant $c < \frac{1}{2}$, no polynomial-time $(c \log n)$-approximation algorithm for the minimum MEG-set problem exists, unless $\mathsf{P} = \mathsf{NP}$.

Partial differential equation (PDE) models for infectious disease have received renewed interest in recent years. Most models of this type extend classical compartmental formulations with additional terms accounting for spatial dynamics, with Fickian diffusion being the most common such term. However, while diffusion may be appropriate for modeling vector-borne diseases, or diseases among plants or wildlife, the spatial propagation of airborne diseases in human populations is heavily dependent on human contact and mobility patterns, which are not necessarily well-described by diffusion. By including an additional chemotaxis-inspired term, in which the infection is propagated along the positive gradient of the susceptible population (from regions of low- to high-density of susceptibles), one may provide a more suitable description of these dynamics. This article introduces and analyzes a mathematical model of infectious disease incorporating a modified chemotaxis-type term. The model is analyzed mathematically and the well-posedness of the resulting PDE system is demonstrated. A series of numerical simulations are provided, demonstrating the ability of the model to naturally capture important phenomena not easily observed in standard diffusion models, including propagation over long spatial distances over short time scales and the emergence of localized infection hotspots

In this paper we study the optimal control of a parabolic initial-boundary value problem of viscous Cahn-Hilliard type with zero Neumann boundary conditions. Phase field systems of this type govern the evolution of diffusive phase transition processes with conserved order parameter. It is assumed that the nonlinear functions driving the physical processes within the spatial domain are double-well potentials of logarithmic type whose derivatives become singular at the boundary of their respective domains of definition. For such systems, optimal control problems have been studied in the past. We focus here on the situation when the cost functional of the optimal control problem contains a nondifferentiable term like the L1-norm, which leads to sparsity of optimal controls. For such cases, we establish first-order necessary and second-order sufficient optimality conditions for locally optimal controls. In the approach to second-order sufficient conditions, the main novelty of this paper, we adapt a technique introduced by E. Casas, C. Ryll and F. Tr\"oltzsch in the paper [SIAM J. Control Optim. 53 (2015), 2168-2202]. In this paper, we show that this method can also be successfully applied to systems of viscous Cahn-Hilliard type with logarithmic nonlinearity. Since the Cahn-Hilliard system corresponds to a fourth-order partial differential equation in contrast to the second-order systems investigated before, additional technical difficulties have to be overcome., Comment: 43 pages. The keywords are viscous Cahn-Hilliard equation, singular potential, optimal control, first-order and second-order optimality conditions

This paper deals with a nonlocal model for a hyperbolic phase field system coupling the standard energy balance equation for temperature with a dynamic for the phase variable: the latter includes an inertial term and a nonlocal convolution-type operator where the family of kernels depends on a small parameter. We rigorously study the asymptotic convergence of the system as the approximating parameter tends to zero and we obtain at the limit the local system with the elliptic laplacian operator acting on the phase variable. Our analysis is based on some asymptotic properties on nonlocal-to-local convergence that have been recently and successfully applied to families of Cahn--Hilliard models.

This work investigates the well-posedness and optimal control of a sixth-order Cahn-Hilliard equation, a higher-order variant of the celebrated and well-established Cahn-Hilliard equation. The equation is endowed with a source term, where the control variable enters as a distributed mass regulator. The inclusion of additional spatial derivatives in the sixth-order formulation enables the model to capture curvature effects, leading to a more accurate depiction of isothermal phase separation dynamics in complex materials systems. We provide a well-posedness result for the aforementioned system when the corresponding nonlinearity of double-well shape is regular and then analyze a corresponding optimal control problem. For the latter, existence of optimal controls is established, and the first-order necessary optimality conditions are characterized via a suitable variational inequality. These results aim at contributing to improve the understanding of the mathematical properties and control aspects of the sixth-order Cahn-Hilliard equation, offering potential applications in the design and optimization of materials with tailored microstructures and properties., Comment: Key words: Sixth-order Cahn--Hilliard equation, functionalized Cahn--Hilliard equation, Willmore regularization, curvature effects, well-posedness, optimal control, first-order necessary optimality conditions

We investigate a new diffuse-interface model that describes creeping two-phase flows (i.e., flows exhibiting a low Reynolds number), especially flows that permeate a porous medium. The system of equations consists of a Brinkman equation for the volume averaged velocity field as well as a convective Cahn--Hilliard equation with dynamic boundary conditions for the phase-field, which describes the location of the two fluids within the domain. The dynamic boundary conditions are incorporated to model the interaction of the fluids with the wall of the container more precisely. In particular, they allow for a dynamic evolution of the contact angle between the interface separating the fluids and the boundary, and also for a convection-induced motion of the corresponding contact line. For our model, we first prove the existence of global-in-time weak solutions in the case where regular potentials are used in the Cahn--Hilliard subsystem. In this case, we can further show the uniqueness of the weak solution under suitable additional assumptions. Moreover, we further prove the existence of weak solutions in the case of singular potentials. Therefore, we regularize such singular potentials by a Yosida approximation, such that the results for regular potentials can be applied, and eventually pass to the limit in this approximation scheme.

Two-scale topology optimization, combined with the design of microstructure families with a broad range of effective material parameters, is increasingly widely used in many fabrication applications to achieve a target deformation behavior for a variety of objects. The main idea of this approach is to optimize the distribution of material properties in the object partitioned into relatively coarse cells, and then replace each cell with microstructure geometry that mimics these material properties. In this paper, we focus on adapting this approach to complex shapes in situations when preserving the shape's surface is important. Our approach extends any regular (i.e. defined on a regular lattice grid) microstructure family to complex shapes, by enriching it with individually optimized cut-cell tiles adapted to the geometry of the cut-cell. We propose an automated and robust pipeline based on this approach, and we show that the performance of the regular microstructure family is only minimally affected by our extension while allowing its use on 2D and 3D shapes of high complexity.

In this paper, we explore the solvability and the optimal control problem for a compartmental model based on reaction-diffusion partial differential equations describing a transmissible disease. The nonlinear model takes into account the disease spreading due to the human social diffusion, under a dynamic heterogeneity in infection risk. The analysis of the resulting system provides the existence proof for a global solution and determines the conditions of optimality to reduce the concentration of the infected population in certain spatial areas., Comment: Keywords: existence and uniqueness of solutions, optimal control in coefficients, partial differential equations, reaction-diffusion system, epidemic models, COVID- 19

In this paper we derive a new model for visco-elasticity with large deformations where the independent variables are the stretch and the rotation tensors which intervene with second gradients terms accounting for physical properties in the principle of virtual power. Another basic feature of our model is that there is conditional compatibility, entering the model as kinematic constraints and depending on the magnitude of an internal force associated to dislocations. Moreover, due to the kinematic constraints, the virtual velocities depend on the solutions of the problem. As a consequence, the variational formulation of the problem and the related mathematical analysis are neither standard nor straightforward. We adopt the strategy to invert the kinematic constraints through Green propagators, obtaining a system of integro-differential coupled equations. As a first mathematical step, we develop the analysis of the model in a simplified setting, i.e. considering the quasi-stationary version of the full system where we neglect inertia. In this context, we prove the existence of a global in time strong solution in three space dimensions for the system, employing techniques from PDEs and convex analysis, thus obtaining a novel breakthrough in the field of three-dimensional finite visco-elasticity described in terms of the stretch and rotation variables. We also study a limit problem, letting the magnitude of the internal force associated to dislocations tend to zero, in which case the deformation becomes incompatible and the equations takes the form of a coupled system of PDEs. For the limit problem we obtain global existence, uniqueness and continuous dependence from data in three space dimensions.

This paper is concerned with the well-posedness and optimal control problem of a reaction-diffusion system for an epidemic Susceptible-Infected-Recovered-Susceptible (SIRS) mathematical model in which the dynamics develops in a spatially heterogeneous environment. Using as control variables the transmission rates $u_{i}$ and $u_{e}$ of contagion resulting from the contact with both asymptomatic and symptomatic persons, respectively, we optimize the number of exposed and infected individuals at a final time $T$ of the controlled evolution of the system. More precisely, we search for the optimal $u_{i}$ and $u_{e}$ such that the number of infected plus exposed does not exceed at the final time a threshold value $\Lambda$, fixed a priori. We prove here the existence of optimal controls in a proper functional framework and we derive the first-order necessary optimality conditions in terms of the adjoint variables., Comment: Keywords: COVID-19, partial differential equations, reaction-diffusion system, epidemic models, existence of solutions, uniqueness, optimal control

In this note, we study the optimal control of a nonisothermal phase field system of Cahn-Hilliard type that constitutes an extension of the classical Caginalp model for nonisothermal phase transitions with a conserved order parameter. It couples a Cahn-Hilliard type equation with source term for the order parameter with the universal balance law of internal energy. In place of the standard Fourier form, the constitutive law of the heat flux is assumed in the form given by the theory developed by Green and Naghdi, which accounts for a possible thermal memory of the evolution. This has the consequence that the balance law of internal energy becomes a second-order in time equation for the thermal displacement or freezing index, that is, a primitive with respect to time of the temperature. Another particular feature of our system is the presence of the source term in the equation for the order parameter, which entails further mathematical difficulties because the mass conservation of the order parameter is no longer satisfied. We study the case that the double-well potential driving the evolution of the phase transition is given by the nondifferentiable double obstacle potential, thereby complementing recent results obtained for the differentiable cases of regular and logarithmic potentials. We derive first-order necessary optimality conditions for the control problem. The analysis is carried out by employing the so-called deep quench approximation in which the nondifferentiable double obstacle potential is approximated by a family of logarithmic potentials for which meaningful first-order necessary optimality conditions are available. Since the results for the logarithmic potentials crucially depend on the validity of the so-called strict separation property which is only available in the spatially two-dimensional situation, our whole analysis is restricted to the two-dimensional case., Comment: Keywords: optimal control, nonisothermal Cahn-Hilliard equation, thermal memory, Cahn-Hilliard equation with source term, Cahn-Hilliard-Oono equation. arXiv admin note: text overlap with arXiv:2207.00375, arXiv:2303.00488

In this note, we study the optimal control of a nonisothermal phase field system of Cahn-Hilliard type that constitutes an extension of the classical Caginalp model for nonisothermal phase transitions with a conserved order parameter. The system couples a Cahn-Hilliard type equation with source term for the order parameter with the universal balance law of internal energy. In place of the standard Fourier form, the constitutive law of the heat flux is assumed in the form given by the theory developed by Green and Naghdi, which accounts for a possible thermal memory of the evolution. This has the consequence that the balance law of internal energy becomes a second-order in time equation for the thermal displacement or freezing index, that is, a primitive with respect to time of the temperature. Another particular feature of our system is the presence of the source term in the equation for the order parameter, which entails additional mathematical difficulties because the mass conservation of the order parameter, typical of the classic Cahn-Hilliard equation, is no longer satisfied. In this paper, we analyze the case that the double-well potential driving the evolution of the phase transition is differentiable, either (in the regular case) on the whole set of reals or (in the singular logarithmic case) on a finite open interval; nondifferentiable cases like the double obstacle potential are excluded from the analysis. We prove the Fr\'echet differentiability of the control-to-state operator between suitable Banach spaces for both the regular and the logarithmic cases and establish the solvability of the corresponding adjoint systems in order to derive the associated first-order necessary optimality conditions for the optimal control problem.

Subgradient methods are the natural extension to the non-smooth case of the classical gradient descent for regular convex optimization problems. However, in general, they are characterized by slow convergence rates, and they require decreasing step-sizes to converge. In this paper we propose a subgradient method with constant step-size for composite convex objectives with $\ell_1$-regularization. If the smooth term is strongly convex, we can establish a linear convergence result for the function values. This fact relies on an accurate choice of the element of the subdifferential used for the update, and on proper actions adopted when non-differentiability regions are crossed. Then, we propose an accelerated version of the algorithm, based on conservative inertial dynamics and on an adaptive restart strategy, that is guaranteed to achieve a linear convergence rate in the strongly convex case. Finally, we test the performances of our algorithms on some strongly and non-strongly convex examples., Comment: 18 pages, 3 figures. Minor changes, extended bibliographical references, new example in Section 4

We propose a new Cahn-Hilliard phase field model coupled to incompressible viscoelasticity at large strains, obtained from a diffuse interface mixture model and formulated in the Eulerian configuration. A new kind of diffusive regularization, of Allen-Cahn type, is introduced in the transport equation for the deformation gradient, together with a regularizing interface term depending on the gradient of the deformation gradient in the free energy of the system. We study the global existence of a weak solution for the model. While standard diffusive regularizations of the transport equation for the deformation gradient presented in literature allows the existence study only for simplified cases, i.e. in two space dimensions and for convex elastic free energy densities of Neo-Hookean type which are independent from the phase field variable, the present regularization allows to study more general cases. In particular, we obtain the global existence of a weak solution in three space dimensions and for generic nonlinear elastic energy densities with polynomial growth. Our analysis considers elastic free energy densities which depend on the phase field variable and which can possibly degenerate for some values of the phase field variable. By means of an iterative argument based on elliptic regularity bootstrap steps, we find the maximum allowed polynomial growths of the Cahn-Hilliard potential and the elastic energy density which guarantee the existence of a solution in three space dimensions. We propose two unconditionally energy stable finite element approximations of the model, based on convex splitting ideas and on the use of a scalar auxiliary variable, proving the existence and stability of discrete solutions. We finally report numerical results for different test cases with shape memory alloy type free energy with pure phases characterized by different elastic properties.