Source author record

Elisabetta Rocca

Elisabetta Rocca appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

41works
12topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

41 published item(s)

preprint2023arXiv

A Cahn-Hilliard phase field model coupled to an Allen-Cahn model of viscoelasticity at large strains

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.

preprint2022arXiv

A Cahn-Hilliard model coupled to viscoelasticity with large deformations

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.

preprint2022arXiv

On a Cahn-Hilliard-Keller-Segel model with generalized logistic source describing tumor growth

We propose a new type of diffuse interface model describing the evolution of a tumor mass under the effects of a chemical substance (e.g., a nutrient or a drug). The process is described by utilizing the variables $φ$, an order parameter representing the local proportion of tumor cells, and $σ$, representing the concentration of the chemical. The order parameter $φ$ is assumed to satisfy a suitable form of the Cahn-Hilliard equation with mass source and logarithmic potential of Flory-Huggins type (or generalizations of it). The chemical concentration $σ$ satisfies a reaction-diffusion equation where the cross-diffusion term has the same expression as in the celebrated Keller-Segel model. In this respect, the model we propose represents a new coupling between the Cahn-Hilliard equation and a subsystem of the Keller-Segel model. We believe that, compared to other models, this choice is more effective in capturing the chemotactic effects that may occur in tumor growth dynamics (chemically induced tumor evolution and consumption of nutrient/drug by tumor cells). Note that, in order to prevent finite time blowup of $σ$, we assume a chemical source term of logistic type. Our main mathematical result is devoted to proving existence of weak solutions in a rather general setting that covers both the two- and three- dimensional cases. Under more restrictive assumptions on coefficients and data, and in some cases on the spatial dimension, we prove various regularity results. Finally, in a proper class of smooth solutions we show uniqueness and continuous dependence on the initial data in a number of significant cases.

preprint2022arXiv

Well-posedness and optimal control for a Cahn-Hilliard-Oono system with control in the mass term

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 φ(t) = \bar φ_0 e^{-t} + \int_0^t e^{-(t-s)} \bar u(s) \, ds$ for every $t\in[0,T]$, and $-M - (\barφ_0)^- , M + (\barφ_0)^+$ belong to the interior of $D(β)$, respectively. The formula just after (2.21) becomes $-M - (\barφ_0)^- \leq \bar φ(t) \leq M + (\barφ_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φ_0)^- , (M+ R) + (\barφ_0)^+$ belong to the interior of $D(β)$. Similar changes must be done in the interval written just two lines before (4.18) and in (4.18) itself.

preprint2021arXiv

Mixed variational formulations for structural topology optimization based on the phase-field approach

We propose a variational principle combining a phase-field functional for structural topology optimization with a mixed (three-field) Hu-Washizu functional, then including directly in the formulation equilibrium, constitutive, and compatibility equations. The resulting mixed variational functional is then specialized to derive a classical topology optimization formulation (where the amount of material to be distributed is an \emph{a priori} assigned quantity acting as a global constraint for the problem) as well as a novel topology optimization formulation (where the amount of material to be distributed is minimized, hence with no pre imposed constraint for the problem). Both formulations are numerically solved by implementing a mixed finite element scheme, with the second approach avoiding the introduction of a global constraint, hence respecting the convenient local nature of the finite element discretization. Furthermore, within the proposed approach it is possible to obtain guidelines for settings proper values of phase-field-related simulation parameters and, thanks to the combined phase-field and Hu-Washizu rationale, a monolithic algorithm solution scheme can be easily adopted. An insightful and extensive numerical investigation results in a detailed convergence study and a discussion on the obtained final designs. The numerical results clearly highlight differences between the two formulations as well as advantages related to the monolithic solution strategy; numerical investigations address both two-dimensional and three-dimensional applications.

preprint2020arXiv

Optimal control of cytotoxic and antiangiogenic therapies on prostate cancer growth

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ï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.

preprint2019arXiv

Optimal control of stochastic phase-field models related to tumor growth

We study a stochastic phase-field model for tumor growth dynamics coupling a stochastic Cahn-Hilliard equation for the tumor phase parameter with a stochastic reaction-diffusion equation governing the nutrient proportion. We prove strong well-posedness of the system in a general framework through monotonicity and stochastic compactness arguments. We introduce then suitable controls representing the concentration of cytotoxic drugs administered in medical treatment and we analyze a related optimal control problem. We derive existence of an optimal strategy and deduce first-order necessary optimality conditions by studying the corresponding linearized system and the backward adjoint system.

