BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Ronald D. Haynes

Ronald D. Haynes contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
math.NAmath.OCphysics.comp-phComputationComputational Engineering, Finance, and Sciencemath.APNumerical Analysisphysics.flu-dynphysics.geo-ph

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
10works
0followers
9topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

10 published item(s)

preprint2018arXiv

Moving Mesh simulation of contact sets in two dimensional models of elastic-electrostatic deflection problems

Numerical and analytical methods are developed for the investigation of contact sets in electrostatic-elastic deflections modeling micro-electro mechanical systems. The model for the membrane deflection is a fourth-order semi-linear partial differential equation and the contact events occur in this system as finite time singularities. Primary research interest is in the dependence of the contact set on model parameters and the geometry of the domain. An adaptive numerical strategy is developed based on a moving mesh partial differential equation to dynamically relocate a fixed number of mesh points to increase density where the solution has fine scale detail, particularly in the vicinity of forming singularities. To complement this computational tool, a singular perturbation analysis is used to develop a geometric theory for predicting the possible contact sets. The validity of these two approaches are demonstrated with a variety of test cases.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Stochastic domain decomposition for the solution of the two-dimensional magnetotelluric problem

Stochastic domain decomposition is proposed as a novel method for solving the two-dimensional Maxwell's equations as used in the magnetotelluric method. The stochastic form of the exact solution of Maxwell's equations is evaluated using Monte-Carlo methods taking into consideration that the domain may be divided into neighboring sub-domains. These sub-domains can be naturally chosen by splitting the sub-surface domain into regions of constant (or at least continuous) conductivity. The solution over each sub-domain is obtained by solving Maxwell's equations in the strong form. The sub-domain solver used for this purpose is a meshless method resting on radial basis function based finite differences. The method is demonstrated by solving a number of classical magnetotelluric problems, including the quarter-space problem, the block-in-half-space problem and the triangle-in-half-space problem.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A Closer Look At Differential Evolution For The Optimal Well Placement Problem

Energy demand has increased considerably with the growth of world population, increasing the interest in the hydrocarbon reservoir management problem. Companies are concerned with maximizing oil recovery while minimizing capital investment and operational costs. A first step in solving this problem is to consider optimal well placement. In this work, we investigate the Differential Evolution (DE) optimization method, using distinct configurations with respect to population size, mutation factor, crossover probability, and mutation strategy, to solve the well placement problem. By assuming a bare control procedure, one optimizes the parameters representing positions of injection and production wells. The Tenth SPE Comparative Solution Project and MATLAB Reservoir Simulation Toolbox (MRST) are the benchmark dataset and simulator used, respectively. The goal is to evaluate the performance of DE in solving this important real-world problem. We show that DE can find high-quality solutions, when compared with a reference from the literature, and a preliminary analysis on the results of multiple experiments gives useful information on how DE configuration impacts its performance.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Stochastic domain decomposition for time dependent adaptive mesh generation

The efficient generation of meshes is an important component in the numerical solution of problems in physics and engineering. Of interest are situations where global mesh quality and a tight coupling to the solution of the physical partial differential equation (PDE) is important. We consider parabolic PDE mesh generation and present a method for the construction of adaptive meshes in two spatial dimensions using stochastic domain decomposition that is suitable for an implementation in a multi- or many-core environment. Methods for mesh generation on periodic domains are also provided. The mesh generator is coupled to a time dependent physical PDE and the system is evolved using an alternating solution procedure. The method uses the stochastic representation of the exact solution of a parabolic linear mesh generator to find the location of an adaptive mesh along the (artificial) subdomain interfaces. The deterministic evaluation of the mesh over each subdomain can then be obtained completely independently using the probabilistically computed solutions as boundary conditions. The parallel performance of this general stochastic domain decomposition approach has previously been shown. We demonstrate the approach numerically for the mesh generation context and compare the mesh obtained with the corresponding single domain mesh using a representative mesh quality measure.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Well Control Optimization using Derivative-Free Algorithms and a Multiscale Approach

In this paper, we use numerical optimization algorithms and a multiscale approach in order to find an optimal well management strategy over the life of the reservoir. The large number of well rates for each control step make the optimization problem more difficult and at a high risk of achieving a suboptimal solution. Moreover, the optimal number of adjustments is not known a priori. Adjusting well controls too frequently will increase unnecessary well management and operation cost, and an excessively low number of control adjustments may not be enough to obtain a good yield. We investigate three derivative-free optimization algorithms, chosen for their robust and parallel nature, to determine optimal well control strategies. The algorithms chosen include generalized pattern search (GPS), particle swarm optimization (PSO) and covariance matrix adaptation evolution strategy (CMA-ES). These three algorithms encompass the breadth of available black-box optimization strategies: deterministic local search, stochastic global search and stochastic local search. In addition, we hybridize the three derivative-free algorithms with a multiscale regularization approach. Starting with a reasonably small number of control steps, the control intervals are subsequently refined during the optimization. Results for experiments studied indicate that CMA-ES performs best among the three algorithms in solving both small and large scale problems. When hybridized with a multiscale regularization approach, the ability to find the optimal solution is further enhanced, with the performance of GPS improving the most. Topics affecting the performance of the multiscale approach are discussed in this paper, including the effect of control frequency on the well control problem. The parameter settings for GPS, PSO, and CMA-ES, within the multiscale approach are considered.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

