Search

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.

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

This paper is concerned with the well-posedness of a diffusion-reaction system for a Susceptible-Exposed-Infected-Recovered (SEIR) mathematical model. This model is written in terms of four nonlinear partial differential equations with nonlinear diffusions, depending on the total amount of the SEIR populations. The model aims at describing the spatio-temporal spread of the COVID-19 pandemic and is a variation of the one recently introduced, discussed and tested in [A. Viguerie et al, Diffusion-reaction compartmental models formulated in a continuum mechanics framework: application to COVID-19, mathematical analysis, and numerical study, Comput. Mech. 66 (2020) 1131-1152]. Here, we deal with the mathematical analysis of the resulting Cauchy-Neumann problem: the existence of solutions is proved in a rather general setting and a suitable time discretization procedure is employed. It is worth mentioning that the uniform boundedness of the discrete solution is shown by carefully exploiting the structure of the system. Uniform estimates and passage to the limit with respect to the time step allow to complete the existence proof. Then, two uniqueness theorems are offered, one in the case of a constant diffusion coefficient and the other for more regular data, in combination with a regularity result for the solutions., Comment: Key words: COVID-19, compartmental model, partial differential equations, diffusion-reaction system, initial-boundary value problem, existence of solutions, uniqueness

We propose a new class of phase field models coupled to viscoelasticity with large deformations, obtained from a diffuse interface mixture model composed by a phase with elastic properties and a liquid phase. The model is formulated in the Eulerian configuration and it is derived by imposing the mass balance for the mixture components and the momentum balance that comes from a generalized form of the principle of virtual powers. The latter considers the presence of a system of microforces and microstresses associated to the microscopic interactions between the mixture's constituents together with a system of macroforces and macrostresses associated to their viscoelastic behavior, taking into account also the friction between the phases. The free energy density of the system is given as the sum of a Cahn-Hilliard term and an elastic polyconvex term, with a coupling between the phase field variable and the elastic deformation gradient in the elastic contribution. General constitutive assumptions complying with a mechanical version of the second law of thermodynamics in isothermal situations are taken. We study the global existence of a weak solution for a simplified and regularized version of the general model, which considers an incompressible elastic free energy of Neo-Hookean type with elastic coefficients depending on the phase field variable. The regularization is properly designed to deal with the coupling between the phase field variable and the elastic deformation gradient in the elastic energy density. The analysis is made both in two and three space dimensions.

We propose a diffuse interface model to describe tumor as a multicomponent deformable porous medium. We include mechanical effects in the model by coupling the mass balance equations for the tumor species and the nutrient dynamics to a mechanical equilibrium equation with phase-dependent elasticity coefficients. The resulting PDE system couples two Cahn--Hilliard type equations for the tumor phase and the healthy phase with a PDE linking the evolution of the interstitial fluid to the pressure of the system, a reaction-diffusion type equation for the nutrient proportion, and a quasistatic momentum balance. We prove here that the corresponding initial-boundary value problem has a solution in appropriate function spaces.

The paper treats the problem of optimal distributed control of a Cahn-Hilliard-Oono system in $\mathbb{R}^d$, $1\leq d\leq 3$, with the control located in the mass term and admitting general potentials that include both the case of a regular potential and the case of some singular potential. The first part of the paper is concerned with the dependence of the phase variable on the control variable. In the case of a logarithmic potential, we need to prove an ad hoc strict separation property, and for this reason we have to restrict ourselves to the case $d = 2$. In the rest of the work, we study the necessary first-order optimality conditions, which are proved under suitable compatibility conditions on the initial datum of the phase variable and the time derivative of the control, at least in case of potentials having unbounded domain. PLEASE NOTE: A revised version of this paper has been published in Discrete Contin. Dyn. Syst. Ser. S 15 (2022), 2135-2172: the authors point out that both the published version and the arXiv preprint contain a trivial mistake in formula (2.21), which leads to suitable changes in the assumption (2.10) and in some occurrences of it. Namely, (2.21) and (2.10) become $ \bar \varphi(t) = \bar \varphi_0 e^{-t} + \int_0^t e^{-(t-s)} \bar u(s) \, ds$ for every $t\in[0,T]$, and $-M - (\bar\varphi_0)^- , M + (\bar\varphi_0)^+$ belong to the interior of $D(\beta)$, respectively. The formula just after (2.21) becomes $-M - (\bar\varphi_0)^- \leq \bar \varphi(t) \leq M + (\bar\varphi_0)^+$ and a similar change must be done in (4.20). The first line of (2.24) is now implied by the new (2.21). Next, (2.39) becomes $-(M+ R) - (\bar\varphi_0)^- , (M+ R) + (\bar\varphi_0)^+$ belong to the interior of $D(\beta)$. Similar changes must be done in the interval written just two lines before (4.18) and in (4.18) itself., Comment: Key words: Cahn-Hilliard equation, logarithmic potential, regular potential, optimal control problem, strict separation property, necessary optimality conditions

