Researcher profile

Annalisa Quaini

Annalisa Quaini contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
13works
0followers
9topics
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

13 published item(s)

preprint2022arXiv

A data-driven Reduced Order Method for parametric optimal blood flow control: application to coronary bypass graft

We consider an optimal flow control problem in a patient-specific coronary artery bypass graft with the aim of matching the blood flow velocity with given measurements as the Reynolds number varies in a physiological range. Blood flow is modelled with the steady incompressible Navier-Stokes equations. The geometry consists in a stenosed left anterior descending artery where a single bypass is performed with the right internal thoracic artery. The control variable is the unknown value of the normal stress at the outlet boundary, which is need for a correct set-up of the outlet boundary condition. For the numerical solution of the parametric optimal flow control problem, we develop a data-driven reduced order method that combines proper orthogonal decomposition (POD) with neural networks. We present numerical results showing that our data-driven approach leads to a substantial speed-up with respect to a more classical POD-Galerkin strategy proposed in [59], while having comparable accuracy.

preprint2022arXiv

A novel Large Eddy Simulation model for the Quasi-Geostrophic Equations in a Finite Volume setting

We present a Large Eddy Simulation (LES) approach based on a nonlinear differential low-pass filter for the simulation of two-dimensional barotropic flows with under-refined meshes. For the implementation of such model, we choose a segregated three-step algorithm combined with a computationally efficient Finite Volume method. We assess the performance of our approach on the classical double-gyre wind forcing benchmark. The numerical experiments we present demonstrate that our nonlinear filter is an improvement over a linear filter since it is able to recover the four-gyre pattern of the time-averaged stream function even with extremely coarse meshes. In addition, our LES approach provides an average kinetic energy that compares well with the one computed with a Direct Numerical Simulation.

preprint2022arXiv

A POD-Galerkin reduced order model for the Navier-Stokes equations in stream function-vorticity formulation

We develop a Proper Orthogonal Decomposition (POD)-Galerkin based Reduced Order Model (ROM) for the efficient numerical simulation of the parametric Navier-Stokes equations in the stream function-vorticity formulation. Unlike previous works, we choose different reduced coefficients for the vorticity and stream function fields. In addition, for parametric studies we use a global POD basis space obtained from a database of time dependent full order snapshots related to sample points in the parameter space. We test the performance of our ROM strategy with the vortex merger benchmark. Accuracy and efficiency are assessed for both time reconstruction and physical parametrization.

preprint2022arXiv

Data-Driven Enhanced Model Reduction for Bifurcating Models in Computational Fluid Dynamics

We investigate various data-driven methods to enhance projection-based model reduction techniques with the aim of capturing bifurcating solutions. To show the effectiveness of the data-driven enhancements, we focus on the incompressible Navier-Stokes equations and different types of bifurcations. To recover solutions past a Hopf bifurcation, we propose an approach that combines proper orthogonal decomposition with Hankel dynamic mode decomposition. To approximate solutions close to a pitchfork bifurcation, we combine localized reduced models with artificial neural networks. Several numerical examples are shown to demonstrate the feasibility of the presented approaches.

preprint2021arXiv

A Comparison of Reduced-Order Modeling Approaches Using Artificial Neural Networks for PDEs with Bifurcating Solutions

This paper focuses on reduced-order models (ROMs) built for the efficient treatment of PDEs having solutions that bifurcate as the values of multiple input parameters change. First, we consider a method called local ROM that uses k-means algorithm to cluster snapshots and construct local POD bases, one for each cluster. We investigate one key ingredient of this approach: the local basis selection criterion. Several criteria are compared and it is found that a criterion based on a regression artificial neural network (ANN) provides the most accurate results for a channel flow problem exhibiting a supercritical pitchfork bifurcation. The same benchmark test is then used to compare the local ROM approach with the regression ANN selection criterion to an established global projection-based ROM and a recently proposed ANN based method called POD-NN. We show that our local ROM approach gains more than an order of magnitude in accuracy over the global projection-based ROM. However, the POD-NN provides consistently more accurate approximations than the local projection-based ROM.