preprint2016arXiv

Optimal control of treatment time in a diffuse interface model of tumor growth

We consider an optimal control problem for a diffuse interface model of tumor growth. The state equations couples a Cahn-Hilliard equation and a reaction-diffusion equation, which models the growth of a tumor in the presence of a nutrient and surrounded by host tissue. The introduction of cytotoxic drugs into the system serves to eliminate the tumor cells and in this setting the concentration of the cytotoxic drugs will act as the control variable. Furthermore, we allow the objective functional to depend on a free time variable, which represents the unknown treatment time to be optimized. As a result, we obtain first order necessary optimality conditions for both the cytotoxic concentration and the treatment time.

preprint2016arXiv

Solvability of an unsaturated porous media flow problem with thermomechanical interaction

A PDE system consisting of the momentum balance, mass balance, and energy balance equations for displacement, capillary pressure, and temperature as a model for unsaturated fluid flow in a porous viscoelastoplastic solid is shown to admit a solution under appropriate assumptions on the constitutive behavior. The problem involves two hysteresis operators accounting for plastic and capillary hysteresis.

preprint2015arXiv

A rate-independent gradient system in damage coupled with plasticity via structured strains

This contribution deals with a class of models combining isotropic damage with plasticity. We are inspired by It has been inspired by a work by Freddi and Royer-Carfagni, including the case where the inelastic part of the strain only evolves in regions where the material is damaged. The evolution both of the damage and of the plastic variable is assumed to be rate-independent. Existence of solutions is established in the abstract energetic framework elaborated by Mielke and coworkers.

preprint2015arXiv

Analysis of a diffuse interface model of multispecies tumor growth

We consider a diffuse interface model for tumor growth recently proposed in [Y. Chen, S.M. Wise, V.B. Shenoy, J.S. Lowengrub, A stable scheme for a nonlinear, multiphase tumor growth model with an elastic membrane, Int. J. Numer. Methods Biomed. Eng., 30 (2014), 726-754]. In this new approach sharp interfaces are replaced by narrow transition layers arising due to adhesive forces among the cell species. Hence, a continuum thermodynamically consistent model is introduced. The resulting PDE system couples four different types of equations: a Cahn-Hilliard type equation for the tumor cells (which include proliferating and dead cells), a Darcy law for the tissue velocity field, whose divergence may be different from 0 and depend on the other variables, a transport equation for the proliferating (viable) tumor cells, and a quasi-static reaction diffusion equation for the nutrient concentration. We establish existence of weak solutions for the PDE system coupled with suitable initial and boundary conditions. In particular, the proliferation function at the boundary is supposed to be nonnegative on the set where the velocity ${\bf u}$ satisfies ${\bf u}\cdotν>0$, where $ν$ is the outer normal to the boundary of the domain. We also study a singular limit as the diffuse interface coefficient tends to zero.

preprint2015arXiv

Asymptotic analyses and error estimates for a Cahn-Hilliard type phase field system modelling tumor growth

This paper is concerned with a phase field system of Cahn-Hilliard type that is related to a tumor growth model and consists of three equations in terms of the variables order parameter, chemical potential and nutrient concentration. This system has been investigated in the recent contributions arXiv:1401.5943 [math.AP] and arXiv:1501.07057 [math.AP] from the viewpoint of well-posedness, long time behavior and asymptotic convergence as two positive viscosity coefficients tend to zero at the same time. Here, we continue the analysis performed in arXiv:1501.07057 [math.AP] by showing two independent sets of results as just one of the coefficients tends to zero, the other remaining fixed. We prove convergence results, uniqueness of solutions to the two resulting limit problems, and suitable error estimates.

preprint2015arXiv

Generalized gradient flow structure of internal energy driven phase field systems

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.

preprint2015arXiv

On the strongly damped wave equation with constraint

A weak formulation for the so-called "semilinear strongly damped wave equation with constraint" is introduced and a corresponding notion of solution is defined. The main idea in this approach consists in the use of duality techniques in Sobolev-Bochner spaces, aimed at providing a suitable "relaxation" of the constraint term. A global in time existence result is proved under the natural condition that the initial data have finite "physical" energy.

