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

Anne E. Staples

Anne E. Staples contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
physics.flu-dyncond-mat.stat-mechhep-thphysics.ao-ph

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
9works
0followers
4topics
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

9 published item(s)

preprint2014arXiv

Explicit Filtering in Large Eddy Simulation of Barotropic Turbulence in Spectral Space

Explicit filtering in large eddy simulation (LES) of a turbulent barotropic flow on the sphere in spectral space is studied and compared to implicit filtering. Here, a smooth filter is applied to the nondivergent barotropic vorticity equation (BVE) on the rotating sphere to divide the flow field into resolved scale (RS) and subfilter scale (SFS) motions. A portion of the SFS motions are reconstructed theoretically using the approximate deconvolution model (ADM). While the unreconstructed portion consists of the subgrid scale (SGS) motions and needs to be modeled separately. In order to investigate the effects of the explicit filtering alone no SGS model is used. It is shown that the explicit filtering accurately tracks the evolution of the coherent structures in two-dimensional turbulent flow on the rotating sphere, whereas the implicit filtering does not. It is also shown that explicit filtering improves the results of the temporal variation of the total kinetic energy and the total enstrophy and the variation of the energy spectrum with wavenumber compare to the implicit filtering. Although explicit filtering is more expensive than implicit filtering it increases the accuracy of the computations and improves the results, particularly where the location of coherent structures is concerned, a topic of particular importance in LES of atmospheric flows for climate and weather applications.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

An Efficient Coarse Grid Projection Method for Quasigeostrophic Models of Large-Scale Ocean Circulation

This paper puts forth a coarse grid projection (CGP) multiscale method to accelerate computations of quasigeostrophic (QG) models for large scale ocean circulation. These models require solving an elliptic sub-problem at each time step, which takes the bulk of the computational time. The method we propose here is a modular approach that facilitates data transfer with simple interpolations and uses black-box solvers for solving the elliptic sub-problem and potential vorticity equations in the QG flow solvers. After solving the elliptic sub-problem on a coarsened grid, an interpolation scheme is used to obtain the fine data for subsequent time stepping on the full grid. The potential vorticity field is then updated on the fine grid with savings in computational time due to the reduced number of grid points for the elliptic solver. The method is applied to both single layer barotropic and two-layer stratified QG ocean models for mid-latitude oceanic basins in the beta plane, which are standard prototypes of more realistic ocean dynamics. The method is found to accelerate these computations while retaining the same level of accuracy in the fine-resolution field. A linear acceleration rate is obtained for all the cases we consider due to the efficient linear-cost fast Fourier transform based elliptic solver used. We expect the speed-up of the CGP method to increase dramatically for versions of the method that use other, suboptimal, elliptic solvers, which are generally quadratic cost. It is also demonstrated that numerical oscillations due to lower grid resolutions, in which the Munk scales are not resolved adequately, are effectively eliminated with CGP method.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

A coarse-grid projection method for accelerating incompressible flow computations

We present a coarse-grid projection (CGP) method for accelerating incompressible flow computations, which is applicable to methods involving Poisson equations as incompressibility constraints. The CGP methodology is a modular approach that facilitates data transfer with simple interpolations and uses black-box solvers for the Poisson and advection-diffusion equations in the flow solver. After solving the Poisson equation on a coarsened grid, an interpolation scheme is used to obtain the fine data for subsequent time stepping on the full grid. A particular version of the method is applied here to the vorticity-stream function, primitive variable, and vorticity-velocity formulations of incompressible Navier-Stokes equations. We compute several benchmark flow problems on two-dimensional Cartesian and non-Cartesian grids, as well as a three-dimensional flow problem. The method is found to accelerate these computations while retaining a level of accuracy close to that of the fine resolution field, which is significantly better than the accuracy obtained for a similar computation performed solely using a coarse grid. A linear acceleration rate is obtained for all the cases we consider due to the linear-cost elliptic Poisson solver used, with reduction factors in computational time between 2 and 42. The computational savings are larger when a suboptimal Poisson solver is used. We also find that the computational savings increase with increasing distortion ratio on non-Cartesian grids, making the CGP method a useful tool for accelerating generalized curvilinear incompressible flow solvers.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

An improved model for reduced-order physiological fluid flows

An improved one-dimensional mathematical model based on Pulsed Flow Equations (PFE) is derived by integrating the axial component of the momentum equation over the transient Womersley velocity profile, providing a dynamic momentum equation whose coefficients are smoothly varying functions of the spatial variable. The resulting momentum equation along with the continuity equation and pressure-area relation form our reduced-order model for physiological fluid flows in one dimension, and are aimed at providing accurate and fast-to-compute global models for physiological systems represented as networks of quasi one-dimensional fluid flows. The consequent nonlinear coupled system of equations is solved by the Lax-Wendroff scheme and is then applied to an open model arterial network of the human vascular system containing the largest fifty-five arteries. The proposed model with functional coefficients is compared with current classical one-dimensional theories which assume steady state Hagen-Poiseuille velocity profiles, either parabolic or plug-like, throughout the whole arterial tree. The effects of the nonlinear term in the momentum equation and different strategies for bifurcation points in the network, as well as the various lumped parameter outflow boundary conditions for distal terminal points are also analyzed. The results show that the proposed model can be used as an efficient tool for investigating the dynamics of reduced-order models of flows in physiological systems and would, in particular, be a good candidate for the one-dimensional, system-level component of geometric multiscale models of physiological systems.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Dynamics of pulsatile flows through elastic microtubes

