Quaderni MOX
Pubblicazioni
del Laboratorio di Modellistica e Calcolo Scientifico MOX. I lavori riguardano prevalentemente il campo dell'analisi numerica, della statistica e della modellistica matematica applicata a problemi di interesse ingegneristico. Il sito del Laboratorio MOX è raggiungibile
all'indirizzo mox.polimi.it
Trovati 1349 prodotti
-
23/2018 - 10/04/2018
Benacchio, T.;Bonaventura,L.
A seamless extension of DG methods for hyperbolic problems to unbounded domains | Abstract | | We consider spectral discretizations of hyperbolic problems on unbounded domains using Laguerre basis functions. Taking as model problem the scalar advection equation, we perform a comprehensive stability analysis that includes strong collocation formulations, nodal and modal weak formulations, with either inflow or outflow boundary conditions, using either Gauss - Laguerre or Gauss - Laguerre - Radau quadrature nodes and based on either scaled Laguerre functions or scaled Laguerre polynomials. We show that some of these combinations give rise to intrinsically unstable schemes, while the combination of scaled Laguerre functions with Gauss - Laguerre - Radau nodes appears to be stable for both strong and weak formulations. We then show how a modal discretization approach for hyperbolic systems on an unbounded domain can be naturally and seamlessly coupled to a discontinuous finite element discretization on a finite domain. Examples of one dimensional hyperbolic systems are solved with the proposed domain decomposition technique. The errors obtained with the proposed approach are found to be small, enabling the use of the coupled scheme for the simulation of Rayleigh damping layers in the semi-infinite part. Energy errors and reflection ratios of the scheme in absorbing
wavetrains and single Gaussian signals show that a small number of nodes in the semi-infinite domain are sufficient to damp the waves. The theoretical insight and numerical results corroborate previous findings by the authors and establish the scaled Laguerre functions-based discretization as a flexible and efficient tool for absorbing layers as well as for the accurate simulation of waves in unbounded regions. |
-
22/2018 - 10/04/2018
Pegolotti, L.; Dede', L.; Quarteroni, A.
Isogeometric Analysis of the electrophysiology in the human heart: numerical simulation of the bidomain equations on the atria | Abstract | | We consider Isogeometric Analysis (IGA) for the numerical solution of the electrophysiology of the atria, which in this work is modeled by means of the bidomain equations on thin surfaces. First, we consider the bidomain equations coupled with the Roger-McCulloch ionic model on simple slabs. Here, our goal is to evaluate the effects of the spatial discretization by IGA and the use of different B-spline basis functions on the accuracy of the approximation, in particular regarding the accuracy of the front
velocity and the dispersion error. Specifically, we consider basis functions with high polynomial degree, p, and global high order continuity, Cp?1, in the computational domain: our results show that the use of such basis functions is beneficial to the accurate approximation of the solution. Then, we consider a realistic application of the bidomain equations coupled with the Courtemanche-Ramirez-Nattel ionic model on the two human atria, which are represented by means of two NURBS surfaces. |
-
21/2018 - 10/04/2018
Gervasio, P.; Dede', L.; Chanon, O.; Quarteroni, A.
Comparing Isogeometric Analysis and Spectral Element Methods: accuracy and spectral properties | Abstract | | In this paper, we carry out a systematic comparison between the theoretical properties of Spectral Element Methods (SEM) and NURBS-based Isogeometric Analysis (IGA) in the framework of the Galerkin method for the approximation of the Poisson problem. Our focus is on their convergence properties, the algebraic structure and the spectral properties of the corresponding discrete arrays (mass and stiffness matrices). We review the available theoretical results for these methods and verify them numerically by performing error analysis on the solution of the Poisson problem. Where theory is lacking, we use numerical investigation of the results to draw conjectures on the behavior of the corresponding theoretical laws in terms of the design parameters, such as the (mesh) element size, the local polynomial degree, the smoothness of the NURBS basis functions, the space dimension, and the total number of degrees of freedom involved in the computations. |
-
20/2018 - 25/03/2018
Bassi, C. ; Abbà, A.; Bonaventura L.; Valdettaro, L.
A priori tests of a novel LES approach to compressible variable density turbulence | Abstract | | We assess the viability of a recently proposed novel approach to LES for compressible variable density flows by means of a priori tests. The a priori tests have been carried out filtering a two-dimensional
DNS database of the classic lock-exchange benchmark. The tests confirm that additional terms should be accounted for in subgrid scale modeling of variable density flows, with respect to the terms usually considered in the traditional approach. Several alternatives for the modeling of these terms are assessed and discussed. |
-
19/2018 - 01/03/2018
Menghini, F.; Dede', L.; Quarteroni, A.
Variational Multiscale LES modeling of blood flow in an idealized left human heart | Abstract | | In this work we build a realistic, although idealized, computational model of the left human heart for the study of the blood flow dynamics. We prescribe the left heart wall displacement based on physiological data and we take into account the presence of both the mitral and aortic valves through a resistive method. We simulate the left heart hemodynamics by means of the Finite Element method and compare different numerical stabilization techniques to account for the transitional and nearly turbulent nature of the blood flow in a physiological regime. In particular, we apply a Variational Multiscale Large Eddy Simulation (LES) model and a Streamwise Upwind Petrov-Galerkin method and we critically analyze the corresponding numerical results. |
-
18/2018 - 26/02/2018
Antonietti, P.F.; Bonaldi, F.; Mazzieri, I.
A high-order discontinuous Galerkin approach to the elasto-acoustic problem | Abstract | | We address the spatial discretization of an evolution problem arising from the coupling of viscoelastic and acoustic wave propagation phenomena by employing a discontinuous Galerkin scheme on polygonal and polyhedral meshes. The coupled nature of the problem is ascribed to suitable transmission conditions imposed at the interface between the solid (elastic) domain and the fluid (acoustic) domain. We state and prove a well-posedness result for the strong formulation of the problem, present a stability analysis for the semi-discrete formulation, and finally prove an a priori hp-version error estimate for the resulting formulation in a suitable (mesh-dependent) energy norm. The conver- gence results are validated by numerical experiments carried out in a two-dimensional setting. |
-
17/2018 - 24/02/2018
Agosti, A.; Giverso, C.; Faggiano, E.;Stamm,A.; Ciarletta, P.
A personalized mathematical tool for neuro-oncology: a clinical case study | Abstract | | This work evaluates the predictive ability of a novel personalised computational tool for simulating the growth of brain tumours using the neuroimaging data collected during one clinical case study.
The mathematical model consists in an evolutionary fourth-order partial differential equation with degenerate motility, in which the spreading dynamics of the multiphase tumour is coupled through a growth term with a parabolic equation determining the diffusing oxygen within the brain. The model also includes a reaction term describing the effects of radiotherapy, that is simulated in accordance to the clinical schedule.
We collect Magnetic Resonance (MRI) and Diffusion Tensor (DTI) imaging data for one patient at given times of key clinical interest, from the first diagnosis of a giant glioblastoma to its surgical removal and the subsequent radiation therapies.
These neuroimaging data allow reconstructing the patient-specific brain geometry in a finite element virtual environment, that is used for simulating the tumour recurrence pattern after the surgical resection. In particular, we characterize the different brain tissues and the tumour location from MRI data, whilst we extrapolate the heterogeneous nutrient diffusion parameters and cellular mobility from DTI data.
The numerical results of the simulated tumour are found in good qualitative and quantitative agreement with the volume and the boundaries observed in MRI data. Moreover, the simulations point out a consistent regression of the tumour mass in correspondence to the application of radiotherapy, with an average growth rate which is of the same order as the one calculated from the neuroimaging data. Remarkably, our results display the highest Jaccard index of the tumour region reported in the biomathematical literature.
In conclusion, this work represents an important proof-of-concept of the ability of this mathematical framework to predict the tumour recurrence and its response to therapies in a patient-specific manner. |
-
16/2018 - 24/02/2018
Calissano, A.; Vantini, S.; Arnaboldi, M.
An elephant in the room: Twitter samplingmethodology. | Abstract | | The usage of social media data is spreading among the broad scientific community: 30000 papers dealing with this type of data are indexed in Scopus in the last decade. On the one hand, this data are very appealing, creating a rich bucket of information. On the other one, gathering them through a repeatable sampling strategy is increasing in complexity (or maybe becoming impossible?). The aim of this paper is to map the scientific community awareness about the sampling strategies used to download on-line data, focusing on the most studied social media: Twitter. This review unveils two unexpected results: the downloaded data are typically far from being randomly sampled, and around 99% of papers does not explicitly declare the sampling strategy used to download the data. These two facts pose some worrisome doubts about the trustworthiness of all the results presented in this stream of literature. |
|