preprint2015arXiv

Vanishing viscosities and error estimate for a Cahn-Hilliard type phase field system related to tumor growth

In this paper we perform an asymptotic analysis for two different vanishing viscosity coefficients occurring in a phase field system of Cahn-Hilliard type that was recently introduced in order to approximate a tumor growth model. In particular, we extend some recent results obtained in the preprint arXiv:1401.5943, letting the two positive viscosity parameters tend to zero independently from each other and weakening the conditions on the initial data in such a way as to maintain the nonlinearities of the PDE system as general as possible. Finally, under proper growth conditions on the interaction potential, we prove an error estimate leading also to the uniqueness result for the limit system.

preprint2014arXiv

"Entropic" solutions to a thermodynamically consistent PDE system for phase transitions and damage

In this paper we analyze a PDE system modelling (non-isothermal) phase transitions and damage phenomena in thermoviscoelastic materials. The model is thermodynamically consistent: in particular, no {\em small perturbation assumption} is adopted, which results in the presence of quadratic terms on the right-hand side of the temperature equation, only estimated in $L^1$. The whole system has a highly nonlinear character. We address the existence of a weak notion of solution, referred to as "entropic", where the temperature equation is formulated with the aid of an entropy inequality, and of a total energy inequality. This solvability concept reflects the basic principles of thermomechanics as well as the thermodynamical consistency of the model. It allows us to obtain \emph{global-in-time} existence theorems without imposing any restriction on the size of the initial data. We prove our results by passing to the limit in a time discretization scheme, carefully tailored to the nonlinear features of the PDE system (with its "entropic" formulation), and of the a priori estimates performed on it. Our time-discrete analysis could be useful towards the numerical study of this model.

preprint2014arXiv

A diffuse interface model for two-phase incompressible flows with nonlocal interactions and nonconstant mobility

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.

preprint2014arXiv

Damage processes in thermoviscoelastic materials with damage-dependent thermal expansion coefficients

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.

preprint2014arXiv

Existence of solutions to a two-dimensional model for nonisothermal two-phase flows of incompressible fluids

We consider a thermodynamically consistent diffuse interface model describing two-phase flows of incompressible fluids in a non-isothermal setting. This model was recently introduced in a previous paper of ours, where we proved existence of weak solutions in three space dimensions. Here, we aim at studying the mathematical properties of the model in the two-dimensional case. In particular, we can show existence of global in time strong solutions. Moreover, we can admit slightly more general conditions on some material coefficients of the system.

preprint2014arXiv

On a diffuse interface model of tumor growth

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 $φ$ nonlinearly coupled with a reaction-diffusion equation for $ψ$, which represents the nutrient-rich extracellular water volume fraction. The coupling is expressed through a suitable proliferation function $p(φ)$ multiplied by the differences of the chemical potentials for $φ$ and $ψ$. The system is equipped with no-flux boundary conditions which entails the conservation of the total mass, that is, the spatial average of $φ+ψ$. 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.

preprint2014arXiv

On a non-isothermal diffuse interface model for two-phase flows of incompressible fluids

We introduce a diffuse interface model describing the evolution of a mixture of two different viscous incompressible fluids of equal density. The main novelty of the present contribution consists in the fact that the effects of temperature on the flow are taken into account. In the mathematical model, the evolution of the macroscopic velocity is ruled by the Navier-Stokes system with temperature-dependent viscosity, while the order parameter representing the concentration of one of the components of the fluid is assumed to satisfy a convective Cahn-Hilliard equation. The effects of the temperature are prescribed by a suitable form of the heat equation. However, due to quadratic forcing terms, this equation is replaced, in the weak formulation, by an equality representing energy conservation complemented with a differential inequality describing production of entropy. The main advantage of introducing this notion of solution is that, while the thermodynamical consistency is preserved, at the same time the energy-entropy formulation is more tractable mathematically. Indeed, global-in-time existence for the initial-boundary value problem associated to the weak formulation of the model is proved by deriving suitable a-priori estimates and showing weak sequential stability of families of approximating solutions.

preprint2014arXiv

Optimal control for a phase field system with a possibly singular potential

In this paper we study a distributed control problem for a phase field system of Caginalp type with logarithmic potential. The main aim of this work would be to force the location of the diffuse interface to be as close as possible to a prescribed set. However, due to the discontinuous character of the cost functional, we have to approximate it by a regular one and, in this case, we solve the associated control problem and derive the related first order necessary optimality conditions.

