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

Charles S. Peskin

Charles S. Peskin contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
math.NANumerical Analysisphysics.flu-dynTissues and Organsmath.DSMolecular Networksquant-phSubcellular Processes

Trust snapshot

Quick read

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

7 published item(s)

preprint2026arXiv

State change via one-dimensional scattering in quantum mechanics

This study aims to address the nature of state change, measurement, and probabilistic outcomes in non-relativistic quantum mechanics. We consider a pair of particles that interact in a one-dimensional setting via a delta-function potential. One of the particles is confined to a one-dimensional box, and the other particle is free. The free particle is incident from the left with specified energy, and it may cause changes in state of the confined particle before flying away to the left or to the right. We present a formulation and computational scheme that avoids the use of perturbation theory and determines the probability of any such outcome as a function of the initial state of the confined particle and the energy of the incident particle. As demonstrated by a direct comparison, this presented method holds multiple advantages over a standard perturbative method. The problem formulation and corresponding computational scheme may have applications in physical settings which admit one-dimensional scattering, e.g., in the study of quantum wires or quantum dots.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Optimal fenestration of the Fontan circulation

In this paper, we develop a pulsatile compartmental model of the Fontan circulation and use it to explore the effects of a fenestration added to this physiology. A fenestration is a shunt between the systemic and pulmonary veins that is added either at the time of Fontan conversion or at a later time for the treatment of complications. This shunt increases cardiac output and decreases systemic venous pressure. However, these hemodynamic benefits are achieved at the expense of a decrease in the arterial oxygen saturation. The model developed this paper incorporates fenestration size as a parameter and describes both blood flow and oxygen transport. It is calibrated to clinical data from Fontan patients, and we use it to study the impact of a fenestration on several hemodynamic variables. In certain scenarios corresponding to high-risk Fontan physiology, we demonstrate the existence of an optimal fenestration size that maximizes oxygen delivery to the systemic tissues.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

The hydrodynamics of a twisting, bending, inextensible fiber in Stokes flow

In swimming microorganisms and the cell cytoskeleton, inextensible fibers resist bending and twisting, and interact with the surrounding fluid to cause or resist large-scale fluid motion. In this paper, we develop a novel numerical method for the simulation of cylindrical fibers by extending our previous work on inextensible bending fibers [Maxian et al., Phys. Rev. Fluids 6 (1), 014102] to fibers with twist elasticity. In our "Euler" model, twist is a scalar function that measures the deviation of the fiber cross section relative to a twist-free frame, the fiber exerts only torque parallel to the centerline on the fluid, and the perpendicular components of the rotational fluid velocity are discarded in favor of the translational velocity. In the first part of this paper, we justify this model by comparing it to another commonly-used "Kirchhoff" formulation where the fiber exerts both perpendicular and parallel torque on the fluid, and the perpendicular angular fluid velocity is required to be consistent with the translational fluid velocity. We then develop a spectral numerical method for the hydrodynamics of the Euler model. We define hydrodynamic mobility operators using integrals of the Rotne-Prager-Yamakawa tensor, and evaluate these integrals through a novel slender-body quadrature, which requires on the order of 10 points along the fiber to obtain several digits of accuracy. We demonstrate that this choice of mobility removes the unphysical negative eigenvalues in the translation-translation mobility associated with asymptotic slender body theories, and ensures strong convergence of the fiber velocity and weak convergence of the fiber constraint forces. We pair the spatial discretization with a semi-implicit temporal integrator to confirm the negligible contribution of twist elasticity to the relaxation dynamics of a bent fiber and study the instability of a twirling fiber.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

An Analysis of the Numerical Stability of the Immersed Boundary Method