Prostate cancer can be lethal in advanced stages, for which chemotherapy may become the only viable therapeutic option. While there is no clear clinical management strategy fitting all patients, cytotoxic chemotherapy with docetaxel is currently regarded as the gold standard. However, tumors may regain activity after treatment conclusion and become resistant to docetaxel. This situation calls for new delivery strategies and drug compounds enabling an improved therapeutic outcome. Combination of docetaxel with antiangiogenic therapy has been considered a promising strategy. Bevacizumab is the most common antiangiogenic drug, but clinical studies have not revealed a clear benefit from its combination with docetaxel. Here, we capitalize on our prior work on mathematical modeling of prostate cancer growth subjected to combined cytotoxic and antiangiogenic therapies, and propose an optimal control framework to robustly compute the drug-independent cytotoxic and antiangiogenic effects enabling an optimal therapeutic control of tumor dynamics. We describe the formulation of the optimal control problem, for which we prove the existence of at least a solution and determine the necessary first order optimality conditions. We then present numerical algorithms based on isogeometric analysis to run a preliminary simulation study over a single cycle of combined therapy. Our results suggest that only cytotoxic chemotherapy is required to optimize therapeutic performance and we show that our framework can produce superior solutions to combined therapy with docetaxel and bevacizumab. We also illustrate how the optimal drug-na\"{i}ve cytotoxic effects computed in these simulations may be successfully leveraged to guide drug production and delivery strategies by running a nonlinear least-square fit of protocols involving docetaxel and a new design drug.

In this paper a phase-field approach for structural topology optimization for a 3D-printing process which includes stress constraint and potentially multiple materials or multiscales is analyzed. First order necessary optimality conditions are rigorously derived and a numerical algorithm which implements the method is presented. A sensitivity study with respect to some parameters is conducted for a two-dimensional cantilever beam problem. Finally, a possible workflow to obtain a 3D-printed object from the numerical solutions is described and the final structure is printed using a fused deposition modeling (FDM) 3D printer.

We present a functionally graded material design (FGMD) approach relying on at opology optimization procedure based on asymptotic homogenization and phase-field method. We also introduce a complete framework from which numerical results lead to 3D printed structures. We firstly present numerical experiments to verify the proposed methodology and, subsequently, we discuss experimental measurements comparing optimized FGMD with a constant density structure having the same weight. Experimental evidence shows the effectiveness of the proposed methodology in improving the overall stiffness of optimized structures.

In the present work we introduce a novel graded-material design for additive manufacturing based on phase-field and topology optimization. The main novelty of this work comes from the introduction of an additional phase-field variable in the classical single-material phase-field topology optimization algorithm. This new variable is used to grade the material properties in a continuous fashion. Two different numerical examples are discussed, in both of them we perform sensitivity studies to asses the effects of different model parameters onto the resulting structure. From the presented results we can observe that the proposed algorithm adds additional freedom in the design, exploiting the higher flexibility coming from additive manufacturing technology.

In the present paper, a continuum model is introduced for fluid flow in a deformable porous medium, where the fluid may undergo phase transitions. Typically, such problems arise in modeling liquid-solid phase transformations in groundwater flows. The system of equations is derived here from the conservation principles for mass, momentum, and energy and from the Clausius-Duhem inequality for entropy. It couples the evolution of the displacement in the matrix material, of the capillary pressure, of the absolute temperature, and of the phase fraction. Mathematical results are proved under the additional hypothesis that inertia effects and shear stresses can be neglected. For the resulting highly nonlinear system of two PDEs, one ODE and one ordinary differential inclusion with natural initial and boundary conditions, existence of global in time solutions is proved by means of cut-off techniques and suitable Moser-type estimates.

In this paper we introduce a general abstract formulation of a variational thermomechanical model, by means of a unified derivation via a generalization of the principle of virtual powers for all the variables of the system, including the thermal one. In particular, choosing as thermal variable the entropy of the system, and as driving functional the internal energy, we get a gradient flow structure (in a suitable abstract setting) for the whole nonlinear PDE system. We prove a global in time existence of (weak) solutions result for the Cauchy problem associated to the abstract PDE system as well as uniqueness in case of suitable smoothness assumptions on the functionals.

