skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.

Attention:

The NSF Public Access Repository (PAR) system and access will be unavailable from 10:00 PM to 12:00 PM ET on Tuesday, March 25 due to maintenance. We apologize for the inconvenience.

Explore Research Products in the PAR
It may take a few hours for recently added research products to appear in PAR search results.


Title: Local sensitivity analysis of the “membrane shape equation” derived from the Helfrich energy
The Helfrich energy is commonly used to model the elastic bending energy of lipid bilayers in membrane mechanics. The governing differential equations for certain geometric characteristics of the shape of the membrane can be obtained by applying variational methods (minimization principles) to the Helfrich energy functional and are well studied in the axisymmetric framework. However, the Helfrich energy functional and the resulting differential equations involve a number of parameters, and there is little explanation of the choice of parameters in the literature, particularly with respect to the choice of the “spontaneous curvature” term that appears in the functional. In this paper, we present a careful analytical and numerical study of certain aspects of parametric sensitivity of Helfrich’s model. Using simulations of specific model systems, we demonstrate the application of our scheme to the formation of spherical buds and pearled shapes in membrane vesicles.  more » « less
Award ID(s):
2012857
PAR ID:
10546996
Author(s) / Creator(s):
 ;  ;  
Publisher / Repository:
SAGE Publications
Date Published:
Journal Name:
Mathematics and Mechanics of Solids
Volume:
26
Issue:
3
ISSN:
1081-2865
Format(s):
Medium: X Size: p. 356-385
Size(s):
p. 356-385
Sponsoring Org:
National Science Foundation
More Like this
  1. ; ; (, Journal of Fluid Mechanics)
    null (Ed.)
    Cellular membranes are elastic lipid bilayers that contain a variety of proteins, including ion channels, receptors and scaffolding proteins. These proteins are known to diffuse in the plane of the membrane and to influence the bending of the membrane. Experiments have shown that lipid flow in the plane of the membrane is closely coupled with the diffusion of proteins. Thus, there is a need for a comprehensive framework that accounts for the interplay between these processes. Here, we present a theory for the coupled in-plane viscous flow of lipids, diffusion of transmembrane proteins and elastic deformation of lipid bilayers. The proteins in the membrane are modelled such that they influence membrane bending by inducing a spontaneous curvature. We formulate the free energy of the membrane with a Helfrich-like curvature elastic energy density function modified to account for the chemical potential energy of proteins. We derive the conservation laws and equations of motion for this system. Finally, we present results from dimensional analysis and numerical simulations and demonstrate the effect of coupled transport processes in governing the dynamics of membrane bending and protein diffusion. 
    more » « less
  2. ; ; ; ; ; (, Mathematical Medicine and Biology: A Journal of the IMA)
    null (Ed.)
    Abstract 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 (Randerson et al., 1982, Artif. Organs, 6, 43–49). We also show that the algebraic differential 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. 
    more » « less
  3. ; ; (, Digital Discovery)
    Rheology-informed neural networks (RhINNs) have recently been popularized as data-driven platforms for solving rheologically relevant differential equations. While RhINNs can be employed to solve different constitutive equations of interest in a forward or inverse manner, their ability to do so strictly depends on the type of data and the choice of models embedded within their structure. Here, the applicability of RhINNs in general, and the interplay between the choice of models, parameters of the neural network itself, and the type of data at hand are studied. To do so, a RhINN is informed by a series of thixotropic elasto-visco-plastic (TEVP) constitutive models, and its ability to accurately recover model parameters from stress growth and oscillatory shear flow protocols is investigated. We observed that by simplifying the constitutive model, RhINN convergence is improved in terms of parameter recovery accuracy and computation speed while over-simplifying the model is detrimental to accuracy. Moreover, several hyperparameters, e.g., the learning rate, activation function, initial conditions for the fitting parameters, and error heuristics, should be at the top of the checklist when aiming to improve parameter recovery using RhINNs. Finally, the given data form plays a pivotal role, and no convergence is observed when one set of experiments is used as the given data for either of the flow protocols. The range of parameters is also a limiting factor when employing RhINNs for parameter recovery, and ad hoc modifications to the constitutive model can be trivial remedies to guarantee convergence when recovering fitting parameters with large values. 
    more » « less
  4. ; ; ; ; (, AIChE Journal)
    Abstract Bi‐continuous jammed emulsion (bijel) membrane reactors, integrating simultaneous reaction and separation, offer a promising avenue for enhancing membrane reactor processes. In this study, we present a comprehensive macroscopic‐scale physicochemical model for tubular bijel membrane reactors and a numerical solution strategy for solving the governing partial differential equations. The model captures the co‐continuous network of two immiscible phases stabilized by nanoparticles at the liquid–liquid interface. We present the derivation of model equations and an efficient numerical solution strategy. The model is validated with experimental results from a conventional enzymatic biphasic membrane reactor for oleuropein hydrolysis, already reported in the literature. Simulation results indicate accurate prediction of reactor behavior, highlighting the potential superiority of bijel membrane reactors over current technologies. This research contributes a valuable tool for scale‐up, design, and optimization of bijel membrane reactors, filling a critical gap in this emerging field. 
    more » « less
  5. ; ; (, Inverse Problems)
    Abstract Bayesian particle filters (PFs) are a viable alternative to sampling methods such as Markov chain Monte Carlo methods to estimate model parameters and related uncertainties when the forward model is a dynamical system, and the data are time series that depend on the state vector. PF techniques are particularly attractive when the dimensionality of the state space is large and the numerical solution of the dynamical system over the time interval corresponding to the data is time consuming. Moreover, information contained in the PF solution can be used to infer on the sensitivity of the unknown parameters to different temporal segments of the data. This, in turn, can guide the design of more efficient and effective data collection procedures. In this article the PF method is applied to the problem of estimating cell membrane permeability to gases from pH measurements on or near the cell membrane. The forward model in this case comprises a spatially distributed system of coupled reaction–diffusion differential equations. The high dimensionality of the state space and the need to account for the micro-environment created by the pH electrode measurement device are additional challenges that are addressed by the solution method. 
    more » « less
  • Citation Formats
  • MLA
  • APA
  • Chicago
  • BibTeX
  • Save / Share this Record
  • Send to Email