122 research outputs found
Acoustic modelling of exhaust devices with nonconforming finite element meshes and transfer matrices
[EN] Transfer matrices are commonly considered in the numerical modelling of the acoustic behaviour associated with exhaust devices in the breathing system of internal combustion engines, such as catalytic converters, particulate filters, perforated mufflers and charge air coolers. In a multidimensional finite element approach, a transfer matrix provides a relationship between the acoustic fields of the nodes located at both sides of a particular region. This approach can be useful, for example, when one-dimensional propagation takes place within the region substituted by the transfer matrix. As shown in recent investigations, the sound attenuation of catalytic converters can be properly predicted if the monolith is replaced by a plane wave four-pole matrix. The finite element discretization is retained for the inlet/outlet and tapered ducts, where multidimensional acoustic fields can exist. In this case, only plane waves are present within the capillary ducts, and three-dimensional propagation is possible in the rest of the catalyst subcomponents. Also, in the acoustic modelling of perforated mufflers using the finite element method, the central passage can be replaced by a transfer matrix relating the pressure difference between both sides of the perforated surface with the acoustic velocity through the perforations. The approaches in the literature that accommodate transfer matrices and finite element models consider conforming meshes at connecting interfaces, therefore leading to a straightforward evaluation of the coupling integrals. With a view to gaining flexibility during the mesh generation process, it is worth developing a more general procedure. This has to be valid for the connection of acoustic subdomains by transfer matrices when the discretizations are nonconforming at the connecting interfaces. In this work, an integration algorithm similar to those considered in the mortar finite element method, is implemented for nonmatching grids in combination with acoustic transfer matrices. A number of numerical test problems related to some relevant exhaust devices are then presented to assess the accuracy and convergence performance of the proposed procedure.Authors gratefully acknowledge the financial support of Ministerio de Ciencia e Innovacion and the European Regional Development Fund by means of the Projects DPI2007-62635 and DPI2010-15412.Denia, F.; MartĂnez-Casas, J.; Baeza, L.; Fuenmayor, F. (2012). Acoustic modelling of exhaust devices with nonconforming finite element meshes and transfer matrices. Applied Acoustics. 73(8):713-722. https://doi.org/10.1016/j.apacoust.2012.02.003S71372273
Space-time domain decomposition for advection-diffusion problems in mixed formulations
This paper is concerned with the numerical solution of porous-media flow and
transport problems , i. e. heterogeneous, advection-diffusion problems. Its aim
is to investigate numerical schemes for these problems in which different time
steps can be used in different parts of the domain. Global-in-time,
non-overlapping domain-decomposition methods are coupled with operator
splitting making possible the different treatment of the advection and
diffusion terms. Two domain-decomposition methods are considered: one uses the
time-dependent Steklov--Poincar{\'e} operator and the other uses optimized
Schwarz waveform relaxation (OSWR) based on Robin transmission conditions. For
each method, a mixed formulation of an interface problem on the space-time
interface is derived, and different time grids are employed to adapt to
different time scales in the subdomains. A generalized Neumann-Neumann
preconditioner is proposed for the first method. To illustrate the two methods
numerical results for two-dimensional problems with strong heterogeneities are
presented. These include both academic problems and more realistic prototypes
for simulations for the underground storage of nuclear waste
Infrastructure for the Coupling of Dune Grids
We describe an abstract interface for the geometric coupling of finite element grids. The scope of the interface encompasses a wide range of domain decomposition techniques in use today, including nonconforming grids and grids of different dimensions. The couplings are described as sets of remote intersections, which encapsulate the relationships between pairs of elements on the coupling interface.
The abstract interface is realized in a module dune-grid-glue for the software framework dune. Several implementations of this interface exist, including one for general nonconforming couplings and a special efficient implementation for conforming interfaces. We present two numerical examples to show the flexibility of the approach
Robin Schwarz algorithm for the NICEM Method: the Pq finite element case
In Gander et al. [2004] we proposed a new non-conforming domain decomposition
paradigm, the New Interface Cement Equilibrated Mortar (NICEM) method, based on
Schwarz type methods that allows for the use of Robin interface conditions on
non-conforming grids. The error analysis was done for P1 finite elements, in 2D
and 3D. In this paper, we provide new numerical analysis results that allow to
extend this error analysis in 2D for piecewise polynomials of higher order and
also prove the convergence of the iterative algorithm in all these cases.Comment: arXiv admin note: substantial text overlap with arXiv:0705.028
A partition of unity approach to fluid mechanics and fluid-structure interaction
For problems involving large deformations of thin structures, simulating
fluid-structure interaction (FSI) remains challenging largely due to the need
to balance computational feasibility, efficiency, and solution accuracy.
Overlapping domain techniques have been introduced as a way to combine the
fluid-solid mesh conformity, seen in moving-mesh methods, without the need for
mesh smoothing or re-meshing, which is a core characteristic of fixed mesh
approaches. In this work, we introduce a novel overlapping domain method based
on a partition of unity approach. Unified function spaces are defined as a
weighted sum of fields given on two overlapping meshes. The method is shown to
achieve optimal convergence rates and to be stable for steady-state Stokes,
Navier-Stokes, and ALE Navier-Stokes problems. Finally, we present results for
FSI in the case of a 2D mock aortic valve simulation. These initial results
point to the potential applicability of the method to a wide range of FSI
applications, enabling boundary layer refinement and large deformations without
the need for re-meshing or user-defined stabilization.Comment: 34 pages, 15 figur
- …