1
|
Xu F, Zhang H, Fan Z, Zheng Z, Sun A, Zhang J. Investigating biomechanical alterations and emptying patterns after various gastrojejunostomy strategy. Food Res Int 2024; 176:113815. [PMID: 38163719 DOI: 10.1016/j.foodres.2023.113815] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/02/2023] [Revised: 11/22/2023] [Accepted: 12/02/2023] [Indexed: 01/03/2024]
Abstract
Gastrojejunostomy is a prominent approach in managing distal gastric cancer that is unresectable due to gastric outlet obstruction (GOO). Research has demonstrated that stomach-partitioning gastrojejunostomy (SPGJ) exhibits superior clinical efficacy compared to conventional gastrojejunostomy (CGJ), however, the underlying mechanism of this phenomenon remains elusive. This study constructed 3D models of the SPGJ and CGJ based on the computed tomography (CT) images obtained from a patient diagnosed with distal gastric cancer. The biomechanical patterns of these procedures in the digestive system were subsequently compared through numerical simulations and in vitro experiments. The results of the numerical simulation demonstrated that the model following SPGJ promoted the discharge of food through the anastomotic orifice and into the lower jejunum. Furthermore, a decrease in passage size after partitioning, the low-level velocity of esophageal, and an increase in contents viscosity effectively inhibited the flow through the passage to the pylorus, ultimately reducing stimulation to tumor. The study also revealed that favorable gastric emptying is associated with a smaller passage and faster inlet velocity, and that lower contents viscosity. The experimental findings conducted in vitro demonstrated that SPGJ exhibited superior efficacy in obstructing the flow near the pylorus in comparison to CGJ. Moreover, a decrease in passage size correlates with a reduction in fluid flow towards the pylorus. These results provide the foundation of theory and practice for the surgical management of patients with GOO resulting from unresectable distal gastric cancer, and have potential implications for clinical interventions.
Collapse
Affiliation(s)
- Fengyan Xu
- School of Mechanical Engineering, Jiangsu University of Technology, Changzhou, Jiangsu, China
| | - Haiqiao Zhang
- Department of General Surgery, Beijing Friendship Hospital, Capital Medical University, Beijing, China
| | - Zhenmin Fan
- School of Mechanical Engineering, Jiangsu University of Technology, Changzhou, Jiangsu, China.
| | - Zhi Zheng
- Department of General Surgery, Beijing Friendship Hospital, Capital Medical University, Beijing, China
| | - Anqiang Sun
- Key Laboratory for Biomechanics and Mechanobiology of Ministry of Education, School of Biological Science and Medical Engineering, Beihang University, Beijing 100191, China
| | - Jun Zhang
- Department of General Surgery, Beijing Friendship Hospital, Capital Medical University, Beijing, China.
| |
Collapse
|
2
|
Elisha G, Halder S, Carlson DA, Kahrilas PJ, Pandolfino JE, Patankar NA. A mechanics-based perspective on the pressure-cross-sectional area loop within the esophageal body. Front Physiol 2023; 13:1066351. [PMID: 36699676 PMCID: PMC9868904 DOI: 10.3389/fphys.2022.1066351] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/10/2022] [Accepted: 12/12/2022] [Indexed: 01/12/2023] Open
Abstract
Introduction: Plotting the pressure-cross-sectional area (P-CSA) hysteresis loops within the esophagus during a contraction cycle can provide mechanistic insights into esophageal motor function. Pressure and cross-sectional area during secondary peristalsis can be obtained from the functional lumen imaging probe (FLIP). The pressure-cross-sectional area plots at a location within the esophageal body (but away from the sphincter) reveal a horizontal loop shape. The horizontal loop shape has phases that appear similar to those in cardiovascular analyses, whichinclude isometric and isotonic contractions followed by isometric and isotonic relaxations. The aim of this study is to explain the various phases of the pressurecross-sectional area hysteresis loops within the esophageal body. Materials and Methods: We simulate flow inside a FLIP device placed inside the esophagus lumen. We focus on three scenarios: long functional lumen imaging probe bag placed insidethe esophagus but not passing through the lower esophageal sphincter, long functional lumen imaging probe bag that crosses the lower esophageal sphincter, and a short functional lumen imaging probe bag placed in the esophagus body that does not pass through the lower esophageal sphincter. Results and Discussion: Horizontal P-CSA area loop pattern is robust and is reproduced in all three cases with only small differences. The results indicate that the horizontal loop pattern is primarily a product of mechanical conditions rather than any inherently different function of the muscle itself. Thus, the distinct phases of the loop can be explained solely based on mechanics.
Collapse
Affiliation(s)
- Guy Elisha
- Department of Mechanical Engineering, McCormick School of Engineering, Northwestern University, Evanston, IL, United States
| | - Sourav Halder
- Theoretical and Applied Mechanics Program, McCormick School of Engineering, Northwestern University, Evanston, IL, United States
| | - Dustin A. Carlson
- Division of Gastroenterology and Hepatology, Feinberg School of Medicine, Northwestern University, Chicago, IL, United States
| | - Peter J. Kahrilas
- Division of Gastroenterology and Hepatology, Feinberg School of Medicine, Northwestern University, Chicago, IL, United States
| | - John E. Pandolfino
- Division of Gastroenterology and Hepatology, Feinberg School of Medicine, Northwestern University, Chicago, IL, United States
| | - Neelesh A. Patankar
- Department of Mechanical Engineering, McCormick School of Engineering, Northwestern University, Evanston, IL, United States,Theoretical and Applied Mechanics Program, McCormick School of Engineering, Northwestern University, Evanston, IL, United States,*Correspondence: Neelesh A. Patankar,
| |
Collapse
|
3
|
Lee JH, Griffith BE. On the Lagrangian-Eulerian Coupling in the Immersed Finite Element/Difference Method. JOURNAL OF COMPUTATIONAL PHYSICS 2022; 457:111042. [PMID: 35300097 PMCID: PMC8923617 DOI: 10.1016/j.jcp.2022.111042] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/17/2023]
Abstract
The immersed boundary (IB) method is a non-body conforming approach to fluid-structure interaction (FSI) that uses an Eulerian description of the momentum, viscosity, and incompressibility of a coupled fluid-structure system and a Lagrangian description of the deformations, stresses, and resultant forces of the immersed structure. Integral transforms with Dirac delta function kernels couple the Eulerian and Lagrangian variables, and in practice, discretizations of these integral transforms use regularized delta function kernels. Many different kernel functions have been proposed, but prior numerical work investigating the impact of the choice of kernel function on the accuracy of the methodology has often been limited to simplified test cases or Stokes flow conditions that may not reflect the method's performance in applications, particularly at intermediate-to-high Reynolds numbers, or under different loading conditions. This work systematically studies the effect of the choice of regularized delta function in several fluid-structure interaction benchmark tests using the immersed finite element/difference (IFED) method, which is an extension of the IB method that uses a finite element structural discretization combined with a Cartesian grid finite difference method for the incompressible Navier-Stokes equations. Whereas the conventional IB method spreads forces from the nodes of the structural mesh and interpolates velocities to those nodes, the IFED formulation evaluates the regularized delta function on a collection of interaction points that can be chosen to be denser than the nodes of the Lagrangian mesh. This opens the possibility of using structural discretizations with wide node spacings that would produce gaps in the Eulerian force in nodally coupled schemes (e.g., if the node spacing is comparable to or broader than the support of the regularized delta functions). Earlier work with this methodology suggested that such coarse structural meshes can yield improved accuracy for shear-dominated cases and, further, found that accuracy improves when the structural mesh spacing is increased. However, these results were limited to simple test cases that did not include substantial pressure loading on the structure. This study investigates the effect of varying the relative mesh widths of the Lagrangian and Eulerian discretizations in a broader range of tests. Our results indicate that kernels satisfying a commonly imposed even-odd condition require higher resolution to achieve similar accuracy as kernels that do not satisfy this condition. We also find that narrower kernels are more robust, in the sense that they yield results that are less sensitive to relative changes in the Eulerian and Lagrangian mesh spacings, and that structural meshes that are substantially coarser than the Cartesian grid can yield high accuracy for shear-dominated cases but not for cases with large normal forces. We verify our results in a large-scale FSI model of a bovine pericardial bioprosthetic heart valve in a pulse duplicator.
Collapse
Affiliation(s)
- Jae H Lee
- Department of Mathematics, University of North Carolina, Chapel Hill, NC, USA
- Department of Mechanical Engineering, Johns Hopkins University, Baltimore, MD, USA
- Institute for Computational Medicine, Johns Hopkins University, Baltimore, MD, USA
| | - Boyce E Griffith
- Departments of Mathematics, Applied Physical Sciences, and Biomedical Engineering, University of North Carolina, Chapel Hill, NC, USA
- Carolina Center for Interdisciplinary Applied Mathematics, University of North Carolina, Chapel Hill, NC, USA
- Computational Medicine Program, University of North Carolina School of Medicine, Chapel Hill, NC, USA
- McAllister Heart Institute, University of North Carolina School of Medicine, Chapel Hill, NC, USA
| |
Collapse
|
4
|
Halder S, Acharya S, Kou W, Campagna RAJ, Triggs JR, Carlson DA, Aadam AA, Hungness ES, Kahrilas PJ, Pandolfino JE, Patankar NA. Myotomy technique and esophageal contractility impact blown-out myotomy formation in achalasia: an in silico investigation. Am J Physiol Gastrointest Liver Physiol 2022; 322:G500-G512. [PMID: 35170365 PMCID: PMC8993593 DOI: 10.1152/ajpgi.00281.2021] [Citation(s) in RCA: 8] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 09/01/2021] [Revised: 01/11/2022] [Accepted: 02/10/2022] [Indexed: 01/31/2023]
Abstract
We used in silico models to investigate the impact of the dimensions of myotomy, contraction pattern, the tone of the esophagogastric junction (EGJ), and musculature at the myotomy site on esophageal wall stresses potentially leading to the formation of a blown-out myotomy (BOM). We performed three sets of simulations with an in silico esophagus model, wherein the myotomy-influenced region was modeled as an elliptical section devoid of muscle fibers. These sets investigated the effects of the dimensions of myotomy, differing esophageal contraction types, and differing esophagogastric junction (EGJ) tone and wall stiffness at the myotomy affected region on esophageal wall stresses potentially leading to BOM. Longer myotomy was found to be accompanied by a higher bolus volume accumulated at the myotomy site. With respect to esophageal contractions, deformation at the myotomy site was greatest with propagated peristalsis, followed by combined peristalsis and spasm, and pan-esophageal pressurization. Stronger EGJ tone with respect to the wall stiffness at the myotomy site was found to aid in increasing deformation at the myotomy site. In addition, we found that an esophagus with a shorter myotomy performed better at emptying the bolus than that with a longer myotomy. Shorter myotomies decrease the chance of BOM formation. Propagated peristalsis with EGJ outflow obstruction has the highest chance of BOM formation. We also found that abnormal residual EGJ tone may be a co-factor in the development of BOM, whereas remnant muscle fibers at myotomy site reduce the risk of BOM formation.NEW & NOTEWORTHY Blown-out myotomy (BOM) is a complication observed after myotomy, which is performed to treat achalasia. In silico simulations were performed to identify the factors leading to BOM formation. We found that a short myotomy that is not transmural and has some structural architecture intact reduces the risk of BOM formation. In addition, we found that high esophagogastric junction tone due to fundoplication is found to increase the risk of BOM formation.
Collapse
Affiliation(s)
- Sourav Halder
- Theoretical and Applied Mechanics Program, McCormick School of Engineering, Northwestern University, Evanston, Illinois
| | - Shashank Acharya
- Department of Mechanical Engineering, McCormick School of Engineering, Northwestern University, Evanston, Illinois
| | - Wenjun Kou
- Division of Gastroenterology and Hepatology, Department of Medicine, Feinberg School of Medicine, Northwestern University, Chicago, Illinois
| | - Ryan A J Campagna
- Department of Surgery, Feinberg School of Medicine, Northwestern University, Chicago, Illinois
| | - Joseph R Triggs
- Division of Gastroenterology, Perelman School of Medicine, University of Pennsylvania, Philadelphia, Pennsylvania
| | - Dustin A Carlson
- Division of Gastroenterology and Hepatology, Department of Medicine, Feinberg School of Medicine, Northwestern University, Chicago, Illinois
| | - Abdul Aziz Aadam
- Division of Gastroenterology and Hepatology, Department of Medicine, Feinberg School of Medicine, Northwestern University, Chicago, Illinois
| | - Eric S Hungness
- Department of Surgery, Feinberg School of Medicine, Northwestern University, Chicago, Illinois
| | - Peter J Kahrilas
- Division of Gastroenterology and Hepatology, Department of Medicine, Feinberg School of Medicine, Northwestern University, Chicago, Illinois
| | - John E Pandolfino
- Division of Gastroenterology and Hepatology, Department of Medicine, Feinberg School of Medicine, Northwestern University, Chicago, Illinois
| | - Neelesh A Patankar
- Theoretical and Applied Mechanics Program, McCormick School of Engineering, Northwestern University, Evanston, Illinois
- Department of Mechanical Engineering, McCormick School of Engineering, Northwestern University, Evanston, Illinois
| |
Collapse
|
5
|
Acharya S, Halder S, Kou W, Kahrilas PJ, Pandolfino JE, Patankar NA. A fully resolved multiphysics model of gastric peristalsis and bolus emptying in the upper gastrointestinal tract. Comput Biol Med 2022; 143:104948. [PMID: 35091365 PMCID: PMC9014465 DOI: 10.1016/j.compbiomed.2021.104948] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/30/2021] [Revised: 10/09/2021] [Accepted: 10/13/2021] [Indexed: 11/22/2022]
Abstract
Over the past few decades, in silico modeling of organ systems has significantly furthered our understanding of their physiology and biomechanical function. In spite of the relative importance of the digestive system in normal functioning of the human body, there is a scarcity of high-fidelity models for the upper gastrointestinal tract including the esophagus and the stomach. In this work, we present a detailed numerical model of the upper gastrointestinal tract that not only accounts for the fiber architecture of the muscle walls, but also the multiphasic components they help transport during normal digestive function. Construction details for 3D models of representative stomach geometry are presented along with a simple strategy for assigning circular and longitudinal muscle fiber orientations for each layer. We developed a fully resolved model of the stomach to simulate gastric peristalsis by systematically activating muscle fibers embedded in the stomach. Following this, for the first time, we simulate gravity-driven bolus emptying into the stomach due to density differences between ingested contents and fluid contents of the stomach. Finally, we present a case of retrograde flow of fluid from the stomach into the esophagus, resembling the phenomenon of acid reflux. This detailed computational model of the upper gastrointestinal tract provides a foundation for future models to investigate the biomechanics of acid reflux and probe various strategies for gastric bypass surgeries to address the growing problem of obesity.
Collapse
Affiliation(s)
- Shashank Acharya
- Department of Mechanical Engineering, Northwestern University, Evanston, IL, 60208, USA
| | - Sourav Halder
- Theoretical and Applied Mechanics Program, Northwestern University, Evanston, IL, 60208, USA
| | - Wenjun Kou
- Division of Gastroenterology and Hepatology, Feinberg School of Medicine, Northwestern University, Chicago, IL, 60611, USA
| | - Peter J Kahrilas
- Division of Gastroenterology and Hepatology, Feinberg School of Medicine, Northwestern University, Chicago, IL, 60611, USA
| | - John E Pandolfino
- Division of Gastroenterology and Hepatology, Feinberg School of Medicine, Northwestern University, Chicago, IL, 60611, USA
| | - Neelesh A Patankar
- Department of Mechanical Engineering, Northwestern University, Evanston, IL, 60208, USA; Theoretical and Applied Mechanics Program, Northwestern University, Evanston, IL, 60208, USA.
| |
Collapse
|
6
|
Kou W, Carlson DA, Baumann AJ, Donnan E, Luo Y, Pandolfino JE, Etemadi M. A deep-learning-based unsupervised model on esophageal manometry using variational autoencoder. Artif Intell Med 2021; 112:102006. [PMID: 33581826 DOI: 10.1016/j.artmed.2020.102006] [Citation(s) in RCA: 15] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/12/2020] [Revised: 10/19/2020] [Accepted: 12/28/2020] [Indexed: 12/27/2022]
Abstract
High-resolution manometry (HRM) is the primary method for diagnosing esophageal motility disorders and its interpretation and classification are based on variables (features) from data of each swallow. Modeling and learning the semantics directly from raw swallow data could not only help automate the feature extraction, but also alleviate the bias from pre-defined features. With more than 32-thousand raw swallow data, a generative model using the approach of variational auto-encoder (VAE) was developed, which, to our knowledge, is the first deep-learning-based unsupervised model on raw esophageal manometry data. The VAE model was reformulated to include different types of loss motivated by domain knowledge and tuned with different hyper-parameters. Training of the VAE model was found sensitive on the learning rate and hence the evidence lower bound objective (ELBO) was further scaled by the data dimension. Case studies showed that the dimensionality of latent space have a big impact on the learned semantics. In particular, cases with 4-dimensional latent variables were found to encode various physiologically meaningful contraction patterns, including strength, propagation pattern as well as sphincter relaxation. Cases with so-called hybrid L2 loss seemed to better capture the coherence of contraction/relaxation transition. Discriminating capability was further evaluated using simple linear discriminative analysis (LDA) on predicting swallow type and swallow pressurization, which yields clustering patterns consistent with clinical impression. The current work on modeling and understanding swallow-level data will guide the development of study-level models for automatic diagnosis as the next stage.
Collapse
Affiliation(s)
- Wenjun Kou
- Department of Medicine, Feinberg School of Medicine, Northwestern University, 676 North Saint Clair Street, 14th Floor, Chicago, IL 60611, USA.
| | - Dustin A Carlson
- Department of Medicine, Feinberg School of Medicine, Northwestern University, 676 North Saint Clair Street, 14th Floor, Chicago, IL 60611, USA
| | - Alexandra J Baumann
- Department of Medicine, Feinberg School of Medicine, Northwestern University, 676 North Saint Clair Street, 14th Floor, Chicago, IL 60611, USA
| | - Erica Donnan
- Department of Medicine, Feinberg School of Medicine, Northwestern University, 676 North Saint Clair Street, 14th Floor, Chicago, IL 60611, USA
| | - Yuan Luo
- Department of Preventive Medicine, Feinberg School of Medicine, Northwestern University, 750 North Lake Shore Drive, 11th Floor, Chicago, IL 60611, USA
| | - John E Pandolfino
- Department of Medicine, Feinberg School of Medicine, Northwestern University, 676 North Saint Clair Street, 14th Floor, Chicago, IL 60611, USA
| | - Mozziyar Etemadi
- Department of Anesthesiology, Feinberg School of Medicine, Northwestern University, Chicago, IL 60611, USA; Department of Biomedical Engineering, McCormick School of Engineering, Northwestern University, Evanston, IL 60201, USA
| |
Collapse
|
7
|
Kolahdouz EM, Bhalla APS, Craven BA, Griffith BE. An Immersed Interface Method for Discrete Surfaces. JOURNAL OF COMPUTATIONAL PHYSICS 2020; 400:108854. [PMID: 31802781 PMCID: PMC6892596 DOI: 10.1016/j.jcp.2019.07.052] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/10/2023]
Abstract
Fluid-structure systems occur in a range of scientific and engineering applications. The immersed boundary (IB) method is a widely recognized and effective modeling paradigm for simulating fluid-structure interaction (FSI) in such systems, but a difficulty of the IB formulation of these problems is that the pressure and viscous stress are generally discontinuous at fluid-solid interfaces. The conventional IB method regularizes these discontinuities, which typically yields low-order accuracy at these interfaces. The immersed interface method (IIM) is an IB-like approach to FSI that sharply imposes stress jump conditions, enabling higher-order accuracy, but prior applications of the IIM have been largely restricted to numerical methods that rely on smooth representations of the interface geometry. This paper introduces an immersed interface formulation that uses only a C 0 representation of the immersed interface, such as those provided by standard nodal Lagrangian finite element methods. Verification examples for models with prescribed interface motion demonstrate that the method sharply resolves stress discontinuities along immersed boundaries while avoiding the need for analytic information about the interface geometry. Our results also demonstrate that only the lowest-order jump conditions for the pressure and velocity gradient are required to realize global second-order accuracy. Specifically, we demonstrate second-order global convergence rates along with nearly second-order local convergence in the Eulerian velocity field, and between first- and second-order global convergence rates along with approximately first-order local convergence for the Eulerian pressure field. We also demonstrate approximately second-order local convergence in the interfacial displacement and velocity along with first-order local convergence in the fluid traction along the interface. As a demonstration of the method's ability to tackle more complex geometries, the present approach is also used to simulate flow in a patient-averaged anatomical model of the inferior vena cava, which is the large vein that carries deoxygenated blood from the lower extremities back to the heart. Comparisons of the general hemodynamics and wall shear stress obtained by the present IIM and a body-fitted discretization approach show that the present method yields results that are in good agreement with those obtained by the body-fitted approach.
Collapse
Affiliation(s)
- Ebrahim M Kolahdouz
- Department of Mathematics, University of North Carolina, Chapel Hill, NC, USA
- Division of Applied Mechanics, Office of Science and Engineering Laboratories, Center for Devices and Radiological Health, United States Food and Drug Administration, Silver Spring, MD, USA
| | | | - Brent A Craven
- Division of Applied Mechanics, Office of Science and Engineering Laboratories, Center for Devices and Radiological Health, United States Food and Drug Administration, Silver Spring, MD, USA
| | - Boyce E Griffith
- Departments of Mathematics, Applied Physical Sciences, and Biomedical Engineering, University of North Carolina, Chapel Hill, NC, USA
- Carolina Center for Interdisciplinary Applied Mathematics, University of North Carolina, Chapel Hill, NC, USA
- McAllister Heart Institute, University of North Carolina, Chapel Hill, NC, USA
| |
Collapse
|
8
|
Griffith BE, Patankar NA. Immersed Methods for Fluid-Structure Interaction. ANNUAL REVIEW OF FLUID MECHANICS 2019; 52:421-448. [PMID: 33012877 PMCID: PMC7531444 DOI: 10.1146/annurev-fluid-010719-060228] [Citation(s) in RCA: 38] [Impact Index Per Article: 7.6] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/24/2023]
Abstract
Fluid-structure interaction is ubiquitous in nature and occurs at all biological scales. Immersed methods provide mathematical and computational frameworks for modeling fluid-structure systems. These methods, which typically use an Eulerian description of the fluid and a Lagrangian description of the structure, can treat thin immersed boundaries and volumetric bodies, and they can model structures that are flexible or rigid or that move with prescribed deformational kinematics. Immersed formulations do not require body-fitted discretizations and thereby avoid the frequent grid regeneration that can otherwise be required for models involving large deformations and displacements. This article reviews immersed methods for both elastic structures and structures with prescribed kinematics. It considers formulations using integral operators to connect the Eulerian and Lagrangian frames and methods that directly apply jump conditions along fluid-structure interfaces. Benchmark problems demonstrate the effectiveness of these methods, and selected applications at Reynolds numbers up to approximately 20,000 highlight their impact in biological and biomedical modeling and simulation.
Collapse
Affiliation(s)
- Boyce E Griffith
- Departments of Mathematics, Applied Physical Sciences, and Biomedical Engineering, University of North Carolina, Chapel Hill, North Carolina 27599, USA
| | - Neelesh A Patankar
- Department of Mechanical Engineering, Northwestern University, Evanston, Illinois 60208, USA
| |
Collapse
|