preprint2014arXiv

Optimal distributed control of a nonlocal Cahn-Hilliard/Navier-Stokes system in 2D

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.

preprint2014arXiv

Optimal distributed control of a nonlocal convective Cahn-Hilliard equation by the velocity in 3D

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.

preprint2014arXiv

Sharp interface control in a Penrose-Fife model

In this paper we study a singular control problem for a system of PDEs describing a phase-field model of Penrose-Fife type. The main novelty of this contribution consists in the idea of forcing a sharp interface separation between the states of the system by using heat sources distributed in the domain and at the boundary. We approximate the singular cost functional with a regular one, which is based on the Legendre-Fenchel relations. Then, we obtain a regularized control problem for which we compute the first order optimality conditions using an adapted penalization technique. The proof of some convergence results and the passage to the limit in these optimality conditions lead to the characterization of the desired optimal controller.

preprint2013arXiv

A degenerating PDE system for phase transitions and damage

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.

preprint2013arXiv

Global weak solution and blow-up criterion of the general Ericksen-Leslie system for nematic liquid crystal flows

In this paper we investigate the three dimensional general Ericksen-Leslie (E--L) system with Ginzburg-Landau type approximation modeling nematic liquid crystal flows. First, by overcoming the difficulties from lack of maximum principle for the director equation and high order nonlinearities for the stress tensor, we prove existence of global-in-time weak solutions under physically meaningful boundary conditions and suitable assumptions on the Leslie coefficients, which ensures that the total energy of the E--L system is dissipated. Moreover, for the E--L system with periodic boundary conditions, we prove the local well-posedness of classical solutions under the so-called Parodi's relation and establish a blow-up criterion in terms of the temporal integral of both the maximum norm of the curl of the velocity field and the maximum norm of the gradient of the liquid crystal director field.

preprint2013arXiv

Nonisothermal nematic liquid crystal flows with the Ball-Majumdar free energy

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.

preprint2013arXiv

The Penrose-Fife phase-field model with dynamic boundary conditions

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.

preprint2012arXiv

Evolution of non-isothermal Landau-de Gennes nematic liquid crystals flows with singular potential

We discuss a 3D model describing the time evolution of nematic liquid crystals in the framework of Landau-de Gennes theory, where the natural physical constraints are enforced by a singular free energy bulk potential proposed by J.M. Ball and A. Majumdar. The thermal effects are present through the component of the free energy that accounts for intermolecular interactions. The model is consistent with the general principle of thermodynamics and mathematically tractable. We identify the a priori estimates for the associated system of evolutionary partial differential equations and construct global-in-time weak solutions for arbitrary physically relevant initial data.

preprint2012arXiv

On a 3D isothermal model for nematic liquid crystals accounting for stretching terms

The present contribution investigates the well-posedness of a PDE system describing the evolution of a nematic liquid crystal flow under kinematic transports for molecules of different shapes. More in particular, the evolution of the {\em velocity field} $\ub$ is ruled by the Navier-Stokes incompressible system with a stress tensor exhibiting a special coupling between the transport and the induced terms. The dynamic of the {\em director field} $\bd$ is described by a variation of a parabolic Ginzburg-Landau equation with a suitable penalization of the physical constraint $|\bd|=1$. Such equation accounts for both the kinematic transport by the flow field and the internal relaxation due to the elastic energy. The main aim of this contribution is to overcome the lack of a maximum principle for the director equation and prove (without any restriction on the data and on the physical constants of the problem) the existence of global in time weak solutions under physically meaningful boundary conditions on $\bd$ and $\ub$.

preprint2012arXiv

Trajectory attractors for the Sun-Liu model for nematic liquid crystals in 3D

In this paper we prove the existence of a trajectory attractor (in the sense of V.V. Chepyzhov and M.I. Vishik) for a nonlinear PDE system coming from a 3D liquid crystal model accounting for stretching effects. The system couples a nonlinear evolution equation for the director d (introduced in order to describe the preferred orientation of the molecules) with an incompressible Navier-Stokes equation for the evolution of the velocity field u. The technique is based on the introduction of a suitable trajectory space and of a metric accounting for the double-well type nonlinearity contained in the director equation. Finally, a dissipative estimate is obtained by using a proper integrated energy inequality. Both the cases of (homogeneous) Neumann and (non-homogeneous) Dirichlet boundary conditions for d are considered.