We investigate the dynamics of pressure driven transient flows of incompressible Newtonian fluids through circular microtubes having thin elastic walls under the long-wavelength and small deformation assumptions, which are valid for many industrial and biological processes. An analytical solution of the coupled fluid and solid equations is found using Navier slip boundary conditions and is shown to include some existing Womersley solutions as limiting cases. The effect of the slip length at the fluid-solid interface is analyzed for oscillatory pressure gradients using a range of slip ratio and frequency parameters. The solutions for elastic and rigid walls are compared for the cases with and without slip boundary conditions for a broad range of the relevant parameters. It is shown that the elastic behavior of the microtube couples nonlinearly with the slip velocity, which greatly enhances the achievable flow rate and pumping efficacy compared to the inelastic case. In addition, it is observed that increasing the slip length produces less shear stress, which is consistent with the nearly frictionless interfaces observed in many microscale experiments.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

High-order methods for decaying two-dimensional homogeneous isotropic turbulence

Numerical schemes used for the integration of complex flow simulations should provide accurate solutions for the long time integrations these flows require. To this end, the performance of various high-order accurate numerical schemes is investigated for direct numerical simulations (DNS) of homogeneous isotropic two-dimensional decaying turbulent flows. The numerical accuracy of compact difference, explicit central difference, Arakawa, and dispersion-relation-preserving schemes are analyzed and compared with the Fourier-Galerkin pseudospectral scheme. In addition, several explicit Runge-Kutta schemes for time integration are investigated. We demonstrate that the centered schemes suffer from spurious Nyquist signals that are generated almost instantaneously and propagate into much of the field when the numerical resolution is insufficient. We further show that the order of the scheme becomes increasingly important for increasing cell Reynolds number. Surprisingly, the sixth-order schemes are found to be in perfect agreement with the pseudospectral method. Considerable reduction in computational time compared to the pseudospectral method is also reported in favor of the finite difference schemes. Among the fourth-order schemes, the compact scheme provides better accuracy than the others for fully resolved computations. The fourth-order Arakawa scheme provides more accurate results for under-resolved computations, however, due to its conservation properties. Our results show that, contrary to conventional wisdom, difference methods demonstrate superior performance in terms of accuracy and efficiency for fully resolved DNS computations of the complex flows considered here. For under-resolved simulations, however, the choice of difference method should be made with care.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Stationary two-dimensional turbulence statistics using a Markovian forcing scheme

In this study we investigate the statistics of two-dimensional stationary turbulence using a Markovian forcing scheme, which correlates the forcing process in the current time step to the previous time step according to a defined memory coefficient. In addition to the Markovian forcing mechanism, the hyperviscous dissipation mechanism for small scales and the Ekman friction type of linear damping mechanism for the large scales are included in the model. We examine the effects of various dissipation and forcing parameters on the turbulence statistics in both wave space and physical space. Our analysis includes the effects of the effective forcing scale, the bandwidth of the forcing, the memory correlation coefficient, and the forcing amplitude, along with the large scale friction and small scale dissipation coefficients. Scaling exponents of structure functions and energy spectra are calculated, and the role of the parameters associated with the Markovian forcing is discussed. We found that the scaling exponents are approximately invariant and show a universal behavior for the various forms of forcing schemes used. We found, however, that the final states strongly depend on the large scale friction mechanism considered. When the large scale friction mechanism is included in the model with a high friction coefficient, we demonstrate that the behavior is no longer universal. Our analysis also shows that the second-order vorticity structure function has an asymptotic scaling exponent for larger dissipation. Additionally, we confirmed that vorticity behaves as a passive scalar when the dissipation mechanism becomes less effective. Finally, although turbulence is not believed to have a separation of time scales in the dynamics of the velocity field, we conjectured that a separation of time scales exists in the dynamics of the energy spectrum.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Approximate deconvolution large eddy simulation of a barotropic ocean circulation model

This paper puts forth a new large eddy simulation closure modeling strategy for two-dimensional turbulent geophysical flows. This closure modeling approach utilizes approximate deconvolution, which is based solely on mathematical approximations and does not employ additional phenomenological arguments to the model. The new approximate deconvolution model is tested in the numerical simulation of the wind-driven circulation in a shallow ocean basin, a standard prototype of more realistic ocean dynamics. The model employs the barotropic vorticity equation driven by a symmetric double-gyre wind forcing, which yields a four-gyre circulation in the time mean. The approximate deconvolution model yields the correct four-gyre circulation structure predicted by a direct numerical simulation, on a coarser mesh but at a fraction of the computational cost. This first step in the numerical assessment of the new model shows that approximate deconvolution could represent a viable tool for under-resolved computations in the large eddy simulation of more realistic turbulent geophysical flows.

0 citations0 reviewsNo code linkedOpen access