Researcher profile

Elisabetta Rocca

Elisabetta Rocca contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
7works
0followers
7topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

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