995 research outputs found
A stabilized finite element method for the mixed wave equation in an ALE framework with application to diphthong production
The archived file is not the final published version of the article.
© (2016) S. Hirzel Verlag/European Acoustics Association
The definitive publisher-authenticated version is available online at http://www.ingentaconnect.com/contentone/dav/aaua/2016/00000102/00000001/art00012
Readers must contact the publisher for reprint or permission to use the material in any form.Working with the wave equation in mixed rather than irreducible form allows one to directly account for both, the acoustic pressure field and the acoustic particle velocity field. Indeed, this becomes the natural option in many problems, such as those involving waves propagating in moving domains, because the equations can easily be set in an arbitrary Lagrangian-Eulerian (ALE) frame of reference. Yet, when attempting a standard Galerkin finite element solution (FEM) for them, it turns out that an inf-sup compatibility constraint has to be satisfied, which prevents from using equal interpolations for the approximated acoustic pressure and velocity fields. In this work it is proposed to resort to a subgrid scale stabilization strategy to circumvent this condition and thus facilitate code implementation. As a possible application, we address the generation of diphthongs in voice production.Peer ReviewedPostprint (author's final draft
Residual-based stabilization of the finite element approximation to the acoustic perturbation equations for low Mach number aeroacoustics
This is the peer reviewed version of the following article: [Guasch, O., Sánchez-MartÃn, P., Pont, A., Baiges, J., and Codina, R. (2016) Residual-based stabilization of the finite element approximation to the acoustic perturbation equations for low Mach number aeroacoustics. Int. J. Numer. Meth. Fluids, 82: 839–857. doi: 10.1002/fld.4243], which has been published in final form at http://onlinelibrary.wiley.com/doi/10.1002/fld.4243/abstract. This article may be used for non-commercial purposes in accordance with Wiley Terms and Conditions for Self-Archiving.The acoustic perturbation equations (APE) are suitable to predict aerodynamic noise in the presence of a non-uniform mean flow. As for any hybrid computational aeroacoustics approach, a first computational fluid dynamics simulation is carried out from which the mean flow characteristics and acoustic sources are obtained. In a second step, the APE are solved to get the acoustic pressure and particle velocity fields. However, resorting to the finite element method (FEM) for that purpose is not straightforward. Whereas mixed finite elements satisfying an appropriate inf–sup compatibility condition can be built in the case of no mean flow, that is, for the standard wave equation in mixed form, these are difficult to implement and their good performance is yet to be checked for more complex wave operators. As a consequence, strong simplifying assumptions are usually considered when solving the APE with FEM. It is possible to avoid them by resorting to stabilized formulations. In this work, a residual-based stabilized FEM is presented for the APE at low Mach numbers, which allows one to deal with the APE convective and reaction terms in its full extent. The key of the approach resides in the design of the matrix of stabilization parameters. The performance of the formulation and the contributions of the different terms in the equations are tested for an acoustic pulse propagating in sheared-solenoidal mean flow, and for the aeolian tone generated by flow past a two-dimensional cylinder.Peer ReviewedPostprint (author's final draft
Adaptive multiscale model reduction with Generalized Multiscale Finite Element Methods
In this paper, we discuss a general multiscale model reduction framework
based on multiscale finite element methods. We give a brief overview of related
multiscale methods. Due to page limitations, the overview focuses on a few
related methods and is not intended to be comprehensive. We present a general
adaptive multiscale model reduction framework, the Generalized Multiscale
Finite Element Method. Besides the method's basic outline, we discuss some
important ingredients needed for the method's success. We also discuss several
applications. The proposed method allows performing local model reduction in
the presence of high contrast and no scale separation
Nonlinear nonlocal multicontinua upscaling framework and its applications
In this paper, we discuss multiscale methods for nonlinear problems. The main
idea of these approaches is to use local constraints and solve problems in
oversampled regions for constructing macroscopic equations. These techniques
are intended for problems without scale separation and high contrast, which
often occur in applications. For linear problems, the local solutions with
constraints are used as basis functions. This technique is called Constraint
Energy Minimizing Generalized Multiscale Finite Element Method (CEM-GMsFEM).
GMsFEM identifies macroscopic quantities based on rigorous analysis. In
corresponding upscaling methods, the multiscale basis functions are selected
such that the degrees of freedom have physical meanings, such as averages of
the solution on each continuum.
This paper extends the linear concepts to nonlinear problems, where the local
problems are nonlinear. The main concept consists of: (1) identifying
macroscopic quantities; (2) constructing appropriate oversampled local problems
with coarse-grid constraints; (3) formulating macroscopic equations. We
consider two types of approaches. In the first approach, the solutions of local
problems are used as basis functions (in a linear fashion) to solve nonlinear
problems. This approach is simple to implement; however, it lacks the nonlinear
interpolation, which we present in our second approach. In this approach, the
local solutions are used as a nonlinear forward map from local averages
(constraints) of the solution in oversampling region. This local fine-grid
solution is further used to formulate the coarse-grid problem. Both approaches
are discussed on several examples and applied to single-phase and two-phase
flow problems, which are challenging because of convection-dominated nature of
the concentration equation
Self-Evaluation Applied Mathematics 2003-2008 University of Twente
This report contains the self-study for the research assessment of the Department of Applied Mathematics (AM) of the Faculty of Electrical Engineering, Mathematics and Computer Science (EEMCS) at the University of Twente (UT). The report provides the information for the Research Assessment Committee for Applied Mathematics, dealing with mathematical sciences at the three universities of technology in the Netherlands. It describes the state of affairs pertaining to the period 1 January 2003 to 31 December 2008
Variational multiscale stabilization of finite and spectral elements for dry and moist atmospheric problems
In this thesis the finite and spectral element methods (FEM and SEM, respectively) applied to
problems in atmospheric simulations are explored through the common thread of Variational
Multiscale Stabilization (VMS). This effort is justified by three main reasons. (i) the recognized
need for new solvers that can efficiently execute on massively parallel architectures ¿a spreading
framework in most fields of computational physics in which numerical weather prediction
(NWP) occupies a prominent position. Element-based methods (e.g. FEM, SEM, discontinuous
Galerkin) have important advantages in parallel code development; (ii) the inherent flexibility of
these methods with respect to the geometry of the grid makes them a great candidate for dynamically
adaptive atmospheric codes; and (iii) the localized diffusion provided by VMS represents
an improvement in the accurate solution of multi-physics problems where artificial diffusion may
fail. Its application to atmospheric simulations is a novel approach within a field of research
that is still open. First, FEM and VMS are described and derived for the solution of stratified
low Mach number flows in the context of dry atmospheric dynamics. The validity of the method
to simulate stratified flows is assessed using standard two- and three-dimensional benchmarks
accepted by NWP practitioners. The problems include thermal and gravity driven simulations.
It will be shown that stability is retained in the regimes of interest and a numerical comparison
against results from the the literature will be discussed. Second, the ability of VMS to stabilize
the FEM solution of advection-dominated problems (i.e. Euler and transport equations) is taken
further by the implementation of VMS as a stabilizing tool for high-order spectral elements with
advection-diffusion problems. To the author¿s knowledge, this is an original contribution to the
literature of high order spectral elements involved with transport in the atmosphere. The problem
of monotonicity-preserving high order methods is addressed by combining VMS-stabilized
SEM with a discontinuity capturing technique. This is an alternative to classical filters to treat
the Gibbs oscillations that characterize high-order schemes. To conclude, a microphysics scheme
is implemented within the finite element Euler solver, as a first step toward realistic atmospheric
simulations. Kessler microphysics is used to simulate the formation of warm, precipitating clouds.
This last part combines the solution of the Euler equations for stratified flows with the solution
of a system of transport equations for three classes of water: water vapor, cloud water, and rain.
The method is verified using idealized two- and three-dimensional storm simulations.En esta tesis los métodos de elementos finitos y espectrales (FEM - finite element method y SEM- spectral element method, respectivamente), aplicados a los problemas de simulaciones atmosféricas, se exploran a través del método de estabilización conocidocomo Variational Multiscale Stabilization (VMS). Tres razones fundamentales justifican este esfuerzo: (i) la necesidad de tener nuevos métodos de solución de las ecuaciones diferenciales a las derivadas parciales usando máquinas paralelas de gran escala –un entorno en expansión en muchos campos de la mecánica computacional, dentro de la cual la predicción numérica de la dinámica atmosférica (NWP-numerical weather prediction)representa una aplicación importante. Métodos del tipo basado en elementos(por ejemplo, FEM, SEM, Galerkin discontinuo) presentan grandes ventajas en el desarrollo de códigos paralelos; (ii) la flexibilidad intrÃnseca de tales métodos respecto a lageometrÃa de la malla computacional hace que esos métodos sean los candidatos ideales para códigos atmosféricos con mallas adaptativas; y (iii) la difusión localizada que VMSintroduce representa una mejora en las soluciones de problemas con fÃsica compleja en los cuales la difusión artificial clásica no funcionarÃa. La aplicación de FEM o SEM con VMS a problemas de simulaciones atmosféricas es una estrategia innovadora en un campo de investigación abierto. En primera instancia, FEM y VMS vienen descritos y derivados para la solución de flujos estratificados a bajo número de Mach en el contexto de la dinámica atmosférica. La validez del método para simular flujos estratificados es verificada por medio de test estándar aceptado por la comunidad dentro del campo deNWP. Los test incluyen simulaciones de flujos térmicos con efectos de gravedad. Se demostrará que la estabilidad del método numérico se preserva dentro de los regÃmenesde interés y se discutirá una comparación numérica de los resultados frente a aquellos hallados en la literatura. En segunda instancia, la capacidad de VMS para estabilizarmétodos FEM en problemas de advección dominante (i.e. ecuaciones de Euler y ecuaciones de transporte) se implementa además en la solución a elementos espectrales de alto orden en problemas de advección-difusión. Hasta donde el autor sabe, esta es una contribución original a la literatura de métodos basados en elementos espectrales en problemas de transporte atmosférico. El problema de monotonicidad con métodos de alto orden es tratado mediante la combinación de SEM+VMS con una técnica de shockcapturing para un mejor tratamiento de las discontinuidades. Esta es una alternativa a los filtros que normalmente se aplican a SEM para eilminar las oscilaciones de Gibbsque caracterizan las soluciones de alto orden. Como último punto, se implementa un esquema de humedad acoplado con el núcleo en elementos finitos; este es un primer paso hacia simulaciones atmosféricas más realistas. La microfÃsica de Kessler se emplea para simular la formación de nubes y tormentas cálidas (warm clouds: no permite la formación de hielo). Esta última parte combina la solución de las ecuaciones de Eulerpara atmósferas estratificadas con la solución de un sistema de ecuaciones de transporte de tres estados de agua: vapor, nubes y lluvia. La calidad del método es verificadautilizando simulaciones de tormenta en dos y tres dimensiones
A Review of Element-Based Galerkin Methods for Numerical Weather Prediction: Finite Elements, Spectral Elements, and Discontinuous Galerkin
Numerical weather prediction (NWP) is in a period of transition. As resolutions increase, global models are moving towards fully nonhydrostatic dynamical cores, with the local and global models using the same governing equations; therefore we have reached a point where it will be necessary to use a single model for both applications. The new dynamical cores at the heart of these unified models are designed to scale efficiently on clusters with hundreds of thousands or even millions of CPU cores and GPUs. Operational and research NWP codes currently use a wide range of numerical methods: finite differences, spectral transform, finite volumes and, increasingly, finite/spectral elements and discontinuous Galerkin, which constitute element-based Galerkin (EBG) methods.Due to their important role in this transition, will EBGs be the dominant power behind NWP in the next 10 years, or will they just be one of many methods to choose from? One decade after the review of numerical methods for atmospheric modeling by Steppeler et al. (Meteorol Atmos Phys 82:287–301, 2003), this review discusses EBG methods as a viable numerical approach for the next-generation NWP models. One well-known weakness of EBG methods is the generation of unphysical oscillations in advection-dominated flows; special attention is hence devoted to dissipation-based stabilization methods. Since EBGs are geometrically flexible and allow both conforming and non-conforming meshes, as well as grid adaptivity, this review is concluded with a short overview of how mesh generation and dynamic mesh refinement are becoming as important for atmospheric modeling as they have been for engineering applications for many years.The authors would like to thank Prof. Eugenio Oñate (U. Politècnica de Catalunya) for his invitation to submit this review article. They are also thankful to Prof. Dale Durran (U. Washington), Dr. Tommaso Benacchio (Met Office), and Dr. Matias Avila (BSC-CNS) for their comments and corrections, as well as
insightful discussion with Sam Watson, Consulting Software Engineer (Exa Corp.) Most of the contribution to this article by the first author stems from his Ph.D. thesis carried out at the Barcelona Supercomputing Center (BSCCNS) and Universitat Politècnica de Catalunya, Spain, supported by a BSC-CNS student grant, by Iberdrola EnergÃas Renovables, and by grant N62909-09-1-4083 of the Office of Naval Research Global. At NPS, SM, AM, MK, and FXG were supported by the Office of Naval Research through program element PE-0602435N, the Air Force Office of Scientific Research through the Computational Mathematics program, and the National Science Foundation (Division of Mathematical Sciences) through program element 121670. The scalability studies of the
atmospheric model NUMA that are presented in this paper used resources of the Argonne Leadership Computing Facility, which is a DOE Office of Science User Facility supported under Contract DE-AC02-06CH11357. SM, MK, and AM are grateful to the National Research Council of the National Academies.Peer ReviewedPostprint (author's final draft
Computational Engineering
The focus of this Computational Engineering Workshop was on the mathematical foundation of state-of-the-art and emerging finite element methods in engineering analysis. The 52 participants included mathematicians and engineers with shared interest on discontinuous Galerkin or Petrov-Galerkin methods and other generalized nonconforming or mixed finite element methods
- …