2020-10-24T21:05:05Z
https://cimec.org.ar/ojs/index.php/cmm/oai
oai:ojs.www.cimec.org.ar:article/2379
2009-12-04T16:20:34Z
cmm:ART
nmb a2200000Iu 4500
dc
Approximating optimization problems over convex functions
Aguilera, Néstor E.
Morin, Pedro
Many problems of theoretical and practical interest involve finding an optimum over a family of convex functions. For instance, finding the projection on the convex functions in $H^k(\Omega)$, and optimizing functionals arising from some problems in economics.
In the continuous setting and assuming smoothness, the convexity constraints may be given locally by asking the Hessian matrix to be positive semidefinite, but in making discrete approximations two difficulties arise: the continuous solutions may be not smooth, and functions with positive semidefinite discrete Hessian need not be convex in a discrete sense.
Previous work has concentrated on non-local descriptions of convexity, making the number of constraints to grow super-linearly with the number of nodes even in dimension 2, and these descriptions are very difficult to extend to higher dimensions.
In this paper we propose a finite difference approximation using positive semidefinite programs and discrete Hessians, and prove convergence under very general conditions, even when the continuous solution is not smooth, working on any dimension, and requiring a linear number of constraints in the number of nodes.
Using semidefinite programming codes, we show concrete examples of approximations to problems in two and three dimensions.
Published: Numerische Mathematik 111 (2008), 1-34; DOI: 10.1007/s00211-008-0176-4.
Centro Internacional de Métodos Numéricos en Ingeniería
2008-04-09 10:25:49
Articles
application/pdf
https://cimec.org.ar/ojs/index.php/cmm/article/view/2379
Cuadernos de Matemática y Mecánica; 2008
eng
Copyright (c)
oai:ojs.www.cimec.org.ar:article/2695
2010-06-08T22:34:30Z
cmm:ART
nmb a2200000Iu 4500
dc
Regularity of the Schrödinger equation for the harmonic oscillator
Bongioanni, Bruno
Rogers, Keith M.
We consider the Schrödinger equation for the harmonic oscillator $i \partial_t u = Hu$, where $H = -\Delta + |x|^2$, with initial data in the Hermite-Sobolev space $H^{-s/2} L^2(\real^n)$. We obtain smoothing and maximal estimates and apply these to perturbations of the equation and almost everywhhere convergence problems.
Accepted: Arkiv för matematik
Centro Internacional de Métodos Numéricos en Ingeniería
2009-04-08 20:22:29
Articles
application/pdf
https://cimec.org.ar/ojs/index.php/cmm/article/view/2695
Cuadernos de Matemática y Mecánica; 2009
eng
Copyright (c)
oai:ojs.www.cimec.org.ar:article/2994
2011-09-26T01:17:03Z
cmm:ART
nmb a2200000Iu 4500
dc
End-Point Estimates for Iterated Commutators of Multilinear Singular Integrals
Pérez, Carlos
Pradolini, Gladis
Torres, Rodolfo H.
Trujillo-González, Rodrigo
Iterated commutators of multilinear Calderón-Zygmund operators and pointwise multiplication with functions in BMO are studied in products of Lebesgue spaces. Both strong type and weak end-point estimates are obtained, including weighted results involving the vectors weights of the multilinear Calderón-Zygmund
theory recently introduced in the literature. Some better than expected estimates for certain multilinear operators are presented too.
Centro Internacional de Métodos Numéricos en Ingeniería
2010-06-07 08:59:51
Articles
application/pdf
https://cimec.org.ar/ojs/index.php/cmm/article/view/2994
Cuadernos de Matemática y Mecánica; 2010
eng
Copyright (c)
oai:ojs.www.cimec.org.ar:article/457
2009-07-31T22:36:23Z
cmm:ART
nmb a2200000Iu 4500
"090731 2009 eng "
dc
A Well Behaved Quasi-distance for Spaces of Homogeneous Type
Macías, Roberto A.
Segovia, Carlos
For any space of homogeneous type a quasi-distance equivalent to the original one is obtained satisfying that, if B and B' are balls such the center of the B' belongs to B and the radius of B' is smaller than the radius of B then, the measure of the intersection of B and B' is smaller than a constant fraction of the measure of the B'. An application to weighted norm inequalities for Hardy-Littlewood maximal function, which extends a result of A. P. Calderón, is given.
Centro Internacional de Métodos Numéricos en Ingeniería
2012-02-21 18:12:01
Articles
application/pdf
https://cimec.org.ar/ojs/index.php/cmm/article/view/457
Cuadernos de Matemática y Mecánica; 1980
eng
Copyright (c)
oai:ojs.www.cimec.org.ar:article/2533
2009-09-07T19:30:44Z
cmm:ART
nmb a2200000Iu 4500
dc
Convergence of finite elements adapted for weaker norms
Morin, Pedro
Siebert, Kunibert G.
Veeser, Andreas
We consider finite elements that are adapted to a (semi)norm that is weaker than the one of the trial space. We establish convergence of the finite element solutions to the exact one under the following conditions: refinement relies on unique quasi-regular element subdivisions and generates locally quasi-uniform grids; the finite element spaces are conforming, nested, and satisfy the inf-sup condition; the error estimator is reliable and appropriately locally efficient; the indicator of a non-marked element is bounded by the estimator contribution associated with the marked elements, and each marked element is subdivided at least once. This abstract convergence result is illustrated by two examples.
Keywords: Adaptivity, conforming finite elements, convergence, weaker (semi)norms, mesh-dependent norms
Published: Applied and Industrial Matematics in Italy II, Selected Contributions from the 8th SIMAI Conference, Vincenzo Cutello, Giorgio Fotia, Luigia Puccio, eds. Series on Advances in Mathematics for Applied Sciences - Vol. 75. World Scientific, 2007.
Centro Internacional de Métodos Numéricos en Ingeniería
2007-03-06 10:29:57
Articles
application/pdf
https://cimec.org.ar/ojs/index.php/cmm/article/view/2533
Cuadernos de Matemática y Mecánica; 2007
eng
Copyright (c)
oai:ojs.www.cimec.org.ar:article/2978
2011-09-26T01:17:03Z
cmm:ART
nmb a2200000Iu 4500
dc
Simulation of free surface flows by a finite element interface capturing technique
Battaglia, Laura
Storti, Mario Alberto
D'Elía, Jorge
Transient free surface flows are numerically simulated by a finite element interface capturing method based on a level set approach. The methodology consists of the solution of two-fluid viscous incompressible flows for a single domain, where the liquid phase is identified by positive values of the level set function, the gaseous one by negative ones, and the free surface by the zero level set. The numerical solution at each time step is performed in three stages: (i) a two-fluid Navier--Stokes stage, (ii) an advection stage for the transport of the level set function, and (iii) a bounded reinitialization with continuous penalization stage for keeping smoothness of the level set function. The proposed procedure, and particularly the renormalization stage, are evaluated in two typical two and three-dimensional problems. [Submitted to Int. J Comp Fluid Dynamics, ISSN 1061-8562]
Centro Internacional de Métodos Numéricos en Ingeniería
2010-06-07 08:59:51
Articles
application/pdf
https://cimec.org.ar/ojs/index.php/cmm/article/view/2978
Cuadernos de Matemática y Mecánica; 2010
eng
Copyright (c)
oai:ojs.www.cimec.org.ar:article/230
2009-08-14T21:35:03Z
cmm:ART
nmb a2200000Iu 4500
dc
A finite element method for surface diffusion: the parametric case
Morin, Pedro
Bänsch, Eberhard
Nochetto, Ricardo H.
Abstract: Surface diffusion is a (4th order highly nonlinear) geometric driven motion of a surface with normal velocity proportional to the surface Laplacian of mean curvature. We present a novel variational formulation for parametric surfaces with or without boundaries. The method is semi-implicit, requires no explicit parametrization, and yields a linear system of elliptic PDE to solve at each time step. We next develop a finite element method, propose a Schur complement approach to solve the resulting linear systems, and show several significant simulations, some with pinch-off in finite time. We introduce a mesh regularization algorithm, which helps prevent mesh distortion, and discuss the use of time and space adaptivity to increase accuracy while reducing complexity.
Keywords: Surface diffusion, fourth-order parabolic problem, finite elements, Schur complement, smoothing effect, pinch-off.
AMS Subject Classifications: 35K55, 65M12, 65M15, 65M60, 65Z05.
Published: Journal of Computational Physics 203 (2005) 321--343.
Published: Journal of Computational Physics 203 (2005) 321--343.
Centro Internacional de Métodos Numéricos en Ingeniería
2006-08-12 12:06:41
Articles
application/pdf
https://cimec.org.ar/ojs/index.php/cmm/article/view/230
Cuadernos de Matemática y Mecánica; 2003
eng
Copyright (c)
oai:ojs.www.cimec.org.ar:article/1391
2010-06-08T22:34:30Z
cmm:ART
nmb a2200000Iu 4500
"090408 2009 eng "
dc
Aeroelastic study of the start-up of a rocket engine nozzle
Garelli, Luciano
Paz, Rodrigo Rafael
Storti, Mario Alberto
The aim of this paper is to analyze the aeroelastic processes developed during the starting phase of a rocket engine via a coupling fluid/structure code. This analysis gives a better understanding of the behavior of the structure as the shock waves propagate inside the engine nozzle. The gasdynamics Euler equations are solved for the fluid, and constitutive linear elastic solid under the assumption of large displacements and rotations with no material damping is adopted for the structure. The coupling of each subproblem is carried out with a Gauß-Seidel algorithm over the fluid and structure states. For the fluid problem an ALE (Arbitrary Lagrangian-Eulerian Formulation) formulation is used. It allows to define a reference system following the moving boundaries while the structure is deformed. The code is validated with a study of the flutter phenomena that may occur when a supersonic compressible fluid flows over a flat solid plate. Regarding the rocket engine ignition problem, a modal analysis of the structure is performed in order to analyze the eigenfrequency shifts when considering the coupling with the fluid flow. [Submitted to Computers and Fluids 2009]
Centro Internacional de Métodos Numéricos en Ingeniería
2009-04-08 20:22:29
Articles
application/pdf
https://cimec.org.ar/ojs/index.php/cmm/article/view/1391
Cuadernos de Matemática y Mecánica; 2009
eng
Copyright (c)
oai:ojs.www.cimec.org.ar:article/2393
2010-06-08T22:34:30Z
cmm:ART
nmb a2200000Iu 4500
dc
LDR a package for likelihood-based sufficient dimension reduction
Cook, Dennis
Forzani, Liliana
Tomassi, Diego
We introduce a software package running under Matlab that implements several recently proposed likelihood-based methods for sufficient dimension reduction. Current capabilities include estimation of reduced subspaces with a fixed dimension d, as well as estimation of d by use of likelihood-ratio testing, permutation testing and information criteria. The methods are suitable for preprocessing data for both regression and classification. Implementations of related estimators are also available. Although the software is more oriented to command-line operation, a graphical user interface is also provided for prototype computations.
Centro Internacional de Métodos Numéricos en Ingeniería
2009-04-08 20:22:29
Articles
application/pdf
https://cimec.org.ar/ojs/index.php/cmm/article/view/2393
Cuadernos de Matemática y Mecánica; 2009
eng
Copyright (c)
oai:ojs.www.cimec.org.ar:article/2963
2010-06-08T22:34:30Z
cmm:ART
nmb a2200000Iu 4500
dc
Weighted local BMO spaces and the local Hardy-Littlewood maximal operator
Chicco Ruiz, Aníbal
Harboure, Eleonor O.
Centro Internacional de Métodos Numéricos en Ingeniería
2009-04-08 20:22:29
Articles
application/pdf
https://cimec.org.ar/ojs/index.php/cmm/article/view/2963
Cuadernos de Matemática y Mecánica; 2009
eng
Copyright (c)
oai:ojs.www.cimec.org.ar:article/3
2007-02-21T01:20:37Z
cmm:ART
nmb a2200000Iu 4500
"060624 2006 eng "
dc
Sloshing in a Multi-Physics Parallel Programming Paradigm
Battaglia, Laura
D'Elía, Jorge
Storti, Mario Alberto
Nigro, Norberto Marcelo
In this work, transient free surface flows of a viscous incompressible fluid are numerically solved through parallel computation. Transient free surface flows are boundary-value problems of moving type that involve geometrical non-linearities. In contrast to more conventional CFD problems, the computational flow domain is partially bounded by a free surface which is not known a priori, since its shape must be computed as part of the solution. In steady-flow the free surface is obtained by an iterative process, but when the free surface evolves with time the problem is more difficult, generating large distortions in the computational flow domain. The incompressible Navier-Stokes numerical solver is based on the finite element method with equal order elements for pressure and velocity (linear elements), and it uses a Streamline Upwind/Petrov-Galerkin (SUPG) scheme combined with a Pressure-Stabilizing/Petrov-Galerkin (PSPG) one. At each time step, the fluid equations are solved with constant pressure and null viscous traction conditions at the free surface and the velocities obtained in this way are used for updating the positions of the surface nodes. Then, a pseudo elastic problem is solved in the fluid domain in order to relocate the interior nodes so as to keep mesh distortion controlled. This has been implemented in PETSc-FEM by running two parallel instances of the code and exchanging information between them. Some numerical examples are presented. [Journal of Applied Mechanics (JAM-ASME) vol 73, pp. 1017-1025 (2006)]
Centro Internacional de Métodos Numéricos en Ingeniería
2006-06-24 09:54:16
Articles
application/pdf
https://cimec.org.ar/ojs/index.php/cmm/article/view/3
Cuadernos de Matemática y Mecánica; 2006
eng
Copyright (c)
oai:ojs.www.cimec.org.ar:article/1318
2009-12-04T16:20:34Z
cmm:ART
nmb a2200000Iu 4500
dc
An h-Adaptive Unstructured Mesh Refinement Strategy for Unsteady Problems
Ríos Rodríguez, Gustavo
Nigro, Norberto Marcelo
Storti, Mario Alberto
An h-adaptive unstructured mesh refinement strategy to solve unsteady problems by the finite element method is presented. The maximum level of refinement for the mesh is prescribed beforehand. The core operation of the strategy, namely the elements refinement, is described in detail. It is shown through numerical tests that one of the advantages of the chosen refinement procedure is to keep bounded the decrease of the mesh's quality. The type of element is not changed and no transition templates are used, therefore hanging nodes appear in the adapted mesh. The 1-irregular nodes refinement constraint is applied and the refinement process driven by this criterion is recursive. Both the strength and weakness of the adaptivity algorithm are mentioned, based on clock time measures and implementation issues. To show the proper working of the strategy, an axisymmetric, compressible non-viscous starting flow in a bell-shaped nozzle is solved over an unstructured mesh of hexaedra. [Submitted for publication to Latin American Applied Research]
Centro Internacional de Métodos Numéricos en Ingeniería
2008-04-09 10:25:49
Articles
application/pdf
https://cimec.org.ar/ojs/index.php/cmm/article/view/1318
Cuadernos de Matemática y Mecánica; 2008
eng
Copyright (c)
oai:ojs.www.cimec.org.ar:article/2384
2009-09-07T19:30:44Z
cmm:ART
nmb a2200000Iu 4500
dc
Maximal function estimates for the parabolic mean value kernel
Aimar, Hugo A.
Gómez, Ivana
Iaffei, Bibiana
We obtain parabolic and one-sided maximal function estimates for
nonisotropic dilations of the mean value kernel for the heat
equation.
Published: Revista Matemática Complutense 21 (2008), no. 2, 519–527.
Links:
http://www.mat.ucm.es/serv/revista/
http://revistas.ucm.es/mat/11391138/articulos/REMA0808220519A.PDF
Centro Internacional de Métodos Numéricos en Ingeniería
2007-03-06 10:29:57
Articles
application/pdf
https://cimec.org.ar/ojs/index.php/cmm/article/view/2384
Cuadernos de Matemática y Mecánica; 2007
eng
Copyright (c)
oai:ojs.www.cimec.org.ar:article/2700
2009-10-02T12:27:30Z
cmm:ART
nmb a2200000Iu 4500
dc
Geometry of Robinson joint consistency in Lukasiewicz logic
Busaniche, Manuela
Mundici, D.
We establish the Robinson joint consistency theorem for the infinitevalued propositional logic of Lukasiewicz. As a corollary we easily obtain the amalgamation property for MV-algebras the algebras of Lukasiewicz logic: all pre-existing proofs of this latter result make essential use of Pierce amalgamation theorem for abelian lattice-ordered groups (with strong unit) together with the categorical equivalence $\Gamma$ between these groups and MV-algebras. Our main tools are elementary and geometric.
Published: Annals of Pure and Applied Logic, 147 (2007), 1-22.
Centro Internacional de Métodos Numéricos en Ingeniería
2007-03-06 10:29:57
Articles
text/plain
https://cimec.org.ar/ojs/index.php/cmm/article/view/2700
Cuadernos de Matemática y Mecánica; 2007
eng
Copyright (c)
oai:ojs.www.cimec.org.ar:article/3000
2012-02-21T00:04:35Z
cmm:ART
nmb a2200000Iu 4500
"110215 2011 eng "
dc
On Haar Bases for Generalized Dyadic Hardy Spaces
Aimar, Hugo
Bernardis, Ana
Nowak, Luis
In this note we prove that Haar type systems are unconditional basis in the generalized dyadic Hardy space in the setting of spaces of homogeneous type. As a consequence, we obtain an alternative proof of the
unconditionality of such basis in Lebesgue spaces on spaces of homogeneous type.
Centro Internacional de Métodos Numéricos en Ingeniería
2012-02-20 21:06:46
Articles
application/pdf
https://cimec.org.ar/ojs/index.php/cmm/article/view/3000
Cuadernos de Matemática y Mecánica; 2011
eng
Copyright (c)
oai:ojs.www.cimec.org.ar:article/877
2007-04-27T03:25:57Z
cmm:ART
nmb a2200000Iu 4500
dc
Numerical simulations of axisymmetric inertial waves in a rotating sphere by finite elements
D'Elia, Jorge
Nigro, Noberto
Storti, Mario Alberto
Axisymmetric inertial waves of a viscous fluid that fill a perturbed
rotating spherical container are numerically simulated by finite
elements. A laminar flow of an incompressible viscous fluid of
Newtonian type is assumed in the numerical simulations. A monolithic
computational code is employed, which is based on stabilized finite
elements by means of a Streamline Upwind Petrov Galerkin (SUPG) and
Pressure Stabilized Petrov Galerkin
(PSPG) composed scheme. The Reynolds number is fixed as 50 000, while the ranges of the Rossby and Ekman
numbers are 0.2 Ro 1 and 2 × 10-5 Ek 10-4 , respectively. Some flow visualizations are performed. The
pressure coefficient spectrum at the centre of the sphere is plotted as a function of the frequency ratio and some
resonant frequencies are identified. The position of these resonant frequencies are in good agreement with previous
experimental and analytical ones in the inviscid limit.
[To appear in International Journal of Computational Fluid Dynamics.Author Posting. (c) Taylor & Francis, 2007. ]
Centro Internacional de Métodos Numéricos en Ingeniería
2006-06-24 09:54:16
Articles
application/pdf
https://cimec.org.ar/ojs/index.php/cmm/article/view/877
Cuadernos de Matemática y Mecánica; 2006
eng
Copyright (c)
oai:ojs.www.cimec.org.ar:article/2373
2010-06-08T22:34:30Z
cmm:ART
nmb a2200000Iu 4500
dc
AFEM for the Laplace-Beltrami Operator on Graphs: Design and Conditional Contraction Property
Mekchay, Khamron
Morin, Pedro
Nochetto, Ricardo H.
We present an adaptive finite element method (AFEM) of any polynomial degree for the Laplace-Beltrami operator on $C^1$ graphs $\Gamma$ in $\R^d$, $(d\ge2)$. We first derive residual-type a posteriori error estimates that account for the interaction of both the energy error in $H^1(\Gamma)$ and the surface error in $W^1_\infty(\Gamma)$ due to approximation of $\Gamma$. We devise a marking strategy to reduce the total error estimator, namely a suitably scaled sum of the energy, geometric, and inconsistency error estimators. We prove a conditional contraction property for the sum of the energy error and the total estimator; the conditional statement encodes resolution of $\Gamma$ in $W^1_\infty$. We conclude with one numerical experiment that illustrates the theory.
Centro Internacional de Métodos Numéricos en Ingeniería
2009-04-08 20:22:29
Articles
application/pdf
https://cimec.org.ar/ojs/index.php/cmm/article/view/2373
Cuadernos de Matemática y Mecánica; 2009
eng
Copyright (c)
oai:ojs.www.cimec.org.ar:article/2538
2009-08-15T13:46:05Z
cmm:ART
nmb a2200000Iu 4500
"090815 2009 eng "
dc
Local Problems on Stars: A Posteriori Error Estimators, Convergence, and Performance
Morin, Pedro
Nochetto, Ricardo H.
Siebert, Kunibert G.
A new computable a posteriori error estimator is introduced, which relies on the solution of small discrete problems on stars. It exhibits built-in flux equilibration and is equivalent to the energy error up to data oscillation without any saturation assumption. A simple adaptive strategy is designed, which simultaneously reduces error and data oscillation, and is shown to converge without mesh preadaptation nor explicit knowledge of constants. Numerical experiments reveal a competitive performance, show extremely good effectivity indices, and yield quasi-optimal meshes.
Keywords: A posteriori error estimators, local problems, stars, data oscillation, adaptivity, convergence, performance
AMS Subject Classifications: 65N12, 65N15, 65N30, 65N50, 65Y20
Published: Mathematics of Computation 72 (2003), 1067-1097.
Centro Internacional de Métodos Numéricos en Ingeniería
2006-06-24 11:50:44
Articles
application/pdf
https://cimec.org.ar/ojs/index.php/cmm/article/view/2538
Cuadernos de Matemática y Mecánica; 2001
eng
Copyright (c)
oai:ojs.www.cimec.org.ar:article/2985
2011-09-26T01:17:03Z
cmm:ART
nmb a2200000Iu 4500
dc
Global Saturation of Regularization Methods for Inverse Ill-Posed Problems
Herdman, Terry
Spies, Ruben Daniel
Temperini, Karina Guadalupe
In this article the concept of saturation of an arbitrary regularization method is formalized based upon the original idea of saturation for spectral regularization methods introduced by Neubauer in 1994 (A. Neubauer, "On converse and saturation results for regularization methods", In BeitrÄage zur angewandten Analysis und Informatik, pp. 262-270. Shaker, Aachen, 1994). Necessary and sufficient conditions for a regularization method to have global saturation are provided. It is shown that for a method to have global saturation the total error must be optimal in two senses, namely as optimal order of convergence over a certain set which at the same time, must be optimal (in a very precise sense) with respect to the error. Finally, two converse results are proved and the theory is applied to find suffcient conditions which ensure the existence of global saturation for spectral methods with classical qualifcation of finite positive order and for methods with maximal qualifcation. Finally, several examples of regularization methods possessing global saturation are shown.
Centro Internacional de Métodos Numéricos en Ingeniería
2010-06-07 08:59:51
Articles
application/pdf
https://cimec.org.ar/ojs/index.php/cmm/article/view/2985
Cuadernos de Matemática y Mecánica; 2010
eng
Copyright (c)
oai:ojs.www.cimec.org.ar:article/456
2009-07-30T17:24:12Z
cmm:ART
nmb a2200000Iu 4500
dc
ON Riesz Transforms and Maximal Functions in the Context of Gaussian Harmonic Analysis
Aimar, Hugo A.
Forzani, Liliana
Scotto, Roberto
The purpose of this paper is twofold. We introduce a general maximal function on the Gaussian setting which dominates the Ornstein-Uhlenbeck maximal operator and prove its weak type (1; 1) by using a covering lemma which is halfway between Besicovitch and Wiener. On the other hand, by taking as starting point the generalized Cauchy-Riemann equations, we introduce a new class of Gaussian Riesz Transforms. We prove, using the maximal function de¯ned in the ¯rst part of the paper, that unlike the ones already studied these new Riesz Transforms are weak type (1; 1) independently of their orders.
Centro Internacional de Métodos Numéricos en Ingeniería
2006-06-24 12:11:10
Articles
application/pdf
https://cimec.org.ar/ojs/index.php/cmm/article/view/456
Cuadernos de Matemática y Mecánica; 2005
eng
Copyright (c)
oai:ojs.www.cimec.org.ar:article/2399
2009-12-04T16:20:34Z
cmm:ART
nmb a2200000Iu 4500
dc
The Mahler measure of linear forms as special values of solutions of algebraic differential equations
Toledano, Ricardo D.
We prove that for each $n\geq 4$ there is an analytic function
$F_n(x)$ satisfying an algebraic differential equation of degree
$n+1$ such that the logarithmic Mahler measure of the linear form
$\lf_n=x_1+\cdots + x_n$ can be essentially computed as the evaluation of $F_n(z)$ at $z=n^{-1}$. We
show that the coefficients of the series representing $F_n(z)$ can
be computed recursively using the $n$-th. symmetric power of a
second order linear algebraic differential equation and we give
an estimate on the growth of these coefficients.
Published: Rocky Mountain Journal of Mathematics. Vol.39, 4, 2009, 1323-1338.
Link: http://rmmc.asu.edu/rmj/rmj.html
Centro Internacional de Métodos Numéricos en Ingeniería
2008-04-09 10:25:49
Articles
application/pdf
https://cimec.org.ar/ojs/index.php/cmm/article/view/2399
Cuadernos de Matemática y Mecánica; 2008
eng
Copyright (c)
oai:ojs.www.cimec.org.ar:article/2971
2010-06-08T22:34:30Z
cmm:ART
nmb a2200000Iu 4500
dc
An Algebraic Composite Finite Element Mesh Method
Sarraf, Sofía
López, Ezequiel Javier
Sonzogni, Victorio Enrique
Bergallo, Marta Beatriz
The composite ﬁnite element mesh method is useful for discretization error estimation and, in addition, for solution improvement with no increment in the computational cost. The technique consists in redeﬁne over a given mesh the linear operator that arises from the discretization of a partial diﬀerential equation. This operator is modiﬁed according to an appropriate linear combination between the operators of the given mesh and of a coarse mesh, which must be a coarsening of the ﬁrst one. In this work a novel algebraic composite mesh technique is proposed. The technique uses some tools from the Algebraic Multigrid method for the deﬁnition of the coarse mesh and the discrete space associated with it. Mesh coarsening is based on the fusion of elements in macroelements, with a new deﬁnition of the grid topology and basis functions. The agglomeration of elements is made in order to reduce the mesh anisotropy, which is of importance in the discretization of convection-diﬀusion-reaction problems. The discrete operator for the coarser mesh is obtained by the Galerkin Coarse Approximation, where inter-grid transfer operators are obtained using the graph of the coarse mesh. Several test problems with diﬀerent boundary conditions are presented.
Centro Internacional de Métodos Numéricos en Ingeniería
2009-04-08 20:22:29
Articles
application/pdf
https://cimec.org.ar/ojs/index.php/cmm/article/view/2971
Cuadernos de Matemática y Mecánica; 2009
eng
Copyright (c)
oai:ojs.www.cimec.org.ar:article/8
2009-07-30T17:24:12Z
cmm:ART
nmb a2200000Iu 4500
dc
MPI for Python
Dalcín, Lisandro Daniel
Paz, Rodrigo Rafael
Storti, Mario Alberto
MPI for Python provides bindings of the Message Passing Interface (MPI) standard for the Python programming language and allows any Python program to exploit multiple processors. This package is constructed on top of the MPI1 specification and defines an ob ject oriented interface which closely follows MPI2 C++ bin ings. It supports point-to-point (sends, receives) and collective (broadcasts, scatters, gathers) communications of general Python ob jects. Efficiency has been tested in a Beowulf class cluster and satisfying results were obtained. MPI for Python is open source and available for download on the web (http://www.cimec.org.ar/python) [Journal of Parallel and Distributed Computing, vol 65/9, pp. 1108-1115 (2005)]
Centro Internacional de Métodos Numéricos en Ingeniería
2006-06-24 12:11:10
Articles
application/pdf
https://cimec.org.ar/ojs/index.php/cmm/article/view/8
Cuadernos de Matemática y Mecánica; 2005
eng
Copyright (c)
oai:ojs.www.cimec.org.ar:article/1386
2009-12-04T16:20:34Z
cmm:ART
nmb a2200000Iu 4500
dc
Doubling Property on Hutchinson Orbits of Some Families of Contractive Similitudes
Aimar, Hugo
Carena, Marilina
Iaffei, Bibiana
We are interested in the behavior of the dynamical system generated by successive applications of Hutchinson similitudes starting from a metric measure space. We prove that, in the case of families of similitudes with the same contraction ratio, even when no point of the orbit is a doubling space, a gradual doubling property is taking place and the limit point recovers the homogeneity property.
Centro Internacional de Métodos Numéricos en Ingeniería
2008-04-09 10:25:49
Articles
application/pdf
https://cimec.org.ar/ojs/index.php/cmm/article/view/1386
Cuadernos de Matemática y Mecánica; 2008
eng
Copyright (c)
oai:ojs.www.cimec.org.ar:article/2389
2010-06-08T22:34:30Z
cmm:ART
nmb a2200000Iu 4500
dc
Boundedness of fractional operators in weighted variable exponent spaces with non doubling measures
Gorosito, Osvaldo
Pradolini, Gladis
Salinas, Oscar
In the context of variable exponent Lebesgue spaces equipped with a lower Ahlfors measure we obtain
weighted norm inequalities over bounded domains for the
centered fractional maximal function
and the fractional integral operator.
Centro Internacional de Métodos Numéricos en Ingeniería
2009-04-08 20:22:29
Articles
application/pdf
https://cimec.org.ar/ojs/index.php/cmm/article/view/2389
Cuadernos de Matemática y Mecánica; 2009
eng
Copyright (c)
oai:ojs.www.cimec.org.ar:article/2706
2010-06-08T22:34:30Z
cmm:ART
nmb a2200000Iu 4500
dc
Study of canonicity in subvarieties of BL-algebras
Busaniche, Manuela
Cabrer, Leonardo
We prove that every subvariety of BL-algebras which is not finitely generated is not $\sigma$-canonical. We also prove $\pi$-canonicity for an infinite family of subvarieties of BL-algebras that are not finitely generated. To do so we study the behavior of canonical extensions of ordered sums of posets.
Accepted: Algebra Universalis.
Centro Internacional de Métodos Numéricos en Ingeniería
2009-04-08 20:22:29
Articles
text/plain
https://cimec.org.ar/ojs/index.php/cmm/article/view/2706
Cuadernos de Matemática y Mecánica; 2009
eng
Copyright (c)
oai:ojs.www.cimec.org.ar:article/4035
2012-02-28T20:24:07Z
cmm:ART
nmb a2200000Iu 4500
"120228 2012 eng "
dc
Smoothness improvement for temperatures in terms of the Besov regularity of initial and Dirichlet data
Aimar, Hugo A.
Gómez, Ivana
Jerison and Kenig in J. Funct. Anal. 130 (1995), no.1, 161-219, gave a precise region $\mathcal{R}$ in the square $[0,1]^2$ for the pairs $(s,\tfrac{1}{p})$ for which every harmonic function in the Lipschitz domain $D$, with Dirichlet data in $B^s_p(\partial D)$, belongs to $B^{s+\tfrac{1}{p}}_p(D)$. We prove that every temperature $u$ in $\Omega=D\times (0,T)$ belongs to $\mathbb{B}^{\alpha}_{\tau}(\Omega)$ with $\tfrac{1}{\tau}=\tfrac{1}{p}+\tfrac{\alpha}{d}$, $0
Centro Internacional de Métodos Numéricos en Ingeniería
2012-02-28 17:07:11
Articles
application/pdf
https://cimec.org.ar/ojs/index.php/cmm/article/view/4035
Cuadernos de Matemática y Mecánica; 2012
eng
Copyright (c)
oai:ojs.www.cimec.org.ar:article/1089
2009-09-07T19:30:44Z
cmm:ART
nmb a2200000Iu 4500
dc
Approximation of spaces of homogeneous type
Aimar, Hugo
Carena, Marilina
Iaffei, Bibiana
In this note we combine the dyadic families introduced by M. Christ and the discrete partitions introduced by J. M. Wu in order to get approximation of a compact space of homogeneous type by a uniform sequence of finite spaces of homogeneous type. The convergence holds in the sense of a metric built on the Hausdorff distance between sets and on the Kantorovich-Rubinshtein metric between measures.
Centro Internacional de Métodos Numéricos en Ingeniería
2007-03-06 10:29:57
Articles
application/pdf
https://cimec.org.ar/ojs/index.php/cmm/article/view/1089
Cuadernos de Matemática y Mecánica; 2007
eng
Copyright (c)
oai:ojs.www.cimec.org.ar:article/2380
2009-12-04T16:20:34Z
cmm:ART
nmb a2200000Iu 4500
dc
Convergence rates for adaptive finite elements
Gaspoz, Fernando D.
Morin, Pedro
In this article we prove that it is possible to construct, using newest-vertex bisection, meshes that equidistribute the error in $H^1$-norm, whenever the function to approximate can be decomposed as a sum of a regular part plus a singular part with singularities around a finite number of points. This decomposition is usual in regularity results of Partial Differential Equations (PDE). As a consequence, the meshes turn out to be quasi-optimal, and convergence rates for adaptive finite element methods (AFEM) using Lagrange finite elements of any polynomial degree are obtained.
Published: IMA Journal of Numerical Analysis 2008; doi: 10.1093/imanum/drn039
Link: http://imajna.oxfordjournals.org/cgi/reprint/drn039?ijkey=aULHrx3AxOoNGTm&keytype=ref
Centro Internacional de Métodos Numéricos en Ingeniería
2008-04-09 10:25:49
Articles
application/pdf
https://cimec.org.ar/ojs/index.php/cmm/article/view/2380
Cuadernos de Matemática y Mecánica; 2008
eng
Copyright (c)
oai:ojs.www.cimec.org.ar:article/2696
2009-10-02T12:13:06Z
cmm:ART
nmb a2200000Iu 4500
dc
Free algebras in varieties of BL-algebras generated by a chain
Busaniche, Manuela
In this paper we describe finitely generated free algebras in varieties of BL-algebras generated by one BL-chain which is an ordinal sum of a finite MV-chain and a generalized BL-chain. We also give some particular examples of these free algebras.
Published: Algebra Universalis 50 (2003), 259-277.
Centro Internacional de Métodos Numéricos en Ingeniería
2006-08-12 12:06:41
Articles
https://cimec.org.ar/ojs/index.php/cmm/article/view/2696
Cuadernos de Matemática y Mecánica; 2003
eng
Copyright (c)
oai:ojs.www.cimec.org.ar:article/2996
2011-09-26T01:17:03Z
cmm:ART
nmb a2200000Iu 4500
dc
A DGCL Compliant FEM Formulation Based on Averaged Jacobians
Storti, Mario Alberto
Garelli, Luciano
Paz, Rodrigo Rafael
In this article a new methodology for developing DGCL (for Discrete Geometric Conservation Law) compliant formulations is presented. It is carried out in the context of the Finite Element Method (FEM) for general advective-diffusive systems on moving domains using an Arbitrary Lagrangian Eulerian (ALE) scheme. There is an extensive literature about the impact of DGCL compliance on the stability and precision of time integration methods. In those articles it has been proved that satisfying the DGCL is a necessary and sufﬁcient condition for any ALE scheme to maintain on moving grids the nonlinear stability properties of its ﬁxed-grid counterpart. However, only a few works propose a methodology for obtaining a compliant scheme. In this work, a DGCL compliant scheme based on an Averaged ALE Jacobians Formulation (AJF) is obtained. This new formulation is applied to the theta-family of time integration methods. In addition, an extension to the three-point Backward Difference Formula (BDF) is given. With the aim to validate the AJF formulation a set of numerical tests are performed. These tests include 2D and 3D diffusion problems with different mesh movements, and the 2D Euler ﬂow over a pitching NACA0012 airfoil.
Centro Internacional de Métodos Numéricos en Ingeniería
2010-06-07 08:59:51
Articles
application/pdf
https://cimec.org.ar/ojs/index.php/cmm/article/view/2996
Cuadernos de Matemática y Mecánica; 2010
eng
Copyright (c)
oai:ojs.www.cimec.org.ar:article/462
2009-07-30T17:24:12Z
cmm:ART
nmb a2200000Iu 4500
dc
Oscillation and variation for the Gaussian Riesz transforms and Poisson integral
Harboure, Eleonor
Macías, Roberto A.
Menárguez, Trinidad
Torrea, José L.
For the family of truncations of the Gaussian Riesz transforms and Poisson integral we study their rate of convergence through the oscillation and variation operators.
More precisely, we search for their L^p(dgamma)-boundedness properties, being dgamma the Gauss
measure. We achieve our results by looking at the oscillation and variation operators from a vector valued point of view.
Centro Internacional de Métodos Numéricos en Ingeniería
2006-06-24 12:11:10
Articles
application/pdf
https://cimec.org.ar/ojs/index.php/cmm/article/view/462
Cuadernos de Matemática y Mecánica; 2005
eng
Copyright (c)
oai:ojs.www.cimec.org.ar:article/1658
2009-04-27T18:39:38Z
cmm:ART
nmb a2200000Iu 4500
dc
A Lagrangian Panel Method in the Time Domain for Moving Free-surface Potential Flows
D'Elía, Jorge
Storti, Mario Alberto
Oñate, Eugenio
Idelsohn, Sergio R.
A Lagrangian-type panel method, in the time domain, is proposed for potential flows with a moving free surface. After a spatial semi.discretization, with a low-order scheme, the instantaneous velocity-potential and normal displacement on the moving free surface, are obtained by eans of a time-marching scheme. The kinematic and dynamic boundary conditions at the free surface are non-linear restrictions over the related Ordinary Differential Equation system and, in order to handle them an alternative Stekhlov-Poincaré operator technique is proposed. The method is applied to sloshing like flow problems.
Centro Internacional de Métodos Numéricos en Ingeniería
2006-08-10 14:21:22
Articles
application/pdf
https://cimec.org.ar/ojs/index.php/cmm/article/view/1658
Cuadernos de Matemática y Mecánica; 2002
eng
Copyright (c)
oai:ojs.www.cimec.org.ar:article/2534
2009-09-07T19:30:44Z
cmm:ART
nmb a2200000Iu 4500
dc
A basic convergence result for conforming adaptive finite element methods
Morin, Pedro
Siebert, Kunibert G.
Veeser, Andreas
We consider the approximate solution with adaptive finite elements of a class of linear boundary value problems, which includes problems of `saddle point' type. For the adaptive algorithm we suppose the following framework: refinement relies on unique quasi-regular element subdivisions and generates locally quasi-uniform grids, the finite element spaces are conforming, nested, and satisfy the inf-sup conditions, the error estimator is reliable as well as locally and discretely efficient, and marked elements are subdivided at least once. Under these assumptions, we give a sufficient and essentially necessary condition on marking for the convergence of the finite element solutions to the exact one. This condition is not only satisfied by Dörfler's strategy, but also by the maximum strategy and the equidistribution strategy.
Keywords: Adaptivity, conforming finite elements, convergence
Published: Mathematical Models and Methods in Applied Sciences (M3AS) 18 (2008) 707--737.
Centro Internacional de Métodos Numéricos en Ingeniería
2007-03-06 10:29:57
Articles
application/pdf
https://cimec.org.ar/ojs/index.php/cmm/article/view/2534
Cuadernos de Matemática y Mecánica; 2007
eng
Copyright (c)
oai:ojs.www.cimec.org.ar:article/2979
2011-09-26T01:17:03Z
cmm:ART
nmb a2200000Iu 4500
dc
A Semi-Analytical Computation of the Kelvin Kernel for Potential Flows with a Free Surface
D'Elía, Jorge
Storti, Mario Alberto
Battaglia, Laura
Abstract. A semi-analytical computation of the three dimen sional Green function for seakeeping ﬂow problems is proposed. A potential ﬂow model is assumed with an harmonic dependence in time and a linearized free-surface boundary condition. The multiplicative Green function is expressed as the product of a time and a spatial parts. The spatial part is known as the Kelvin kernel, which is the sum of two Rankine sources and a wave-like kernel, being the last one written using the Haskind- Havelock representation. Numerical eﬃciency is improved by an analytical integration of the two Rankine kernels and the use of a singularity subtractive technique for the Haskind-Havelock integral, where a globally adaptive quadrature is performed for the regular part and an analytic integration is used for the singular one. The proposed computation is employed in a low order panel method with ﬂat triangular elements. As a numerical example, an oscillating ﬂoating unit hemisphere in heave and surge modes is considered, where analytical and semi-analytical solutions are taken as a reference.[Submitted to Computational and Applied Mathematics ISSN: 0101-8205]
Centro Internacional de Métodos Numéricos en Ingeniería
2010-06-07 08:59:51
Articles
application/pdf
https://cimec.org.ar/ojs/index.php/cmm/article/view/2979
Cuadernos de Matemática y Mecánica; 2010
eng
Copyright (c)
oai:ojs.www.cimec.org.ar:article/231
2009-08-14T21:33:03Z
cmm:ART
nmb a2200000Iu 4500
"060817 2006 eng "
dc
Discrete Gradient Flows for Shape Optimization and Applications
Morin, Pedro
Dogan, Gunay
Nochetto, Ricardo H.
Verani, Marco
We present a variational framework for shape optimization problems that establishes clear and explicit connections among the continuous formulation, its full discretization and the resulting linear algebraic systems. Our approach hinges on the following essential features: shape differential calculus, a semi-implicit time discretization and a finite element method for space discretization. We use shape differential calculus to express variations of bulk and surface energies with respect to domain changes. The semi-implicit time discretization allows us to track the domain boundary without an explicit parametrization, and has the flexibility to choose different descent directions by varying the scalar product used for the computation of normal velocity. We propose a Schur complement approach to solve the resulting linear systems efficiently. We discuss applications of this framework to image segmentation, optimal shape design for PDE, and surface diffusion, along with the choice of suitable scalar products in each case. We illustrate the method with several numerical experiments, some developing pinch-off and topological changes in finite time.
Keywords: Shape optimization, scalar product, gradient flow, semi-implicit discretization, finite elements, surface diffusion, image segmentation.
Published: Computer Methods in Applied Mechanics and Engineering 196 (2007), 3898--3914.
Centro Internacional de Métodos Numéricos en Ingeniería
2006-06-24 12:11:10
Articles
application/pdf
https://cimec.org.ar/ojs/index.php/cmm/article/view/231
Cuadernos de Matemática y Mecánica; 2005
eng
Copyright (c)
oai:ojs.www.cimec.org.ar:article/2394
2009-12-04T16:20:34Z
cmm:ART
nmb a2200000Iu 4500
dc
Covariance reducing models: An alternative to spectral modelling of covariance matrices
Cook, Dennis
Forzani, Liliana
We introduce a software package running under Matlab that implements several recently proposed likelihood-based methods for sufficient dimension reduction. Current capabilities include estimation of reduced subspaces with a fixed dimension d, as well as estimation of d by use of likelihood-ratio testing, permutation testing and information criteria. The methods are suitable for preprocessing data for both regression and classification. Implementations of related estimators are also available. Although the software is more oriented to command-line operation, a graphical user interface is also provided for prototype computations.
Published: Biometrika 95 (2008), no. 4, 799--812
Centro Internacional de Métodos Numéricos en Ingeniería
2008-04-09 10:25:49
Articles
application/pdf
https://cimec.org.ar/ojs/index.php/cmm/article/view/2394
Cuadernos de Matemática y Mecánica; 2008
eng
Copyright (c)
oai:ojs.www.cimec.org.ar:article/2964
2010-06-08T22:34:30Z
cmm:ART
nmb a2200000Iu 4500
dc
Composition of Fractional Orlicz maximal operators and A1-weights on spaces of homogeneous type
Bernardis, Ana L.
Lorente, María
Pradolini, Gladis G.
Riveros, María Silvina
For a Young function $\Theta$ and $0\leq \alpha
Centro Internacional de Métodos Numéricos en Ingeniería
2009-04-08 20:22:29
Articles
application/pdf
https://cimec.org.ar/ojs/index.php/cmm/article/view/2964
Cuadernos de Matemática y Mecánica; 2009
eng
Copyright (c)
oai:ojs.www.cimec.org.ar:article/4
2007-02-21T01:01:17Z
cmm:ART
nmb a2200000Iu 4500
dc
A Minimal Element Distortion Strategy for Computational Mesh Dynamics
López, Ezequiel
Nigro, Norberto Marcelo
Storti, Mario Alberto
Toth, Jorge
Mesh motion strategy is one of the key points in many
fluid-structure interaction (FSI) problems. Due to the increasing
application of FSI to solve the current challenging engineering
problems this topic has deserved a highlight interest. There are
several different strategies to solve this problem, some of them using
a discrete and lumped spring-mass system to propagate the boundary
motion into the volume mesh and many others using an elastostatic
problem to deform the mesh. In all these strategies there is always a
risk of producing an invalid mesh, a mesh with some elements inverted.
Normally this condition is irreversible and once an invalid mesh is
obtained it is difficult to continue. In this paper the mesh motion
strategy is defined as an optimization problem. By its definition this
strategy may be classified as a particular case of an elastostatic
problem where the material constitutive law is defined in terms of the
minimization of certain energy functional that takes into account the
degree of element distortion. Some advantages of this strategy is its
natural tendency to high quality meshes, its robustness and its
straightforward extension to 3D problems. Several examples included in
this paper show these capabilities. Even though this strategy seems
to be very robust it is not able to recover a valid mesh starting from
an invalid one. This improvement is left for future work.[International Journal for Numerical Methods in Engineering
(accepted)]
Centro Internacional de Métodos Numéricos en Ingeniería
2006-06-24 09:54:16
Articles
application/pdf
https://cimec.org.ar/ojs/index.php/cmm/article/view/4
Cuadernos de Matemática y Mecánica; 2006
eng
Copyright (c)
oai:ojs.www.cimec.org.ar:article/1319
2009-12-04T16:20:34Z
cmm:ART
nmb a2200000Iu 4500
dc
A CVBEM Formulation for Multiple Profiles and Cascades
D'Elía, Jorge
Storti, Mario Alberto
Idelsohn, Sergio
A numerical algorithm based on the Complex Variable Boundary Element Method (CVBEM) for plane incompress- ible potential flow around aerofoils and cascades is described. The method is based on the representation of the complex disturbance velocity by means of a Cauchy-type integral around the foil. The Cauchy density function is approximated piecewise linearly and a linear system on the nodal values is obtained by collocation at the nodes. The Kutta condition is imposed via a Lagrange multiplier, in contrast with the least-squares formulation used in a previous work. For cascades, the problem is conformally mapped by a simple hyperbolic function (exponential or hyperbolic tangent) to a related problem with only one profile and one or two poles. Thus, the cascade problem is accurately solved with minor modifications to the single profile code and at the same cost of a single profile computation. Finally, several numerical examples are shown: single Joukowski and NACA profiles, interference coefficients for the flat plate cascade and a plane cascade at the external cylindrical section of an industrial fan. [Submitted to Journal of Applied Mechanics ASME ]
Centro Internacional de Métodos Numéricos en Ingeniería
2008-04-09 10:25:49
Articles
application/pdf
https://cimec.org.ar/ojs/index.php/cmm/article/view/1319
Cuadernos de Matemática y Mecánica; 2008
eng
Copyright (c)
oai:ojs.www.cimec.org.ar:article/2385
2009-09-07T19:30:44Z
cmm:ART
nmb a2200000Iu 4500
dc
Parabolic mean values and maximal estimates for gradients of temperatures
Aimar, Hugo A.
Gómez, Ivana
Iaffei, Bibiana
The main result of this paper is a pointwise estimate for
the space time gradients of a temperature on a cylindrical
domain in terms of an iteration of two maximal operators. The
result is an extension to the parabolic setting of the elliptic
inequalities proved by D.Jerison and C.Kenig and also by S.Dahlke and
R.DeVore. After an improvement of the parabolic mean value formula and the
analysis of the kernel and the operator that provides the space
derivatives
of temperatures, we obtain a pointwise estimate for space
gradients weighted by powers of the distance to the parabolic
boundary in terms of an iteration of two maximal operators which
are well known in harmonic analysis: the one-sided maximal
Hardy-Littlewood operator in the time variable and the
Calderón maximal operator in the space variable.
Published: Journal of Functional Analysis. Volume 255, Issue 8, 2008, Pages 1939-1956
Link: doi:10.1016/j.jfa.2008.06.006
Centro Internacional de Métodos Numéricos en Ingeniería
2007-03-06 10:29:57
Articles
application/pdf
https://cimec.org.ar/ojs/index.php/cmm/article/view/2385
Cuadernos de Matemática y Mecánica; 2007
eng
Copyright (c)
oai:ojs.www.cimec.org.ar:article/2701
2009-10-02T12:33:31Z
cmm:ART
nmb a2200000Iu 4500
dc
Spectral duality for finitely generated nilpotent minimum algebras, with applications
Aguzzoli, S.
Busaniche, Manuela
Marra, V.
We establish a categorical duality for the finitely generated Lindenbaum-Tarski algebras of propositional nilpotent minimum logic. The latter's conjunction is semantically interpreted by a left-continuous (but not continuous) triangular norm; implication is obtained through residuation. Our duality allows one to transfer to nilpotent minimum logic several known results about inutitionistic logic with the prelineari\-ty axiom (also called G\"{o}del-Dummett logic), mutatis mutandis. We give several such applications.
Published: Journal of Logic and Computation 17 (2007), 749-765.
Centro Internacional de Métodos Numéricos en Ingeniería
2007-03-06 10:29:57
Articles
text/plain
https://cimec.org.ar/ojs/index.php/cmm/article/view/2701
Cuadernos de Matemática y Mecánica; 2007
eng
Copyright (c)
oai:ojs.www.cimec.org.ar:article/3001
2012-02-21T00:04:57Z
cmm:ART
nmb a2200000Iu 4500
"110215 2011 eng "
dc
Dyadic Fefferman-Stein Inequalities and the Equivalence of Haar Bases on Weighted Lebesgue Spaces
Aimar, Hugo
Bernardis, Ana
Nowak, Luis
In this note we give conditions on two dyadic systems to obtain the equivalence of corresponding Haar systems on dyadic weighted Lebesgue spaces on spaces of homogeneous type. In order to obtain these result we prove a Fefferman-Stein weighted inequality for vector valued dyadic Hardy-
Littlewood maximal operators with dyadic weights in this general setting.
Centro Internacional de Métodos Numéricos en Ingeniería
2012-02-20 21:06:46
Articles
application/pdf
https://cimec.org.ar/ojs/index.php/cmm/article/view/3001
Cuadernos de Matemática y Mecánica; 2011
eng
Copyright (c)
oai:ojs.www.cimec.org.ar:article/913
2009-09-07T19:30:44Z
cmm:ART
nmb a2200000Iu 4500
dc
A Preconditioning Mass Matrix to Avoid the Ill-Posed Two-Fluid Model
Zanotti, Angel L.
Mendez, Carlos G.
Nigro, Norberto M.
Storti, Mario Alberto
Two-fluid models are central to the simulation of transport
processes in two-phase homogenized systems. Even though this physical
model has been widely accepted, an inherently non-hyperbolic and
non-conservative ill-posed problem arises from the mathematical point
of view. It has been demonstrated that this drawback occurs even for a
very simplified model, i.e., an inviscid model with no interfacial
terms. Lots of efforts have been made to remedy this anomaly and in
the literature two different types of approaches can be found. On one
hand, extra terms with physical origin are added to model the
interphase interaction, but even though this methodology seems to be
realistic, several extra parameters arise from each added term with
the associated difficulty in their estimation. On the other hand,
mathematical based-work has been done to find the way to remove the
complex eigenvalues obtained with two-fluid model
equations. Preconditioned systems, characterized as a projection of
the complex eigenvalues over the real axis, may be one of the choices.
The aim of this paper is to introduce a simple and novel mathematical
strategy based on the application of a preconditioning mass matrix
that circumvents the drawback caused by the non-hyperbolic behavior of
the original model. Although the mass and momentum conservation
equations are modified, the target of this methodology is to present
another way to reach a steady state solution (using a time marching
scheme), greatly valued by researchers in industrial process
design. Attaining this goal is possible because only the temporal term
is affected by the preconditioner. The obtained matrix has two
parameters that correct the non-hyperbolic behavior of the model: the
first one modifies the eigenvalues removing their imaginary part and
the second one recovers the real part of the original
eigenvalues. Besides the theoretical development of the
preconditioning matrix, several numerical results are presented to
show the validity of the method. [To appear in Journal of Applied Mechanics]
Centro Internacional de Métodos Numéricos en Ingeniería
2007-03-06 10:29:57
Articles
application/pdf
https://cimec.org.ar/ojs/index.php/cmm/article/view/913
Cuadernos de Matemática y Mecánica; 2007
eng
Copyright (c)
oai:ojs.www.cimec.org.ar:article/2375
2010-06-08T22:34:30Z
cmm:ART
nmb a2200000Iu 4500
dc
Convergence and quasi-optimality of adaptive FEM for Steklov eigenvalue problems
Garau, Eduardo M.
Morin, Pedro
In this article we prove convergence of adaptive finite element methods for Steklov eigenvalue problems under very general assumptions for simple as well as multiple eigenvalues starting from any initial triangulation. We also prove the optimality of the approximations assuming Dörfler's Strategy for marking, when we consider simple eigenvalues.
Centro Internacional de Métodos Numéricos en Ingeniería
2009-04-08 20:22:29
Articles
application/pdf
https://cimec.org.ar/ojs/index.php/cmm/article/view/2375
Cuadernos de Matemática y Mecánica; 2009
eng
Copyright (c)
oai:ojs.www.cimec.org.ar:article/2539
2009-08-15T13:46:05Z
cmm:ART
nmb a2200000Iu 4500
"090815 2009 eng "
dc
Data Oscillation and Convergence of Adaptive FEM
Morin, Pedro
Nochetto, Ricardo H.
Siebert, Kunibert G.
Data oscillation is intrinsic information missed by the averaging process associated with finite element methods (FEM) regardless of quadrature. Ensuring a reduction rate of data oscillation, together with an error reduction based on a posteriori error estimators, we construct a simple and efficient adaptive FEM for elliptic PDE with linear rate of convergence without any preliminary mesh adaptation nor explicit knowledge of constants. Any prescribed error tolerance is thus achieved in a finite number of steps. A number of numerical experiments in 2d and 3d yield quasi-optimal meshes along with a competitive performance.
Keywords: A posteriori error estimators, data oscillation, adaptive mesh refinement, convergence, performance, quasi-optimal meshes
AMS Subject Classifications: 65N12, 65N15, 65N30, 65N50, 65Y20
Published: SIAM Journal on Numerical Analysis, Volume 38, Number 2 (2000), 466-488.
Centro Internacional de Métodos Numéricos en Ingeniería
2006-06-24 11:50:44
Articles
application/postscript
https://cimec.org.ar/ojs/index.php/cmm/article/view/2539
Cuadernos de Matemática y Mecánica; 2001
eng
Copyright (c)
oai:ojs.www.cimec.org.ar:article/2986
2011-09-26T01:17:03Z
cmm:ART
nmb a2200000Iu 4500
dc
Generalized Qualification and Qualification Levels for Spectral Regularization Methods
Herdman, Terry
Spies, Ruben D.
Temperini, Karina G.
The concept of qualifcation for spectral regularization methods (SRM) for inverse ill-posed problems is strongly associated to the optimal order of convergence of the regularization error. In this article, the definition of
qualification is extended and three different levels are introduced: weak, strong and optimal. It is shown that the weak qualifcation extends the definition introduced
by Mathé and Pereverzev in 2003 (Mathé, P. and Pereverzev, S. V.: "Geometry of linear ill-posed problems in variable Hilbert scales", Inverse Problems, 19(3), 789-803 (2003)), mainly in the sense that the functions associated to orders of convergence and source sets need not be the same. It is shown that certain methods possessing infinite classical qualification, e.g. truncated singular value decomposition (TSVD), Landweber's method and Showalter's method, also
have generalized qualification leading to an optimal order of convergence of the regularization error. Sufficient conditions for a SRM to have weak qualification are provided and necessary and sufficient conditions for a given order of convergence to be strong or optimal qualification are found. Examples of all three qualification levels are provided and the relationships between them as well as with the classical concept of qualification and the qualification introduced by Mathé and Pereverzev are shown. In particular, SRMs
having extended qualification in each one of the three levels and having zero or infinite classical qualification are presented. Finally several implications of this theory in the
context of orders of convergence, converse results and maximal source sets for inverse ill-posed problems, are shown.
Centro Internacional de Métodos Numéricos en Ingeniería
2010-06-07 08:59:51
Articles
application/pdf
https://cimec.org.ar/ojs/index.php/cmm/article/view/2986
Cuadernos de Matemática y Mecánica; 2010
eng
Copyright (c)
oai:ojs.www.cimec.org.ar:article/458
2009-09-07T19:32:50Z
cmm:ART
nmb a2200000Iu 4500
dc
Dimension Functions of Cantor Sets
García, Ignacio
Molter, Ursula
Scotto, Roberto
We estimate the packing measure of Cantor sets associated to nonincreasing sequences through their decay. This result, dual to one obtained by Besicovitch and Taylor, allows us to characterize the dimension functions recently found by Cabrelli et al for these sets.
Published: Proceedings of the American Mathematical Society, Vol 135 (2007), 3151-3161.
Centro Internacional de Métodos Numéricos en Ingeniería
2006-06-24 12:11:10
Articles
application/pdf
https://cimec.org.ar/ojs/index.php/cmm/article/view/458
Cuadernos de Matemática y Mecánica; 2005
eng
Copyright (c)
oai:ojs.www.cimec.org.ar:article/2400
2009-12-04T16:20:34Z
cmm:ART
nmb a2200000Iu 4500
dc
Polynomials with all their roots on $S^1$
Toledano, Ricardo D.
We give a characterization of monic polynomials with coefficients in the ring of integers of a Galois number field having all of their roots on the unit circle. Such a characterization is given in terms of finitely many sum of powers of the roots of the considered polynomials.
Published: Lithuanian Mathematical Journal. Vol. 49, Nro 3, (2009), 331-336.
Link: http://www.springer.com/math/journal/10986
Centro Internacional de Métodos Numéricos en Ingeniería
2008-04-09 10:25:49
Articles
application/pdf
https://cimec.org.ar/ojs/index.php/cmm/article/view/2400
Cuadernos de Matemática y Mecánica; 2008
eng
Copyright (c)
oai:ojs.www.cimec.org.ar:article/2972
2010-06-08T22:34:30Z
cmm:ART
nmb a2200000Iu 4500
dc
Measuring the level sets of generalized homogeneous functions
Aimar, Hugo
Gómez, Ivana
In this note we prove a formula for the volume of level sets of generalized homogeneous functions in terms of measures supported on the level surfaces. We relate the results to some well known mean value formulas for solutions of PDE.
Centro Internacional de Métodos Numéricos en Ingeniería
2009-04-08 20:22:29
Articles
application/pdf
https://cimec.org.ar/ojs/index.php/cmm/article/view/2972
Cuadernos de Matemática y Mecánica; 2009
eng
Copyright (c)
oai:ojs.www.cimec.org.ar:article/9
2009-07-30T17:24:12Z
cmm:ART
nmb a2200000Iu 4500
dc
An interface strip preconditioner for domain decomposition methods: Application to hydrology
Paz, Rodrigo Rafael
Storti, Mario Alberto
In this paper, the efficiency of a parallelizable preconditioner for Domain Decomposition Methods in the context of the solution of non-symmetric linear equations arising from discretization of the Saint-Venant equations, is investigated. The proposed Interface Strip Preconditioner (IS) is based on solving a problem in a narrow strip around the interface. It requires much less memory and computing time than classical Neumann-Neumann preconditioner, and handles correctly the flux splitting among sub-domains that share the interface. The performance of this preconditioner is assessed with an analytical study of Schur complement matrix eigenvalues and numerical experiments conducted in a 2 parallel computational environment (consisting of a Beowulf cluster of twenty-nodes). [International Journal for Numerical Methods in Engineering 62(13), pp. 1873-1894 (2005)]
Centro Internacional de Métodos Numéricos en Ingeniería
2006-06-24 12:11:10
Articles
application/pdf
https://cimec.org.ar/ojs/index.php/cmm/article/view/9
Cuadernos de Matemática y Mecánica; 2005
eng
Copyright (c)
oai:ojs.www.cimec.org.ar:article/1387
2009-12-04T16:20:34Z
cmm:ART
nmb a2200000Iu 4500
dc
On the Uniform Doubling of Hutchinson Orbits of Contractive Mappings
Aimar, Hugo
Carena, Marilina
Iaffei, Bibiana
We are interested in the preservation of doubling properties along the Hutchinson orbit generated by successive applications of contraction mappings on a metric measure space. We construct some elementary examples, built on Muckenhoupt weights, showing that even when the initial and the limit points of the orbit are doubling, no iteration of the IFS remains doubling. We also obtain positive results under some quantitative assumptions on the separation of the images through the IFS. We also explore the completeness, in the Hausdorff-Kantorovich metric, of a version of the doubling property which is suitable for the application of Hutchinson type contractions.
Centro Internacional de Métodos Numéricos en Ingeniería
2008-04-09 10:25:49
Articles
application/pdf
https://cimec.org.ar/ojs/index.php/cmm/article/view/1387
Cuadernos de Matemática y Mecánica; 2008
eng
Copyright (c)
oai:ojs.www.cimec.org.ar:article/2390
2009-12-04T16:20:34Z
cmm:ART
nmb a2200000Iu 4500
dc
What is a Sobolev space for the Laguerre function systems?
Bongioanni, Bruno
Torrea, José L.
We discuss the concept of Sobolev space associated to the Laguerre operator $ L_\al = - y\,\frac{d^2}{dy^2} - \frac{d}{dy} + \frac{y}{4} + \frac{\al^2}{4y},\quad y\in (0,\infty).$ We show that the natural definition does not fit with the concept of potential space, defined via the potentials $ (L_\al)^{-s}.$ An appropriate Laguerre-Sobolev spaces are defined in order to have the mentioned coincidence. An application is given to the almost everywhere convergence of solutions of the Schr\"odinger equation. Other Laguerre operators are also considered.
Published: Studia Math. 192 (2009), no. 2, 147--172.
Centro Internacional de Métodos Numéricos en Ingeniería
2008-04-09 10:25:49
Articles
application/pdf
https://cimec.org.ar/ojs/index.php/cmm/article/view/2390
Cuadernos de Matemática y Mecánica; 2008
eng
Copyright (c)
oai:ojs.www.cimec.org.ar:article/2707
2010-06-08T22:34:30Z
cmm:ART
nmb a2200000Iu 4500
dc
Classification of finitely generated lattice-ordered groups with order-unit
Busaniche, Manuela
Cabrer, Leonardo
Mundici, D.
A unital $\ell$-group $(G,u)$ is an abelian group $G$ equipped with a translation-invariant lattice-order and a distinguished element $u$, called order-unit, whose positive integer multiples eventually dominate each element of $G$. We classify finitely generated unital $\ell$-groups by sequences $\mathcal W = (W_{0},W_{1},\ldots)$ of weighted abstract simplicial complexes, where $W_{t+1}$ is obtained from $W_{t}$ either by the classical Alexander binary stellar operation, or by deleting a maximal simplex of $W_{t}$. A simple criterion is given to recognize when two such sequences classify isomorphic unital $\ell$-groups. Many properties of the unital $\ell$-group $(G,u)$ can be directly read off from its associated sequence: for instance, the properties of being totally ordered, archimedean, finitely presented, simplicial, free.
Centro Internacional de Métodos Numéricos en Ingeniería
2009-04-08 20:22:29
Articles
text/plain
https://cimec.org.ar/ojs/index.php/cmm/article/view/2707
Cuadernos de Matemática y Mecánica; 2009
eng
Copyright (c)
oai:ojs.www.cimec.org.ar:article/4335
2013-01-18T20:35:47Z
cmm:ART
nmb a2200000Iu 4500
dc
A FFT Preconditioning Technique for the Solution of Incompressible Flow on GPUs
Storti, Mario Alberto
Paz, Rodrigo Rafael
Dalcin, Lisandro Daniel
Costarelli, Santiago Daniel
Idelsohn, Sergio Rodolfo
Graphic Processing Units have received much attention in last years. Compute-intensive algorithms operating on multidimensional arrays that have nearest neighbor dependency and/or exploit data locality can achieve massive speedups. Simulation of problems modeled by time-dependent Partial Differential Equations by using explicit time-stepping methods on structured grids is an instance of such GPU-friendly algorithms. Solvers for transient incompressible fluid flow cannot be developed in a fully explicit manner due to the incompressibility constraint. Segregated algorithms like the fractional step method require the solution of a Poisson problem for the pressure field at each time level. This stage is usually the most time-consuming one. This work discuss a solver for the pressure problem in applications using immersed boundary techniques in order to account for moving solid bodies. This solver is based on standard Conjugate Gradients iterations and depends on the availability of a fast Poisson solver on the whole domain to define a preconditioner. We provide a theoretical and numerical evidence on the advantages of our approach versus classical techniques based on fixed point iterations such as the Iterated Orthogonal Projection method.
Centro Internacional de Métodos Numéricos en Ingeniería
2012-02-28 17:07:11
Articles
application/pdf
https://cimec.org.ar/ojs/index.php/cmm/article/view/4335
Cuadernos de Matemática y Mecánica; 2012
eng
Copyright (c)
oai:ojs.www.cimec.org.ar:article/1090
2009-09-07T19:30:44Z
cmm:ART
nmb a2200000Iu 4500
dc
Numerical simulations of the flow around a spinning projectile in subsonic regime
Garibaldi, Javier
Storti, Mario Alberto
Battaglia, Laura
D'Elía, Jorge
The unsteady flow around a 155 mm projectile governed by the Navier-Stokes (NS) equations is numerically solved with a Large Eddy Simulation (LES) scheme, together with the Sub- Grid Scale (SGS) solved by a Smagorinsky model and the van Driest near-wall damping. The computed results are obtained in the subsonic flow regime for a viscous and incompressible Newtonian fluid in order to determine the axial drag coefficient, and they are validated against experimental data. The problem was solved by a monolithic finite element code for parallel computing on a Beowulf cluster.
Centro Internacional de Métodos Numéricos en Ingeniería
2007-03-06 10:29:57
Articles
application/pdf
https://cimec.org.ar/ojs/index.php/cmm/article/view/1090
Cuadernos de Matemática y Mecánica; 2007
eng
Copyright (c)
oai:ojs.www.cimec.org.ar:article/2381
2009-12-04T16:20:34Z
cmm:ART
nmb a2200000Iu 4500
dc
Convergence of adaptive finite element methods for eigenvalue problems
Garau, Eduardo M.
Morin, Pedro
Zuppa, Carlos
In this article we prove convergence of adaptive finite element methods for second order elliptic eigenvalue problems. We consider Lagrange finite elements of any degree and prove convergence for simple as well as multiple eigenvalues under a minimal refinement of marked elements, for all reasonable marking strategies, and starting from any initial triangulation.
Published: Mathematical Models and Methods in Applied Sciences (M3AS), Vol. 19, No. 5 (2009), 721--747.
Centro Internacional de Métodos Numéricos en Ingeniería
2008-04-09 10:25:49
Articles
application/pdf
https://cimec.org.ar/ojs/index.php/cmm/article/view/2381
Cuadernos de Matemática y Mecánica; 2008
eng
Copyright (c)
oai:ojs.www.cimec.org.ar:article/2697
2009-10-02T12:15:27Z
cmm:ART
nmb a2200000Iu 4500
dc
Decomposition of BL-chains
Busaniche, Manuela
A simple and self contained proof of decomposition of BL-chains into ordinal sums of Wajsberg hoops is given.
Published: Algebra Universalis, 52 (2004), 519-525.
Centro Internacional de Métodos Numéricos en Ingeniería
2006-08-12 12:07:07
Articles
text/plain
https://cimec.org.ar/ojs/index.php/cmm/article/view/2697
Cuadernos de Matemática y Mecánica; 2004
eng
Copyright (c)
oai:ojs.www.cimec.org.ar:article/2997
2011-09-26T01:17:03Z
cmm:ART
nmb a2200000Iu 4500
dc
Galerkin Boundary Elements for Exterior Stokes Flows
D'Elía, Jorge
Battaglia, Laura
Cardona, Alberto
Storti, Mario Alberto
Franck, Gerardo
An indirect boundary integral equation for steady Stokes ﬂow around a rigid body in the three-dimensional space is proposed, and is numerically solved by using collocation and Galerkin weight- ing procedures. The resulting double surface integrals of the Galerkin technique that express the pairwise interaction among all boundary elements, which are quadruple integrals, are computed on ﬂat simplex triangles using a regularized quadrature scheme. Numerical examples include the steady creeping ﬂow around the sphere of unit radius and the cube of unit edge length, covering issues on the convergence under mesh reﬁnement and stability under small mesh perturbations.
Centro Internacional de Métodos Numéricos en Ingeniería
2010-06-07 08:59:51
Articles
application/pdf
https://cimec.org.ar/ojs/index.php/cmm/article/view/2997
Cuadernos de Matemática y Mecánica; 2010
eng
Copyright (c)
oai:ojs.www.cimec.org.ar:article/874
2007-04-27T16:48:12Z
cmm:ART
nmb a2200000Iu 4500
dc
Dynamic boundary conditions in CFD
Storti, Mario Alberto
Nigro, Norberto Marcelo
Paz, Rodrigo Rafael
Dalcín, Lisandro Daniel
The number and type of boundary conditions to be used in the numerical modeling of fluid mechanics problems is normally chosen according to a simplified analysis of the character- istics, and also from the experience of the modeler. The problem is harder at input/output boundaries which are, in most cases, artificial boundaries, so that a bad decision about the boundary conditions to be imposed may affect the precision and stability of the whole computation. For inviscid flows, the analysis of the sense of propagation in the normal di- rection to the boundaries gives the number of conditions to be imposed and, in addition, the conditions that are "absorbing" for the waves impinging normal to the boundary. In practice, it amounts to counting the number of positive and negative eigenvalues of the advective flux Jacobian projected onto the normal. The problem is still harder when the number of incoming characteristics varies during the computation, and to correctly treat these cases poses both mathematical and practical problems. One example considered here is compressible flow where the flow regime at a certain part of an inlet/outlet boundary can change from subsonic to supersonic and the flow can revert. In this work the technique for dynamically imposing the correct number of boundary conditions along the computation, using Lagrange multipliers and penalization is discussed, and several numerical examples are presented. [Submitted to Computer Methods in Applied Mechanics and Engineering (2007)]
Centro Internacional de Métodos Numéricos en Ingeniería
2006-06-24 09:54:16
Articles
application/pdf
https://cimec.org.ar/ojs/index.php/cmm/article/view/874
Cuadernos de Matemática y Mecánica; 2006
eng
Copyright (c)
oai:ojs.www.cimec.org.ar:article/1659
2009-04-28T02:33:33Z
cmm:ART
nmb a2200000Iu 4500
dc
Iterative solution of panel method discretizations for 3D potential flow problems. The modal multipolar preconditioning
D'Elía, Jorge
Storti, Mario Alberto
Idelsohn, Sergio R.
The iterative solution of linear systems arising from panel method
discretization of three-dimensional (3D) exterior potential problems
coming mainly from aero-hydrodynamic engineering problems, is
discussed. We propose an original peconditioning based on an
approximate eigenspace decomposition, that corrects bad conditioning
arising from pair of surfaces that are very close from each other,
which is a very common situation in slender wings and other
aerodynamoc profiles. This preconditioning has been tested with the
standard Bi-Conjugate Gradient (BCG) and Conjugate Gradient Squared
(CGS) iterative methods.
Centro Internacional de Métodos Numéricos en Ingeniería
2006-06-24 11:50:44
Articles
application/pdf
https://cimec.org.ar/ojs/index.php/cmm/article/view/1659
Cuadernos de Matemática y Mecánica; 2001
eng
Copyright (c)
oai:ojs.www.cimec.org.ar:article/2535
2009-08-14T21:45:16Z
cmm:ART
nmb a2200000Iu 4500
dc
Surface Diffusion of Graphs: Variational Formulation, Error Analysis and Simulation
Bänsch, Eberhard
Morin, Pedro
Nochetto, Ricardo H.
Abstract: Surface diffusion is a (4th order highly nonlinear) geometric driven motion of a surface with normal velocity proportional to the surface Laplacian of mean curvature. We present a novel variational formulation for graphs and derive a priori error estimates for a time-continuous finite element discretization. We also introduce a semi-implicit time discretization and a Schur complement approach to solve the resulting fully discrete, linear systems. After computational verification of the orders of convergence for polynomial degrees 1 and 2, we show several simulations in 1d and 2d with and without forcing which explore the smoothing effect of surface diffusion as well as the onset of singularities in finite time, such as infinite slopes and cracks.
Keywords: Surface diffusion, fourth-order parabolic problem, finite elements, a priori error estimates, Schur complement, smoothing effect.
AMS Subject Classifications: 35K55, 65M12, 65M15, 65M60, 65Z05.
Published: SIAM Journal on Numerical Analysis, Volume 42, Number 2 (2004), 773--799.
Centro Internacional de Métodos Numéricos en Ingeniería
2006-08-12 12:06:41
Articles
application/pdf
https://cimec.org.ar/ojs/index.php/cmm/article/view/2535
Cuadernos de Matemática y Mecánica; 2003
eng
Copyright (c)
oai:ojs.www.cimec.org.ar:article/2980
2011-09-26T01:17:03Z
cmm:ART
nmb a2200000Iu 4500
dc
Flow Study and Wetting Efficiency of a Perforated-Plate Tray Distributor in a Trickle Bed Reactor
Ramajo, Damián
Márquez-Damián, Santiago
Raviculé, Marcela
Monsalvo, María
Storti, Mario Alberto
Nigro, Norberto Marcelo
Trickle bed reactors (TBRs) are massively employed in petrochemical and chemical plants. In general these reactors have one or more beds filled up with catalyst particles. Efficient catalyst utilization relies on a good liquid charge distribution across the catalyst beds. However, normally the distribution is not perfect and some parts of the beds will get less liquid reactants while others will get more than the average. In zones where there is a maldistribution of reactants the reaction will progress in undesirable way leading to deactivation of the catalyst and towards low conversions. Bad tray efficiency due to non-uniform liquid distribution will result in low reactor efficiency and shorten the catalyst's cycle time. The TBR analyzed here is a hydrogenation one that processes C4 (liquid) and hydrogen (gas) to produce butene 1 (also named α butylene). The two-phase charge is introduced through the upper side of the TBR and the liquid phase accumulates on the tray to a certain level swamping the perforated-plate tray. The liquid phase flows down through 68 small holes while the gas phase circulates through 7 gas chimneys. There is another ceramic-ball bed above the catalyst bed with the aim to obtain a better distribution of the charge. In this work a computational fluid dynamics analysis (CFD) employing the Eulerian two-fluid model was carried out with the aim to understand the fluid dynamics of the distribution process and to determine the wetting efficiency of the tray distributor under different operating conditions. The small tray holes were modeled by sinks (drains) and sources, firstly employing numerical and experimental models to obtain the flow rate versus liquid height response. Because of the scarce liquid sloshing above the tray, little differences on the liquid discharge through the holes were found. Due to the low gas fraction of the charge the liquid flows only by gravity following an almost vertical trajectory from the holes to the ceramic ball bed. So, the extension of the wetted zone at the top of the ceramic ball bed is small. A suitable correlation to estimate liquid diffusion inside the ceramic-ball bed was employed, showing that the overall catalyst bed surface is wetted but significant differences on liquid fraction are found. Moreover, a possible additional cause of the low TBR efficiency could be the well known fouling vulnerability of this kind of tray distributors. In this sense, two simple geometric modifications were proposed to enhance tray performance; firstly reducing the amount of gas chimneys to only one, thus adding additional drip points, secondly replacing the holes by short risers in order to reduce the vulnerability to plugging. [Submitted to Int J Chem React Engng ISSN: 1542-6580]
Centro Internacional de Métodos Numéricos en Ingeniería
2010-06-07 08:59:51
Articles
application/pdf
https://cimec.org.ar/ojs/index.php/cmm/article/view/2980
Cuadernos de Matemática y Mecánica; 2010
eng
Copyright (c)
oai:ojs.www.cimec.org.ar:article/266
2006-09-13T12:22:48Z
cmm:ART
nmb a2200000Iu 4500
dc
Comparison of Hardy-Littlewood and Dyadic Maximal Functions on Spaces of Homogeneous Type
Aimar, Hugo A.
Bernardis, Ana L.
Iaffei, Bibiana
We obtain a comparison of the level sets for two maximal functions on a space of homogeneous type: the Hardy-Littlewood maximal function of mean values over balls and the dyadic maximal function of mean values over the
dyadic sets introduced by M. Christ in [4]. As applications to the theory of Ap weights on this setting, we compare the standard and the dyadic Muckenhoupt
classes and we give an alternative proof of reverse Hölder type inequalities.
Centro Internacional de Métodos Numéricos en Ingeniería
2006-06-24 09:54:16
Articles
application/pdf
https://cimec.org.ar/ojs/index.php/cmm/article/view/266
Cuadernos de Matemática y Mecánica; 2006
eng
Copyright (c)
oai:ojs.www.cimec.org.ar:article/2395
2009-09-07T19:30:44Z
cmm:ART
nmb a2200000Iu 4500
dc
Principal fitted components for dimension reduction in regression
Cook, Dennis
Forzani, Liliana
We provide a remedy for two concerns that have dogged the use of principal components in regression: (i) principal components are computed from the predictors alone and do not make apparent use of the response, and (ii) principal components are not invariant or equivariant under full rank linear transformation of the predictors. The development begins with principal fitted components [Cook, R. D. (2007). Fisher lecture: Dimension reduction in regression (with discussion). Statist. Sci. 22 126] and uses normal models for the inverse regression of the predictors on the response to gain reductive information for the forward regression of interest. This approach includes methodology for testing hypotheses about the number of components and about conditional independencies among the predictors.
Published: Statistical Science 2008, Vol. 23, No. 4, 485–501.
DOI: 10.1214/08-STS275
Centro Internacional de Métodos Numéricos en Ingeniería
2007-03-06 10:29:57
Articles
application/pdf
https://cimec.org.ar/ojs/index.php/cmm/article/view/2395
Cuadernos de Matemática y Mecánica; 2007
eng
Copyright (c)
oai:ojs.www.cimec.org.ar:article/2969
2010-06-08T22:34:30Z
cmm:ART
nmb a2200000Iu 4500
dc
Odd BMO(R) functions and Carleson measures in the Bessel setting
Chicco Ruiz, Aníbal
Fariña, J.C.
Rodríguez-Mesa, L.
In this paper we characterize the odd functions in BMO(R) by using Carleson measures associated with Poisson and heat semigroups for Bessel operators.
Accepted: Integral Equations and Operator Theory
Centro Internacional de Métodos Numéricos en Ingeniería
2009-04-08 20:22:29
Articles
application/pdf
https://cimec.org.ar/ojs/index.php/cmm/article/view/2969
Cuadernos de Matemática y Mecánica; 2009
eng
Copyright (c)
oai:ojs.www.cimec.org.ar:article/5
2007-02-21T19:04:40Z
cmm:ART
nmb a2200000Iu 4500
dc
Strong coupling strategy for fluid structure interaction problems in supersonic regime via fixed point iteration
Storti, Mario Alberto
Nigro, Norberto Marcelo
Paz, Rodrigo Rafael
In this paper we present some results on the time integration stability when solving fluid/structure interaction problems with strong coupled staged strategies. The flutter of a flat plate faced with a fluid at Mach between 1.8 and 3 is taken as a benchmark. The precision of differ- ent predictor strategies and the influence of the partitioned strong coupling on stability is discussed. [Submitted to Journal of Sound and Vibration]
Centro Internacional de Métodos Numéricos en Ingeniería
2006-06-24 09:54:16
Articles
application/pdf
https://cimec.org.ar/ojs/index.php/cmm/article/view/5
Cuadernos de Matemática y Mecánica; 2006
eng
Copyright (c)
oai:ojs.www.cimec.org.ar:article/1380
2009-12-04T16:20:34Z
cmm:ART
nmb a2200000Iu 4500
dc
Full Numerical Quadrature in Galerkin Boundary Element Methods
D'Elía, Jorge
Storti, Mario Alberto
Battaglia, Laura
When a Galerkin discretization of a boundary integral equation with a weakly singular kernel is performed over triangles, a double surface integral must be evaluated for each pair of them. If these pairs are not contiguous nor the same, the kernel is regular and a Gauss-Legendre quadrature can be employed. But, when they have a common edge or a common vertex, then edge and vertex weak sin- gularities appear, while autointegrals have both facets coincidents and the whole integration domain is weakly singular. Taylor (IEEE Trans. on Antennas and Propagation, 51(7):1630–1637 (2003)) proposed a systematic evaluation, based on a reordering and partition of the integration domain, together a use of the Duffy transformations in order to remove the singularities, in such a way that a Gauss-Legendre quadrature was performed on three coordinates with an analytic one in the fourth one. Since this scheme is a bit restrictive because it was designed for electromagnetic wave propagation kernels, a full numer- ical quadrature on the four coordinates is proposed here in order to handle other kernels, like those of creeping flows. A numerical test is also proposed based on slight modification of the Wang-Atalla one (Comm. in Num. Meth. Eng., 13(0):1–7 (1997)). [Submitted for publicaciont to Communications in Numerical Methods in Engineering]
Centro Internacional de Métodos Numéricos en Ingeniería
2008-04-09 10:25:49
Articles
application/pdf
https://cimec.org.ar/ojs/index.php/cmm/article/view/1380
Cuadernos de Matemática y Mecánica; 2008
eng
Copyright (c)
oai:ojs.www.cimec.org.ar:article/2386
2010-06-08T22:34:30Z
cmm:ART
nmb a2200000Iu 4500
dc
Weighted inequalities and pointwise estimates for the multilinear fractional integral and maximal operators
Pradolini, Gladis G.
In this article we prove weighted norm inequalities and pointwise estimates between the multilinear fractional integral operator and the multilinear fractional maximal. As a consequence of these estimations we obtain weighted weak and strong inequalities for the multilinear fractional integral operator. In particular, we extend some results given in \cite{CPSS} to the multilinear context. On the other hand we prove weighted pointwise estimates between the multilinear fractional maximal operator ${\cal M}_{\alpha,B}$ associated to a Young function $B$ and the multilinear maximal operators ${\cal M}_{\psi}={\cal M}_{0,\psi}$, $\psi(t)=B(t^{1-\alpha/(nm)})^{{nm}/{(nm-\alpha)}}$. As an application of these estimate we obtain a direct proof of the $L^p-L^q$ boundedness results of ${\cal M}_{\alpha,B}$ for the case $B(t)=t$ and $B_k(t)=t(1+\log^+t)^k$ when $1/q=1/p-\alpha/n$. We also give sufficient conditions on the weights involved in the boundedness results of ${\cal M}_{\alpha,B}$ that generalizes those given in \cite{M} for $B(t)=t$. Finally, we prove some boundedness results in Banach function spaces for a generalized version of the multilinear fractional maximal operator.
Centro Internacional de Métodos Numéricos en Ingeniería
2009-04-08 20:22:29
Articles
application/pdf
https://cimec.org.ar/ojs/index.php/cmm/article/view/2386
Cuadernos de Matemática y Mecánica; 2009
eng
Copyright (c)
oai:ojs.www.cimec.org.ar:article/2702
2009-12-04T16:20:34Z
cmm:ART
nmb a2200000Iu 4500
dc
Free MVn-algebras
Busaniche, Manuela
Cignoli, Roberto
In this note we characterize free algebras in varieties of MV-algebras generated by a finite chain ${\bf L}_n$ as algebras of continuous functions from the spectrum of the Boolean skeleton of the free algebra into ${\bf L}_n$.
Published: Algebra Universalis 58 (2008), 335-340.
Centro Internacional de Métodos Numéricos en Ingeniería
2008-04-09 10:25:49
Articles
text/plain
https://cimec.org.ar/ojs/index.php/cmm/article/view/2702
Cuadernos de Matemática y Mecánica; 2008
eng
Copyright (c)
oai:ojs.www.cimec.org.ar:article/3002
2012-02-21T00:05:24Z
cmm:ART
nmb a2200000Iu 4500
"110215 2011 eng "
dc
Haar Type Bases in Lorentz Spaces Via Extrapolation
Aimar, Hugo
Bernardis, Ana
Nowak, Luis
In this note we consider Haar systems as unconditional bases for Lorentz spaces defined on spaces of homogeneous type. We also give characterizations of these spaces in terms of the Haar coefficients. The basic tools are the Rubio de Francia extrapolation technique and the characterization of weighted Lebesgue spaces with Haar bases.
Centro Internacional de Métodos Numéricos en Ingeniería
2012-02-20 21:06:46
Articles
application/pdf
https://cimec.org.ar/ojs/index.php/cmm/article/view/3002
Cuadernos de Matemática y Mecánica; 2011
eng
Copyright (c)
oai:ojs.www.cimec.org.ar:article/1017
2009-07-30T17:24:12Z
cmm:ART
nmb a2200000Iu 4500
dc
An Interface Strip Preconditioner for Domain Decomposition Methods Application to Hydrology
Paz, Rodrigo R.
Storti, Mario Alberto
In this paper, the efficiency of a parallelizable preconditioner for
domain decomposition methods in the context of the solution of
non-symmetric linear equations arising from discretization of the
Saint-Venant equations, is investigated. The proposed interface strip
preconditioner (IS) is based on solving a problem in a narrow strip
around the interface. It requires much less memory and computing time
than classical NeumannNeumann preconditioner, and handles
correctly the flux splitting among sub-domains that share the
interface. The performance of this preconditioner is assessed with an
analytical study of Schur complement matrix eigenvalues and numerical
experiments conducted in a parallel computational environment
(consisting of a Beowulf cluster of 20 nodes). [To appear in Int. J Num. Meth. Engng. ]
Centro Internacional de Métodos Numéricos en Ingeniería
2006-06-24 12:11:10
Articles
application/pdf
https://cimec.org.ar/ojs/index.php/cmm/article/view/1017
Cuadernos de Matemática y Mecánica; 2005
eng
Copyright (c)
oai:ojs.www.cimec.org.ar:article/2376
2010-06-08T22:34:30Z
cmm:ART
nmb a2200000Iu 4500
dc
On Besov regularity of temperatures
Aimar, Hugo A.
Gómez, Ivana
Iaffei, Bibiana
We prove space-time parabolic Besov regularity in terms
of integrability of Besov norms in the space variable for solutions of the
heat equation on cylindrical regions based on Lipschitz domains.
Accepted: The Journal of Fourier Analysis and Applications. January 2010.
Centro Internacional de Métodos Numéricos en Ingeniería
2009-04-08 20:22:29
Articles
application/pdf
https://cimec.org.ar/ojs/index.php/cmm/article/view/2376
Cuadernos de Matemática y Mecánica; 2009
eng
Copyright (c)
oai:ojs.www.cimec.org.ar:article/2540
2010-06-08T22:34:30Z
cmm:ART
nmb a2200000Iu 4500
dc
Validation of A 0D/1D Computational Code for the Design of Several Kind of Internal Combustion Engines
López, Ezequiel José
Nigro, Norberto Marcelo
A code for computational simulation of internal combustion engines is presented. One-dimensional gas dynamics equations are used for model the ﬂow through pipes and manifolds, and the remaining components in the engine (cylinders, valves, etc.) are modeled by using thermodynamic or 0D models. The numerical code developed is able to simulate spark-ignition and compressionignition, two-stroke and four-stroke, multicylinder and multi-valve engines, naturally aspirated or turbo-charged, and diﬀerent geometries of the combustion chamber. The code was implemented in the scripting language Python, which is a dynamic object-oriented programming language that oﬀers strong support for integration with other languages and tools. The numerical methods used in the discretization of the equations and implementation details are presented. Several test cases are included in order to show the performance of the code.
Centro Internacional de Métodos Numéricos en Ingeniería
2009-04-08 20:22:29
Articles
application/pdf
https://cimec.org.ar/ojs/index.php/cmm/article/view/2540
Cuadernos de Matemática y Mecánica; 2009
eng
Copyright (c)
oai:ojs.www.cimec.org.ar:article/2987
2011-09-26T01:17:03Z
cmm:ART
nmb a2200000Iu 4500
dc
Arbitrary Divergence Speed of the Least-Squares Method in Infinite-Dimensional Inverse Ill-Posed Problems
Spies, Ruben D.
Temperini, Karina G.
A standard engineering procedure for approximating the solutions of an infinite-dimensional inverse problem of the form Ax = y, where A is a given compact linear operator on a Hilbert space X and y is the given data, is to find a sequence
{XN} of finite-dimensional approximating subspaces of X whose union is dense in X and to construct the sequence {xN} of least squares solutions of the problem in XN.
In 1980, Seidman (Nonconvergence Results for the Application of Least-Squares Estimation to Ill-Posed Problems, Journal of Optimization Theory and Applications, 30, 4, 535-547, 1980) showed that if the problem is ill-posed, then, without any additional assumptions on the exact solution or on the sequence of approximating subspaces {XN}, it cannot be guaranteed that the sequence fxNg will converge to the exact solution. In this article this result is extended in the following sense: it is shown that if X is separable, then for any y in X, y not equal 0, and for any arbitrarily given function s defined over the natural numbers with values in R+, there exists an injective, compact linear operator A and an increasing sequence of finite-dimensional subspaces XN contained in X such that ||xN - A^{-1}y||>= s(N) for all N, where xN is the least squares solution of Ax = y in XN.
Centro Internacional de Métodos Numéricos en Ingeniería
2010-06-07 08:59:51
Articles
application/pdf
https://cimec.org.ar/ojs/index.php/cmm/article/view/2987
Cuadernos de Matemática y Mecánica; 2010
eng
Copyright (c)
oai:ojs.www.cimec.org.ar:article/459
2007-04-27T12:12:42Z
cmm:ART
nmb a2200000Iu 4500
dc
On the Composition of the Integral and Derivative Operators of Functional Order
Hartzstein, Silvia I.
Viviani, Beatriz E.
The Integral, I_phi, and Derivative, D_phi, operators of order phi, with phi a function of positive lower type and upper type lower than 1, were defined in [HV2] in the setting of spaces of homogeneous-type. These definitions generalize those of the fractional integral and derivative operators of order alpha, where phi(t) = t^alpha, given in [GSV].
In this work we show that the composition T_phi = D_phi I_phi is a singular integral operator. This result in addition with the results obtained in [HV2] of boundedness of I_phi and D_phi or the T1-theorems proved in [HV1] yield the fact that T_phi is a Calderón-Zygmund operator bounded on the generalized Besov, ˙B^{psi,q}_p 1 leq p, q < infty, and Triebel-Lizorkin spaces, ˙F^{psi,q}_p , 1 < p, q < infty, of order psi = psi_ 1/psi_ 2, where psi_1 and psi_2 are two quasi-increasing functions of adequate upper types s_1 and s_2, respectively.
Centro Internacional de Métodos Numéricos en Ingeniería
2006-06-24 11:50:44
Articles
application/pdf
https://cimec.org.ar/ojs/index.php/cmm/article/view/459
Cuadernos de Matemática y Mecánica; 2001
eng
Copyright (c)
oai:ojs.www.cimec.org.ar:article/2401
2010-06-08T22:34:30Z
cmm:ART
nmb a2200000Iu 4500
dc
Boundedness of the fractional maximal operator on variable exponent Lebesgue spaces: a short proof
Gorosito, Osvaldo P.
Pradolini, Gladis G.
Salinas, Oscar
We give a simple proof of the boundedness of the fractional maximal operator providing in this way an alternative approach to the one given by C. Capone, D. Cruz Uribe and A. Fiorenza in The fractional maximal operator and fractional integrals on variable L^p spaces, 743--770.
Centro Internacional de Métodos Numéricos en Ingeniería
2009-04-08 20:22:29
Articles
application/pdf
https://cimec.org.ar/ojs/index.php/cmm/article/view/2401
Cuadernos de Matemática y Mecánica; 2009
eng
Copyright (c)
oai:ojs.www.cimec.org.ar:article/2981
2010-06-08T22:40:05Z
cmm:ART
nmb a2200000Iu 4500
dc
On local times, density estimation and classification rules for functional data
Llop, Pamela Nerina
Forzani, Liliana
Fraiman, Ricardo
In this paper we define a $\sqrt{n}$-consistent nonparametric density estimator for functional data. Under mild conditions we obtain strong consistency, strong orders of convergence and derive the asymptotic distribution of the estimator. We propose a nonparametric classification rule based on local times (occupation measure) and include some simulations studies
Centro Internacional de Métodos Numéricos en Ingeniería
2009-04-08 20:22:29
Articles
application/pdf
https://cimec.org.ar/ojs/index.php/cmm/article/view/2981
Cuadernos de Matemática y Mecánica; 2009
eng
Copyright (c)
oai:ojs.www.cimec.org.ar:article/26
2009-07-30T17:24:12Z
cmm:ART
nmb a2200000Iu 4500
"060704 2006 eng "
dc
A Preconditioner for the Schur Complement Matrix
Storti, Mario Alberto
Dalcín, Lisandro Daniel
Paz, Rodrigo R.
Yommi, Andrea
Sonzogni, Victorio E.
Nigro, Norberto M.
A preconditioner for iterative solution of the interface problem in Schur Complement Domain Decomposition Methods is presented. This preconditioner is based on solving a problem in a narrow strip around the interface. It requires much less memory and computing time than classical Neumann-Neumann preconditioner and its variants, and handles correctly the flux splitting among subdomains that share the interface. The performance of this preconditioner is assessed with an analytical study of Schur complement matrix eigenvalues and several numerical experiments conducted in a sequential computational environment. Results in a production parallel finite element code are given in a companion paper. [Advances in Engineering Software 37, pp. 754-762 (2006)]
Centro Internacional de Métodos Numéricos en Ingeniería
2006-06-24 12:11:10
Articles
application/pdf
https://cimec.org.ar/ojs/index.php/cmm/article/view/26
Cuadernos de Matemática y Mecánica; 2005
eng
Copyright (c)
oai:ojs.www.cimec.org.ar:article/1390
2010-11-18T22:44:30Z
cmm:ART
nmb a2200000Iu 4500
"090408 2009 eng "
dc
Local Absorbent Boundary Condition for Nonlinear Hyperbolic Problems with Unknown Riemann Invariants
Paz, Rodrigo Rafael
Garelli, Luciano
Storti, Mario Alberto
A general methodology for developing absorbing boundary conditions for general nonlinear hyperbolic advective-diffusive equations with unknown Riemann invariants is presented. In problems where the Riemann invariants (RI) are known (e.g. the flow in a shallow rectangular channel, the gas flow equations), the imposition of non-reflective boundary conditions is straightforward. In problems where Riemann invariants are unknown (e.g. the flow in a non-rectangular channels, the stratified 2D shallow water flows) it is possible to impose that kind of conditions analyzing the projection of the Jacobians of advective flux functions onto normal directions to fictitious surfaces or boundaries. The advantage of the method is that it is very easy to implement in a finite element code and is only based on computing the advective flux functions (and the their Jacobian projections), then, imposing non-linear constraints via Lagrange Multipliers or Penalty Methods. The application of the dynamic absorbing boundary conditions to typical wave propagation problems with unknown Riemann invariants, like non-linear Saint-Venant system of conservation laws for non-rectangular and non-prismatic 1D channels and stratified 1D/2D shallow water equations, is presented. Also, the new absorbent/dynamic condition can handle automatically the change of Jacobians structure when the flow regime changes from subcritical to supercritical and viceversa, or when recirculating zones are present in regions near fictitious walls.
Centro Internacional de Métodos Numéricos en Ingeniería
2009-04-08 20:22:29
Articles
application/pdf
https://cimec.org.ar/ojs/index.php/cmm/article/view/1390
Cuadernos de Matemática y Mecánica; 2009
eng
Copyright (c)
oai:ojs.www.cimec.org.ar:article/2391
2009-09-07T19:30:44Z
cmm:ART
nmb a2200000Iu 4500
dc
Composition of operators in Orlicz spaces
Bongioanni, Bruno
Harboure, Eleonor
In this work we find sharp conditions for boundedness on Orlicz spaces of the composition of $j$ operators, each one being of restricted weak type $(p,p)$ for some $p>1$, and of strong type $(\infty,\infty)$. Particularly, we find necessary and sufficient conditions to obtain modular inequalities for the $j$-times composition of the Cesàro Maximal function of order $\alpha$. With this approach we treat a kind of strong maximal function related to Cesàro averages over $n$-dimensional rectangles.
Published: Rocky Mountain J. Math. 38 (2008), no. 1, 41--59.
Centro Internacional de Métodos Numéricos en Ingeniería
2007-03-06 10:29:57
Articles
application/pdf
https://cimec.org.ar/ojs/index.php/cmm/article/view/2391
Cuadernos de Matemática y Mecánica; 2007
eng
Copyright (c)
oai:ojs.www.cimec.org.ar:article/2961
2010-06-08T22:34:30Z
cmm:ART
nmb a2200000Iu 4500
dc
An h-Adaptive Solution of the Spherical Blast Wave Problem
Ríos Rodríguez, Gustavo
Storti, Mario Alberto
López, Ezequiel Javier
Sarraf, Sofía
Compressible ﬂow problems usually exhibit features in the ﬂow ﬁeld that need a very reﬁned mesh in order to achieve an acceptable accuracy in the numerical solution. Uniform reﬁnement of the whole mesh quickly becomes prohibitive in three dimensions. In these situations, the use of an adaptive technique shows advantageous. Here, a h-reﬁnement/coarsening strategy for unstructured ﬁnite element meshes is described and applied to solve the spherical blast wave problem driven by a point-like intense explosion. The resulting meshes are non-conformal. A 3-D reﬁnement constraint is introduced to enforce a smooth size distribution among neighbor elements in the mesh and a gradient error indicator based on the ﬂow variables is used to track discontinuities. It is shown that the adaption strategy allows to improve the accuracy and to reduce the spurious oscillations in the numerical solution, requiring just a small fraction of the overall computing time.
Centro Internacional de Métodos Numéricos en Ingeniería
2009-04-08 20:22:29
Articles
application/pdf
https://cimec.org.ar/ojs/index.php/cmm/article/view/2961
Cuadernos de Matemática y Mecánica; 2009
eng
Copyright (c)
oai:ojs.www.cimec.org.ar:article/4337
2013-01-20T13:27:33Z
cmm:ART
nmb a2200000Iu 4500
dc
Fastmat: A C++ Library for Multi-Index Array Computations
Paz, Rodrigo Rafael
Storti, Mario Alberto
Dalcin, Lisandro Daniel
Castro, Hugo Guillermo
Kler, Pablo
In this paper we introduce and describe an efficient thread-safe
matrix library for computing element/cell residuals and Jacobians in
Finite Elements and Finite Volumes-like codes. The library provides
a wide range of multi-index tensor operations that are normally used
in scientific numerical computations. The library implements an
algorithm for choosing the optimal computation order when a product
of several tensors is performed (i.e., the so-called `multi-product'
operation). Another key-point of the FastMat approach is that
some computations (for instance the optimal order in the
multi-product operation mentioned before) are computed in the first
iteration of the loop body and stored in a cache object, so
that in the second and subsequent executions these computations are
retrieved from the cache, and then not recomputed. The library is
open source and freely available within the multi-physics parallel
FEM code PETSc-FEM http://www.cimec.org.ar/petscfem and it can
be exploited on distributed and shared memory architectures as well
as in hybrid approaches. We studied the performance of the library
in the context of typical FEM tensor contractions.
Centro Internacional de Métodos Numéricos en Ingeniería
2012-02-28 17:07:11
Articles
application/pdf
https://cimec.org.ar/ojs/index.php/cmm/article/view/4337
Cuadernos de Matemática y Mecánica; 2012
eng
Copyright (c)
oai:ojs.www.cimec.org.ar:article/1091
2009-09-07T19:30:44Z
cmm:ART
nmb a2200000Iu 4500
dc
CFD analysis of a polybutene reactor to diagnose causes of polymer adhesion at inner walls
Ramajo, Damián
Raviculé, Marcela
Horowitz, Gabriel
Storti, Mario Alberto
Nigro, Norberto Marcelo
In this work the flow patterns inside a polybutene reactor were studied by CFD (Computational Fluid Dynamics) using finite elements. Research was carried out with the aim to find possible causes of excessive adhesion of polymer and catalyst particles at the reactor inner walls. The multi-fluid formulation for a three-phase system formed by liquid reactor mixture, solid catalyst particles and small gas bubbles generated by the reaction was applied. Deposition of solid particles and a non-homogeneous flow distribution over the lower reactor walls were founded. Based on the hypothesis that adhesion phenomena is related to a combination of catalyst-particle precipitation at walls and locally low shear stresses, several operative and constructive modifications were proposed in order to reduce this phenomenon. [Submitted to LAAR]
Centro Internacional de Métodos Numéricos en Ingeniería
2007-03-06 10:29:57
Articles
application/pdf
https://cimec.org.ar/ojs/index.php/cmm/article/view/1091
Cuadernos de Matemática y Mecánica; 2007
eng
Copyright (c)
oai:ojs.www.cimec.org.ar:article/2382
2010-06-08T22:34:30Z
cmm:ART
nmb a2200000Iu 4500
dc
Pointwise estimates for gradients of temperatures in terms of maximal functions
Aimar, Hugo A.
Gómez, Ivana
Iaffei, Bibiana
We give a detailed proof, in the case of one space dimension, of a pointwise upper estimate for the space gradient of a temperature. The operators involved are a one-sided Hardy-Littlewood maximal in time and the Calderón sharp maximal operator in space.
Published: Revista de la Unión Matemática Argentina, vol. 50, no. 2, 105-114, (2009).
Centro Internacional de Métodos Numéricos en Ingeniería
2009-04-08 20:22:29
Articles
application/pdf
https://cimec.org.ar/ojs/index.php/cmm/article/view/2382
Cuadernos de Matemática y Mecánica; 2009
eng
Copyright (c)
oai:ojs.www.cimec.org.ar:article/2698
2009-10-02T12:17:58Z
cmm:ART
nmb a2200000Iu 4500
dc
Free algebras in varieties of BL-algebras generated by BLn-chains
Busaniche, Manuela
Cignoli, Roberto
Free algebras with an arbitrary number of free generators in varieties of BL-algebras generated by one BL-chain which is an ordinal sum of a finite MV-chain ${\bf L}_n$ and a generalized BL-chain ${\bf B}$ are described in terms of weak Boolean products of BL-algebras that are ordinal sums of subalgebras of ${\bf L}_n$ and free algebras in the variety of basic hoops generated by ${\bf B}.$ The Boolean products are taken over the Stone spaces of the Boolean subalgebras of idempotents of free algebras in the variety of MV-algebras generated by ${\bf L}_n.$
Published: Journal of the Australian Mathematical Society, 80 (2006), 419-439.
Centro Internacional de Métodos Numéricos en Ingeniería
2006-06-24 09:54:16
Articles
text/plain
https://cimec.org.ar/ojs/index.php/cmm/article/view/2698
Cuadernos de Matemática y Mecánica; 2006
eng
Copyright (c)
oai:ojs.www.cimec.org.ar:article/3003
2012-02-21T00:05:49Z
cmm:ART
nmb a2200000Iu 4500
dc
Modelling and High Performance Simulation of Electrophoretic Techniques in Microfluidic Chips
Kler, Pablo
Berli, Claudio
Gusrnieri, Fabio Ariel
Electrophoretic separations comprise a group of analytical tech-
niques such as capillary zone electrophoresis, isoelectric focusing, isota-
chophoresis and free ﬂow electrophoresis. These techniques have been
miniaturized in the last years and now represent one of the most im-
portant applications of the lab-on-a-chip technology. A 3D and time-
dependent numerical model of electrophoresis on microﬂuidic devices is
presented. The model is based on the set of equations that governs elec-
trical phenomena, ﬂuid dynamics, mass transport and chemical reactions.
The relationship between the buﬀer characteristics (ionic strength, pH)
and surface potential of channel walls is taken into consideration. Numer-
ical calculations were performed by using PETSc-FEM, in a Python en-
vironment, employing high performance parallel computing. The method
includes a set of last generation preconditioners and solvers, especially
addressed to 3D microﬂuidic problems, which signiﬁcantly improve the
numerical eﬃciency in comparison with typical commercial software for
multiphysics. In this work, after discussing two validation examples, the
numerical prototyping of a microﬂuidic chip for two-dimensional elec-
trophoresis is presented.
Centro Internacional de Métodos Numéricos en Ingeniería
2012-02-20 21:06:46
Articles
application/pdf
https://cimec.org.ar/ojs/index.php/cmm/article/view/3003
Cuadernos de Matemática y Mecánica; 2011
eng
Copyright (c)
oai:ojs.www.cimec.org.ar:article/875
2007-02-21T17:23:26Z
cmm:ART
nmb a2200000Iu 4500
dc
MPI for Python: Performance Improvements and MPI2 Extensions
Dalcin, Lisandro Daniel
Storti, Mario Alberto
Paz, Rodrigo Rafael
D'Elia, Jorge
MPI for Python provides bindings of the Message Passing Interface (MPI) standard for the Python programming language and allows any Python program to exploit multiple processors. In its first release, MPI for Python was constructed on top of the MPI-1 speci- fication defining an object oriented interface that closely followed the MPI-2 C++ bindings, and provided support for communications of general Python objects. In the latest release, this package is improved to enable direct blocking/nonblocking communication of numeric arrays, and to support almost all MPI-2 features. Improvements in communication performance have been tested in a Beowulf class cluster. Results showed a negligible overhead in comparison to compiled C code. MPI for Python is open source and available for download on the web (http: //www.cimec.org.ar/python) [submitted to Journal of Parallel and Distributed Computing (2006)]
Centro Internacional de Métodos Numéricos en Ingeniería
2006-06-24 09:54:16
Articles
application/pdf
https://cimec.org.ar/ojs/index.php/cmm/article/view/875
Cuadernos de Matemática y Mecánica; 2006
eng
Copyright (c)
oai:ojs.www.cimec.org.ar:article/1660
2009-04-28T14:55:05Z
cmm:ART
nmb a2200000Iu 4500
dc
Computing Ship Wave Resistance from Wave Amplitude with a Non-Local Absorbing Boundary Condition
Storti, Mario Alberto
D'Elía, Jorge
Idelsohn, Sergio
A method for computing ship wave resistance from a momentum flux balance is presented. It is based on computing the momentum flux carried by the gravity waves that exit the computational domain through the outlet plane. It can be shown that this method ensures a non-negative wave-resistance, in contrast with straightforward integration of the normal pressure forces. However, this calculation should be performed on a transverse plane located far behind the ship. Traditional Dawson-like methods add a numerical viscosity that dampens the wave pattern so that some amount of momentum flux is lost, and resulting in an error in the momentum balance. The flow field is computed, then, with a centred scheme with absorbing boundary conditions.
Centro Internacional de Métodos Numéricos en Ingeniería
2006-06-24 11:50:44
Articles
application/pdf
https://cimec.org.ar/ojs/index.php/cmm/article/view/1660
Cuadernos de Matemática y Mecánica; 2001
eng
Copyright (c)
oai:ojs.www.cimec.org.ar:article/2536
2009-08-14T21:49:12Z
cmm:ART
nmb a2200000Iu 4500
dc
Convergence of Adaptive Finite Element Methods
Morin, Pedro
Nochetto, Ricardo H.
Siebert, Kunibert G.
Adaptive finite element methods (FEM) have been widely used in applications for over 20 years now. In practice, they converge starting from coarse grids, although no mathematical theory has been able to prove this assertion. Ensuring an error reduction rate based on a posteriori error estimators, together with a reduction rate of data oscillation (information missed by the underlying averaging process), we construct a simple and efficient adaptive FEM for elliptic partial differential equations. We prove that this algorithm converges with linear rate without any preliminary mesh adaptation nor explicit knowledge of constants. Any prescribed error tolerance is thus achieved in a finite number of steps. A number of numerical experiments in two and three dimensions yield quasi-optimal meshes along with a competitive performance. Extensions to higher order elements and applications to saddle point problems are discussed as well.
Keywords: A posteriori error estimators, data oscillation, adaptive mesh refinement, convergence, Stokes, Uzawa
AMS Subject Classifications: 65N12, 65N15, 65N30, 65N50, 65Y20
Published: SIAM Review, 44 (2002) 631--658.
Centro Internacional de Métodos Numéricos en Ingeniería
2006-06-24 11:50:44
Articles
application/pdf
https://cimec.org.ar/ojs/index.php/cmm/article/view/2536
Cuadernos de Matemática y Mecánica; 2001
eng
Copyright (c)
oai:ojs.www.cimec.org.ar:article/2982
2011-09-26T01:17:03Z
cmm:ART
nmb a2200000Iu 4500
dc
On convex regression estimators
Aguilera, Néstor E.
Forzani, Liliana
Morin, Pedro
A new nonparametric estimator of a convex regression function in any dimension is proposed and its convergence properties are studied. We start by using any estimator of the regression function and we convexify it by taking the convex envelope of a sample of the approximation obtained.
We prove that the uniform rate of convergence of the estimator is maintained after the convexification is applied.
The finite sample properties of the new estimator are investigated by means of a simulation study and the application of the new method is demonstrated in examples.
Keywords: approximation, convex regression, convexity, data-smoothing, nonparametric regression
Submitted: mayo 2010.
Centro Internacional de Métodos Numéricos en Ingeniería
2010-06-07 08:59:51
Articles
application/pdf
https://cimec.org.ar/ojs/index.php/cmm/article/view/2982
Cuadernos de Matemática y Mecánica; 2010
eng
Copyright (c)
oai:ojs.www.cimec.org.ar:article/267
2006-09-13T12:17:20Z
cmm:ART
nmb a2200000Iu 4500
dc
Multiresolution Approximations and Unconditional Bases on Weighted Lebesgue Spaces on Spaces of Homogeneous Type
Aimar, Hugo A.
Bernardis, Ana L.
Iaffei, Bibiana
Starting from a slight modification of the dyadic sets introduced by M. Christ in [M. Christ, A T(b) theorem with remarks on analytic capacity and the Cauchy integral, Colloq. Math. 60/61 (1990) 601–628] on a space of homogeneous type (X, d, µ), an MRA type structure and a Haar system H controled by the quasi distance d, can be constructed in this general setting in such a way that H is an orthonormal basis for L2(dµ). This paper is devoted
to explore under which conditions on the measure measure nu, the system H is also an unconditional basis for the Lebesgue spaces Lp(dnu). As a consequence, we obtain a characterization of these spaces in terms of the H–coefficients.
Centro Internacional de Métodos Numéricos en Ingeniería
2006-06-24 09:54:16
Articles
application/pdf
https://cimec.org.ar/ojs/index.php/cmm/article/view/267
Cuadernos de Matemática y Mecánica; 2006
eng
Copyright (c)
oai:ojs.www.cimec.org.ar:article/2396
2009-12-04T16:20:34Z
cmm:ART
nmb a2200000Iu 4500
dc
Likelihood-based sufficient dimension reduction
Cook, Dennis
Forzani, Liliana
We obtain the maximum likelihood estimator of the central subspace under conditional normality of the predictors given the response. Analytically and in simulations we found that our new estimator can preform much better than sliced inverse regression, sliced average variance estimation and directional regression, and that it seems quite robust to deviations from normality.
Published: J. Amer. Statist. Assoc. 104 (2009), no. 485, 197--208.
Centro Internacional de Métodos Numéricos en Ingeniería
2008-04-09 10:25:49
Articles
application/pdf
https://cimec.org.ar/ojs/index.php/cmm/article/view/2396
Cuadernos de Matemática y Mecánica; 2008
eng
Copyright (c)
oai:ojs.www.cimec.org.ar:article/2968
2009-12-04T16:24:25Z
cmm:ART
nmb a2200000Iu 4500
dc
Maximal operators, Riesz transforms and Littlewood–Paley functions associated with Bessel operators on BMO
Betancor, J.J.
Chicco Ruiz, Aníbal
Fariña, J.C.
Rodríguez-Mesa, L.
In this paper we study boundedness properties of certain harmonic analysis operators (maximal operators for heat and Poisson semigroups, Riesz transform and Littlewood-Paley g-functions) associated with Bessel operators, on the space BMO_o (R) that consists of the odd functions with bounded mean oscillation on R
Published: Journal of Mathematical Analysis and Applications 363 (1), 310-326.
Centro Internacional de Métodos Numéricos en Ingeniería
2008-04-09 10:25:49
Articles
application/pdf
https://cimec.org.ar/ojs/index.php/cmm/article/view/2968
Cuadernos de Matemática y Mecánica; 2008
eng
Copyright (c)
oai:ojs.www.cimec.org.ar:article/6
2007-04-27T12:12:42Z
cmm:ART
nmb a2200000Iu 4500
dc
Hot-pressing process modeling for medium density fiberboard (MDF)
Nigro, Norberto Marcelo
Storti, Mario Alberto
In this paper we present a numerical solution for the mathematical modeling of the hot-pressing process applied to medium density fiber- board. The model is based in the work of Humphrey [1982], Humphrey and Bolton [1989] and Carvalho and Costa [1998], with some modifica- tions and extensions in order to take into account mainly the convec- tive effects on the phase change term and also a conservative numerical treatment of the resulting system of partial differential equations. [International Journal of Mathematics and Mathematical Sciences 26, pp. 713-729 (2001)]
Centro Internacional de Métodos Numéricos en Ingeniería
2006-06-24 11:50:44
Articles
application/pdf
https://cimec.org.ar/ojs/index.php/cmm/article/view/6
Cuadernos de Matemática y Mecánica; 2001
eng
Copyright (c)
oai:ojs.www.cimec.org.ar:article/1383
2009-12-04T16:20:34Z
cmm:ART
nmb a2200000Iu 4500
dc
High Performance Simulations of Electrokinetic Flow and Transport in Microfluidic Chips
Kler, Pablo
Dalcín, Lisandro
Guarnieri, Fabio Ariel
Storti, Mario Alberto
This article discusses high performance numerical simulations of electroki- netic flow and transport phenomena in microfluidic chips. Modeling grounds on conservation equations of mass, momentum and electric charge in the framework of continuum mechanics. Two examples of interest in microfluidics are consid- ered as study cases. Three dimension effects and whole chip geometries are taking into account. All numerical simulations presented are performed with PETSc-FEM within a Python programming environment employing parralel computing. Computation time and parallel efficiency are measured in order to study additive Schwarz method performance as domain decomposition technique in solving common ill-conditioned microfluidics problems. [Submitted to Comp. Method in Applied Mech Engng]
Centro Internacional de Métodos Numéricos en Ingeniería
2008-04-09 10:25:49
Articles
application/pdf
https://cimec.org.ar/ojs/index.php/cmm/article/view/1383
Cuadernos de Matemática y Mecánica; 2008
eng
Copyright (c)
oai:ojs.www.cimec.org.ar:article/2387
2009-12-04T16:20:34Z
cmm:ART
nmb a2200000Iu 4500
dc
Weighted inequalities for negative powers of Schrödinger operators
Bongioanni, Bruno
Harboure, Eleonor
Salinas, Oscar
In this article we obtain boundedness of the operator $(-\Delta~+~V)^{-\alpha/2}$ from $L^{p,\infty}(w)$ into weighted bounded mean oscillation type spaces $BMO_\LL^{\beta}(w)$ under appropriate conditions on the weight $w$. We also show that these weighted spaces also have a point-wise description for $0 < \beta < 1$. Finally, we study the behaviour of the operator $(-\Delta~+~V)^{-\alpha/2}$ when acting on $BMO_\LL^{\beta}(w)$.
Published: J. Math. Anal. Appl. 348 (2008), no. 1, 12--27.
Centro Internacional de Métodos Numéricos en Ingeniería
2008-04-09 10:25:49
Articles
application/pdf
https://cimec.org.ar/ojs/index.php/cmm/article/view/2387
Cuadernos de Matemática y Mecánica; 2008
eng
Copyright (c)
oai:ojs.www.cimec.org.ar:article/2704
2010-06-08T22:34:30Z
cmm:ART
nmb a2200000Iu 4500
dc
Constructive logic with strong negation as a substructural logic
Busaniche, Manuela
Cignoli, Roberto
Spinks and Veroff have shown that constructive logic with strong negation (CLSN for short), can be considered as a substructural logic. We use algebraic tools developed to study substructural logics to investigate some axiomatic extensions of CLSN. For instance we prove that Nilpotent Minimum Logic is the extension of CLSN by the prelinearity axiom. This generalizes the well known result by Monteiro and Vakarelov that three-valued \Lukasiewicz logic is an extension of CLSN. A Glivenko-like theorem relating CLSN and three-valued \Lukasiewicz logic is proved.
Accepted: Journal of Logic and Computation ; doi: 10.1093/logcom/exn081.
Centro Internacional de Métodos Numéricos en Ingeniería
2009-04-08 20:22:29
Articles
text/plain
https://cimec.org.ar/ojs/index.php/cmm/article/view/2704
Cuadernos de Matemática y Mecánica; 2009
eng
Copyright (c)
oai:ojs.www.cimec.org.ar:article/3729
2012-02-21T00:06:08Z
cmm:ART
nmb a2200000Iu 4500
"110215 2011 eng "
dc
Parabolic Besov Regularity for the Heat Equation
Aimar, Hugo
Gómez, Ivana
We obtain parabolic Besov smoothness improvement for temperatures on cylindrical regions based on Lipschitz domains. The results extend those for harmonic functions obtained by S. Dahlke and R. DeVore in “Besov regularity for elliptic boundary value problems” published in Communications in Partial Differential Equations 22, no.1-2, 1-16, 1997.
Centro Internacional de Métodos Numéricos en Ingeniería
2012-02-20 21:06:46
Articles
application/pdf
https://cimec.org.ar/ojs/index.php/cmm/article/view/3729
Cuadernos de Matemática y Mecánica; 2011
eng
Copyright (c)
oai:ojs.www.cimec.org.ar:article/1018
2009-09-07T19:30:44Z
cmm:ART
nmb a2200000Iu 4500
dc
Numerical simulation of the Ahmed vehicle model near-wake
Franck, Gerardo
Nigro, Norberto M.
Storti, Mario Alberto
D'Elia, Jorge
The near-wake structure of the flow around the Ahmed vehicle model
is numerically achieved by a time-averaged procedure of the unsteady
flow modeled by the Navier-Stokes equations with a Large Eddy
Simulation (LES) model for the turbulence . The Reynolds (based on
the model length) and the Mach numbers are fixed in 4.25 million and
0.18, respectively. A viscous and incompressible fluid model of
Newtonian type is assumed. A LES technique together with the
Smagorinsky model as Subgrid Scale Modeling (SGM) and the van
Driest near-wall damping is employed. The coherent macro structures in
the near-wake are estimated through the second invariant of the
velocity gradient (Q-criterion) applied on the time-average
flow. A monolithic computational code is employed, which is based on
the finite element method with equal order basis functions (linear)
for pressure and velocity and uses a Streamline Upwind Petrov-Galerkin
(SUPG) scheme combined with a Pressure Stabilizing
Petrov-Galerkin(PSPG) one. Parallel computing on a Beowulf
cluster with a domain decomposition technique for solving the
algebraic system is used. [Submitted to Int J Num Meth Fluids]
Centro Internacional de Métodos Numéricos en Ingeniería
2007-03-06 10:29:57
Articles
application/pdf
https://cimec.org.ar/ojs/index.php/cmm/article/view/1018
Cuadernos de Matemática y Mecánica; 2007
eng
Copyright (c)
252c519e67c29ad54097341e9c092edc