9,174 research outputs found
A Generalized Multiscale Finite Element Method for the Brinkman Equation
In this paper we consider the numerical upscaling of the Brinkman equation in
the presence of high-contrast permeability fields. We develop and analyze a
robust and efficient Generalized Multiscale Finite Element Method (GMsFEM) for
the Brinkman model. In the fine grid, we use mixed finite element method with
the velocity and pressure being continuous piecewise quadratic and piecewise
constant finite element spaces, respectively. Using the GMsFEM framework we
construct suitable coarse-scale spaces for the velocity and pressure that yield
a robust mixed GMsFEM. We develop a novel approach to construct a coarse
approximation for the velocity snapshot space and a robust small offline space
for the velocity space. The stability of the mixed GMsFEM and a priori error
estimates are derived. A variety of two-dimensional numerical examples are
presented to illustrate the effectiveness of the algorithm.Comment: 22 page
Multilevel Markov Chain Monte Carlo Method for High-Contrast Single-Phase Flow Problems
In this paper we propose a general framework for the uncertainty
quantification of quantities of interest for high-contrast single-phase flow
problems. It is based on the generalized multiscale finite element method
(GMsFEM) and multilevel Monte Carlo (MLMC) methods. The former provides a
hierarchy of approximations of different resolution, whereas the latter gives
an efficient way to estimate quantities of interest using samples on different
levels. The number of basis functions in the online GMsFEM stage can be varied
to determine the solution resolution and the computational cost, and to
efficiently generate samples at different levels. In particular, it is cheap to
generate samples on coarse grids but with low resolution, and it is expensive
to generate samples on fine grids with high accuracy. By suitably choosing the
number of samples at different levels, one can leverage the expensive
computation in larger fine-grid spaces toward smaller coarse-grid spaces, while
retaining the accuracy of the final Monte Carlo estimate. Further, we describe
a multilevel Markov chain Monte Carlo method, which sequentially screens the
proposal with different levels of approximations and reduces the number of
evaluations required on fine grids, while combining the samples at different
levels to arrive at an accurate estimate. The framework seamlessly integrates
the multiscale features of the GMsFEM with the multilevel feature of the MLMC
methods following the work in \cite{ketelson2013}, and our numerical
experiments illustrate its efficiency and accuracy in comparison with standard
Monte Carlo estimates.Comment: 29 pages, 6 figure
A review of surrogate models and their application to groundwater modeling
The spatially and temporally variable parameters and inputs to complex groundwater models typically result in long runtimes which hinder comprehensive calibration, sensitivity, and uncertainty analysis. Surrogate modeling aims to provide a simpler, and hence faster, model which emulates the specified output of a more complex model in function of its inputs and parameters. In this review paper, we summarize surrogate modeling techniques in three categories: data-driven, projection, and hierarchical-based approaches. Data-driven surrogates approximate a groundwater model through an empirical model that captures the input-output mapping of the original model. Projection-based models reduce the dimensionality of the parameter space by projecting the governing equations onto a basis of orthonormal vectors. In hierarchical or multifidelity methods the surrogate is created by simplifying the representation of the physical system, such as by ignoring certain processes, or reducing the numerical resolution. In discussing the application to groundwater modeling of these methods, we note several imbalances in the existing literature: a large body of work on data-driven approaches seemingly ignores major drawbacks to the methods; only a fraction of the literature focuses on creating surrogates to reproduce outputs of fully distributed groundwater models, despite these being ubiquitous in practice; and a number of the more advanced surrogate modeling methods are yet to be fully applied in a groundwater modeling context
Generalized multiscale finite element methods (GMsFEM)
In this paper, we propose a general approach called Generalized Multiscale Finite Element Method (GMsFEM) for performing multiscale simulations for problems without scale separation over a complex input space. As in multiscale finite element methods (MsFEMs), the main idea of the proposed approach is to construct a small dimensional local solution space that can be used to generate an efficient and accurate approximation to the multiscale solution with a potentially high dimensional input parameter space. In the proposed approach, we present a general procedure to construct the offline space that is used for a systematic enrichment of the coarse solution space in the online stage. The enrichment in the online stage is performed based on a spectral decomposition of the offline space. In the online stage, for any input parameter, a multiscale space is constructed to solve the global problem on a coarse grid. The online space is constructed via a spectral decomposition of the offline space and by choosing the eigenvectors corresponding to the largest eigenvalues. The computational saving is due to the fact that the construction of the online multiscale space for any input parameter is fast and this space can be re-used for solving the forward problem with any forcing and boundary condition. Compared with the other approaches where global snapshots are used, the local approach that we present in this paper allows us to eliminate unnecessary degrees of freedom on a coarse-grid level. We present various examples in the paper and some numerical results to demonstrate the effectiveness of our method
Mini-Workshop: Numerical Upscaling for Media with Deterministic and Stochastic Heterogeneity
This minisymposium was third in series of similar events, after two very successful meetings in 2005 and 2009. The aim was to provide a forum for an extensive discussion on the theoretical aspects and on the areas of application and validity of numerical upscaling approaches for heterogeneous problems with deterministic and stochastic coefficients. The intensive discussions during the meeting contributed to a better understanding of upscaling approaches for multiscale problems with stochastic coefficients, and for synergy between scientists coming to this topic from the area of deterministic multiscale problems on one hand, and those coming from the area of SPDE on the other hand. Recent advanced results on upscaling approaches for deterministic multiscale problems were presented, well mixed with strong presentations on SDE and SPDE. The open problems in these areas were discussed, with emphasis on the case of stochastic coefficients brainstorming numerous numerical upscaling approaches. A number of young researchers, very actively working in these areas, were involved in the workshop discussing the links between scales., thus ensuring the continuity between the generations of researchers
Anomalous transport in the crowded world of biological cells
A ubiquitous observation in cell biology is that diffusion of macromolecules
and organelles is anomalous, and a description simply based on the conventional
diffusion equation with diffusion constants measured in dilute solution fails.
This is commonly attributed to macromolecular crowding in the interior of cells
and in cellular membranes, summarising their densely packed and heterogeneous
structures. The most familiar phenomenon is a power-law increase of the MSD,
but there are other manifestations like strongly reduced and time-dependent
diffusion coefficients, persistent correlations, non-gaussian distributions of
the displacements, heterogeneous diffusion, and immobile particles. After a
general introduction to the statistical description of slow, anomalous
transport, we summarise some widely used theoretical models: gaussian models
like FBM and Langevin equations for visco-elastic media, the CTRW model, and
the Lorentz model describing obstructed transport in a heterogeneous
environment. Emphasis is put on the spatio-temporal properties of the transport
in terms of 2-point correlation functions, dynamic scaling behaviour, and how
the models are distinguished by their propagators even for identical MSDs.
Then, we review the theory underlying common experimental techniques in the
presence of anomalous transport: single-particle tracking, FCS, and FRAP. We
report on the large body of recent experimental evidence for anomalous
transport in crowded biological media: in cyto- and nucleoplasm as well as in
cellular membranes, complemented by in vitro experiments where model systems
mimic physiological crowding conditions. Finally, computer simulations play an
important role in testing the theoretical models and corroborating the
experimental findings. The review is completed by a synthesis of the
theoretical and experimental progress identifying open questions for future
investigation.Comment: review article, to appear in Rep. Prog. Phy
- …