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

John M. Stockie

John M. Stockie appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
physics.flu-dynmath.NAmath.APBiological Physicsphysics.comp-phApplicationsComputational Engineering, Finance, and ScienceDistributed, Parallel, and Cluster Computingmath.FAmath.HOphysics.ao-ph

Catalog footprint

What is connected

18works
11topics
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

18 published item(s)

preprint2016arXiv

Bayesian estimation of airborne fugitive emissions using a Gaussian plume model

A new method is proposed for estimating the rate of fugitive emissions of particulate matter from multiple time-dependent sources via measurements of deposition and concentration. We cast this source inversion problem within the Bayesian framework, and use a forward model based on a Gaussian plume solution. We present three alternate models for constructing the prior distribution on the emission rates as functions of time. Next, we present an industrial case study in which our framework is applied to estimate the rate of fugitive emissions of lead particulates from a smelter in Trail, British Columbia, Canada. The Bayesian framework not only provides an approximate solution to the inverse problem, but also quantifies the uncertainty in the solution. Using this information we perform an uncertainty propagation study in order to assess the impact of the estimated sources on the area surrounding the industrial site.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

An immersed boundary model of the cochlea with parametric forcing

The cochlea or inner ear has a remarkable ability to amplify sound signals. This is understood to derive at least in part from some active process that magnifies vibrations of the basilar membrane (BM) and the cochlear partition in which it is embedded, to the extent that it overcomes the effect of viscous damping from the surrounding cochlear fluid. Many authors have associated this amplification ability to some type of mechanical resonance within the cochlea, however there is still no consensus regarding the precise cause of amplification. Our work is inspired by experiments showing that the outer hair cells within the cochlear partition change their lengths when stimulated, which can in turn cause periodic distortions of the BM and other structures in the cochlea. This paper investigates a novel fluid-mechanical resonance mechanism that derives from hydrodynamic interactions between an oscillating BM and the surrounding cochlear fluid. We present a model of the cochlea based on the immersed boundary method, in which a small-amplitude periodic internal forcing due to outer hair cells can induce parametric resonance. A Floquet stability analysis of the linearized equations demonstrates the existence of resonant (unstable) solutions within the range of physical parameters corresponding to the human auditory system. Numerical simulations of the immersed boundary equations support the analytical results and clearly demonstrate the existence of resonant solution modes. These results are then used to illustrate the influence of parametric resonance on wave propagation along the BM and explicit comparisons are drawn with results from another two-dimensional cochlea model.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Mathematics For Industry: A Personal Perspective

"I am an industrial mathematician." When asked to identify my profession or academic field of study, this is the most concise answer I can provide. However, this seemingly straightforward statement is commonly greeted by a blank stare or an uncomfortable silence, regardless of whether I am speaking to a fellow mathematician or a non-mathematician. I usually follow up with the clarification: "I am an applied mathematician who derives much of my inspiration from the study of industrial problems that I encounter through collaborations with companies." This dispels some confusion, but unfortunately still leaves a great deal open to interpretation owing to the vagueness of the words mathematics, industry and company, each of which covers an extremely broad range of scientific or socio-economic activity. To those academics who actually work in the field of industrial mathematics (and whose "perspective" referred to in the title is the focus of this article) this ambiguity is familiar and untroubling. However, for anyone less acquainted with the work of industrial mathematicians, some clarification is desirable especially for anyone who might be considering entering the field. This essay therefore aims to shed light upon the nature of research being done at the interface between mathematics and industry, paying particular attention to the following questions: What is industrial mathematics? Where is industrial mathematics? How does one do industrial mathematics? Why (or more precisely, what value is there in doing) industrial mathematics? I will attempt to answer these questions by means of several case studies drawn from my own experience in tackling mathematical problems from industry.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Multi-layer asymptotic solution for wetting fronts in porous media with exponential moisture diffusivity

