Search results for "Discontinuous Galerkin"
showing 10 items of 27 documents
High-fidelity analysis of multilayered shells with cut-outs via the discontinuous Galerkin method
2021
Abstract A novel numerical method for the analysis of multilayered shells with cut-outs is presented. In the proposed approach, the shell geometry is represented via either analytical functions or NURBS parametrizations , while generally-shaped cut-outs are defined implicitly within the shell modelling domain via a level set function . The multilayered shell problem is addressed via the Equivalent-Single-Layer approach whereby high-order polynomial functions are employed to approximate the covariant components of the displacement field throughout the shell thickness. The shell governing equations are then derived from the Principle of Virtual Displacements of three-dimensional elasticity an…
Analysis of blood flow in one dimensional elastic artery using Navier-Stokes conservation laws
2017
En los últimos años, la simulación computacional en ámbitos médicos ha aumentado notablemente en múltiples ramas de la ciencia, desde modelización a métodos numéricos, pasando por informática. Los principales objetivos de esta disciplina incipiente son comprobar hipótesis antes de una intervención, o ver qué efecto podría tener un medicamento antes de tomarlo, entre otros. En este trabajo deduciremos desde los principios físicos más básicos un modelo unidimensional para la simulación del flujo sanguíneo en arterias elásticas. Proporcionaremos un marco histórico, así como una revisión de este tipo de modelos. Estudiaremos también desde el punto de vista del análisis matemático las ecuaciones…
Modeling wave propagation in elastic solids via high-order accurate implicit-mesh discontinuous Galerkin methods
2022
A high-order accurate implicit-mesh discontinuous Galerkin framework for wave propagation in single-phase and bi-phase solids is presented. The framework belongs to the embedded-boundary techniques and its novelty regards the spatial discretization, which enables boundary and interface conditions to be enforced with high-order accuracy on curved embedded geometries. High-order accuracy is achieved via high-order quadrature rules for implicitly-defined domains and boundaries, whilst a cell-merging strategy addresses the presence of small cut cells. The framework is used to discretize the governing equations of elastodynamics, written using a first-order hyperbolic momentum-strain formulation…
Shear-Thinning in Oligomer Melts—Molecular Origins and Applications
2021
We investigate the molecular origin of shear-thinning in melts of flexible, semiflexible and rigid oligomers with coarse-grained simulations of a sheared melt. Entanglements, alignment, stretching and tumbling modes or suppression of the latter all contribute to understanding how macroscopic flow properties emerge from the molecular level. In particular, we identify the rise and decline of entanglements with increasing chain stiffness as the major cause for the non-monotonic behaviour of the viscosity in equilibrium and at low shear rates, even for rather small oligomeric systems. At higher shear rates, chains align and disentangle, contributing to shear-thinning. By performing simulations …
Comparison of continuous and discontinuous Galerkin approaches for variable-viscosity Stokes flow
2015
We describe a Discontinuous Galerkin (DG) scheme for variable-viscosity Stokes flow which is a crucial aspect of many geophysical modelling applications and conduct numerical experiments with different elements comparing the DG approach to the standard Finite Element Method (FEM). We compare the divergence-conforming lowest-order Raviart-Thomas (RT0P0) and Brezzi-Douglas-Marini (BDM1P0) element in the DG scheme with the bilinear Q1P0 and biquadratic Q2P1 elements for velocity and their matching piecewise constant/linear elements for pressure in the standard continuous Galerkin (CG) scheme with respect to accuracy and memory usage in 2D benchmark setups. We find that for the chosen geodynami…
Discrete Maximum Principle for Galerkin Finite Element Solutions to Parabolic Problems on Rectangular Meshes
2004
One of the most important problems in numerical simulation is the preservation of qualitative properties of solutions of mathematical models. For problems of parabolic type, one of such properties is the maximum principle. In [5], Fujii analyzed the discrete analogue of the (continuous) maximum principle for the linear parabolic problems, and derived sufficient conditions guaranteeing its validity for the Galerkin finite element approximations built on simplicial meshes. In our paper, we present the sufficient conditions for the validity of the discrete maximum principle for the case of bilinear finite element space approximations on rectangular meshes.
A discontinuous Galerkin formulation for nonlinear analysis of multilayered shells refined theories
2023
A novel pure penalty discontinuous Galerkin method is proposed for the geometrically nonlinear analysis of multilayered composite plates and shells, modelled via high-order refined theories. The approach allows to build different two-dimensional equivalent single layer structural models, which are obtained by expressing the covariant components of the displacement field through-the-thickness via Taylor’s polynomial expansion of different order. The problem governing equations are deduced starting from the geometrically nonlinear principle of virtual displacements in a total Lagrangian formulation. They are addressed with a pure penalty discontinuous Galerkin method using Legendre polynomial…
Comparison between adaptive and uniform discontinuous Galerkin simulations in dry 2D bubble experiments
2013
Accepted by the Journal of Computational Physics Adaptive mesh refinement generally aims to increase computational efficiency without compromising the accuracy of the numerical solution. However it is an open question in which regions the spatial resolution can actually be coarsened without affecting the accuracy of the result. This question is investigated for a specific example of dry atmospheric convection, namely the simulation of warm air bubbles. For this purpose a novel numerical model is developed that is tailored towards this specific application. The compressible Euler equations are solved with a Discontinuous Galerkin method. Time integration is done with an IMEXmethod and the dy…
Equivalent-Single-Layer discontinuous Galerkin methods for static analysis of multilayered shells
2021
Abstract An original formulation for the elastic analysis of multilayered shells is presented in this work. The key features of the formulation are: the representation of the shell mean surface via a generic system of curvilinear coordinates; the unified treatment of general shell theories via an Equivalent-Single-Layer approach based on the through-the-thickness expansion of the covariant components of the displacement field; and an Interior Penalty discontinuous Galerkin scheme for the solution of the set of governing equations. The combined use of these features enables a high-order solution of the multilayered shell problem. Several numerical tests are presented for isotropic, orthotrop…
High-order simulation scheme for active particles driven by stress boundary conditions
2020
Abstract We study the dynamics and interactions of elliptic active particles in a two dimensional solvent. The particles are self-propelled through prescribing a fluid stress at one half of the fluid-particle boundary. The fluid is treated explicitly solving the Stokes equation through a discontinuous Galerkin scheme, which allows to simulate strictly incompressible fluids. We present numerical results for a single particle and give an outlook on how to treat suspensions of interacting active particles.