We prove existence, uniqueness, regularity and separation properties for a nonlocal Cahn-Hilliard equation with a reaction term. We deal here with the case of logarithmic potential and degenerate mobility as well an uniformly lipschitz in $u$ reaction term $g(x,t,u).$

In the present contribution the sliding mode control (SMC) problem for a phase-field model of Caginalp type is considered. First we prove the well-posedness and some regularity results for the phase-field type state systems modified by the state- feedback control laws. Then, we show that the chosen SMC laws force the system to reach within finite time the sliding manifold (that we chose in order that one of the physical variables or a combination of them remains constant in time). We study three different types of feedback control laws: the first one appears in the internal energy balance and forces a linear combination of the temperature and the phase to reach a given (space dependent) value, while the second and third ones are added in the phase relation and lead the phase onto a prescribed target $phi^*$. While the control law is non-local in space for the first two problems, it is local in the third one, i.e., its value at any point and any time just depends on the value of the state.

We consider a full Navier-Stokes and $Q$-tensor system for incompressible liquid crystal flows of nematic type. In the two dimensional periodic case, we prove the existence and uniqueness of global strong solutions that are uniformly bounded in time. This result is obtained without any smallness assumption on the physical parameter $\xi$ that measures the ratio between tumbling and aligning effects of a shear flow exerting over the liquid crystal directors. Moreover, we show the uniqueness of asymptotic limit for each global strong solution as time goes to infinity and provide an uniform estimate on the convergence rate.

We consider a diffuse interface model for incompressible isothermal mixtures of two immiscible fluids with matched constant densities. This model consists of the Navier-Stokes system coupled with a convective nonlocal Cahn-Hilliard equation with non-constant mobility. We first prove the existence of a global weak solution in the case of non-degenerate mobilities and regular potentials of polynomial growth. Then we extend the result to degenerate mobilities and singular (e.g. logarithmic) potentials. In the latter case we also establish the existence of the global attractor in dimension two. Using a similar technique, we show that there is a global attractor for the convective nonlocal Cahn-Hilliard equation with degenerate mobility and singular potential in dimension three., Comment: 47 pages

We study a diffuse interface model for incompressible isothermal mixtures of two immiscible fluids coupling the Navier--Stokes system with a convective nonlocal Cahn--Hilliard equation in two dimensions of space. We apply recently proved well-posedness and regularity results in order to establish existence of optimal controls as well as first-order necessary optimality conditions for an associated optimal control problem in which a distributed control is applied to the fluid flow., Comment: 32 pages

In this paper we prove existence of global in time weak solutions for a highly nonlinear PDE system arising in the context of damage phenomena in thermoviscoelastic materials. The main novelty of the present contribution with respect to the ones already present in the literature consists in the possibility of taking into account a damage-dependent thermal expansion coefficient. This term implies the presence of nonlinear couplings in the PDE system, which make the analysis more challenging.

In this paper we study a distributed optimal control problem for a nonlocal convective Cahn-Hilliard equation with degenerate mobility and singular potential in three dimensions of space. While the cost functional is of standard tracking type, the control problem under investigation cannot easily be treated via standard techniques for two reasons: The state system is a highly nonlinear system of PDEs containing singular and degenerating terms, and the control variable, which is given by the velocity of the motion occurring in the convective term, is nonlinearly coupled to the state variable. The latter fact makes it necessary to state rather special regularity assumptions for the admissible controls, which, while looking a bit nonstandard, are however quite natural in the corresponding analytical framework. In fact, they are indispensable prerequisites to guarantee the well-posedness of the associated state system. In this contribution, we employ recently proved existence, uniqueness and regularity results for the solution to the associated state system in order to establish the existence of optimal controls and appropriate first-order necessary optimality conditions for the optimal control problem.

We consider a diffuse interface model of tumor growth proposed by A.~Hawkins-Daruud et al. This model consists of the Cahn-Hilliard equation for the tumor cell fraction $\varphi$ nonlinearly coupled with a reaction-diffusion equation for $\psi$, which represents the nutrient-rich extracellular water volume fraction. The coupling is expressed through a suitable proliferation function $p(\varphi)$ multiplied by the differences of the chemical potentials for $\varphi$ and $\psi$. The system is equipped with no-flux boundary conditions which entails the conservation of the total mass, that is, the spatial average of $\varphi+\psi$. Here we prove the existence of a weak solution to the associated Cauchy problem, provided that the potential $F$ and $p$ satisfy sufficiently general conditions. Then we show that the weak solution is unique and continuously depends on the initial data, provided that $p$ satisfies slightly stronger growth restrictions. Also, we demonstrate the existence of a strong solution and that any weak solution regularizes in finite time. Finally, we prove the existence of the global attractor in a phase space characterized by an a priori bounded energy., Comment: 31 pages

