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 1153 products
-
11/2007 - 07/26/2007
Badia, Santiago; Nobile, Fabio; Vergara, Christian
Fluid-structure partitioned procedures based on Robin transmission conditions | Abstract | | In this article we design new partitioned procedures for fluid-structure
interaction problems, based on Robin-type transmission conditions. The
choice of the coefficient in the Robin conditions is justified via simplified
models. The strategy is effective whenever an incompressible fluid interacts
with a relatively thin membrane, as in haemodynamics applications. We
analyze theoretically the new iterative procedures on a model problem,
which represents a simplified blood-vessel system. In particular, the Robin-
Neumann scheme exhibits enhanced convergence properties with respect to
the existing partitioned procedures. The theoretical results are checked
using numerical experimentation.
|
-
10/2007 - 04/28/2007
Parolini, Nicola; Quarteroni, Alfio
Modelling and Numerical Simulation for Yacht Engineering | Abstract | | In the past few years, Computational Fluid Dynamics (CFD) has become an essential tool in the design and optimization of racing sailboats and in particular America s Cup yachts.
The prevalent role of CFD in the design
process is demonstrated by the number of numerical simulations on different boat features, ranging from hull and appendage design to sail optimization, that each America s Cup syndicate carries on during the boat design and its further development.
In this work, we report some of the numerical results obtained in the framework of the research partnership between the Ecole Polytechnique Fédérale de Lausanne (EPFL) and the Alinghi Team, in preparation to the 32nd edition of the America s Cup.
A particular attention is devoted to the innovative aspects of the numerical models that have been recently developed. |
-
09/2007 - 04/21/2007
May, Caterina; Flournoy, Nancy
Asymptotics in response-adaptive designs generated by a two-colors, randomly reinforced urn | Abstract | | This paper illustrates asymptotic properties for a response-adaptive
design generated by a two-color, randomly reinforced urn model. The
design considered is optimal in the sense that it assigns patients to
the best treatment with probability converging to one. An approach
to show the joint asymptotic normality of the estimators of the mean
responses to the treatments is provided in spite of the fact that allocation
proportions converge to zero and one. Results on the rate of
convergence of the number of patients assigned to each treatment are
also obtained. Finally, we study the asymptotic behavior of a suitable
test statistic.
|
-
08/2007 - 04/14/2007
Formaggia, Luca; Miglio, Edie; Mola, Andrea; Parolini, Nicola
Fluid-structure interaction problems in free surface flows: application to boat dynamics | Abstract | | We present some recent studies on fluid-structure
interaction problems in the presence of free surface flow. We consider
the dynamics of boats simulated as rigid bodies. Several hydrodynamic
models are presented, ranging from full Reynolds Averaged
Navier-Stokes equations down to reduced models based on potential flow
theory.
|
-
07/2007 - 03/22/2007
Ern, Alexandre; Stephansen, Annette F.; Zunino, Paolo
A Discontinuous Galerkin method with weighted averages for advection-diffusion equations with locally vanishing and anisotropic diffusivity | Abstract | | We consider Discontinuous Galerkin approximations of advection-diffusion equations with anisotropic and discontinuous diffusivity, and propose the symmetric weighted interior penalty (SWIP) method for better coping with locally vanishing diffusivity. The analysis yields convergence results for the natural energy norm that are optimal (with respect to mesh-size) and robust (fully independent of the diffusivity). The convergence results for the advective derivative are optimal with respect to mesh-size and robust for isotropic diffusivity, as well as for anisotropic diffusivity in the dominant advection regime. In the dominant diffusivity regime, an optimal convergence result for the the $L^2$-norm is also recovered.
Numerical results are presented to illustrate the performance of the scheme.
|
-
06/2007 - 02/20/2007
Deponti, Alberto; Bonaventura, Luca; Fraccarollo, Luigi; Miglio, Edie; Rosatti, Giorgio
Analysis of Hyperbolic Systems for Mobile--Bed, Free--Surface Flow Modelling in Arbitrary Cross Sections | Abstract | | A model for mobile--bed river hydraulics based on the conservation equations of liquid mass, solid sediment mass and momentum is presented and analysed. The equations are reduced to one dimension by averaging over the cross section. These equations differ from the classical one--dimensional equations for mobile--bed, free--surface flows, since they take into account the effect of non--uniformities in velocity distribution along cross sections of arbitrary shape. By using appropriate closure formulae for sediment transport and for bottom friction, a system of three non--linear hyperbolic equations is obtained. Its eigenstructure is studied and its dependency on the closure formulae and on the non--dimensional parameters determining the flow and transport regimes is investigated. The existence and uniqueness of classical solutions of this hyperbolic system are discussed and an energy inequality for the frozen coefficient problem is derived.
|
-
05/2007 - 02/18/2007
Alì, G.; Culpo, M.; Micheletti, S.
Domain decomposition techniques for microelectronic modeling | Abstract | | This paper is meant to be the continuation of the previous work [1] where
a coupled ODE/PDE method for the simulation of semiconductor devices
was introduced. From a strictly mathematical viewpoint, analytical results
on coupled PDE/ODE systems (as arising in integrated circuit simulation)
can be found in [2]. In particular, in the present paper, we numerically
investigate an algorithm of Domain Decomposition type for the simulation
of circuits containing distributed devices (x 1) as well as semiconductors
in which some part is modeled with lumped parameters (x 2). It is worth
noticing the original employment of the Domain Decomposition technique
within the con_nes of a heterogeneous PDE/ODE coupling, versus its
typical use in a homogeneous full-PDE context. The results presented
here have been studied in the seminal work [3], while a more thorough
analysis is ongoing [4].
|
-
04/2007 - 02/15/2007
Nobile, Fabio; Tempone, Raul; Clayton G., Webster
An anisotropic sparse grid stochastic collocation method for elliptic partial differential equations with random input data | Abstract | | This work proposes and analyzes an anisotropic sparse grid stochastic collocation method for solving elliptic partial differential equations with random coefficients and forcing terms (input data of the model). The method consists of a Galerkin approximation in the space variables and
a collocation, in probability space, on sparse tensor product grids utilizing either Clenshaw-Curtis or Gaussian knots. Even in the presence of nonlinearities, the collocation approach leads to the solution of uncoupled
deterministic problems, just as in the Monte Carlo method. This work includes a priori and a posteriori procedures to adapt the anisotropy of
the sparse grids to each given problem. These procedures seem to be very effective for the problems under study. The proposed method combines
the advantages of isotropic sparse collocation with those of anisotropic full tensor product collocation: the first approach is effective for problems depending on random variables which weigh approximately equally in the solution, while the benefits of the latter approach become apparent when solving highly anisotropic problems depending on a relatively small number of random variables, as in the case where input random variables are
Karhunen-Loève truncations of “smooth” random fields. This work also provides a rigorous convergence analysis of the fully discrete problem and demonstrates: (sub)-exponential convergence in the asymptotic regime and algebraic convergence in the pre-asymptotic regime, with respect to the total number of collocation points. Numerical examples illustrate the theoretical results and are used to compare this approach with several others, including the standard Monte Carlo. In particular, for moderately large dimensional problems, the sparse grid approach with a properly chosen anisotropy seems to be very efficient and superior to all examined methods.
|
|