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.

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. (, Applied Sciences)
    Mathematical analysis of the well known model of a piezoelectric energy harvester is presented. The harvester is designed as a cantilever Timoshenko beam with piezoelectric layers attached to its top and bottom faces. Thin, perfectly conductive electrodes are covering the top and bottom faces of the piezoelectric layers. These electrodes are connected to a resistive load. The model is governed by a system of three partial differential equations. The first two of them are the equations of the Timoshenko beam model and the third one represents Kirchhoff’s law for the electric circuit. All equations are coupled due to the piezoelectric effect. We represent the system as a single operator evolution equation in the Hilbert state space of the system. The dynamics generator of this evolution equation is a non-selfadjoint matrix differential operator with compact resolvent. The paper has two main results. Both results are explicit asymptotic formulas for eigenvalues of this operator, i.e., the modal analysis for the electrically loaded system is performed. The first set of the asymptotic formulas has remainder terms of the order O ( 1 n ) , where n is the number of an eigenvalue. These formulas are derived for the model with variable physical parameters. The second set of the asymptotic formulas has remainder terms of the order O ( 1 n 2 ) , and is derived for a less general model with constant parameters. This second set of formulas contains extra term depending on the electromechanical parameters of the model. It is shown that the spectrum asymptotically splits into two disjoint subsets, which we call the α -branch eigenvalues and the θ -branch eigenvalues. These eigenvalues being multiplied by “i” produce the set of the vibrational modes of the system. The α -branch vibrational modes are asymptotically located on certain vertical line in the left half of the complex plane and the θ -branch is asymptotically close to the imaginary axis. By having such spectral and asymptotic results, one can derive the asymptotic representation for the mode shapes and for voltage output. Asymptotics of vibrational modes and mode shapes is instrumental in the analysis of control problems for the harvester. 
    more » « less
  5. ; ; ; ; (, 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
  • Citation Formats
  • MLA
  • APA
  • Chicago
  • BibTeX
  • Save / Share this Record
  • Send to Email