BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Donald L. Brown

Donald L. Brown appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
math.NANumerical AnalysisComputational Engineering, Finance, and Sciencemath-phmath.MPphysics.comp-ph

Catalog footprint

What is connected

7works
6topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

7 published item(s)

preprint2022arXiv

Multiscale Sub-grid Correction Method for Time-Harmonic High-Frequency Elastodynamics with Wavenumber Explicit Bounds

The simulation of the elastodynamics equations at high-frequency suffers from the well known pollution effect. We present a Petrov--Galerkin multiscale sub-grid correction method that remains pollution free in natural resolution and oversampling regimes. This is accomplished by generating corrections to coarse-grid spaces with supports determined by oversampling lengths related to the $\log(k)$, $k$ being the wavenumber. Key to this method are polynomial-in-$k$ bounds for stability constants and related inf-sup constants. To this end, we establish polynomial-in-$k$ bounds for the elastodynamics stability constants in general Lipschitz domains with radiation boundary conditions in $\mathbb{R}^3$. Previous methods relied on variational techniques, Rellich identities, and geometric constraints. In the context of elastodynamics, these suffer from the need to hypothesize a Korn's inequality on the boundary. The methods in this work are based on boundary integral operators and estimation of Green's function's derivatives dependence on $k$ and do not require this extra hypothesis. We also implemented numerical examples in two and three dimensions to show the method eliminates pollution in the natural resolution and oversampling regimes, as well as performs well when compared to standard Lagrange finite elements.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Preconditioning Markov Chain Monte Carlo Method for Geomechanical Subsidence using multiscale method and machine learning technique

In this paper, we consider the numerical solution of the poroelasticity problem with stochastic properties. We present a Two-stage Markov Chain Monte Carlo method for geomechanical subsidence. In this work, we study two techniques of preconditioning: (MS) multiscale method for model order reduction and (ML) machine learning technique. The purpose of preconditioning is the fast sampling, where a new proposal is first testes by a cheap multiscale solver or using fast prediction of the neural network and the full fine grid computations will be conducted only if the proposal passes the first step. To construct a reduced order model, we use the Generalized Multiscale Finite Element Method and present construction of the multiscale basis functions for pressure and displacements in stochastic fields. In order to construct a machine learning based preconditioning, we generate a dataset using a multiscale solver and use it to train neural networks. The Karhunen-Loeve expansion is used to represent the realization of the stochastic field. Numerical results are presented for two- and three-dimensional model examples.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A Generalized Multiscale Finite Element Method for Poroelasticity Problems I: Linear Problems

In this paper, we consider the numerical solution of poroelasticity problems that are of Biot type and develop a general algorithm for solving coupled systems. We discuss the challenges associated with mechanics and flow problems in heterogeneous media. The two primary issues being the multiscale nature of the media and the solutions of the fluid and mechanics variables traditionally developed with separate grids and methods. For the numerical solution we develop and implement a Generalized Multiscale Finite Element Method (GMsFEM) that solves problem on a coarse grid by constructing local multiscale basis functions. The procedure begins with construction of multiscale bases for both displacement and pressure in each coarse block. Using a snapshot space and local spectral problems, we construct a basis of reduced dimension. Finally, after multiplying by a multiscale partitions of unity, the multiscale basis is constructed in the offline phase and the coarse grid problem then can be solved for arbitrary forcing and boundary conditions. We implement this algorithm on two heterogenous media and compute error between the multiscale solution with the fine-scale solutions. Randomized oversampling and forcing strategies are also tested.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A Generalized Multiscale Finite Element Method for Poroelasticity Problems II: Nonlinear Coupling

In this paper, we consider the numerical solution of some nonlinear poroelasticity problems that are of Biot type and develop a general algorithm for solving nonlinear coupled systems. We discuss the difficulties associated with flow and mechanics in heterogenous media with nonlinear coupling. The central issue being how to handle the nonlinearities and the multiscale scale nature of the media. To compute an efficient numerical solution we develop and implement a Generalized Multiscale Finite Element Method (GMsFEM) that solves nonlinear problems on a coarse grid by constructing local multiscale basis functions and treating part of the nonlinearity locally as a parametric value. After linearization with a Picard Iteration, the procedure begins with construction of multiscale bases for both displacement and pressure in each coarse block by treating the staggered nonlinearity as a parametric value. Using a snapshot space and local spectral problems, we construct an offline basis of reduced dimension. From here an online, parametric dependent, space is constructed. Finally, after multiplying by a multiscale partitions of unity, the multiscale basis is constructed and the coarse grid problem then can be solved for arbitrary forcing and boundary conditions. We implement this algorithm on a geometry with a linear and nonlinear pressure dependent permeability field and compute error between the multiscale solution with the fine-scale solutions.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

An energy-stable convex splitting for the phase-field crystal equation

The phase-field crystal equation, a parabolic, sixth-order and nonlinear partial differential equation, has generated considerable interest as a possible solution to problems arising in molecular dynamics. This is because the phase-field crystal model can capture atomic-scale effects at time-scales that are orders of magnitude larger than what molecular dynamics simulations can afford presently. Nonetheless, solving this equation is not a trivial task, as a non-increasing free energy and mass conservation need to be verified for the numerical solution to be valid. This work focuses on these issues, and proposes a novel algorithm that guarantees mass conservation, unconditional energy stability and is second-order accurate in time. This is achieved through a convex-concave splitting of the nonlinearity present in the equation, along with the use of a stabilization term that bounds possible increases in free energy. We present numerical results that validate our mathematical proofs, and show two and three dimensional simulations involving crystal growth that showcase the robustness of the method.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Multiscale Petrov-Galerkin Method for High-Frequency Heterogeneous Helmholtz Equations

This paper presents a multiscale Petrov-Galerkin finite element method for time-harmonic acoustic scattering problems with heterogeneous coefficients in the high-frequency regime. We show that the method is pollution- free also in the case of heterogeneous media provided that the stability bound of the continuous problem grows at most polynomially with the wave number k. By generalizing classical estimates of [Melenk, Ph.D. Thesis 1995] and [Hetmaniuk, Commun. Math. Sci. 5 (2007)] for homogeneous medium, we show that this assumption of polynomially wave number growth holds true for a particular class of smooth heterogeneous material coefficients. Further, we present numerical examples to verify our stability estimates and implement an example in the wider class of discontinuous coefficients to show computational applicability beyond our limited class of coefficients.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

A Multiscale Method for Porous Microstructures

In this paper we develop a multiscale method to solve problems in complicated porous microstructures with Neumann boundary conditions. By using a coarse-grid quasi-interpolation operator to define a fine detail space and local orthogonal decomposition, we construct multiscale corrections to coarse-grid basis functions with microstructure. By truncating the corrector functions we are able to make a computationally efficient scheme. Error results and analysis are presented. A key component of this analysis is the investigation of the Poincaré constants in perforated domains as they may contain micro-structural information. Using a constructive method originally developed for weighted Poincaré inequalities, we are able to obtain estimates on Poincaré constants with respect to scale and separation length of the pores. Finally, two numerical examples are presented to verify our estimates.

0 citations0 reviewsNo code linkedOpen access