Graphic summary
  • Show / hide key
  • Information


Scientific and technological production
  •  

1 to 50 of 718 results
  • Parametric solutions involving geometry: a step towards efficient shape optimization

     Ammar, Amine; Huerta Cerezuela, Antonio; Chinesta, Francisco; Cueto Prendes, Elias; Leygue, Adrien
    Computer methods in applied mechanics and engineering
    Date of publication: 2014-01
    Journal article

    Read the abstract Read the abstract View View Open in new window  Share Reference managers Reference managers Open in new window

    Optimization of manufacturing processes or structures involves the optimal choice of many parameters (process parameters, material parameters or geometrical parameters). Usual strategies proceed by defining a trial choice of those parameters and then solving the resulting model. Then, an appropriate cost function is evaluated and its optimality checked. While the optimum is not reached, the process parameters should be updated by using an appropriate optimization procedure, and then the model must be solved again for the updated process parameters. Thus, a direct numerical solution is needed for each choice of the process parameters, with the subsequent impact on the computing time. In this work we focus on shape optimization that involves the appropriate choice of some parameters defining the problem geometry. The main objective of this work is to describe an original approach for computing an off-line parametric solution. That is, a solution able to include information for different parameter values and also allowing to compute readily the sensitivities. The curse of dimensionality is circumvented by invoking the Proper Generalized Decomposition (PGD) introduced in former works, which is applied here to compute geometrically parametrized solutions. © 2013 Elsevier B.V.

  • Real-time wave monitoring in Harbors

     Huerta Cerezuela, Antonio; Modesto Galende, David; Zlotnik, Sergio
    International Conference on Computational Engineering and Science for Safety and Environmental Problems
    Presentation's date: 2014-04-15
    Presentation of work at congresses

     Share Reference managers Reference managers Open in new window

  • Access to the full text
    Stability of anchored sheet wall in cohesive-frictional soils by FE limit analysis  Open access

     Muñoz Romero, Jose Javier; Lyamin, Andrei; Huerta Cerezuela, Antonio
    International journal for numerical and analytical methods in geomechanics
    Date of publication: 2013-06-03
    Journal article

    Read the abstract Read the abstract Access to the full text Access to the full text Open in new window  Share Reference managers Reference managers Open in new window

    This study extends the limit analysis techniques used for the computation of strict bounds of the load factors in solids to stability problems with interfaces, anchors and joints. The cases considered include the pull-out capacity of multibelled anchors and the stability of retaining walls for multiple conditions at the anchor/soil and wall/soil interfaces. Three types of wall supports are examined: free standing wall, simply supported wall and anchored wall. The results obtained are compared against available experimental and numerical data. The conclusion drawn confirms the validity of numerical limit analysis for the computation of accurate bounds on limit loads and capturing failure modes of structures with multiple inclusions of complex interface and support conditions.

  • Dimensionless analysis of HSDM and application to simulation of breakthrough curves of highly adsorbent porous media

     Perez Foguet, Agusti; Casoni Rero, Eva; Huerta Cerezuela, Antonio
    Journal of environmental engineering (ASCE)
    Date of publication: 2013-05
    Journal article

    Read the abstract Read the abstract View View Open in new window  Share Reference managers Reference managers Open in new window

    The homogeneous surface diffusion model (HSDM) is widely used for adsorption modeling of aqueous solutions. The Biot number is usually used to characterize model behavior. However, some limitations of this characterization have been reported recently, and the Stanton number has been proposed as a complement to be considered. In this work, a detailed dimensionless analysis of HSDM is presented and limit behaviors of the model are characterized, confirming but extending previous results. An accurate and efficient numerical solver is used for these purposes. The intraparticle diffusion equation is reduced to a system of two ordinary differential equations, the transport-reaction equation is discretized by using a discontinuous Galerkin method, and the overall system evolution is integrated with a time-marching scheme. This approach facilitates the simulation of HSDM with a wide range of dimensionless numbers and with a correct treatment of shocks, which appear with nonlinear adsorption isotherms and with large Biot numbers and small surface diffusivity modulus. The approach is applied to simulate the breakthrough curves of granular ferric hydroxide. Published experimental data is adequately simulated.

  • Efficiency of high-order elements for continuous and discontinuous Galerkin methods

     Huerta Cerezuela, Antonio; Angeloski, Aleksandar; Roca Navarro, Xevi; Peraire Guitart, Jaume
    International journal for numerical methods in engineering
    Date of publication: 2013
    Journal article

    View View Open in new window  Share Reference managers Reference managers Open in new window

  • High-order continuous and discontinuous Galerkin methods for wave problems

     Giorgiani, Giorgio; Modesto Galende, David; Fernandez Mendez, Sonia; Huerta Cerezuela, Antonio
    International journal for numerical methods in fluids
    Date of publication: 2013-12-10
    Journal article

    Read the abstract Read the abstract View View Open in new window  Share Reference managers Reference managers Open in new window

    Three Galerkin methods-continuous Galerkin, Compact Discontinuous Galerkin, and hybridizable discontinuous Galerkin-are compared in terms of performance and computational efficiency in 2-D scattering problems for low and high-order polynomial approximations. The total number of DOFs and the total runtime are used for this correlation as well as the corresponding precision. The comparison is carried out through various numerical examples. The superior performance of high-order elements is shown. At the same time, similar capabilities are shown for continuous Galerkin and hybridizable discontinuous Galerkin, when high-order elements are adopted, both of them clearly outperforming compact discontinuous Galerkin.

    hree Galerkin methods—continuous Galerkin, Compact Discontinuous Galerkin, and hybridizable discontinuous Galerkin—are compared in terms of performance and computational efficiency in 2-D scattering problems for low and high-order polynomial approximations. The total number of DOFs and the total runtime are used for this correlation as well as the corresponding precision. The comparison is carried out through various numerical examples. The superior performance of high-order elements is shown. At the same time, similar capabilities are shown for continuous Galerkin and hybridizable discontinuous Galerkin, when high-order elements are adopted, both of them clearly outperforming compact discontinuous Galerkin

  • Hybridizable discontinuous Galerkin p-adaptivity for wave propagation problems

     Giorgiani, Giorgio; Fernandez Mendez, Sonia; Huerta Cerezuela, Antonio
    International journal for numerical methods in fluids
    Date of publication: 2013
    Journal article

    Read the abstract Read the abstract View View Open in new window  Share Reference managers Reference managers Open in new window

    A p-adaptive hybridizable discontinuous Galerkin method for the solution of wave problems is presented in a challenging engineering problem. Moreover, its performance is compared with a high-order continuous Galerkin. The hybridization technique allows to reduce the coupled degrees of freedom to only those on the mesh element boundaries, whereas the particular choice of the numerical fluxes opens the path to a superconvergent postprocessed solution. This superconvergent postprocessed solution is used to construct a simple and inexpensive error estimator. The error estimator is employed to obtain solutions with the prescribed accuracy in the area (or areas) of interest and also drives a proposed iterative mesh adaptation procedure. The proposed method is applied to a nonhomogeneous scattering problem in an unbounded domain. This is a challenging problem because, on the one hand, for high frequencies, numerical difficulties are an important issue because of the loss of the ellipticity and the oscillatory behavior of the solution. And on the other hand, it is applied to real harbor agitation problems. That is, the mild slope equation in frequency domain (Helmholtz equation with nonconstant coefficients) is solved on real geometries with the corresponding perfectly matched layer to damp the diffracted waves. The performance of the method is studied on two practical examples. The adaptive hybridizable discontinuous Galerkin method exhibits better efficiency compared with a high-order continuous Galerkin method using static condensation of the interior nodes.

  • PGD-based ¿Computational Vademecum¿ for efficient design, optimization and control

     Chinesta, Francisco; Leygue, Adrien; Bordeu, F.; Cueto Prendes, Elias; González, David; Alfaro, Icíar; Ammar, Amine; Huerta Cerezuela, Antonio; Diez Mejia, Pedro
    IACM expressions: Bulletin for the international association for computational mechanics
    Date of publication: 2013
    Journal article

    Read the abstract Read the abstract View View Open in new window  Share Reference managers Reference managers Open in new window

    D espite the impressive progresses attained by simulation capabilities and techniques, some challenging problems remain today intractable. These prob- lems, that are common to many branches of science and engineering, are of differ- ent nature. Among them, we can cite those related to high-dimensional models, on which mesh-based approaches fail due to the exponential increase of degrees of freedom. Other challenging scenarios concern problems requiring many direct solutions (optimization, inverse identifica- tion, uncertainty quantification ! ) or those needing very fast solutions (real time simulation, simulation based control ! ).

  • On the natural stabilization of convection dominated problems using high order Bubnov¿Galerkin finite elements

     Cai, Quanji; Kollmannsberger, Stefan; Sala Lardies, Esther; Huerta Cerezuela, Antonio; Rank, Ernst
    Computers & mathematics with applications
    Date of publication: 2013-10-15
    Journal article

    View View Open in new window  Share Reference managers Reference managers Open in new window

  • Computable exact bounds for linear outputs from stabilized solutions of the advection-diffusion-reaction equation

     Parés Mariné, Núria; Diez Mejia, Pedro; Huerta Cerezuela, Antonio
    International journal for numerical methods in engineering
    Date of publication: 2013-02-01
    Journal article

    View View Open in new window  Share Reference managers Reference managers Open in new window

  • Streamline upwind/Petrov¿Galerkin-based stabilization of proper generalized decompositions for high-dimensional advection¿diffusion equations

     González, David; Cueto, Elias; Chinesta, Francisco; Diez Mejia, Pedro; Huerta Cerezuela, Antonio
    International journal for numerical methods in engineering
    Date of publication: 2013-06
    Journal article

    Read the abstract Read the abstract View View Open in new window  Share Reference managers Reference managers Open in new window

    This work is a first attempt to address efficient stabilizations of high dimensional advection-diffusion models encountered in computational physics. When addressing multidimensional models, the use of mesh-based discretization fails because the exponential increase of the number of degrees of freedom related to a multidimensional mesh or grid, and alternative discretization strategies are needed. Separated representations involved in the so-called proper generalized decomposition method are an efficient alternative as proven in our former works; however, the issue related to efficient stabilizations of multidimensional advection-diffusion equations has never been addressed to our knowledge. Thus, this work is aimed at extending some well-experienced stabilization strategies widely used in the solution of 1D, 2D, or 3D advection-diffusion models to models defined in high-dimensional spaces, sometimes involving tens of coordinates

  • PGD-Based computational vademecum for efficient design, optimization and control

     Chinesta, Francisco; Leygue, Adrien; Bordeu, Felipe; Aguado, José Vicente; Cueto, Elias; González, David; Alfaro, Icíar; Ammar, Amine; Huerta Cerezuela, Antonio
    Archives of computational methods in engineering
    Date of publication: 2013
    Journal article

    View View Open in new window  Share Reference managers Reference managers Open in new window

  • Hybridizable discontinuous Galerkin with degree adaptivity for wave and flow problems

     Giorgiani, Giorgio; Fernandez Mendez, Sonia; Huerta Cerezuela, Antonio
    Computational Methods for Coupled Problems in Science and Engineering
    Presentation's date: 2013-06-17
    Presentation of work at congresses

     Share Reference managers Reference managers Open in new window

  • Access to the full text
    Decomposition techniques in computational limit analysis  Open access

     Muñoz Romero, Jose Javier; Rabiei, Syednima; Huerta Cerezuela, Antonio
    International Conference on Computational Plasticity Fundamentals and Applications
    Presentation's date: 2013-09-03
    Presentation of work at congresses

    Read the abstract Read the abstract Access to the full text Access to the full text Open in new window  Share Reference managers Reference managers Open in new window

    Numerical techniques for the computation of strict bounds in limit analyses have been developed for more than thirty years. The efficiency of these techniques have been substantially improved in the last ten years, and have been successfully applied to academic problems, foundations and excavations. We here extend the theoretical background to problems with anchors, interface conditions, and joints. Those extensions are relevant for the analysis of retaining and anchored walls, which we study in this work. The analysis of three-dimensional domains remains as yet very scarce. From the computational standpoint, the memory requirements and CPU time are exceedingly prohibitive when mesh adaptivity is employed. For this reason, we also present here the application of decomposition techniques to the optimisation problem of limit analysis. We discuss the performance of different methodologies adopted in the literature for general optimisation problems, such as primal and dual decomposition, and suggest some strategies that are suitable for the parallelisation of large three-dimensional problems. The proposed decomposition techniques are tested against representative problems.

    Numerical techniques for the computation of strict bounds in limit analyses have been developed for more than thirty years. The efficiency of these techniques have been substantially improved in the last ten years, and have been successfully applied to academic problems, foundations and excavations. We here extend the theoretical background to problems with anchors, interface conditions, and joints. Those extensions are relevant for the analysis of retaining and anchored walls, which we study in this work. The analysis of three-dimensional domains remains as yet very scarce. From the computational standpoint, the memory requirements and CPU time are exceedingly prohibitive when mesh adaptivity is employed. For this reason, we also present here the application of decomposition techniques to the optimisation problem of limit analysis. We discuss the performance of different methodologies adopted in the literature for general optimisation problems, such as primal and dual decomposition, and suggest some strategies that are suitable for the parallelisation of large three-dimensional problems. The proposed decomposition techniques are tested against representative problems.

  • Access to the full text
    Deim-based PGD for parametric nonlinear model order reduction  Open access

     Aguado, José V.; Chinesta, Francisco; Leygue, Adrien; Cueto Prendes, Elias; Huerta Cerezuela, Antonio
    International Conference on Adaptive Modeling and Simulation
    Presentation's date: 2013-06-04
    Presentation of work at congresses

    Read the abstract Read the abstract Access to the full text Access to the full text Open in new window  Share Reference managers Reference managers Open in new window

    Anewtechniqueforefficientlysolvingparametricnonlinearreducedorder modelsintheProperGeneralizedDecomposition(PGD)frameworkispresentedhere. ThistechniqueisbasedontheDiscreteEmpiricalInterpolationMethod(DEIM)[1],and thusthenonlineartermisinterpolatedusingthereducedbasisinsteadofbeingfully evaluated.TheDEIMhasalreadybeendemonstratedtoprovidesatisfactoryresultsin termsofcomputationalcomplexitydecreasewhencombinedwiththeProperOrthogonal Decomposition(POD).However,inthePODcasethereducedbasisis aposteriori known asitcomesfromseveralpre-computedsnapshots.Onthecontrary,thePGDisan apriori modelreductionmethod.ThismakestheDEIM-PGDcouplingratherdelicate,because differentchoicesarepossibleasitisanalyzedinthiswork

  • Improving numerical efficiency with model reduction and High-Order Adaptative Discontinuous Galerkin

     Huerta Cerezuela, Antonio; Modesto Galende, David; Zlotnik, Sergio; Cueto Prendes, Elias; Chinesta, Francisco
    International Conference on Adaptive Modeling and Simulation
    Presentation's date: 2013-06-04
    Presentation of work at congresses

    View View Open in new window  Share Reference managers Reference managers Open in new window

  • Efficiency of high-order elements for continuous and disconrinuous Galerkin

     Angeloski, Aleksandar; Roca Navarro, Xevi; Peraire Guitart, Jaume; Huerta Cerezuela, Antonio
    Congreso de Métodos Numéricos en Ingeniería
    Presentation's date: 2013-06-26
    Presentation of work at congresses

    View View Open in new window  Share Reference managers Reference managers Open in new window

  • X-FEM de alto orden para problemas con inclusiones materiales

     Sala Lardies, Esther; Fernandez Mendez, Sonia; Huerta Cerezuela, Antonio
    Congreso de Métodos Numéricos en Ingeniería
    Presentation's date: 2013-06-26
    Presentation of work at congresses

    View View Open in new window  Share Reference managers Reference managers Open in new window

  • On the numerical efficiency of model reduction and high-order discontinuous Galerkin for waves and fluid problems

     Modesto Galende, David; Giorgiani, Giorgio; Fernandez Mendez, Sonia; Huerta Cerezuela, Antonio
    Congreso de Métodos Numéricos en Ingeniería
    Presentation's date: 2013-06-26
    Presentation of work at congresses

    View View Open in new window  Share Reference managers Reference managers Open in new window

  • Advanced harmonic techniques for solving the transient heat equation

     Aguado, José Vicente; Huerta Cerezuela, Antonio; Chinesta, Francisco; Cueto Prendes, Elias; Leygue, Adrien
    Congreso de Métodos Numéricos en Ingeniería
    Presentation's date: 2013-06-26
    Presentation of work at congresses

    Read the abstract Read the abstract View View Open in new window  Share Reference managers Reference managers Open in new window

    Simulation-based management of processes has always been a big challenge because of the difficulty to merge two apparently opposed idea s: the real time handling demanded by the process, and the traditionally strong time c onsumption linked to the sim- ulation. However, very often processes are monitored or con trolled by performing relatively few measures on some specific points of the domain. Then, one s trategy for accelerating the simulation might be based on trying to compute the respon se only where it is needed. For this reason this work revisits a classical technique, th e harmonic analysis, which is combined with recent model reduction techniques giving pla ce to an advanced harmonic technique for solving linear parabolic equations and, in pa rticular, the transient heat equa- tion. Thanks to fact that Reciprocity principle applies in t he frequency domain, the solu- tion can be computed where it is of interest by performing an i nexpensive postprocessing of a pre-computed transfer function. This transfer functio n, that must be valid for a suffi- ciently large range of frequencies, is efficiently computed b y using the Proper Generalized Decomposition method.

  • Adaptive hybrid discontinuous methods for fluid and wave problems  Open access

     Giorgiani, Giorgio
    Defense's date: 2013-04-02
    Department of Applied Mathematics III, Universitat Politècnica de Catalunya
    Theses

    Read the abstract Read the abstract Access to the full text Access to the full text Open in new window  Share Reference managers Reference managers Open in new window

    Esta tesis propone una técnica de adaptividad p para el método de Galerkin discontinuo híbrido (HDG)El método HDG es un nuevo método de Galerkin discontinuo (DG) con características interesantes. Tiene todas las ventajas de los métodos discontinuos, como la estabilización inherente y las propiedades conservación local, y además permite reducir los grados de libertad acoplados del problema a aquellos de una aproximación de la solución definida sólo en las caras de la malla . Además, las propiedades de convergencia de la solución de HDG permiten postprocesar la solución elemento a elemento, obteniendo una solución superconvergente.Debido al carácter discontinuo de la aproximación en HDG, es posible implementar con simplicidad cálculos con p-variable. En esta tesis el postproceso es superconvergente se utiliza para definir un estimador de error fiable y computacionalmente barato, que se utiliza para guiar un proceso de adaptativo automatico. Se modifica el grado del polinomio en cada elemento con el objetivo de obtener una distribución de error uniforme por debajo de una tolerancia definida por usuario. No se aplican modificas topológicas a la discretización, entonces el proceso adaptativo es muy rapido.Al principio, la tecnica de HDG p-adaptativo se aplica a la solución de problemas de onda. En particular, la ecuación de ¿pendiente suave¿ se utiliza para modelar el problema de propagación de ondas marina en áreas costeras y puertos. El método de HDG se compara con el método de elementos finitos Galerkin continuo (CG), que es hoy en día el método más común utilizado en la práctica de la ingenieril para este tipo de aplicaciones. Experimentos numéricos revelan que la eficiencia de HDG es paracida a la eficiencia de CG para cálculos de grado uniforme, y claramente superior a otros métodos DG como el método de Galerkin discontinuo compacto.Cuando se considere p-adaptatividad, se obtiene un importante ahorro en el coste computacional.Luego, la metodología se aplica a la solución de las ecuaciones de Navier-Stokes para flujo incompresible laminare. Se consideran aplicaciones tanto en estado estacionario como en transitorio. Se presentan varios experimentos numéricos en 2D y 3D, incluyendo ejemplos académicos y aplicaciones de interese ingenieril. A pesar de la simplicidad y bajo coste del estimador de error, se demuestra su alta eficiencia usando ejemplos analíticos. Además, aunque la técnica adaptativa sea basada en una estimación de error para sólo el campo de velocidad, se alcanza alta precisión para todas las variables, y una evaluación precisa de las fuerzas fluidodinamicas. En particular, los altos grados polinomiales son automáticamente situado a lo largo de capas límite, lo que reduce la necesidad de elementos altamente distorsionados. Pruebas numéricas muestran una importante reducción en el coste computacional, comparado con cálculos de grado uniforme, para los cálculos estacionarios y transitorios.

    This PhD thesis proposes a p-adaptive technique for the Hybridizable Discontinuous Galerkin method (HDG). The HDG method is a novel discontinuous Galerkin method (DG) with interesting characteristics. While retaining all the advantages of the common DG methods, such as the inherent stabilization and the local conservation properties, HDG allows to reduce the coupled degrees of freedom of the problem to those of an approximation of the solution de¿ned only on the faces of the mesh. Moreover, the convergence properties of the HDG solution allow to perform an element-by-element postprocess resulting in a superconvergent solution. Due to the discontinuous character of the approximation in HDG, p-variable computations are easily implemented. In this work the superconvergent postprocess is used to de¿ne a reliable and computationally cheap error estimator, that is used to drive an automatic adaptive process. The polynomial degree in each element is automatically adjusted aiming at obtaining a uniform error distribution below a user de¿ned tolerance. Since no topological modi¿cation of the discretization is involved, fast adaptations of the mesh are obtained. First, the p-adaptive HDG is applied to the solution of wave problems. In particular, the Mild Slope equation is used to model the problem of sea wave propagation is coastal areas and harbors. The HDG method is compared with the continuous Galerkin (CG) ¿nite element method, which is nowadays the common method used in the engineering practice for this kind of applications. Numerical experiments reveal that the e¿ciency of HDG is close to CG for uniform degree computations, clearly outperforming other DG methods such as the Compact Discontinuous Galerkin method. When p-adaptivity is considered, an important saving in computational cost is shown. Then, the methodology is applied to the solution of the incompressible Navier-Stokes equations for the simulation of laminar ¿ows. Both steady state and transient applications are considered. Various numerical experiments are presented, in 2D and 3D, including academic examples and more challenging applications of engineering interest. Despite the simplicity and low cost of the error estimator, high e¿ciency is exhibited for analytical examples. Moreover, even though the adaptive technique is based on an error estimate for just the velocity ¿eld, high accuracy is attained for all variables, with sharp resolution of the key features of the ¿ow and accurate evaluation of the ¿uid-dynamic forces. In particular, high degrees are automatically located along boundary layers, reducing the need for highly distorted elements in the computational mesh. Numerical tests show an important reduction in computational cost, compared to uniform degree computations, for both steady and unsteady computations.

  • Proper generalized decomposition based dynamic data-driven control of thermal processes

     Ghnatios, Ch.; Masson, F.; Huerta Cerezuela, Antonio; Leygue, Adrien; Cueto, E.; Chinesta, Francisco
    Computer methods in applied mechanics and engineering
    Date of publication: 2012
    Journal article

    View View Open in new window  Share Reference managers Reference managers Open in new window

  • Are High-order and Hybridizable Discontinuous Galerkin methods competitive?

     Huerta Cerezuela, Antonio; Angeloski, Aleksandar; Roca Navarro, Xevi; Peraire Guitart, Jaume
    Oberwolfach reports
    Date of publication: 2012
    Journal article

    View View Open in new window  Share Reference managers Reference managers Open in new window

  • A note on upper bound formulations in limit analysis

     Muñoz Romero, Jose Javier; Huerta Cerezuela, Antonio; Bonet, Javier; Peraire, Jaume
    International journal for numerical methods in engineering
    Date of publication: 2012-08-24
    Journal article

    Read the abstract Read the abstract View View Open in new window  Share Reference managers Reference managers Open in new window

    In this paper we study some recent formulations for the computation of upper bounds in limit analysis. We show that a previous formulation presented by the authors does not guarantee the strictness of the upper bound, nor does it provide a velocity field that satisfies the normality rule everywhere. We show that these deficiencies are related to the quadrature employed for the evaluation of the dissipation power. We derive a formulation that furnishes a strict upper bound of the load factor, which in fact coincides with a formulation reported in the literature. From the analysis of these formulations we propose a post-process which consists in computing exactly the dissipation power for the optimum upper bound velocity field. This post-process may further reduce the strict upper bound of the load factor in particular situations. Finally, we also determine the quadratures that must be used in the elemental and edge gap contributions so that they are always positive and their addition equals the global bound gap.

  • One-dimensional shock-capturing for high-order discontinuous Galerkin methods

     Casoni Rero, Eva; Peraire Guitart, Jaume; Huerta Cerezuela, Antonio
    International journal for numerical methods in fluids
    Date of publication: 2012-02
    Journal article

    View View Open in new window  Share Reference managers Reference managers Open in new window

  • Interactive simulation on a smart Phone

     Cueto, E.; Huerta Cerezuela, Antonio; Chinesta, Francisco
    Benchmark
    Date of publication: 2012-01
    Journal article

    Read the abstract Read the abstract View View Open in new window  Share Reference managers Reference managers Open in new window

    Traditionally, Simulationbased Engineering Sciences (SBES) made use of static data inputs to perform the simulations. Namely parameters of the model, boundary conditions, etc. were traditionally obtained by experimentation and could not be modified during the course of the simulation. More recently, large efforts have been invested in developing dynamic datadriven application systems (DDDAS): systems in which measurements and simulations are continuously influencing each other in a symbiotic manner. It should be understood that measurements should be incorporated in real time to the simulations, while simulations could eventually control the way in which measurements are done.

  • Estimation of crack opening from a two-dimensional continuum-based finite element computation

     Dufour, F.; Legrain, Dominique; Pijaudier Cabot, Gilles; Huerta Cerezuela, Antonio
    International journal for numerical and analytical methods in geomechanics
    Date of publication: 2012-11
    Journal article

    View View Open in new window  Share Reference managers Reference managers Open in new window

  • Access to the full text
    High-order implicit time integration for unsteady incompressible flows  Open access

     de Villardi de Montlaur, Adeline; Fernandez Mendez, Sonia; Huerta Cerezuela, Antonio
    International journal for numerical methods in fluids
    Date of publication: 2012-10-20
    Journal article

    Read the abstract Read the abstract Access to the full text Access to the full text Open in new window  Share Reference managers Reference managers Open in new window

    The spatial discretization of unsteady incompressible Navier–Stokes equations is stated as a system of differential algebraic equations, corresponding to the conservation of momentum equation plus the constraint due to the incompressibility condition. Asymptotic stability of Runge–Kutta and Rosenbrock methods applied to the solution of the resulting index-2 differential algebraic equations system is analyzed. A critical comparison of Rosenbrock, semi-implicit, and fully implicit Runge–Kutta methods is performed in terms of order of convergence and stability. Numerical examples, considering a discontinuous Galerkin formulation with piecewise solenoidal approximation, demonstrate the applicability of the approaches and compare their performance with classical methods for incompressible flows.

  • Are high-order and hybridizable discontinuous Galerkin methods competitive?

     Huerta Cerezuela, Antonio; Angeloski, Aleksandar; Roca Navarro, Xevi; Peraire Guitart, Jaume
    Oberwolfach Reports
    Presentation's date: 2012-02-14
    Presentation of work at congresses

    Read the abstract Read the abstract View View Open in new window  Share Reference managers Reference managers Open in new window

    The talk covered several issues motivated by a practical engineering wave propagation problem: real-time evaluation of wave agitation in harbors. The first part, presented the application of a reduced order model in the framework of a Helmholtz equation with non-constant coefficients in an unbounded domain. This problem requires large numbers of degrees of freedom (ndof) because relatively high frequencies with small (compared with the domain size) geometric features must be considered.

  • Access to the full text
    Hybridizable discontinuous Galerkin p-adaptivity for wave problems  Open access

     Giorgiani, Giorgio; Fernandez Mendez, Sonia; Huerta Cerezuela, Antonio
    European Community on Computational Methods in Applied Sciences Young Investigators Conference
    Presentation's date: 2012-04-24
    Presentation of work at congresses

    Read the abstract Read the abstract Access to the full text Access to the full text Open in new window  Share Reference managers Reference managers Open in new window

    A p-adaptive Hybridizable Discontinuous Galerkin (HDG) method is presented for the solution of wave problems. The HDG method allows to drastically reduce the coupled degrees of freedom of the computation seeking for an approximation of the solution that is defined only on the edges of the mesh. The particular choice of the numerical fluxes driven by the hybridization technique allows to obtain an optimally converging solution not only for the primal unknown but also for its derivative. This characteristic allows to perform a post-process of the solution that provides a super-convergent solution. The discontinuous character of the solution provides an optimal framework for a p-adaptive technique. The post-processed solution of the HDG method is used to construct a cheap and reliable error estimator that drives an element by element modification of the approximation degree. The proposed p-adaptive HDG method is compared with high-order CG computation with static condensation of the interior nodes. A challenging problem is considered for the comparison: a non-homogeneous scattering problem in an open domain.

  • Efficiency and accuracy of high-order computations and reduced order modelling in coastal engineering wave propagation problems

     Modesto Galende, David; Giorgiani, Giorgio; Zlotnik, Sergio; Huerta Cerezuela, Antonio
    European Community on Computational Methods in Applied Sciences Young Investigators Conference
    Presentation's date: 2012-04-24
    Presentation of work at congresses

    Read the abstract Read the abstract View View Open in new window  Share Reference managers Reference managers Open in new window

    Several numerical issues have to be considered when solving wave propagation problems, between which there are the artificial boundary conditions, the small geometrical features that can be influential or the variable coefficients. Apart from them, two issues are mainly addressed and discussed. Firstly, low order elements need very high wave resolution for capturing the solution in the area of interest, leading to extremely dense meshes. High-order finite elements are proposed to be an efficient and accurate solution for solving the problem. Secondly, the very large number of test cases. When designing harbour models, a huge number of incident waves, in term of wavelengths and directions, have to be studied. The excessive computational cost to carry out all the possible direct problems prevents the whole data evaluation, inducing the lost of important information. Reduced order models may be an alternative if they are computable, efficient and accurate. The applicability of Proper Generalized Decomposition (PGD) is exploited. Unlike previous PGD contributions, which deal with elliptic problems, the present work is focused on a more challenging scenario for the separable representation due to the loss of the elliptic behaviour. The proposed PGD involves a separable representation of the unknown reflected wave in space, wave number and angle of incidence. Such decomposition appears to be really interesting for practical purpose, where goal-oriented results are critical for a wide range of frequencies and incident waves. Moreover, when accuracy and efficiency are of concern, the number of terms in the reduced model are determined by means of an error estimation based on the dual formulation of the problem.

  • The Proper generalized decomosition applied to a parametrized wave propagation problem

     Huerta Cerezuela, Antonio; Modesto Galende, David
    World Congress on Computational Mechanics
    Presentation's date: 2012-07-11
    Presentation of work at congresses

    View View Open in new window  Share Reference managers Reference managers Open in new window

  • The proper generalized decomposition applied to a practical wave propagation problem

     Modesto Galende, David; Zlotnik, Sergio; Huerta Cerezuela, Antonio
    European Congress on Computational Methods in Applied Sciences and Engineering
    Presentation's date: 2012-09-13
    Presentation of work at congresses

    View View Open in new window  Share Reference managers Reference managers Open in new window

  • Efficiency of high-order elements for continuous and discontinuous Galerkin methods

     Huerta Cerezuela, Antonio; Angeloski, Aleksandar; Roca Navarro, Xevi; Peraire Guitart, Jaume
    Spain-Japan Workshop on Computational Mechanics
    Presentation's date: 2012-09-17
    Presentation of work at congresses

    View View Open in new window  Share Reference managers Reference managers Open in new window

  • Hybridizable discontinuous Galerkin p-adaptivity for fluid problems

     Giorgiani, Giorgio; Fernandez Mendez, Sonia; Huerta Cerezuela, Antonio
    European Congress on Computational Methods in Applied Sciences and Engineering
    Presentation's date: 2012-09-12
    Presentation of work at congresses

    View View Open in new window  Share Reference managers Reference managers Open in new window

  • Access to the full text
    Optimally convergent high-order X-FEM for problems with voids and inclusions  Open access

     Sala Lardies, Esther; Fernandez Mendez, Sonia; Huerta Cerezuela, Antonio
    European Congress on Computational Methods in Applied Sciences and Engineering
    Presentation's date: 2012-09-13
    Presentation of work at congresses

    Read the abstract Read the abstract Access to the full text Access to the full text Open in new window  Share Reference managers Reference managers Open in new window

    Solution of multiphase problems shows discontinuities across the material interfaces, which are usually weak. Using the eXtended Finite Element Method (X-FEM), these problems can be solved even for meshes that do not match the geometry. The basic idea is to enrich the interpolation space by means of a ridge function that is able to reproduce the discontinuity inside the elements. This approach yields excellent results for linear elements, but fails to be optimal if high-order interpolations are used. In this work, we propose a formulation that ensures optimal convergence rates for bimaterial problems. The key idea is to enrich the interpolation using a Heaviside function that allows the solution to represent polynomials on both sides of the interface and, provided the interface is accurately approximated, it yields optimal convergence rates. Although the interpolation is discontinuous, the desired continuity of the solution is imposed modifying the weak form. Moreover, in order to ensure optimal convergence, an accurate description of the interface (which also defines an integration rule for the elements cut by the interface) is needed. Here, we comment on different options that have been successfully used to integrate high-order X-FEM elements, and describe a general algorithm based on approximating the interface by piecewise polynomials of the same degree that the interpolation functions.

  • Desarrollo y análisis del métod PGD para procesos de conformado termoplástico

     Diez Mejia, Pedro; Parés Mariné, Núria; Modesto Galende, David; Zlotnik, Sergio; Verdugo Rojano, Francesc; Huerta Cerezuela, Antonio
    Participation in a competitive project

     Share

  • Desarrollo y análisis del método PGD para procesos de conformado termoplástico

     Huerta Cerezuela, Antonio; Diez Mejia, Pedro; Parés Mariné, Núria; Modesto Galende, David; Zlotnik, Sergio; Verdugo Rojano, Francesc
    Participation in a competitive project

     Share

  • Métodos Runge-Kutta implícitos de alto orden para flujos incompresibles

     de Villardi de Montlaur, Adeline; Fernandez Mendez, Sonia; Huerta Cerezuela, Antonio
    Revista internacional de métodos numéricos para cálculo y diseño en ingeniería
    Date of publication: 2011
    Journal article

    View View Open in new window  Share Reference managers Reference managers Open in new window

  • A simple shock-capturing technique for high-order discontinuous Galerkin methods

     Casoni Rero, Eva; Peraire Guitart, Jaume; Huerta Cerezuela, Antonio
    International journal for numerical methods in fluids
    Date of publication: 2011-11-11
    Journal article

    View View Open in new window  Share Reference managers Reference managers Open in new window

  • 3D NURBS-enhanced finite element method (NEFEM)

     Sevilla Cardenas, Ruben; Fernandez Mendez, Sonia; Huerta Cerezuela, Antonio
    International journal for numerical methods in engineering
    Date of publication: 2011-10-14
    Journal article

    Read the abstract Read the abstract View View Open in new window  Share Reference managers Reference managers Open in new window

    This paper presents the extension of the recently proposed NURBS-enhanced finite element method (NEFEM) to 3D domains. NEFEM is able to exactly represent the geometry of the computational domain by means of its CAD boundary representation with non-uniform rational B-splines (NURBS) surfaces. Specific strategies for interpolation and numerical integration are presented for those elements affected by the NURBS boundary representation. For elements not intersecting the boundary, a standard finite element rationale is used, preserving the efficiency of the classical FEM. In 3D NEFEM special attention must be paid to geometric issues that are easily treated in the 2D implementation. Several numerical examples show the performance and benefits of NEFEM compared with standard isoparametric or cartesian finite elements. NEFEM is a powerful strategy to efficiently treat curved boundaries and it avoids excessive mesh refinement to capture small geometric features.

  • Comparison of high-order curved finite elements

     Sevilla Cardenas, Ruben; Fernandez Mendez, Sonia; Huerta Cerezuela, Antonio
    International journal for numerical methods in engineering
    Date of publication: 2011-02-02
    Journal article

    Read the abstract Read the abstract View View Open in new window  Share Reference managers Reference managers Open in new window

    Several finite element techniques used in domains with curved boundaries are discussed and compared, with particular emphasis in two issues: the exact boundary representation of the domain and the consistency of the approximation. The influence of the number of integration points on the accuracy of the computation is also studied. Two-dimensional numerical examples, solved with continuous and discontinuous Galerkin formulations, are used to test and compare all these methodologies. In every example shown, the recently proposed NURBS-enhanced finite element method (NEFEM) provides the maximum accuracy for a given spatial discretization, at least one order of magnitude more accurate than classical isoparametric finite element method (FEM). Moreover, NEFEM outperforms Cartesian FEM and p-FEM, stressing the importance of the geometrical model as well as the relevance of a consistent approximation in finite element simulations.

  • Access to the full text
    NURBS-Enhanced Finite Element Method (NEFEM)  Open access

     Sevilla Cardenas, Ruben; Fernandez Mendez, Sonia; Huerta Cerezuela, Antonio
    Archives of computational methods in engineering
    Date of publication: 2011-11
    Journal article

    Read the abstract Read the abstract Access to the full text Access to the full text Open in new window  Share Reference managers Reference managers Open in new window

    The development of NURBS-Enhanced Finite Element Method (NEFEM) is revisited. This technique allows a seamless integration of the CAD boundary representation of the domain and the finite element method (FEM). The importance of the geometrical model in finite element simulations is addressed and the benefits and potential of NEFEM are discussed and compared with respect to other curved finite element techniques.

  • SHOCK CAPTURING FOR DISCONTINUOUS GALERKIN METHODS  Open access

     Casoni Rero, Eva
    Defense's date: 2011-10-14
    Department of Applied Mathematics III, Universitat Politècnica de Catalunya
    Theses

    Read the abstract Read the abstract Access to the full text Access to the full text Open in new window  Share Reference managers Reference managers Open in new window

    Aquesta tesi doctoral proposa formulacions de Galerkin Discontinu (DG) d’alt ordre per la captura de shocks, obtenint alhora solucions altament precises per problemes de flux compressible. En les últimes dècades, la investigació en els mètodes de DG ha estat en constant creixement. L'èxit dels mètodes DG en problemes hiperbòlics ha conduit el seu desenvolupament en lleis de conservació no lineals i problemes de convecció dominant. Entre els avantatges dels mètodes DG, destaquen la seva estabilitat inherent i les propietats locals de conservació. D'altra banda, els mètodes DG estan especialment dissenyats per l’ús aproximacions d'ordre superior. De fet, en els últims anys s'ha demostrat que la resolució de problemes de convecció dominant ja no es restringeix només a elements d'ordre inferior. De fet, es necessiten models numèrics d'alta precisió per aconseguir prediccions altament fiables dins la dinàmica de fluids computacional (CFD). En aquest context es presenten i discuteixen dos tècniques de captura de shocks. En primer lloc, es presenta una tècnica novedosa i senzilla basada en la introducció d'una nova base de funcions de forma. Aquesta base té la capacitat de canviar a nivell local entre una interpolació contínua o discontínua, depenent de la suavitat de la funció que es vol aproximar. En presència de xocs, les discontinuïtats introduïdes dins l’element permeten incloure l'estabilització necessària gràcies a l’ús dels fluxos numèrics, i alhora exploten les propietats intrínsiques del mètodes DG. En conseqüència, es poden utilitzar malles grolleres amb elements d’ordre superior. Amb aquestes discretitzacions i, utilitzant el mètode proposats, els xocs queden continguts a l’interior de l’element i per tant, és possible evitar l’ús de tècniques de refinament adaptatiu de la malla, alhora que es manté la localitat i compacitat dels esquemes DG. En segon lloc, es proposa una tècnica clàssica i, aparentment simple: la introducció de la viscositat artificial. Primerament es realitza un estudi detallat per al cas unidimensional. S’obté una viscositat d’alta precisió que escala segons el valor hk amb 1 ≤ k ≤ p i essent h la mida de l’element. En conseqüència, s’obté un xoc amb amplitud del mateix ordre. Seguidament, l'estudi de la viscositat unidimensional obtenida s'extén al cas multidimensional per a malles triangulars. L'extensió es basa en la projecció de la viscositat unidimensional en unes determinades direccions espacials dins l’element. Es demostra de manera consistent que la viscositat introduïda és, com a molt, del mateix ordre que la resolució donada per la discretització espacial, és a dir, h/p. El mètode és especialment eficient per aproximacions de Galerkin discontinu d’alt ordre, per exemple p≥ 3. Les dues metodologies es validen mitjançant una àmplia selecció d’exemples numèrics. En alguns exemples, els mètodes proposats permeten una reducció en el nombre de graus de llibertat necessaris per capturar xocs acuradament de fins i tot un ordre de magnitud, en comparació amb mètodes estàndar de refinament adaptatiu amb aproximacions de baix ordre.

    This thesis proposes shock-capturing methods for high-order Discontinuous Galerkin (DG) formulations providing highly accurate solutions for compressible flows. In the last decades, research in DG methods has been very active. The success of DG in hyperbolic problems has driven many studies for nonlinear conservation laws and convection-dominated problems. Among all the advantages of DG, their inherent stability and local conservation properties are relevant. Moreover, DG methods are naturally suited for high-order approximations. Actually, in recent years it has been shown that convection-dominated problems are no longer restricted to low-order elements. In fact, highly accurate numerical models for High-Fidelity predictions in CFD are necessary. Under this rationale, two shock-capturing techniques are presented and discussed. First, a novel and simple technique based on on the introduction of a new basis of shape functions is presented. It has the ability to change locally between a continuous or discontinuous interpolation depending on the smoothness of the approximated function. In the presence of shocks, the new discontinuities inside an element introduce the required stabilization thanks to the numerical fluxes, thus exploiting DG inherent properties. Large high-order elements can therefore be used and shocks are captured within a single element, avoiding adaptive mesh refinement and preserving the locality and compactness of the DG scheme. Second, a classical and, apparently simple, technique is advocated: the introduction of artificial viscosity. First, a one-dimensional study is perfomed. Viscosity of the order O(hk) with 1≤ k≤ p is obtained, hence inducing a shock width of the same order. Second, the study extends the accurate one-dimensional viscosity to triangular multidimensional meshes. The extension is based on the projection of the one-dimensional viscosity into some characteristic spatial directions within the elements. It is consistently shown that the introduced viscosity scales, at most, withthe DG resolutions length scales, h/p. The method is especially reliable for highorder DG approximations, say p≥3. A wide range of different numerical tests validate both methodologies. In some examples the proposed methods allow to reduce by an order of magnitude the number of degrees of freedom necessary to accurately capture the shocks, compared to standard low order h-adaptive approaches.

  • A proper generalized decomposition for the parametrized wave problem

     Modesto Galende, David; Zlotnik, Sergio; Huerta Cerezuela, Antonio
    International Conference on Adaptive Modeling and Simulation
    Presentation's date: 2011-06-07
    Presentation of work at congresses

    View View Open in new window  Share Reference managers Reference managers Open in new window

  • High-order discontinuous Galerkin methods for incompressible flows

     de Villardi de Montlaur, Adeline; Fernandez Mendez, Sonia; Huerta Cerezuela, Antonio
    International Conference on Finite Elements in Flow Problems
    Presentation's date: 2011-04-24
    Presentation of work at congresses

    View View Open in new window  Share Reference managers Reference managers Open in new window

  • Cálculos p-adaptativos de Galerkin discontinuo de alto orden

     Huerta Cerezuela, Antonio; Giorgiani, Giorgio; Modesto Galende, David; Fernandez Mendez, Sonia
    Congreso en Métodos Numéricos em Engenharia
    Presentation's date: 2011-06-24
    Presentation of work at congresses

    View View Open in new window  Share Reference managers Reference managers Open in new window

  • Novel shock-capturing schemes for discontinuous Galerkin methods

     Casoni Rero, Eva; Peraire Guitart, Jaume; Huerta Cerezuela, Antonio
    Numerical Methods for Hyperbolic Equations Theory and Applications
    Presentation's date: 2011-07-05
    Presentation of work at congresses

     Share Reference managers Reference managers Open in new window

  • High-order discontinuous Galerkin methods for compressible flow

     Casoni Rero, Eva; Huerta Cerezuela, Antonio; Peraire, Jaume
    International Conference on Finite Elements in Flow Problems
    Presentation's date: 2011-03-23
    Presentation of work at congresses

    View View Open in new window  Share Reference managers Reference managers Open in new window