preprint2021arXiv

A kinetic theory approach for 2D crowd dynamics with emotional contagion

In this paper, we present a computational modeling approach for the dynamics of human crowds, where the spreading of an emotion (specifically fear) has an influence on the pedestrians' behavior. Our approach is based on the methods of the kinetic theory of active particles. The model allows us to weight between two competing behaviors depending on fear level: the search for less congested areas and the tendency to follow the stream unconsciously (herding). The fear level of each pedestrian influences her walking speed and is influenced by the fear levels of her neighbors. Numerically, we solve our pedestrian model with emotional contagion using an operator splitting scheme. We simulate evacuation scenarios involving two groups of interacting pedestrians to assess how domain geometry and the details of fear propagation impact evacuation dynamics. Further, we reproduce the evacuation dynamics of an experimental study involving distressed ants.

preprint2021arXiv

An unfitted finite element method for two-phase Stokes problems with slip between phases

We present an isoparametric unfitted finite element approach of the CutFEM or Nitsche-XFEM family for the simulation of two-phase Stokes problems with slip between phases. For the unfitted generalized Taylor--Hood finite element pair $\mathbf{P}_{k+1}-P_k$, $k\ge1$, we show an inf-sup stability property with a stability constant that is independent of the viscosity ratio, slip coefficient, position of the interface with respect to the background mesh and, of course, mesh size. In addition, we prove stability and optimal error estimates that follow from this inf-sup property. We provide numerical results in two and three dimensions to corroborate the theoretical findings and demonstrate the robustness of our approach with respect to the contrast in viscosity, slip coefficient value, and position of the interface relative to the fixed computational mesh.

preprint2020arXiv

A microscopic approach to study the onset of a highly infectious disease spreading

We combine a pedestrian dynamics model with a contact tracing method to simulate the initial spreading of a highly infectious airborne disease in a confined environment. We focus on a medium size population (up to 1000 people) with a small number of infectious people (1 or 2) and the rest of the people is divided between immune and susceptible. We adopt a space-continuous model that represents pedestrian dynamics by the forces acting on them, i.e. a microscopic force-based model. Once discretized, the model results in a high-dimensional system of second order ordinary differential equations. Before adding the contact tracing to the pedestrian dynamics model, we calibrate the model parameters, validate the model against empirical data, and show that pedestrian self-organization into lanes can be captured. We consider an explicit approach for contact tracing by introducing a sickness domain around a sick person. A healthy but susceptible person who remains in the sickness domain for a certain amount of time may get infected (with a prescribed probability) and become a so-called secondary contact. As a concrete setting to simulate the onset of disease spreading, we consider terminals in two US airports: Hobby Airport in Houston and the Atlanta International Airport. We consider different scenarios and we quantify the increase in average number of secondary contacts increases as a given terminal becomes more densely populated, the percentage of immune people decreases, the number of primary contacts increases, and areas of high density (such as the boarding buses) are present.

preprint2020arXiv

A Reduced Order Modeling technique to study bifurcating phenomena: application to the Gross-Pitaevskii equation

We propose a computationally efficient framework to treat nonlinear partial differential equations having bifurcating solutions as one or more physical control parameters are varied. Our focus is on steady bifurcations. Plotting a bifurcation diagram entails computing multiple solutions of a parametrized, nonlinear problem, which can be extremely expensive in terms of computational time. In order to reduce these demanding computational costs, our approach combines a continuation technique and Newton's method with a Reduced Order Modeling (ROM) technique, suitably supplemented with a hyper-reduction method. To demonstrate the effectiveness of our ROM approach, we trace the steady solution branches of a nonlinear Schrödinger equation, called Gross-Pitaevskii equation, as one or two physical parameters are varied. In the two parameter study, we show that our approach is 60 times faster in constructing a bifurcation diagram than a standard Full Order Method.

