635 research outputs found
Generalized multiscale finite element method for a strain-limiting nonlinear elasticity model
In this paper, we consider multiscale methods for nonlinear elasticity. In
particular, we investigate the Generalized Multiscale Finite Element Method
(GMsFEM) for a strain-limiting elasticity problem. Being a special case of the
naturally implicit constitutive theory of nonlinear elasticity, strain-limiting
relation has presented an interesting class of material bodies, for which
strains remain bounded (even infinitesimal) while stresses can become
arbitrarily large. The nonlinearity and material heterogeneities can create
multiscale features in the solution, and multiscale methods are therefore
necessary. To handle the resulting nonlinear monotone quasilinear elliptic
equation, we use linearization based on the Picard iteration. We consider two
types of basis functions, offline and online basis functions, following the
general framework of GMsFEM. The offline basis functions depend nonlinearly on
the solution. Thus, we design an indicator function and we will recompute the
offline basis functions when the indicator function predicts that the material
property has significant change during the iterations. On the other hand, we
will use the residual based online basis functions to reduce the error
substantially when updating basis functions is necessary. Our numerical results
show that the above combination of offline and online basis functions is able
to give accurate solutions with only a few basis functions per each coarse
region and updating basis functions in selected iterations.Comment: 19 pages, 2 figures, submitted to Journal of Computational and
Applied Mathematic
Higher-order discontinuous modelling of fracturing in quasi-brittle materials
Quasi-brittle failure is characterised by material degradation, fracturing and potential interaction of fragmented parts. The computational description of this behaviour has presented significant challenges to the mechanics community over the past few decades, driven by the development of technology, the increasing social and economical constraints for safer and more complicated engineering designs and consequently by the increasing requirements for more accurate understanding of macro- and micro-structural processes.
Finite element methods have been pushed to their limits in an attempt to resolve strain localisation and ultimately fracturing in a unified and objective manner, while discrete methods have been utilised by artificial connection of discrete bodies which are identified a priori to act as continua. Neither of these attempts comprises a diritta via for modelling the transition from continuum to discontinuum efficiently and this has led to the investigation of alternative techniques.
Herein, the numerical modelling of quasi-brittle localisation and fracturing is investigated using the Numerical Manifold Method (NMM) as an alternative unifying framework to industry-established techniques such as the Finite Element Method (FEM) and Discrete Element Method (DEM). One of the particularly interesting aspects of NMM is with respect to its potential for modelling both continuum and discontinuum states and providing an efficient framework for modelling the entire transition between continuum to discontinuum, from a continuum point of view, without remeshing. The attractive nature of this capability advocates potential for modelling mechanics of materials such as concrete, rock and masonry, but also a more general class of quasi-brittle materials.
This work investigates and extends NMM primarily with respect to the following characteristics:
1. Discontinuities, such as cracks, are introduced naturally in a discrete manner, but in a continuum setting, without the need for remeshing
2. The approximation is improved globally or locally, for any arbitrary level, without remeshing
3. Integration is undertaken explicitly, for any arbitrary level of local improvement of the approximation
Furthermore, NMM is reformulated using a constrained variational approach for generalised three-dimensional problems. Essential boundary conditions are enforced using Lagrange multipliers and projection matrices and potential higher-order boundary issues are investigated. The developments are implemented algorithmically in MATLAB and higher-order enrichment is demonstrated with the use of adaptivity
Besov regularity for operator equations on patchwise smooth manifolds
We study regularity properties of solutions to operator equations on
patchwise smooth manifolds such as, e.g., boundaries of
polyhedral domains . Using suitable biorthogonal
wavelet bases , we introduce a new class of Besov-type spaces
of functions
. Special attention is paid on the
rate of convergence for best -term wavelet approximation to functions in
these scales since this determines the performance of adaptive numerical
schemes. We show embeddings of (weighted) Sobolev spaces on
into , ,
which lead us to regularity assertions for the equations under consideration.
Finally, we apply our results to a boundary integral equation of the second
kind which arises from the double layer ansatz for Dirichlet problems for
Laplace's equation in .Comment: 42 pages, 3 figures, updated after peer review. Preprint: Bericht
Mathematik Nr. 2013-03 des Fachbereichs Mathematik und Informatik,
Universit\"at Marburg. To appear in J. Found. Comput. Mat
A fast immersed boundary method for external incompressible viscous flows using lattice Green's functions
A new parallel, computationally efficient immersed boundary method for
solving three-dimensional, viscous, incompressible flows on unbounded domains
is presented. Immersed surfaces with prescribed motions are generated using the
interpolation and regularization operators obtained from the discrete delta
function approach of the original (Peskin's) immersed boundary method. Unlike
Peskin's method, boundary forces are regarded as Lagrange multipliers that are
used to satisfy the no-slip condition. The incompressible Navier-Stokes
equations are discretized on an unbounded staggered Cartesian grid and are
solved in a finite number of operations using lattice Green's function
techniques. These techniques are used to automatically enforce the natural
free-space boundary conditions and to implement a novel block-wise adaptive
grid that significantly reduces the run-time cost of solutions by limiting
operations to grid cells in the immediate vicinity and near-wake region of the
immersed surface. These techniques also enable the construction of practical
discrete viscous integrating factors that are used in combination with
specialized half-explicit Runge-Kutta schemes to accurately and efficiently
solve the differential algebraic equations describing the discrete momentum
equation, incompressibility constraint, and no-slip constraint. Linear systems
of equations resulting from the time integration scheme are efficiently solved
using an approximation-free nested projection technique. The algebraic
properties of the discrete operators are used to reduce projection steps to
simple discrete elliptic problems, e.g. discrete Poisson problems, that are
compatible with recent parallel fast multipole methods for difference
equations. Numerical experiments on low-aspect-ratio flat plates and spheres at
Reynolds numbers up to 3,700 are used to verify the accuracy and physical
fidelity of the formulation.Comment: 32 pages, 9 figures; preprint submitted to Journal of Computational
Physic
Adaptive Generalized Multiscale Model Reduction Techniques for Problems in Perforated Domains
Multiscale modeling of complex physical phenomena in many areas, including hydrogeology,
material science, chemistry and biology, consists of solving problems in highly
heterogeneous porous media. In many of these applications, differential equations are formulated
in perforated domains which can be considered as the region outside of inclusions
or connected bodies of various sizes. Due to complicated geometries of these inclusions,
solutions to these problems have multiscale features. Taking into account the uncertainties,
one needs to solve these problems extensively many times. Model reduction techniques
are significant for problems in perforated domains in order to improve the computational
efficiency.
There are some existing approaches for model reduction in perforated domains including
homogenization, heterogeneous multiscale methods and multiscale finite element
methods. These techniques typically consider the case when there is a scale separation or
the perforation distribution is periodic, and assume that the solution space can be approximated
by the solutions of directional cell problems and the effective equations contain a
limited number of effective parameters.
For more complicated problems where the effective properties may be richer, we are
interested in developing systematic local multiscale model reduction techniques to obtain
accurate macroscale representations of the underlying fine-scale problem in highly heterogeneous
perforated domains. In this dissertation, based on the framework of Generalized
Multiscale Finite Element Method, we develop novel methods and algorithms including
(1) development of systematic local model reduction techniques for computing multiscale
basis in perforated domains, (2) numerical analysis and exhaustive simulation utilizing the
proposed basis functions, (3) design of different applicable global coupling frameworks
and (4) applications to various problems with challenging engineering backgrounds. Our
proposed methods can significantly advance the computational efficiency and accuracy for
multiscale problems in perforated media
Development and applications of the Finite Point Method to compressible aerodynamics problems
This work deals with the development and application of the Finite Point Method (FPM) to compressible aerodynamics problems. The research focuses mainly on investigating the capabilities of the meshless technique to address practical problems, one of the most outstanding issues in meshless methods.
The FPM spatial approximation is studied firstly, with emphasis on aspects of the methodology that can be improved to increase its robustness and accuracy. Suitable ranges for setting the relevant approximation parameters and the performance likely to be attained in practice are determined. An automatic procedure to adjust the approximation parameters is also proposed to simplify the application of the method, reducing problem- and user-dependence without affecting the flexibility of the meshless technique.
The discretization of the flow equations is carried out following wellestablished approaches, but drawing on the meshless character of the
methodology. In order to meet the requirements of practical applications, the procedures are designed and implemented placing emphasis on robustness and efficiency (a simplification of the basic FPM technique is proposed to this end). The flow solver is based on an upwind spatial discretization of the convective fluxes (using the approximate Riemann solver of Roe) and an explicit time integration scheme. Two additional artificial diffusion schemes are also proposed to suit those cases of study in which computational cost is a major concern. The performance of the flow solver is evaluated in order to determine the potential of the meshless approach. The accuracy, computational cost and parallel scalability of the method are studied in comparison with a conventional FEM-based technique.
Finally, practical applications and extensions of the flow solution scheme are presented. The examples provided are intended not only to show the
capabilities of the FPM, but also to exploit meshless advantages. Automatic hadaptive procedures, moving domain and fluid-structure interaction problems, as well as a preliminary approach to solve high-Reynolds viscous flows, are a sample of the topics explored.
All in all, the results obtained are satisfactorily accurate and competitive in terms of computational cost (if compared with a similar mesh-based
implementation). This indicates that meshless advantages can be exploited with efficiency and constitutes a good starting point towards more challenging applications.En este trabajo se aborda el desarrollo del Método de Puntos Finitos (MPF) y su aplicación a problemas de aerodinámica de flujos compresibles. El objetivo principal es investigar el potencial de la técnica sin malla para la solución de problemas prácticos, lo cual constituye una de las limitaciones más importantes de los métodos sin malla.
En primer lugar se estudia la aproximación espacial en el MPF, haciendo hincapié en aquéllos aspectos que pueden ser mejorados para incrementar la robustez y exactitud de la metodología. Se determinan rangos adecuados para el ajuste de los parámetros de la aproximación y su comportamiento en situaciones prácticas. Se propone además un procedimiento de ajuste automático de estos parámetros a fin de simplificar la aplicación del método y reducir la dependencia de factores como el tipo de problema y la intervención del usuario, sin afectar la flexibilidad de la técnica sin malla.
A continuación se aborda el esquema de solución de las ecuaciones del flujo. La discretización de las mismas se lleva a cabo siguiendo métodos estándar, pero aprovechando las características de la técnica sin malla. Con el objetivo de abordar problemas prácticos, se pone énfasis en la robustez y eficiencia de la implementación numérica (se propone además una simplificación del procedimiento de solución). El comportamiento del esquema se estudia en detalle para evaluar su potencial y se analiza su exactitud, coste computacional y escalabilidad, todo ello en comparación con un método convencional basado en Elementos Finitos.
Finalmente se presentan distintas aplicaciones y extensiones de la metodología desarrollada. Los ejemplos numéricos pretenden demostrar las
capacidades del método y también aprovechar las ventajas de la metodología sin malla en áreas en que la misma puede ser de especial interés. Los problemas tratados incluyen, entre otras características, el refinamiento automático de la discretización, la presencia de fronteras móviles e
interacción fluido-estructura, como así también una aplicación preliminar a flujos compresibles de alto número de Reynolds. Los resultados obtenidos muestran una exactitud satisfactoria. Además, en comparación con una técnica similar basada en Elementos Finitos, demuestran ser competitivos en términos del coste computacional. Esto indica que las ventajas de la metodología sin malla pueden ser explotadas con eficiencia, lo cual constituye un buen punto de partida para el desarrollo de ulteriores aplicaciones.Postprint (published version
Development and applications of the finite point method to compressible aerodynamics problems
This work deals with the development and application of the Finite Point
Method (FPM) to compressible aerodynamics problems. The research focuses
mainly on investigating the capabilities of the meshless technique to address
practical problems, one of the most outstanding issues in meshless methods.
The FPM spatial approximation is studied firstly, with emphasis on aspects of
the methodology that can be improved to increase its robustness and accuracy.
Suitable ranges for setting the relevant approximation parameters and the
performance likely to be attained in practice are determined. An automatic
procedure to adjust the approximation parameters is also proposed to simplify
the application of the method, reducing problem- and user-dependence
without affecting the flexibility of the meshless technique.
The discretization of the flow equations is carried out following wellestablished
approaches, but drawing on the meshless character of the methodology. In order to meet the requirements of practical applications, the procedures are designed and implemented placing emphasis on robustness and efficiency (a simplification of the basic FPM technique is proposed to this end). The flow solver is based on an upwind spatial discretization of the convective fluxes (using the approximate Riemann solver of Roe) and an explicit time integration scheme. Two additional artificial diffusion schemes are also proposed to suit those cases of study in which computational cost is a major concern. The performance of the flow solver is evaluated in order to determine the potential of the meshless approach. The accuracy, computational cost and parallel scalability of the method are studied in
comparison with a conventional FEM-based technique.
Finally, practical applications and extensions of the flow solution scheme are
presented. The examples provided are intended not only to show the
capabilities of the FPM, but also to exploit meshless advantages. Automatic hadaptive procedures, moving domain and fluid-structure interaction problems,
as well as a preliminary approach to solve high-Reynolds viscous flows, are a
sample of the topics explored.
All in all, the results obtained are satisfactorily accurate and competitive in
terms of computational cost (if compared with a similar mesh-based
implementation). This indicates that meshless advantages can be exploited
with efficiency and constitutes a good starting point towards more challenging
applications
On p-Robust Saturation for hp-AFEM
We consider the standard adaptive finite element loop SOLVE, ESTIMATE, MARK,
REFINE, with ESTIMATE being implemented using the -robust equilibrated flux
estimator, and MARK being D\"orfler marking. As a refinement strategy we employ
-refinement. We investigate the question by which amount the local
polynomial degree on any marked patch has to be increase in order to achieve a
-independent error reduction. The resulting adaptive method can be turned
into an instance optimal -adaptive method by the addition of a coarsening
routine
- …