In this paper we study a distributed optimal control problem for a nonlocal convective Cahn--Hilliard equation with degenerate mobility and singular potential in three dimensions of space. While the cost functional is of standard tracking type, the control problem under investigation cannot easily be treated via standard techniques for two reasons: the state system is a highly nonlinear system of PDEs containing singular and degenerating terms, and the control variable, which is given by the velocity of the motion occurring in the convective term, is nonlinearly coupled to the state variable. The latter fact makes it necessary to state rather special regularity assumptions for the admissible controls, which, while looking a bit nonstandard, are however quite natural in the corresponding analytical framework. In fact, they are indispensable prerequisites to guarantee the well-posedness of the associated state system. In this contribution, we employ recently proved existence, uniqueness and regularity results for the solution to the associated state system in order to establish the existence of optimal controls and appropriate first-order necessary optimality conditions for the optimal control problem.

In this paper we prove the existence of global in time weak solutions for an evolutionary PDE system modelling nonisothermal Landau-de Gennes nematic liquid crystal (LC) flows in three dimensions of space. In our model, the incompressible Navier-Stokes system for the macroscopic velocity $\vu$ is coupled to a nonlinear convective parabolic equation describing the evolution of the Q-tensor $\QQ$, namely a tensor-valued variable representing the normalized second order moments of the probability distribution function of the LC molecules. The effects of the (absolute) temperature $\vt$ are prescribed in the form of an energy balance identity complemented with a global entropy production inequality. Compared to previous contributions, we can consider here the physically realistic singular configuration potential $f$ introduced by Ball and Majumdar. This potential gives rise to severe mathematical difficulties since it introduces, in the Q-tensor equation, a term which is at the same time singular in $\QQ$ and degenerate in $\vt$. To treat it a careful analysis of the properties of $f$, particularly of its blow-up rate, is carried out.

In this paper we derive, starting from the basic principles of Thermodynamics, an extended version of the nonconserved Penrose-Fife phase transition model, in which dynamic boundary conditions are considered in order to take into account interactions with walls. Moreover, we study the well-posedness and the asymptotic behavior of the Cauchy problem for the PDE system associated to the model, allowing the phase configuration of the material to be described by a singular function., Comment: 32 pages

In this paper we present a new thermodynamically consistent phase transition model describing the evolution of a liquid substance, e.g., water, in a rigid container $\Omega$ when we freeze the container. Since the density $\varrho_{2}$ of ice with volume fraction $\beta_{2}$, is lower than the density $\varrho_{1}$ of water with volume fraction $\beta_{1}$, experiments - for instance the freezing of a glass bottle filled with water - show that the water pressure increases up to the rupture of the bottle. When the container is not impermeable, freezing may produce a non-homogeneous material, for instance water ice or sorbet. Here we describe a general class of phase transition processes including this example as particular case. Moreover, we study the resulting nonlinear and singular PDE system from the analytical viewpoint recovering existence of a global (in time) weak solution and also uniqueness for some particular choices of the nonlinear functions involved.

We propose a model for non-isothermal phase transitions with non-conserved order parameter driven by a spatially nonlocal free energy with respect to both the temperature and the order parameter. The resulting system of equations is shown to be thermodynamically consistent and to admit a strong solution.

This work is concerned with the study of an initial and boundary value problem for a nonconserved system of phase field equations arising from the Penrose-Fife approach to phase transitions problems. Several works deal with variations of the same problem coupled with third type boundary conditions for the heat flux. On the contrary, our aim is to consider the case of the nonhomogeneous Neumann boundary condition for the heat flux, to find well-posedness for a weak formulation of this problem, and to prove a regularity result in case of smoother data and a slightly less general heat flux law.

Different systems for modeling tumor growth are presented. We follow the path of diffusive interface models and investigate these tumor models with respect to their well-posedness. Many biological phenomena, such as temporal and spatial nonlocal effects, complex nonlinearities, and mixed-dimensional couplings, are involved in mathematical oncology. As a result, detailed analysis of these complex systems is required, and we provide rigorous proofs for this. Verschiedene Systeme zur Modellierung von Tumorwachstum werden vorgestellt. Wir folgen dem Ansatz der diffusiven Grenzflächenmodelle und untersuchen diese Modelle auf deren mathematische Wohlgestelltheit. Viele biologische Phänomene wie zeitliche und räumliche nichtlokale Effekte, komplexe Nichtlinearitäten und gemischtdimensionale Kopplungen sind in der mathematischen Onkologie involviert. Daher ist eine detaillierte Analysis dieser Systeme erforderlich und wir liefern rigorose Beweise dafür.