preprint2011arXiv

A nonlocal quasilinear multi-phase system with nonconstant specific heat and heat conductivity

In this paper, we prove the existence and global boundedness from above for a solution to an integrodifferential model for nonisothermal multi-phase transitions under nonhomogeneous third type boundary conditions. The system couples a quasilinear internal energy balance ruling the evolution of the absolute temperature with a vectorial integro-differential inclusion governing the (vectorial) phase-parameter dynamics. The specific heat and the heat conductivity are allowed to depend both on the order parameter and on the absolute temperature of the system, and the convex component of the free energy may or may not be singular. Uniqueness and continuous data dependence are also proved under additional assumptions.

preprint2011arXiv

On the long-time behavior of some mathematical models for nematic liquid crystals

A model describing the evolution of a liquid crystal substance in the nematic phase is investigated in terms of two basic state variables: the {\it velocity field} $\bu$ and the {\it director field} $\di$, representing the preferred orientation of molecules in a neighborhood of any point in a reference domain. After recalling a known existence result, we investigate the long-time behavior of weak solutions. In particular, we show that any solution trajectory admits a non-empty $ω$-limit set containing only stationary solutions. Moreover, we give a number of sufficient conditions in order that the $ω$-limit set contains a single point. Our approach improves and generalizes existing results on the same problem.

preprint2011arXiv

Well-posedness of an extended model for water-ice phase transitions

We propose an improved model explaining the occurrence of high stresses due to the difference in specific volumes during phase transitions between water and ice. The unknowns of the resulting evolution problem are the absolute temperature, the volume increment, and the liquid fraction. The main novelty here consists in including the dependence of the specific heat and of the speed of sound upon the phase. These additional nonlinearities bring new mathematical difficulties which require new estimation techniques based on Moser iteration. We establish the existence of a global solution to the corresponding initial-boundary value problem, as well as lower and upper bounds for the absolute temperature. Assuming constant heat conductivity, we also prove uniqueness and continuous data dependence of the solution.

preprint2009arXiv

A bottle in a freezer

We propose here a model for solidification of a liquid contents of an elastic bottle in a freezer. The main goal is to explain the occurrence of high stresses inside the bottle. As a by-product, we derive a formula for the undercooling coefficient in terms of the elasticity constants, latent heat, and the phase expansion coefficient. We investigate the well-posedness of the three-dimensional model: we prove the existence and uniqueness of a solution for the corresponding initial-boundary value problem which couples a PDE with an integrodifferential equation and an ordinary differential inclusion ruling the evolution of the phase parameter. Finally, we prove some results on the long time behavior of solutions.

preprint2009arXiv

A model for shape memory alloys with the possibility of voids

The paper is devoted to the study of a mathematical model for the thermomechanical evolution of metallic shape memory alloys. The main novelty of our approach consists in the fact that we include the possibility for these materials to exhibit voids during the phase change process. Indeed, in the engineering paper has been recently proved that voids may appear when the mixture is produced by the aggregations of powder. Hence, the composition of the mixture varies (under either thermal or mechanical actions) in this way: the martensites and the austenite transform into one another whereas the voids volume fraction evolves. The first goal of this contribution is hence to state a PDE system capturing all these modelling aspects in order then to establish the well-posedness of the associated initial-boundary value problem.

preprint2009arXiv

Phase separation in a gravity field

We prove here well-posedness and convergence to equilibria for the solution trajectories associated to a model for solidification of a liquid content of a rigid container in a gravity field. We observe that the gravity effects, which can be neglected without considerable changes of the process on finite time intervals, have a substantial influence on the long time behavior of the evolution system. Without gravity, we find a temperature interval, in which all phase distributions with a prescribed total liquid contents are admissible equilibria, while, under the influence of gravity, the only equilibrium distribution in a connected container consists in two pure phases separated by one plane interface perpendicular to the gravity force.

preprint2008arXiv

Solid liquid phase changes with different densities

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 $Ω$ when we freeze the container. Since the density $\varrho_{2}$ of ice with volume fraction $β_{2}$, is lower than the density $\varrho_{1}$ of water with volume fraction $β_{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.