We study the asymptotic behaviour of sharp front solutions arising from the nonlinear diffusion equation θ_t = (D(θ)θ_x)_x, where the diffusivity is an exponential function D(θ) = D_o exp(βθ). This problem arises for example in the study of unsaturated flow in porous media where θ represents the liquid saturation. For physical parameters corresponding to actual porous media, the diffusivity at the residual saturation is D(0) = D_o << 1 so that the diffusion problem is nearly degenerate. Such problems are characterised by wetting fronts that sharply delineate regions of saturated and unsaturated flow, and that propagate with a well-defined speed. Using matched asymptotic expansions in the limit of large β, we derive an analytical description of the solution that is uniformly valid throughout the wetting front. This is in contrast with most other related analyses that instead truncate the solution at some specific wetting front location, which is then calculated as part of the solution, and beyond that location the solution is undefined. Our asymptotic analysis demonstrates that the solution has a four-layer structure, and by matching through the adjacent layers we obtain an estimate of the wetting front location in terms of the material parameters describing the porous medium. Using numerical simulations of the original nonlinear diffusion equation, we demonstrate that the first few terms in our series solution provide approximations of physical quantities such as wetting front location and speed of propagation that are more accurate (over a wide range of admissible β values) than other asymptotic approximations reported in the literature.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Multiscale model of a freeze-thaw process for tree sap exudation

Sap transport in trees has long fascinated scientists, and a vast literature exists on experimental and modelling studies of trees during the growing season when large negative stem pressures are generated by transpiration from leaves. Much less attention has been paid to winter months when trees are largely dormant but nonetheless continue to exhibit interesting flow behaviour. A prime example is sap exudation, which refers to the peculiar ability of sugar maple (Acer saccharum) and related species to generate positive stem pressure while in a leafless state. Experiments demonstrate that ambient temperatures must oscillate about the freezing point before significantly heightened stem pressures are observed, but the precise causes of exudation remain unresolved. The prevailing hypothesis attributes exudation to a physical process combining freeze-thaw and osmosis, which has some support from experimental studies but remains a subject of active debate. We address this knowledge gap by developing the first mathematical model for exudation, while also introducing several essential modifications to this hypothesis. We derive a multiscale model consisting of a nonlinear system of differential equations governing phase change and transport within wood cells, coupled to a suitably homogenized equation for temperature on the macroscale. Numerical simulations yield stem pressures that are consistent with experiments and provide convincing evidence that a purely physical mechanism is capable of capturing exudation.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

On regularizations of the delta distribution

In this article we consider regularizations of the Dirac delta distribution with applications to prototypical elliptic and hyperbolic partial differential equations (PDEs). We study the convergence of a sequence of distributions $\mathcal{S}_H$ to a singular term $\mathcal{S}$ as a parameter $H$ (associated with the {support size} of $\mathcal{S}_H$) shrinks to zero. We characterize this convergence in both the weak-$\ast$ topology of distributions, as well as in a weighted Sobolev norm. These notions motivate a framework for constructing regularizations of the delta distribution that includes a large class of existing methods in the literature. This framework allows different regularizations to be compared. The convergence of solutions of PDEs with these regularized source terms is then studied in various topologies such as pointwise convergence on a deleted neighborhood and weighted Sobolev norms. We also examine the lack of symmetry in tensor product regularizations and effects of dissipative error in hyperbolic problems.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Simulating biofilm deformation and detachment with the immersed boundary method

We apply the immersed boundary (or IB) method to simulate deformation and detachment of a periodic array of wall-bounded biofilm colonies in response to a linear shear flow. The biofilm material is represented as a network of Hookean springs that are placed along the edges of a triangulation of the biofilm region. The interfacial shear stress, lift and drag forces acting on the biofilm colony are computed by using fluid stress jump method developed by Williams, Fauci and Gaver [Disc. Contin. Dyn. Sys. B 11(2):519-540, 2009], with a modified version of their exclusion filter. Our detachment criterion is based on the novel concept of an averaged equivalent continuum stress tensor defined at each IB point in the biofilm which is then used to determine a corresponding von Mises yield stress; wherever this yield stress exceeds a given critical threshold the connections to that node are severed, thereby signalling the onset of a detachment event. In order to capture the deformation and detachment behaviour of a biofilm colony at different stages of growth, we consider a family of four biofilm shapes with varying aspect ratio. Our numerical simulations focus on the behaviour of weak biofilms (with relatively low yield stress threshold) and investigate features of the fluid-structure interaction such as locations of maximum shear and increased drag. The most important conclusion of this work is that the commonly employed detachment strategy in biofilm models based only on interfacial shear stress can lead to incorrect or inaccurate results when applied to the study of shear induced detachment of weak biofilms. Our detachment strategy based on equivalent continuum stresses provides a unified and consistent IB framework that handles both sloughing and erosion modes of biofilm detachment, and is consistent with strategies employed in many other continuum based biofilm models.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

An efficient parallel immersed boundary algorithm using a pseudo-compressible fluid solver