We present a numerical stability analysis of the immersed boundary(IB) method for a special case which is constructed so that Fourier analysis is applicable. We examine the stability of the immersed boundary method with the discrete Fourier transforms defined differently on the fluid grid and the boundary grid. This approach gives accurate theoretical results about the stability boundary since it takes the effects of the spreading kernel of the immersed boundary method on the numerical stability into account. In this paper, the spreading kernel is the standard 4-point IB delta function. A three-dimensional incompressible viscous flow and a no-slip planar boundary are considered. The case of a planar elastic membrane is also analyzed using the same analysis framework and it serves as an example of many possible generalizations of our theory. We present some numerical results and show that the observed stability behaviors are consistent with what are predicted by our theory.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

An immersed boundary method with subgrid resolution and improved numerical stability applied to slender bodies in Stokes flow

The immersed boundary method is a numerical and mathematical formulation for solving fluid-structure interaction problems. It relies on solving fluid equations on an Eulerian fluid grid and interpolating the resulting velocity back onto immersed structures. To resolve slender fibers, the grid spacing must be on the order of the fiber radius, and thus the number of required grid points along the filament must be of the same order as the aspect ratio. Simulations of slender bodies using the IB method can therefore be intractable. A technique is presented to address this problem in the context of Stokes flow. The velocity of the structure is split into a component coming from the underlying fluid grid, which is coarser than normally required, and a component proportional to the force (a drag term). The drag coefficient is set so that a single sphere is represented exactly on a grid of arbitrary meshwidth. Implicit treatment of the drag term removes some of the stability restrictions normally associated with the IB method. This comes at a loss of accuracy, although tests are conducted that show 1-2 digits of relative accuracy can be obtained on coarser grids. After its accuracy and stability are tested, the method is applied to two real world examples: fibers in shear flow and a suspension of fibers. These examples show that the method can reproduce existing results and make reasonable predictions about the viscosity of an aligned fiber suspension.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Models for plasma kinetics during simultaneous therapeutic plasma exchange and extracorporeal membrane oxygenation

This paper focuses on the derivation and simulation of mathematical models describing new plasma fraction in blood for patients undergoing simultaneous extracorporeal membrane oxygenation and therapeutic plasma exchange. Models for plasma exchange with either veno-arterial or veno-venous extracorporeal membrane oxygenation are considered. Two classes of models are derived for each case, one in the form of an algebraic delay equation and another in the form of a system of delay differential equations. In special cases, our models reduce to single compartment ones for plasma exchange that have been validated with experimental data. We also show that the algebraic delay equations are forward Euler discretizations of the delay differential equations, with timesteps equal to transit times through model compartments. Numerical simulations are performed to compare different model types, to investigate the impact of plasma device port switching on the efficiency of the exchange process, and to study the sensitivity of the models to their parameters.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The Influence of Spatial Variation in Chromatin Density Determined by X-ray Tomograms on the Time to Find DNA Binding Sites

In this work we examine how volume exclusion caused by regions of high chromatin density might influence the time required for proteins to find specific DNA binding sites. The spatial variation of chromatin density within mouse olfactory sensory neurons is determined from soft X-ray tomography reconstructions of five nuclei. We show that there is a division of the nuclear space into regions of low-density euchromatin and high-density heterochromatin. Volume exclusion experienced by a diffusing protein caused by this varying density of chromatin is modeled by a repulsive potential. The value of the potential at a given point in space is chosen to be proportional to the density of chromatin at that location. The constant of proportionality, called the volume exclusivity, provides a model parameter that determines the strength of volume exclusion. Numerical simulations demonstrate that the mean time for a protein to locate a binding site localized in euchromatin is minimized for a finite, non-zero volume exclusivity. For binding sites in heterochromatin, the mean time is minimized when the volume exclusivity is zero (the protein experiences no volume exclusion). An analytical theory is developed to explain these results. The theory suggests that for binding sites in euchromatin there is an optimal level of volume exclusivity that balances a reduction in the volume searched in finding the binding site, with the height of effective potential barriers the protein must cross during the search process.

0 citations0 reviewsNo code linkedOpen access