MOX Reports
The preprint collection of the Laboratory for Modeling and Scientific Computation MOX. It mainly contains works on numerical
analysis and mathematical modeling applied to engineering problems. MOX web site is mox.polimi.it
Found 1237 products
-
17/2014 - 05/06/2014
Discacciati, M.; Gervasio, P.; Quarteroni, A.
Interface Control Domain Decomposition (ICDD) Method for Stokes-Darcy coupling | Abstract | | We propose the ICDD method to solve the coupling between Stokes and Darcy equations. According to this approach, the problem is formulated as an optimal control problem whose control variables are the traces of the velocity and the pressure on the internal boundaries of the subdomains that provide an overlapping decomposition of the original computational domain. A theoretical analysis is carried out, an efficient solution algorithm is proposed, and several numerical tests are implemented. Our results show the accuracy of the ICDD methos, its computational efficiency and robustness with respect to the different parameters involved (grid-size, polynomial degrees, permeability of the porous domain, thickness of the overlapping region). The ICDD approach is more versatile and easier to implement than the model based on the Beavers, Joseph and Saffman coupling conditions. |
-
16/2014 - 05/05/2014
Dede, L.; Jaggli, C.; Quarteroni, A.
Isogeometric numerical dispersion analysis for elastic wave propagation | Abstract | | In this paper, we carry out a numerical dispersion analysis for the linear elastodynamics equations approximated by means of NURBS-based Isogeometric Analysis in the framwork of the Galerkin method; specifically, we consider the analysis of harminic plane waves in an isotropic and homogeneous elastic medium. We compare and discuss the errors associated to the compressional and shear wave velocities and we provide the anisotropic curves for numerical approximations obtained by considerin B-splines and NURBS basis functions of different regularity, namely globally C^0- and C^(p-1)- continuous, being p the polynomial degree. We conclude our analysis by numerically simulating the seismic wave propagation in a sinusoidal shaped valley with discontinuous elastic parameters across an internal interface. |
-
15/2014 - 05/02/2014
Esfandiar, B.; Porta, G.; Perotto, S.; Guadagnini, A;
Anisotropic mesh and time step adaptivity for solute transport modeling in porous media | Abstract | | We assess the impact of space-time mesh adaptivity on the modeling of solute transport in porous media. This approach allows an automatic selection of both the spatial mesh and the time step on the basis of a suitable recovery-based error estimator. In particular, we deal with an anisotropic control of the spatial mesh. The solver coupled with the adaptive module deals with an advection-dispersion equation to model the transport of dissolved species, which are assumed to be convected by a Darcy flow field. The whole solution-adaptation procedure is assessed through two-dimensional numerical tests. A numerical convergence analysis of the spatial mesh adaptivity is first performed by considering a test-case with analytical solution. Then, we validate the space-time adaptive procedure by reproducing a set of experimental observations associated with solute transport in a homogeneous sand pack. The accuracy and the efficiency of the methodology are discussed and numerical results are compared with those associated with fixed uniform space-time discretizations. This assessment shows that the proposed approach is robust and reliable. In particular, it allows us to obtain a significant improvement of the simulation quality of the early solute arrivals times at the outlet of the medium. |
-
14/2014 - 04/23/2014
Dassi, F.; Formaggia, L.; Zonca, S.
Degenerate Tetrahedra Recovering | Abstract | | Many mesh generation or optimization algorithms could produce a low quality tetrahedral mesh, i.e. a mesh where the tetrahedra have very small solid or dihedral angles. In this paper, we propose a series of operations to recover these degenerate tetrahedra. In particular, we will focus on the standard shapes of these undesired mesh elements (sliver, cap, wedge and spade) and, for each of these configurations, we apply a suitable sequence of classical mesh modification procedures to get a higher quality mesh. The reliability of the proposed mesh optimization algorithm is numerically proved with some examples. |
-
13/2014 - 04/21/2014
Ballarin, F.; Manzoni, A.; Quarteroni, A.; Rozza, G.
Supremizer stabilization of POD-Galerkin approximation of parametrized Navier-Stokes equations | Abstract | | In this work we present a stable proper orthogonal decomposition (POD)-Galerkin approximation for parametrized steady Navier-Stokes equations. The stabilization is guaranteed by the use of supremizers solutions that enrich the reduced velocity space. Numerical results show that an equivalent inf-sup condition is fulfilled, yielding stability for both velocity and pressure. Our stability analysis is first carried out from a theoretical standpoint, then confirmed by numerical tests performed on a parametrized two-dimensional backward facing step flow. |
-
12/2014 - 03/14/2014
Fumagalli, I; Parolini, N.; Verani, M.
Shape optimization for Stokes flow: a reference domain approach | Abstract | | In this paper we analyze a shape optimization problem, with Stokes equations as the state problem, defined on a domain with a part of the boundary that is described as the graph of the control function. The state problem formulation is mapped onto a reference domain, which is independent of the control function, and the analysis is mainly led on such domain. The existence of an optimal control function is proved, and optimality conditions are derived. After the analytical inspection of the problem, finite element discretization is considered for both the control function and the state variables, and a priori convergence error estimates are derived. Numerical experiments assess the validity of the theoretical results. |
-
11/2014 - 02/28/2014
Taddei, T.; Quarteroni, A.; Salsa, S.
An offline-online Riemann solver for one-dimensional systems of conservation laws | Abstract | | In this paper we present a new technique based on two different phases, here called offline and online stages, for the solution of the Riemann problem for one-dimensional hyperbolic systems of conservation laws. After theoretically motivating our offline-online technique, we prove its effectiveness by means of two numerical examples. |
-
10/2014 - 02/17/2014
Antonietti, P.F.; Dedner, A.; Madhavan, P.; Stangalino, S.; Stinner, B.; Verani, M.
High order discontinuous Galerkin methods on surfaces | Abstract | | We derive and analyze high order discontinuous Galerkin methods for
second-order elliptic problems on implicitely defined surfaces in R^3. This is done by carefully adapting the unified discontinuous Galerkin framework of [D.N. Arnold, F. Brezzi, B. Cockburn, and L.D. Marini, SIAM J. Numer. Anal., 2002] on a triangulated surface approximating the smooth surface. We prove optimal error estimates in both a (mesh dependent) energy and L^2 norms. |
|