We propose an efficient algorithm for the immersed boundary method on distributed-memory architectures, with the computational complexity of a completely explicit method and excellent parallel scaling. The algorithm utilizes the pseudo-compressibility method recently proposed by Guermond and Minev [Comptes Rendus Mathematique, 348:581-585, 2010] that uses a directional splitting strategy to discretize the incompressible Navier-Stokes equations, thereby reducing the linear systems to a series of one-dimensional tridiagonal systems. We perform numerical simulations of several fluid-structure interaction problems in two and three dimensions and study the accuracy and convergence rates of the proposed algorithm. For these problems, we compare the proposed algorithm against other second-order projection-based fluid solvers. Lastly, the strong and weak scaling properties of the proposed algorithm are investigated.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Correction to "Parametric Resonance in Immersed Elastic Boundaries"

This note is a correction to a paper of Cortez, Peskin, Stockie & Varela [SIAM J. Appl. Math., 65(2):494-520, 2004], who studied the stability of a parametrically-forced, circular, elastic fiber immersed in an incompressible fluid in 2D, and showed the existence of parametric resonance. The results were represented as plots that separate parameter space into regions where the solution is either stable or unstable. We uncovered two errors in the paper: the first was in the derivation of the eigenvalue problem, and the second was in the code to used to calculate the stability contours.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Homogenization of the Stefan problem, with application to maple sap exudation

The technique of periodic homogenization with two-scale convergence is applied to the analysis of a two-phase Stefan-type problem that arises in the study of a periodic array of melting ice bars. For this "reduced model" we prove results on existence, uniqueness and convergence of the two-scale limit solution in the weak form, which requires solving a macroscale problem for the global temperature field and a reference cell problem at each point in space which captures the underlying phase change process occurring on the microscale. We state a corresponding strong formulation of the limit problem and use it to design an efficient numerical solution algorithm. The same homogenized temperature equations are then applied to solve a much more complicated problem involving multi-phase flow and heat transport in trees, where the sap is present in both frozen and liquid forms and a third gas phase is also present. Our homogenization approach has the advantage that the global temperature field is a solution of the same reduced model equations, while all the remaining physics are relegated to the reference cell problem. Numerical simulations are performed to validate our results and draw conclusions regarding the phenomenon known as sap exudation, which is of great importance in sugar maple trees and few other related species.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Parametric resonance in spherical immersed elastic shells

We perform a stability analysis for a fluid-structure interaction problem in which a spherical elastic shell or membrane is immersed in a 3D viscous, incompressible fluid. The shell is an idealised structure having zero thickness, and has the same fluid lying both inside and outside. The problem is formulated mathematically using the immersed boundary framework in which Dirac delta functions are employed to capture the two-way interaction between fluid and immersed structure. The elastic structure is driven parametrically via a time-periodic modulation of the elastic membrane stiffness. We perform a Floquet stability analysis, considering the case of both a viscous and inviscid fluid, and demonstrate that the forced fluid-membrane system gives rise to parametric resonances in which the solution becomes unbounded even in the presence of viscosity. The analytical results are validated using numerical simulations with a 3D immersed boundary code for a range of wavenumbers and physical parameter values. Finally, potential applications to biological systems are discussed, with a particular focus on the human heart and investigating whether or not fluid-structure interaction-mediated instabilities could play a role in cardiac fluid dynamics.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Simulating flexible fiber suspensions using a scalable immersed boundary algorithm

We present an approach for numerically simulating the dynamics of flexible fibers in a three-dimensional shear flow using a scalable immersed boundary (IB) algorithm based on Guermond and Minev's pseudo-compressible fluid solver. The fibers are treated as one-dimensional Kirchhoff rods that resist stretching, bending, and twisting, within the generalized IB framework. We perform a careful numerical comparison against experiments on single fibers performed by S. G. Mason and co-workers, who categorized the fiber dynamics into several distinct orbit classes. We show that the orbit class may be determined using a single dimensionless parameter for low Reynolds flows. Lastly, we simulate dilute suspensions containing up to hundreds of fibers using a distributed- memory computer cluster. These simulations serve as a stepping stone for studying more complex suspension dynamics including non-dilute suspensions and aggregation of fibers (also known as flocculation).

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Numerical simulations of particle sedimentation using the immersed boundary method

We study the settling of solid particles within a viscous incompressible fluid contained in a two-dimensional channel, where the mass density of the particles is slightly greater than that of the fluid. The fluid-structure interaction problem is simulated numerically using the immersed boundary method, with an added mass term that is incorporated using a Boussinesq approximation. Simulations are performed with a single circular particle, and also with two particles in various initial configurations. The terminal settling velocities for the particles correspond closely with both theoretical and experimental results, and the single-particle dynamics reproduce expected behavior qualitatively. The two-particle simulations exhibit drafting-kissing-tumbling dynamics that is similar to what is observed in other experimental and numerical studies.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