A stochastic domain decomposition method for time dependent mesh generation

We use a time-relaxed linear grid generator of Winslow type to propose a new deterministic-stochastic domain decomposition approach to the generation of adaptive moving meshes. The method uses the probabilistic form of the exact solution of the linear mesh generator to provide the sub-domain interface values that serve as boundary conditions for the domain decomposition. The meshes over the single sub-domains can then be obtained independently of each other in an entirely parallel fashion. We demonstrate numerically the capability of the method to yield smooth, time-evolving meshes.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Joint optimization of well placement and control for nonconventional well types

Optimal well placement and optimal well control are two important areas of study in oilfield development. Although the two problems differ in several respects, both are important considerations in optimizing total oilfield production, and so recent work in the field has considered the problem of addressing both problems jointly. Two general approaches to addressing the joint problem are a simultaneous approach, where all parameters are optimized at the same time, or a sequential approach, where a distinction between placement and control parameters is maintained by separating the optimization problem into two (or more) stages, some of which consider only a subset of the total number of variables. This latter approach divides the problem into smaller ones which are easier to solve, but may not explore search space as fully as a simultaneous approach. In this paper we combine a stochastic global algorithm (Particle Swarm Optimization) and a local search (Mesh Adaptive Direct Search) to compare several simultaneous and sequential approaches to the joint placement and control problem. In particular, we study how increasing the complexity of well models (requiring more variables to describe the well's location and path) affects the respective performances of the two approaches. The results of several experiments with synthetic reservoir models suggest that the sequential approaches are better able to deal with increasingly complex well parameterizations than the simultaneous approaches.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Parallel stochastic methods for PDE based grid generation

The efficient generation of meshes is an important step in the numerical solution of various problems in physics and engineering. We are interested in situations where global mesh quality and tight coupling to the physical solution is important. We consider elliptic PDE based mesh generation and present a method for the construction of adaptive meshes in two spatial dimensions using domain decomposition that is suitable for an implementation on parallel computing architectures. The method uses the stochastic representation of the exact solution of a linear mesh generator of Winslow type to find the points of the adaptive mesh along the subdomain interfaces. The meshes over the single subdomains can then be obtained completely independently of each other using the probabilistically computed solutions along the interfaces as boundary conditions for the linear mesh generator. Further to the previously acknowledged performance characteristics, we demonstrate how the stochastic domain decomposition approach is particularly suited to the problem of grid generation - generating quality meshes efficiently. In addition we show further improvements are possible using interpolation of the subdomain interfaces and smoothing of mesh candidates. An optimal placement strategy is introduced to automatically choose the number and placement of points along the interface using the mesh density function. Various examples of meshes constructed using this stochastic-deterministic domain decomposition technique are shown and compared to the respective single domain solutions using a representative mesh quality measure. A brief performance study is included to show the viability of the stochastic domain decomposition approach and to illustrate the effect of algorithmic choices on the solver's efficiency.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Efficient Optimization of the Likelihood Function in Gaussian Process Modelling

Gaussian Process (GP) models are popular statistical surrogates used for emulating computationally expensive computer simulators. The quality of a GP model fit can be assessed by a goodness of fit measure based on optimized likelihood. Finding the global maximum of the likelihood function for a GP model is typically very challenging as the likelihood surface often has multiple local optima, and an explicit expression for the gradient of the likelihood function is typically unavailable. Previous methods for optimizing the likelihood function (e.g. MacDonald et al. (2013)) have proven to be robust and accurate, though relatively inefficient. We propose several likelihood optimization techniques, including two modified multi-start local search techniques, based on the method implemented by MacDonald et al. (2013), that are equally as reliable, and significantly more efficient. A hybridization of the global search algorithm Dividing Rectangles (DIRECT) with the local optimization algorithm BFGS provides a comparable GP model quality for a fraction of the computational cost, and is the preferred optimization technique when computational resources are limited. We use several test functions and a real application from an oil reservoir simulation to test and compare the performance of the proposed methods with the one implemented by MacDonald et al. (2013) in the R library GPfit. The proposed method is implemented in a Matlab package, GPMfit.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Generating Equidistributed Meshes in 2D via Domain Decomposition

In this paper we consider Schwarz domain decomposition applied to the generation of 2D spatial meshes by a local equidistribution principle. We briefly review the derivation of the local equidistribution principle and the appropriate choice of boundary conditions. We then introduce classical and optimized Schwarz domain decomposition methods to solve the resulting system of nonlinear equations. The implementation of these iterations are discussed, and we conclude with numerical examples to illustrate the performance of the approach.

0 citations0 reviewsNo code linkedOpen access