preprint2020arXiv

Coupling kinetic theory approaches for pedestrian dynamics and disease contagion in a confined environment

The goal of this work is to study an infectious disease spreading in a medium size population occupying a confined environment. For this purpose, we consider a kinetic theory approach to model crowd dynamics in bounded domains and couple it to a kinetic equation to model contagion. The interactions of a person with other pedestrians and the environment are modeled by using tools of game theory. The pedestrian dynamics model allows to weight between two competing behaviors: the search for less congested areas and the tendency to follow the stream unconsciously in a panic situation. Each person in the system has a contagion level that is affected by their neighborhood. For the numerical solution of the coupled problem, we propose a numerical algorithm that at every time step solves one crowd dynamics problem and one contagion problem, i.e. with no subiterations between the two. We test our coupled model on a problem involving a small crowd walking through a corridor.

preprint2020arXiv

Experimental validation of a phase-field model to predict coarsening dynamics of lipid domains in multicomponent membranes

Membrane phase-separation is a mechanism that biological membranes often use to locally concentrate specific lipid species in order to organize diverse membrane processes. Phase separation has also been explored as a tool for the design of liposomes with heterogeneous and spatially organized surfaces. These "patchy" liposomes are promising platforms for delivery purposes, however their design and optimization through experimentation can be expensive and time-consuming. We developed a computationally efficient method based on the surface Cahn-Hilliard phase-field model to complement experimental investigations in the design of patchy liposomes. The method relies on thermodynamic considerations to set the initial state for numerical simulations. We show that our computational approach delivers not only qualitative pictures, but also accurate quantitative information about the dynamics of the membrane organization. In particular, the computational and experimental results are in excellent agreement in terms of raft area fraction, total raft perimeter over time and total number of rafts over time for two different membrane compositions (DOPC:DPPC with a 2:1 molar ratio with 20% Chol and DOPC:DPPC with a 3:1 molar ratio with 20% Chol). Thus, the computational phase-field model informed by experiments has a considerable potential to assist in the design of liposomes with spatially organized surfaces, thereby containing the cost and time required by the design process.

preprint2019arXiv

Numerical modelling of phase separation on dynamic surfaces

The paper presents a model of lateral phase separation in a two component material surface. The resulting fourth order nonlinear PDE can be seen as a Cahn-Hilliard equation posed on a time-dependent surface. Only elementary tangential calculus and the embedding of the surface in $\mathbb{R}^3$ are used to formulate the model, thereby facilitating the development of a fully Eulerian discretization method to solve the problem numerically. A hybrid method, finite difference in time and trace finite element in space, is introduced and stability of its semi-discrete version is proved. The method avoids any triangulation of the surface and uses a surface-independent background mesh to discretize the equation. Thus, the method is capable of solving the Cahn-Hilliard equation numerically on implicitly defined surfaces and surfaces undergoing strong deformations and topological transitions. We assess the approach on a set of test problems and apply it to model spinodal decomposition and pattern formation on colliding surfaces. Finally, we consider the phase separation on a sphere splitting into two droplets.

preprint2019arXiv

Reduced Basis Model Order Reduction for Navier-Stokes equations in domains with walls of varying curvature

We consider the Navier-Stokes equations in a channel with a narrowing and walls of varying curvature. By applying the empirical interpolation method to generate an affine parameter dependency, the offline-online procedure can be used to compute reduced order solutions for parameter variations. The reduced order space is computed from the steady-state snapshot solutions by a standard POD procedure. The model is discretised with high-order spectral element ansatz functions, resulting in 4752 degrees of freedom. The proposed reduced order model produces accurate approximations of steady-state solutions for a wide range of geometries and kinematic viscosity values. The application that motivated the present study is the onset of asymmetries (i.e., symmetry breaking bifurcation) in blood flow through a regurgitant mitral valve, depending on the Reynolds number and the valve shape. Through our computational study, we found that the critical Reynolds number for the symmetry breaking increases as the wall curvature increases.