A mathematical model of sap exudation in maple trees governed by ice melting, gas dissolution and osmosis

We develop a mathematical model for sap exudation in a maple tree that is based on a purely physical mechanism for internal pressure generation in trees in the leafless state. There has been a long-standing controversy in the tree physiology literature over precisely what mechanism drives sap exudation, and we aim to cast light on this issue. Our model is based on the work of Milburn and O'Malley [Can. J. Bot., 62(10):2101-2106, 1984] who hypothesized that elevated sap pressures derive from compressed gas that is trapped within certain wood cells and subsequently released when frozen sap thaws in the spring. We also incorporate the extension of Tyree [in Tree Sap, pp. 37-45, eds. M. Terazawa et al., Hokkaido Univ. Press, 1995] who argued that gas bubbles are prevented from dissolving because of osmotic pressure that derives from differences in sap sugar concentrations and the selective permeability of cell walls. We derive a system of differential-algebraic equations based on conservation principles that is used to test the validity of the Milburn-O'Malley hypothesis and also to determine the extent to which osmosis is required. This work represents the first attempt to derive a detailed mathematical model of sap exudation at the micro-scale.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Riemann solver for a kinematic wave traffic model with discontinuous flux

We investigate a model for traffic flow based on the Lighthill-Whitham-Richards model that consists of a hyperbolic conservation law with a discontinuous, piecewise-linear flux. A mollifier is used to smooth out the discontinuity in the flux function over a small distance epsilon << 1 and then the analytical solution to the corresponding Riemann problem is derived in the limit as epsilon goes to 0. For certain initial data, the Riemann problem can give rise to zero waves that propagate with infinite speed but have zero strength. We propose a Godunov-type numerical scheme that avoids the otherwise severely restrictive CFL constraint that would arise from waves with infinite speed by exchanging information between local Riemann problems and thereby incorporating the effects of zero waves directly into the Riemann solver. Numerical simulations are provided to illustrate the behaviour of zero waves and their impact on the solution. The effectiveness of our approach is demonstrated through a careful convergence study and comparisons to computations using a third-order WENO scheme.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Numerical Simulations of Gravity-Driven Fingering in Unsaturated Porous Media Using a Non-Equilibrium Model

This is a computational study of gravity-driven fingering instabilities in unsaturated porous media. The governing equations and corresponding numerical scheme are based on the work of Nieber et al. [Ch. 23 in Soil Water Repellency, eds. C. J. Ritsema and L. W. Dekker, Elsevier, 2003] in which non-monotonic saturation profiles are obtained by supplementing the Richards equation with a non-equilibrium capillary pressure-saturation relationship, as well as including hysteretic effects. The first part of the study takes an extensive look at the sensitivity of the finger solutions to certain key parameters in the model such as capillary shape parameter, initial saturation, and capillary relaxation coefficient. The second part is a comparison to published experimental results that demonstrates the ability of the model to capture realistic fingering behaviour.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

An inverse Gaussian plume approach for estimating atmospheric pollutant emissions from multiple point sources

A method is developed for estimating the emission rates of contaminants into the atmosphere from multiple point sources using measurements of particulate material deposited at ground level. The approach is based on a Gaussian plume type solution for the advection-diffusion equation with ground-level deposition and given emission sources. This solution to the forward problem is incorporated into an inverse algorithm for estimating the emission rates by means of a linear least squares approach. The results are validated using measured deposition and meteorological data from a large lead-zinc smelting operation in Trail, British Columbia. The algorithm is demonstrated to be robust and capable of generating reasonably accurate estimates of total contaminant emissions over the relatively short distances of interest in this study.

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

A moving mesh method with variable relaxation time

We propose a moving mesh adaptive approach for solving time-dependent partial differential equations. The motion of spatial grid points is governed by a moving mesh PDE (MMPDE) in which a mesh relaxation time τis employed as a regularization parameter. Previously reported results on MMPDEs have invariably employed a constant value of the parameter τ. We extend this standard approach by incorporating a variable relaxation time that is calculated adaptively alongside the solution in order to regularize the mesh appropriately throughout a computation. We focus on singular problems involving self-similar blow-up to demonstrate the advantages of using a variable relaxation ime over a fixed one in terms of accuracy, stability and efficiency.

0 citations0 reviewsNo code linkedOpen access